楊叢彥，趙丹，胡海濤，丁國(guó)良

（上海交通大學(xué)制冷與低溫工程研究所，上海 200240）

多溫區(qū)冰箱用并聯(lián)雙蒸發(fā)器動(dòng)態(tài)性能計(jì)算方法

楊叢彥*，趙丹，胡海濤**，丁國(guó)良

（上海交通大學(xué)制冷與低溫工程研究所，上海 200240）

多溫區(qū)冰箱中采用并聯(lián)雙蒸發(fā)器可以有效地提高各間室的溫度控制精度和冰箱的能效。冰箱并聯(lián)雙蒸發(fā)器運(yùn)行模式復(fù)雜，現(xiàn)有的蒸發(fā)器模型不能直接用于并聯(lián)雙蒸發(fā)器的動(dòng)態(tài)仿真，需要開發(fā)并聯(lián)雙蒸發(fā)器動(dòng)態(tài)性能的算法。本文總結(jié)了并聯(lián)雙蒸發(fā)器所有的運(yùn)行模式以及模式轉(zhuǎn)換方式；分析不同運(yùn)行模式下蒸發(fā)器中制冷劑的分布狀態(tài)，建立不同分布狀態(tài)下質(zhì)量、能量及其他狀態(tài)參數(shù)的計(jì)算公式以及各模式光滑轉(zhuǎn)換準(zhǔn)則；基于以上分析開發(fā)了用于并聯(lián)雙蒸發(fā)器所有運(yùn)行模式下動(dòng)態(tài)性能的計(jì)算方法。最后將本文提出的方法嵌入到成熟的仿真軟件中，仿真計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)相對(duì)偏差均在10%以內(nèi)。

并聯(lián)雙蒸發(fā)器；冰箱；動(dòng)態(tài)性能；快速；穩(wěn)定

0 引言

多溫區(qū)冰箱中采用并聯(lián)雙蒸發(fā)器可以有效地提高各間室的溫度控制精度和冰箱能效[1-5]。

并聯(lián)雙蒸發(fā)器冰箱兩間室所需要的制冷量不同，系統(tǒng)通過(guò)調(diào)節(jié)冷藏和冷凍毛細(xì)管的節(jié)流效應(yīng)，改變兩蒸發(fā)器的換熱溫差以及制冷劑流量從而滿足兩蒸發(fā)器不同換熱量，為了使兩蒸發(fā)器在不同的換熱溫差以及制冷劑流量下?lián)Q熱性能達(dá)到最優(yōu)，需要對(duì)兩蒸發(fā)器的結(jié)構(gòu)進(jìn)行設(shè)計(jì)。在并聯(lián)雙蒸發(fā)器冰箱實(shí)際運(yùn)行過(guò)程中，兩蒸發(fā)器中制冷劑的狀態(tài)參數(shù)時(shí)刻發(fā)生變化，導(dǎo)致其動(dòng)態(tài)性能時(shí)刻變化，因此蒸發(fā)器的結(jié)構(gòu)設(shè)計(jì)需要考慮動(dòng)態(tài)性能。對(duì)于現(xiàn)有的經(jīng)驗(yàn)設(shè)計(jì)方法，需要反復(fù)制作樣機(jī)并進(jìn)行性能測(cè)試，不僅設(shè)計(jì)效率低，而且還會(huì)增加企業(yè)的設(shè)計(jì)成本。而計(jì)算機(jī)仿真方法相較傳統(tǒng)的設(shè)計(jì)方式，具有設(shè)計(jì)周期短、設(shè)計(jì)成本低等優(yōu)點(diǎn)，在制冷裝置動(dòng)態(tài)性能的計(jì)算得到了廣泛應(yīng)用[6-11]。因此，為了對(duì)并聯(lián)雙蒸發(fā)器冰箱進(jìn)行優(yōu)化設(shè)計(jì)，必須開發(fā)冰箱用并聯(lián)雙蒸發(fā)器動(dòng)態(tài)性能計(jì)算方法。

并聯(lián)雙蒸發(fā)器要實(shí)現(xiàn)計(jì)算機(jī)動(dòng)態(tài)仿真需要考慮并聯(lián)雙蒸發(fā)器的實(shí)際運(yùn)行情況。并聯(lián)蒸發(fā)器的運(yùn)行模式復(fù)雜，共有4種運(yùn)行模式包括：1）冷凍蒸發(fā)器開機(jī)、冷藏蒸發(fā)器停機(jī)；2）冷藏蒸發(fā)器開機(jī)、冷凍蒸發(fā)器停機(jī)；3）兩蒸發(fā)器抽真空；4）兩蒸發(fā)器停機(jī)。在冰箱實(shí)際運(yùn)行過(guò)程中，通過(guò)控制4種運(yùn)行模式相互切換，實(shí)現(xiàn)給兩間室制冷的目的。為了避免冷藏蒸發(fā)器開機(jī)過(guò)程中冷凍蒸發(fā)器內(nèi)還殘存制冷劑的問題，冷凍循環(huán)結(jié)束之后通過(guò)抽真空過(guò)程實(shí)現(xiàn)制冷劑的再分配而不能直接切換至冷藏循環(huán)，因此共有11種模式切換方式。制冷劑分布狀態(tài)在蒸發(fā)器運(yùn)行模式切換前后會(huì)發(fā)生動(dòng)態(tài)變化，因此，冰箱用并聯(lián)雙蒸發(fā)器動(dòng)態(tài)性能計(jì)算方法需要能夠描述兩蒸發(fā)器在11種模式切換方式下制冷劑的分布狀態(tài)，并能夠計(jì)算在不同運(yùn)行模式切換下兩蒸發(fā)器的動(dòng)態(tài)性能。

現(xiàn)有的蒸發(fā)器動(dòng)態(tài)仿真計(jì)算方法[12-14]共有3種：分布參數(shù)計(jì)算方法、集中參數(shù)計(jì)算方法和分相集中參數(shù)計(jì)算方法。分布參數(shù)方法通過(guò)將蒸發(fā)器劃分為大量的微元進(jìn)行計(jì)算，在穩(wěn)態(tài)仿真時(shí)具有計(jì)算精度高的優(yōu)點(diǎn)，但存在計(jì)算速度慢的問題；集中參數(shù)計(jì)算方法通過(guò)采用均相假設(shè)，具有計(jì)算速度快、計(jì)算穩(wěn)定等優(yōu)點(diǎn)，但計(jì)算精度不高；分相集中參數(shù)方法依據(jù)制冷劑的分布狀態(tài)將蒸發(fā)器劃分為多個(gè)相區(qū)的組合，在現(xiàn)有的計(jì)算方法中具有綜合優(yōu)勢(shì)，可以快速、穩(wěn)定地仿真單蒸發(fā)器和串聯(lián)雙蒸發(fā)器的開機(jī)制冷以及停機(jī)兩種運(yùn)行模式下的動(dòng)態(tài)性能。上述已有研究主要是針對(duì)單蒸發(fā)器和串聯(lián)雙蒸發(fā)器，到目前為止，關(guān)于并聯(lián)雙蒸發(fā)器計(jì)算機(jī)動(dòng)態(tài)性能仿真的研究非常匱乏。并聯(lián)雙蒸發(fā)器相較單蒸發(fā)器和串聯(lián)雙蒸發(fā)器運(yùn)行模式更加復(fù)雜，同時(shí)在實(shí)際運(yùn)行過(guò)程中，不同運(yùn)行模式還需要相互切換，因此現(xiàn)有的用于單蒸發(fā)器和串聯(lián)雙蒸發(fā)器的動(dòng)態(tài)計(jì)算方法不能直接用于計(jì)算并聯(lián)雙蒸發(fā)器復(fù)雜的動(dòng)態(tài)性能。

并聯(lián)雙蒸發(fā)器動(dòng)態(tài)仿真方法主要有以下難點(diǎn)。1）制冷劑狀態(tài)的確定：在并聯(lián)雙蒸發(fā)器仿真過(guò)程中，制冷劑的動(dòng)態(tài)變化主要發(fā)生在運(yùn)行模式切換的過(guò)程中，因此描述所有運(yùn)行模式下蒸發(fā)器中制冷劑的狀態(tài)存在難點(diǎn)。

2）保證運(yùn)行模式切換中蒸發(fā)器數(shù)學(xué)模型光滑轉(zhuǎn)換：動(dòng)態(tài)性能仿真計(jì)算方法需要描述制冷劑不同分布狀態(tài)下的數(shù)學(xué)模型，在11種模式切換方式中，數(shù)學(xué)模型間的光滑連接是保證系統(tǒng)仿真穩(wěn)定性的關(guān)鍵。

3）保證蒸發(fā)器動(dòng)態(tài)仿真結(jié)果穩(wěn)定的算法：在并聯(lián)雙蒸發(fā)器不同的運(yùn)行過(guò)程中，兩個(gè)蒸發(fā)器之間有質(zhì)量和能量的傳遞，傳遞過(guò)程復(fù)雜，計(jì)算方法需要保證質(zhì)量和能量的嚴(yán)格守恒，這就給雙蒸發(fā)器算法的開發(fā)帶來(lái)困難。

本文為解決以上難點(diǎn)，開發(fā)了冰箱并聯(lián)雙蒸發(fā)器動(dòng)態(tài)性能計(jì)算方法，從而保證快速、精確地實(shí)現(xiàn)并聯(lián)雙蒸發(fā)器冰箱系統(tǒng)的動(dòng)態(tài)性能的仿真計(jì)算。

1 開發(fā)動(dòng)態(tài)性能計(jì)算方法

并聯(lián)雙蒸發(fā)器冰箱系統(tǒng)包含冷藏循環(huán)制冷回路（如圖1中a-b-c-d-e-a）和冷凍循環(huán)制冷回路（如圖1中a-b-c-g-f-a），系統(tǒng)控制策略通過(guò)三通電磁閥控制制冷劑的流向，實(shí)現(xiàn)了兩制冷回路的交替運(yùn)行。

圖1 并聯(lián)雙蒸發(fā)器冰箱的循環(huán)示意圖

并聯(lián)雙蒸發(fā)器的運(yùn)行模式包括：

1）冷凍開機(jī)冷藏停機(jī)：壓縮機(jī)開機(jī)，三通電磁閥的出口1關(guān)閉，制冷劑的流向?yàn)閍-b-c-g-f-a，冰箱給冷凍室進(jìn)行制冷；

2）冷藏開機(jī)冷凍停機(jī)：壓縮機(jī)開機(jī)，三通電磁閥的出口2關(guān)閉，制冷劑的流向?yàn)閍-b-c-d-e-a，冰箱給冷藏室進(jìn)行制冷；

3）兩蒸發(fā)器抽真空：壓縮機(jī)保持開機(jī)，三通電磁閥的兩出口均關(guān)閉，制冷劑通過(guò)壓縮機(jī)，從兩個(gè)蒸發(fā)器進(jìn)入冷凝器；

4）兩蒸發(fā)器停機(jī)：壓縮機(jī)關(guān)機(jī)。

動(dòng)態(tài)性能計(jì)算方法就是要快速、穩(wěn)定地計(jì)算這四種運(yùn)行模式循環(huán)切換下蒸發(fā)器的動(dòng)態(tài)性能。

1.1 計(jì)算思路

為了能夠快速、精確地計(jì)算兩蒸發(fā)器中制冷劑的狀態(tài)參數(shù)，本文針對(duì)并聯(lián)雙蒸發(fā)器的運(yùn)行特點(diǎn)提出如下計(jì)算思路。

1）明確并聯(lián)雙蒸發(fā)器在所有運(yùn)行模式下制冷劑的分布狀態(tài)?？偨Y(jié)并聯(lián)雙蒸發(fā)器所有的運(yùn)行模式以及運(yùn)行模式切換過(guò)程，并分析制冷劑在蒸發(fā)器所有運(yùn)行模式轉(zhuǎn)換前、后的分布狀態(tài)。

2）建立并聯(lián)雙蒸發(fā)器在所有分布狀態(tài)下的計(jì)算模型。開發(fā)并聯(lián)雙蒸發(fā)器中制冷劑不同分布狀態(tài)下具體的計(jì)算公式，實(shí)現(xiàn)所有分布狀態(tài)下制冷劑狀態(tài)參數(shù)計(jì)算。

3）建立并聯(lián)雙蒸發(fā)器不同運(yùn)行模式光滑轉(zhuǎn)換判斷準(zhǔn)則。根據(jù)制冷劑各分布狀態(tài)的計(jì)算模型計(jì)算得到的制冷劑的狀態(tài)參數(shù)應(yīng)滿足光滑轉(zhuǎn)換判斷準(zhǔn)則，從而保證并聯(lián)雙蒸發(fā)器仿真計(jì)算過(guò)程中各運(yùn)行模式光滑轉(zhuǎn)換。

4）開發(fā)實(shí)現(xiàn)不同運(yùn)行模式下并聯(lián)雙蒸發(fā)器質(zhì)量和能量嚴(yán)格守恒的算法，穩(wěn)定地計(jì)算制冷劑在并聯(lián)雙蒸發(fā)器不同運(yùn)行模式切換過(guò)程的質(zhì)量、能量及其他狀態(tài)參數(shù)，保證仿真計(jì)算過(guò)程收斂。

1.2 明確制冷劑分布狀態(tài)

在并聯(lián)雙蒸發(fā)器不同運(yùn)行模式下，蒸發(fā)器中制冷劑分布狀態(tài)有4種基本狀態(tài)，包括SH（全過(guò)熱）、SH+TP（過(guò)熱和兩相共存）、TP1（進(jìn)出口均為兩相）和TP2（進(jìn)口兩相/出口飽和），如圖2所示。

圖2 蒸發(fā)器中制冷劑基本狀態(tài)

在冰箱實(shí)際運(yùn)行過(guò)程中，通過(guò)控制4種運(yùn)行模式相互切換，即任意一種運(yùn)行模式切換至其他3種運(yùn)行模式，實(shí)現(xiàn)給兩箱室制冷的目的。為了避免冷藏蒸發(fā)器開機(jī)過(guò)程中冷凍蒸發(fā)器內(nèi)還殘存制冷劑的問題，冷凍循環(huán)結(jié)束之后要通過(guò)抽真空過(guò)程實(shí)現(xiàn)制冷劑的再分配而不能直接切換至冷藏循環(huán)，因此4種運(yùn)行模式相互轉(zhuǎn)換的共有11種可能情況。在不同運(yùn)行模式切換過(guò)程中，兩蒸發(fā)器中制冷劑的分布狀態(tài)會(huì)發(fā)生改變，下面將分別列出11種模式轉(zhuǎn)換方式下并聯(lián)雙蒸發(fā)器轉(zhuǎn)換前后的運(yùn)行模式，以及模式轉(zhuǎn)換前后兩蒸發(fā)器中制冷劑的分布狀態(tài)，如表1所示。

表1 制冷劑在兩蒸發(fā)器中的狀態(tài)變化

1.3 建立制冷劑不同分布狀態(tài)下質(zhì)量、能量及其它狀態(tài)參數(shù)的計(jì)算公式

在并聯(lián)雙蒸發(fā)器實(shí)際運(yùn)行過(guò)程中，2個(gè)蒸發(fā)器之間有質(zhì)量和能量的傳遞，為保證計(jì)算算法質(zhì)量和能量的嚴(yán)格守恒，將兩蒸發(fā)器的總質(zhì)量和總能量作為狀態(tài)參數(shù)，這樣質(zhì)量和能量守恒方程就作為基本控制方程，如公式(1)和公式(2)所示。

式中：

E——制冷劑熱力學(xué)能，kJ；

M——制冷劑質(zhì)量，kg；

min——蒸發(fā)器進(jìn)口質(zhì)量流量，kg/s；

mout——蒸發(fā)器出口質(zhì)量流量，kg/s；

hin——蒸發(fā)器進(jìn)口焓，kJ/kg；

hout——蒸發(fā)器出口焓，kJ/kg；

Q——制冷劑側(cè)換熱量，kJ；

t——時(shí)間，s。

本文提出了制冷劑不同分布狀態(tài)下的質(zhì)量、能量及其他狀態(tài)參數(shù)的計(jì)算公式。下面將分別列出制冷劑4種分布狀態(tài)下的計(jì)算公式。

1.3.1 制冷劑為過(guò)熱狀態(tài)（SH）

式中：

A——蒸發(fā)器截面積，m2；

L——蒸發(fā)器長(zhǎng)度，m；

p——蒸發(fā)壓力，kPa；

1.3.2 制冷劑為兩相和過(guò)熱共存狀態(tài)（SH+TP）

式中：

C——系數(shù)；

xi，x0——相區(qū)進(jìn)出口干度；

ρg，ρ1——蒸發(fā)器氣態(tài)、液態(tài)密度，kg/m3；

hg，h1——蒸發(fā)器氣態(tài)、液態(tài)焓，kJ/kg；

1.3.3 制冷劑為進(jìn)出口均為兩相（TP1）狀態(tài)或制冷劑為進(jìn)口兩相/出口飽和（TP2）狀態(tài)

各模式的判斷準(zhǔn)則保證了各模式之間的光滑轉(zhuǎn)換。根據(jù)焓值線性分布假設(shè)，具體的各模式的判斷準(zhǔn)則如表2所示。

表2 制冷劑基本狀態(tài)判斷準(zhǔn)則

1.4 算法實(shí)現(xiàn)

確定兩蒸發(fā)器在不同運(yùn)行模式下狀態(tài)參數(shù)的計(jì)算方法如下：

1) 首先輸入邊界條件制冷劑進(jìn)口質(zhì)量流量min、出口質(zhì)量流量mout、進(jìn)口焓hin、出口焓hout以及制冷劑側(cè)換熱量Q；

2) 根據(jù)式(1)～式(2)積分計(jì)算第i時(shí)刻制冷劑的總質(zhì)量Mi和總能量Ei；

3) 輸入制冷劑的總質(zhì)量Mi、總能量Ei；

4) 并聯(lián)雙蒸發(fā)器運(yùn)行模式切換共分為11種，不同切換過(guò)程，兩蒸發(fā)器中制冷劑的狀態(tài)參數(shù)確定流程不同，如果切換模式為模式1～3中的1種，轉(zhuǎn)到步驟5）；如果切換模式為模式4～8中的1種，比較冷藏、冷凍蒸發(fā)器內(nèi)壓力，當(dāng)冷藏蒸發(fā)壓力大于冷凍蒸發(fā)壓力轉(zhuǎn)到步驟6）；否則轉(zhuǎn)到步驟5）；如果切換模式為模式9～10中的1種，比較冷藏、冷凍蒸發(fā)器內(nèi)壓力，當(dāng)冷凍蒸發(fā)壓力大于冷藏蒸發(fā)壓力轉(zhuǎn)到步驟6）；否則轉(zhuǎn)到步驟8）；如果切換模式為模式11，轉(zhuǎn)到步驟8）；

5) 假定壓力，分別計(jì)算冷藏蒸發(fā)器和冷凍蒸發(fā)器中制冷劑的狀態(tài)參數(shù)，如果計(jì)算的兩蒸發(fā)器內(nèi)制冷劑的總質(zhì)量（總能量）和實(shí)際總質(zhì)量（總能量）相等時(shí)，計(jì)算收斂，轉(zhuǎn)到步驟10）；否則調(diào)整壓力假定值，重新計(jì)算；

6) 假設(shè)蒸發(fā)器內(nèi)壓力，計(jì)算冷凍蒸發(fā)器中制冷劑狀態(tài)參數(shù)，如果計(jì)算的總能量與實(shí)際總能量相等時(shí)，計(jì)算收斂，轉(zhuǎn)到步驟7），否則調(diào)整壓力假定值，重新計(jì)算；

7) 假設(shè)蒸發(fā)器內(nèi)壓力，計(jì)算冷藏蒸發(fā)器中制冷劑狀態(tài)參數(shù)，如果計(jì)算的總質(zhì)量與實(shí)際總質(zhì)量相等時(shí)，計(jì)算收斂，轉(zhuǎn)到步驟10）；否則調(diào)整壓力假定值，重新計(jì)算；

8) 假設(shè)壓力，計(jì)算冷凍蒸發(fā)器中制冷劑的狀態(tài)參數(shù)，如果計(jì)算的冷凍蒸發(fā)器內(nèi)制冷劑的總能量與其實(shí)際總能量相等時(shí)，計(jì)算收斂，轉(zhuǎn)到步驟9）；否則調(diào)整壓力假定值，重新計(jì)算；

9) 假設(shè)蒸發(fā)器內(nèi)壓力，計(jì)算冷藏蒸發(fā)器中制冷劑狀態(tài)參數(shù)，如果計(jì)算的總能量與實(shí)際總能量相等時(shí)，計(jì)算收斂，轉(zhuǎn)到步驟10），否則調(diào)整壓力假定值，重新計(jì)算；

10) 輸出蒸發(fā)器中制冷劑的狀態(tài)參數(shù)，計(jì)算i+1時(shí)刻制冷劑的狀態(tài)參數(shù)。

具體算法如圖3所示。

圖3 并聯(lián)雙蒸發(fā)器計(jì)算流程

基于式(3)～式(17)，確定蒸發(fā)器中制冷劑分布狀態(tài)和狀態(tài)參數(shù)的計(jì)算流程如下。

1）輸入已知條件。

8）采用式(11)～式(14)及式(15)～式(17)計(jì)算制冷劑的狀態(tài)參數(shù)，轉(zhuǎn)到9）。

9）輸出制冷劑質(zhì)量、能量及蒸發(fā)壓力等參數(shù)。

具體計(jì)算流程如圖4所示。

圖4 蒸發(fā)器狀態(tài)參數(shù)計(jì)算流程

2 算法驗(yàn)證

為了證明計(jì)算方法的有效性，本文將提出的計(jì)算方法嵌入到成熟的冰箱系統(tǒng)仿真軟件[15]中，仿真軟件在個(gè)人電腦（Intel Pentium D CPU 2.66 GHz；1024 MB RAM）模擬測(cè)試冰箱實(shí)際運(yùn)行24小時(shí)的動(dòng)態(tài)性能，仿真所花費(fèi)的時(shí)間在300秒以內(nèi)，并將仿真結(jié)果中蒸發(fā)器中制冷劑的狀態(tài)參數(shù)與實(shí)驗(yàn)測(cè)試結(jié)果[16]進(jìn)行比較，其中溫度測(cè)量精度為±0.3oC，制冷劑壓力測(cè)量精度為±0.2%。用于測(cè)試的蒸發(fā)器參數(shù)如表3所示。

表3 兩蒸發(fā)器的結(jié)構(gòu)參數(shù)

翅片間距/mm 6.35 8.47

采用本文提出的并聯(lián)雙蒸發(fā)器計(jì)算方法仿真冰箱蒸發(fā)器正常工作時(shí)的性能，仿真結(jié)果和實(shí)驗(yàn)結(jié)果吻合很好，結(jié)果顯示仿真值相對(duì)于實(shí)驗(yàn)值的相對(duì)偏差在均在10%以內(nèi)，如圖5所示。

圖5 兩蒸發(fā)器出口溫度仿真計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果對(duì)比

3 結(jié)論

1）通過(guò)將制冷劑的分布狀態(tài)劃分為SH（全過(guò)熱）、SH+TP（兩相和過(guò)熱共存）、TP1（進(jìn)出口均為兩相）以及TP2（進(jìn)口兩相/出口飽和）4種分布狀態(tài)，給出制冷劑四種分布狀態(tài)下質(zhì)量、能量及其他狀態(tài)參數(shù)的計(jì)算公式，可以實(shí)現(xiàn)并聯(lián)雙蒸發(fā)器所有運(yùn)行模式下動(dòng)態(tài)性能計(jì)算；

3）提出了基于質(zhì)量和能量守恒的并聯(lián)雙蒸發(fā)器動(dòng)態(tài)性能算法，計(jì)算蒸發(fā)器正常工作時(shí)的動(dòng)態(tài)性能與實(shí)驗(yàn)結(jié)果吻合良好，結(jié)果顯示仿真值和實(shí)驗(yàn)值的相對(duì)偏差在均在10%以內(nèi)。在個(gè)人電腦仿真24 h性能所花費(fèi)的時(shí)間仿真小于300 s，滿足設(shè)計(jì)要求。

[1] YOON W J, JUNG H W, CHUNG H J, et al. Performance optimization of a two-circuit cycle with parallel evaporators for a domestic refrigerator–freezer[J]. International Journal of Refrigeration, 2011, 34(1): 216-224.

[2] YOON W J, SEO K, CHUMG H J, et al. Performance optimization of dual-loop cycles using R-600a and hydrocarbon mixtures designed for a domestic refrigerator-freezer[J]. International Journal of Refrigeration, 2012, 35(6): 1657-1667.

[3] WON S, JUNG D, RADERMACHER R. An experimental study of the performance of a dual-loop refrigerator freezer system[J]. International Journal of Refrigeration, 1994, 17(6): 411-416.

[4] 張艷. 三溫室冰箱的節(jié)能方案研究[J]. 制冷技術(shù), 2011, 31(1): 15-17.

[5] 梅勤光. 家用冰箱的節(jié)能減排技術(shù)[J]. 制冷技術(shù), 2008, 28(4): 12-14.

[6] 盧智利, 丁國(guó)良, 張春路. 雙溫雙控冰箱的仿真與應(yīng)用[J]. 上海交通大學(xué)學(xué)報(bào), 2004, 38(2): 268-272.

[7] 丁國(guó)良, 陳芝久, 闕雄才, 等. 制冷空調(diào)裝置設(shè)計(jì)方法現(xiàn)代化與智能化仿真探索[J]. 制冷技術(shù), 1999, 19(1): 10-14.

[8] LU Z L, DING G L, ZHANG C L. Dynamic simulation of natural convection bypass two-circuit cycle refrigerator-freezer and its application Part II: system simulation and application[J]. Applied Thermal Engineering, 2004, 24(10): 1525-1533.

[9] 孫浩然, 任滔, 李智強(qiáng), 等. 結(jié)合用戶數(shù)據(jù)的空調(diào)器仿真平臺(tái)構(gòu)建[J]. 制冷技術(shù), 2014, 34(4): 31-37.

[10] 陸震. 以制冷節(jié)能促進(jìn)可持續(xù)發(fā)展[J]. 制冷技術(shù), 2005, 25(3): 33-37.

[11] 付德鋼, 陸震, 曹衛(wèi)華, 等. 直燃型雙效溴化鋰吸收式制冷機(jī)擾動(dòng)工況的動(dòng)態(tài)仿真[J]. 上海交通大學(xué)學(xué)報(bào), 2004, 38(2): 281-285.

[12] 趙丹, 丁國(guó)良, 胡海濤. 質(zhì)量和能量嚴(yán)格守恒的蒸發(fā)器動(dòng)態(tài)仿真模型[J]. 制冷學(xué)報(bào), 2015, 36(1): 76-83.

[13] GRUHLE W D, ISERMANN R. Modeling and control of a refrigerant evaporator[J]. Journal of Dynamic Systems Measurement and Control, 1985, 107(4): 287-292.

[14] CHI J, DIDION D. A simulation model of the transient performance of a heat pump[J]. International Journal of Refrigeration, 1982, 5(3): 176-184.

[15] 楊叢彥, 趙丹, 丁國(guó)良, 等. 多類型冰箱部件及系統(tǒng)仿真軟件開發(fā)[C]// 中國(guó)家用電器技術(shù)大會(huì)論文集. 北京:《電器》雜志社, 2015: 37-46.

[16] 趙丹. 面向制冷空調(diào)裝置快速穩(wěn)定仿真的部件模型拓展及系統(tǒng)仿真平臺(tái)開發(fā)[D]. 上海: 上海交通大學(xué), 2014.

Calculation Method for Transient Performance of Parallel Dual Evaporators for Multi-temperature Refrigerator

YANG Cong-yan*, ZHAO Dan, HU Hai-tao**, DING Guo-liang
(Institute of Refrigeration & Cryogenics Engineering, Shanghai Jiao Tong University, Shanghai 200240, China)

The adoption of a parallel dual evaporators for multi-temperature refrigerator can improve the accuracy of temperature control and the system energy efficiency. Because of the complicated operation pattern of the parallel dual evaporator, the existing evaporator model cannot apply to the parallel dual evaporator transient simulation directly. The purpose of this study is to develop a fast and stable algorithm for calculating the transient performance of parallel dual evaporators. All phase distributions of the refrigerant in both evaporators are analyzed by summarizing all the operation patterns and the forms of operation pattern switchover of the parallel dual evaporators, and the equations for calculating mass, internal energy and other state parameters of all refrigerant states and the criterion of smooth pattern switchover are established. Based on the above analysis, a calculating method for the transient performance in all operation patterns of parallel dual evaporators is developed. The calculation method is inserted into the simulation software, and the relative deviation between the simulation results and the experimental data is less than 10%.

Parallel dual evaporator; Refrigerator; Transient performance; Fast; Stable

10.3969/j.issn.2095-4468.2016.06.109

*楊叢彥（1990-），女，碩士生。研究方向：冰箱系統(tǒng)仿真及優(yōu)化設(shè)計(jì)。

**胡海濤（通訊作者），男，博士，副教授、博士生導(dǎo)師。研究方向：制冷系統(tǒng)仿真及換熱器優(yōu)化設(shè)計(jì)。聯(lián)系地址：上海市東川路800號(hào)上海交通大學(xué)，郵編：200240。聯(lián)系電話：021-34206295。E-mail：huhaitao2001@sjtu.edu.cn。