Nghiên ứu mô phỏng trong hế độ ổn định máy nén kín dùng ho bơm nhiệt loại nhỏ

Taylor, Bolton Chris 1997, Integration of compressor performance maps and NIST refrigerant database in an air conditioner thermal performance simulation model, 2, pp.. Caillat Jean-Luc,

Máy nén kín dùng cho b

Máy nén pittông

Xi lanh Đỉnh xi lanh Ống tiêu âm Đầu hút Động cơ điện

Hình 1.3 Sơ đồ cấu tạo của máy nén pittông kín

 d d  hút (p2/p1) Khi p2 càng cao thì Vd  2  d 

Hình 1.4 Biểu diễn quá trình nén thực của máy nén pittông trên đồ thị chỉ thị p – V

Ngoài ra, ka sát làm nóng pittông và xilanh, và xilanh  Nên khi 

hút vào xilanh, thành xilanh 

Máy nén scroll

Máy nén scroll          B scroll                nhiên,     nén scroll và m

 thì  [11] Ngày nay, máy nén scroll 

Hình 1.5 Các bộ phận chính của máy nén scroll

Khi máy nén scroll  làm thay

Hình 1.6 Các quá trình nén khí trong máy nén scroll

và êm ái trong quá trình máy nén 

1.2 Tổng quan v mô phề ỏng máy nén kín

1.2.2 Mô hình bán thực nghi m ệ

 Mô hình Hiller và Glicksman 

 , các  cho các môi c Mô hình c phép  

  Tuy nhiên, mô hình trên không 

 Mô hình này cho phép tính toán

 cho máy nén pittông và scroll Li 

 là máy nén pittông, máy nén scroll, máy nén rô x

 k  t o Theo [9] các hàm này có :

H

Ph   u l tt

Ph

Theo [23]  Threlkeld xem quá trình nén và 

Nhập giá trị to, tk, Qo, Vlt,

Tính toán thông số trạng thái tại các điểm của chu trình

Tính lưu lượng khối lƣợng thực tế mtt

Tính hiệu suất thể tích 

Hình 2.4 Sơ đồ tính toán hiệu suất thể tích 

hình 2.4,  2.1o, tk, t qn, tql  , riêng qo trình (2.3o 

Bảng 2.1 Thông số đặc trưng của các máy nén dùng để xác định dạng hàm hiệu suất thể tích và hiệu suất không thuận nghịch

 hình 2.4  bày trong   cho các máy nén khác có trong 

Bảng 2.2 Hiệu suất thể tích của máy nén P2 tk, o C to, o C

Bảng 2.3 Hiệu suất thể tích của máy nén S6 tk, o C to, o C

  av, bv và cv toán hình 2.4 ta

2.2.2 Hiệu su t không thu n nghấ ậ ịch

 tính toán hi u su t không thu n ngh ch c a máy nén        2.5

c hi u su t không thu n ngh ch c a máy nén ng v i m i         ch  làm vic to và tk  s 

  s cho các máy nén khác

Nhập giá trị t o , t k , Q o , V lt , t qn

Tính toán thông số trạng thái tại các điểm của chu trình

 riêng qo và công nén

Tính lưu lượng khối lƣợng thực tế m tt

Tính tổn thất cơ do ma sát N me

Tính công suất nén thực tế N tt

Tính công nén thực tế l r

Tính hiệu suất không thuận nghịch  s

Hình 2.5 Sơ đồ tính toán hiệu suất không thuận nghịch  s

Bảng 2.4 Hiệu suất không thuận nghịch của máy nén P2 tk, o C to, o C

Bảng 2.5 Hiệu suất không thuận nghịch của máy nén S6 tk, o C to, o C

 26) và (2.27)  s,i  1÷ 3) theo 3   máy nén pittông và máy nén scroll : là

 máy nén pittông và máy nén scroll  (2.1) xác

 o , (2.11)  máy nén Ne và (2.13)  tt

  cácmáy nén khác nhau 

Chương 3 XÂY D NG PH N M M MÔ PH NG MÁY NÉN Ự Ầ Ề Ỏ KÍN

3.1 Xây d ng mô hình mô ph ng máy nén kín ự ỏ

3.1.1 L a ch n các thông s u vào ự ọ ố đầ

3.1.2 L a ch n các thông s u ra ự ọ ố đầ

3.1.3 L a ch n các thông s ự ọ ốmô tả

Mô hình mô phỏng máy nén kín

Hình 3.1 Sơ đồ mô hình mô phỏng trong chế độ ổn định máy nén kín

3.2 Xây d ng thu t toán ự ậ chương trình mô phỏng máy nén kín

Tính công nén riêng ls

Vlt, tqn , t ql, me, el

Tính năng suất lạnh riêng q o

Hình 3.2 Lưu đồ thuật toán xác định thông số hoạt động của máy nén

 tiên ta,  o , tk các Ndd,

Vlt, t qn , t ql , el, me  ,

rình con trình bày trên hình 3.4  s e 

 po o  máy nén)  tính theo :

Hình 3.3 Lưu đồ thuật toán chương trình con xác định lưu lượng khối lượng

các đi: to,[i=1÷3], tk,[i=1÷3], Qo,[i=1÷3], Nel,[i=1÷3]

Tính các hệ số a v, bv, cv, as, bs, cs i:=1

3.3 L p trình mô ph ng máy nén ậ ỏ kín

3.3.1 L a ch n công c l p trình mô ph ng máy nén ự ọ ụ ậ ỏ kín

3.3.2 Giới thiệu ph n m m mô ph ng máy nén kín CLD ầ ề ỏ

Hình 3.5 Giao diện phần mềm CLD mô phỏng hoạt động của máy nén kín

Hình 3.6 Giao diện bảng thông số hoạt động của máy nén theo sự thay đổi t k và t o

heo tk và to (hình 3.6) 3.4 Kết luận chương 3

Chương 4 ĐÁNH GIÁ SAI S MÔ PH NG Ố Ỏ

4.1 L a ch n m u máy ự ọ ẫ nén để đánh giá sai số

máy nén scroll Copeland t qn = 10 K, t ql = 0 K và máy nén pittông Danfoss t qn = 11,1 K, t ql = 8,3 K

Bảng 4.1 Thông số đặc trưng của các máy nén

STT Hãng Model Môi chất

Công suất điện định mức, HP

y i , yi 4.2 L a chự ọn chế độ thực nghi m ệ

 này có  áy nén do nhà

 sao cho khi tính toán 

Hình 4.1 Các phương án lựa chọn 3 chế độ thực nghiệm

Bảng 4.2 Giá trị các hệ số của hàm  và  s với các máy nén khác nhau

Ký hiệu máy nén av bv cv as bs cs

4.3 Đánh giá sai số mô phỏng

 o  e   và so sánh Qo và Ne  

Bảng 4.3 So sánh Q o của mô hình mô phỏng với số liệu thực nghiệm (máy nén S7)

Hình 4.2 Đồ thị so sánh Q o của mô hình mô phỏng với số liệu thực nghiệm (máy nén S7)

Bảng 4.4 So sánh Q o của mô hình mô phỏng với số liệu thực nghiệm (máy nén P2)

Hình 4.3 Đồ thị so sánh Q o của mô hình mô phỏng với số liệu thực nghiệm (máy nén P2)

Hình 4.4 Đồ thị so sánh N e của mô hình mô phỏng với số liệu thực nghiệm (máy nén S7)

Bảng 4.6 So sánh N e của mô hình mô phỏng với số liệu thực nghiệm (máy nén P2)

Hình 4.5 Đồ thị so sánh N e của mô hình mô phỏng với số liệu thực nghiệm (máy nén P2)

Bảng 4.7 Tổng hợp sai số giữa kết quả mô phỏng Q o và N e với số liệu thực nghiệm

Sai số tương đối lớn nhất (về phía âm), %

, % lớn nhất (về phía dương)

Nhận xét: máy nén kín dùng cho

Chương 5 ỨNG D NG C A MÔ HÌNH Ụ Ủ MÔ PH NG MÁY NÉN KÍN Ỏ

5.1 Xây dựng đường đặc tích c a máy nén ủ kín

máy nén kín CLD  tính Qo và Ne máy nén kín  Do

Hình 5.1 Đường đặc tính năng suất lạnh Q o của máy nén S6

Hình 5.2 Đường đặc tính công suất điện tiêu thụ N e của máy nén S6

5.2 Ứng dụng mô hình mô phỏng máy nén kín để mô phỏng bơm nhiệt

5.2.1 Mô hình mô phỏngbơm nhiệt

5.2.1.1 Mô hình mô phng máy nén kín

Mô hình máy nén kín 

Mô hình mô phỏng máy nén kín

Thông số đầu vào Thông số đầu ra

Hình 5.3 Sơ đồ mô hình mô phỏng trong chế độ ổn định máy nén kín

5.2.1.2 Mô hình mô phng  a) Gi thi t c th   và n gin hóa mô hình

t bh  tính theo : kk,v kk,r bh kk,v o kk,r o t t t t t ln t t

Thông số đầu ra mkk

Hình 5.4 Sơ đồ mô hình mô phỏng dàn bay hơi

5.2.1.3 Mô hình mô phng thit b  a) Gi thi t c th   n hóa mô hình:

Hình 5.5 Quá trình truyền nhiệt trong thiết bị ngưng tụ a)  b) 

- 61 - nc,r nc,v nt nc,r k nc,v k t t t t t ln t t

Thông số đầu ra mnc

[W/m 2 K] i2r[kJ/kg] tK [ o C] i3 [kJ/kg] tnc,r [ o C]

Hình 5.6 Sơ đồ mô hình mô phỏng thiết bị ngưng tụ

5.2.1.4 Mô hình mô phng thit bti

DÀN BAY HƠI t kk,v [ o C] t o [ o C] m tt [kg/s] i 4 [kJ/kg]

Hình 5.7 Sơ đồ mô hình mô phỏng máy nén và dàn bay hơi

THIẾT BỊ NGƢNG TỤ tnc,v [ o C] i2r [kJ/kg] i3 [kJ/kg] i4[kJ/kg]

Qk[kW] mtt[kg/s] tk [ o C]

Hình 5.8 Sơ đồ mô hình mô phỏng thiết bị ngưng tụ và van tiết lưu

mô hình máy nén và dàn bay

DÀN BAY HƠI tnc,v [ o C] tkk,v [ o C] t o [ o C] t k [ o C] i 2r [kJ/kg] m tt [kg/s] m tt [kg/s] i 3 [kJ/kg] i 4 [kJ/kg]

Hình 5.9 Sơ đồ mô hình mô phỏng bơm nhiệt

5.2.2 Phân tích kết quả mô phỏng bơm nhiệt

Bảng 5.1 Thông số kỹ thuật các bộ phận cấu thành bơm nhiệt

Tên bộ phận Thông số kỹ thuật

Hình 5 .10 Giao diện phần mềm mô phỏng thông số làm việc của bơm nhiệt

Bảng 5.2 Một số thông số làm việc của bơm nhiệt mô phỏng ở điều kiện thiết kế

TT Tên thông số Ký hiệu Đơn vị Giá trị

Trong bài viết này, chúng ta sẽ phân tích các thành phần trong bảng 5.1 và cách xác định liệu máy nén có bị quá tải hay không Để làm điều này, chúng ta sử dụng các thông số làm việc và điều kiện nhiệt độ của không khí, bao gồm cả áp suất và nhiệt độ môi trường Việc hiểu rõ các thông số này là rất quan trọng để đảm bảo hiệu suất tối ưu cho máy nén, đồng thời giúp giảm thiểu tiêu thụ năng lượng Thông qua việc áp dụng các công thức và dữ liệu cụ thể, chúng ta có thể đánh giá chính xác tình trạng hoạt động của thiết bị.

K t qu mô ph  ng h p khi  t nhi t bé nh t và công su  t

n tiêu th c a máy nén l n nh   c trình bày trong các b ng 5.3 và b ng 5.4  

Kết quả phân tích cho thấy sự giảm nhiệt độ môi trường làm giảm hiệu suất hoạt động của máy nén N95E Cụ thể, nhiệt độ trung bình giảm từ 20,7°C xuống 14,6°C, trong khi nhiệt độ tối thiểu giảm từ 8°C xuống 3,8°C Điều này dẫn đến sự suy giảm công suất tiêu thụ của máy nén, từ 6,40 kW xuống 5,89 kW.

 t l nh co (t 15,66 kW xu ng 13,35 kW) Do c   công sun tiêu th c t lnh c  u suy gim dn s suy gi m c  t nhi t c  Q c 20,95 kW xu ng (t 18,14 kW)

Bảng 5.3 Kết quả mô phỏng trong trường hợp năng suất nhiệt đạt bé nhất

TT Tên thông số Ký hiệu Đơn vị Giá trị

Bảng 5.4 Kết quả mô phỏng trong trường hợp công suất điện tiêu thụ của máy nén lớn nhất

TT Tên thông số Ký hiệu Đơn vị Giá trị

Công suất tiêu thụ của máy nén Nét được xác định trong khoảng từ 6,4 kW đến 20,95 kW, với nhiệt độ môi trường tăng từ 20,7°C lên 39,5°C Dữ liệu này được trình bày trên bảng 5.3, cho thấy mối liên hệ giữa công suất tiêu thụ và nhiệt độ môi trường Bảng 5.4 cung cấp thông tin chi tiết về công suất tiêu thụ của máy nén trong các điều kiện khác nhau.

-   cho máy nén bán kín 

            Tạp chí Năng lượng Nhiệt, 114, pp 14-18

Tạp chí Khoa học Giao thông Vận tải

3 Bùi (2007), Tính toán thiết kế thiết bị trao đổi nhiệt 

4  , PTùy (2006), Kỹ thuật lạnh cơ sở  

5 Tạp chí Năng lượng Nhiệt, 61, pp 6- 8.

6 (2005), Nghiên cứu thực nghiệm bơm nhiệt đun nước nóng sử dụng dàn lạnh không khí trong điều kiện Việt nam

7 Switzerland: Heat pump market Available from: http://www.hp- summit.de/en/press/countryreports/country_report_switzerland/

8 http://www.elektronika-sa.com.pl/tclines.php?TID=1

9 Arthur J Howard, Beard J Taylor, Bolton Chris (1997), Integration of compressor performance maps and NIST refrigerant database in an air conditioner thermal performance simulation model, 2, pp 1265-1270

10 Caillat Jean-Luc, Ni Shimao (1988), A computer model for scroll compressors,

International compressor engineering conference, pp 601

11 Carrier (2004), Scroll compressors, Carrier corporation syracuse, New York, U.S.A

12 Carrington C G., Bannister P., Liu Q (1996), Performance of a scroll compressor with R134a at medium temperature heat pump conditions,

International Journal of Energy Research, 20, pp 733-743

13 Duprez Marie-Eve, Dumont Eric, Frère Marc (2007), Modelling of reciprocating and scroll compressors, International Journal of Refrigeration,

14 Emerson (2008), Copeland scroll, Emerson Climate Technologies Inc

15 Fisher S K., Rice C K (1983), The Oak Ridge Heat Pump Models: I A Steady-State Computer Design Model for Air- -to Air Heat Pumps.

16 Haberschill P., et al (1994), Hermetic compressor models determination of parameters from a minimum number of tests, Purdue International Compressor

17 Hanlon Paul C (2001), Compressor handbook, MCGraw-Hill, New York - U.S.A

18 Hiller C C., Glicksman L R (1976), Detailed Modeling and Computer Simulation of Reciprocating Refrigeration Compressors, International Compressor Engineering Conference, pp 162

19 Jọhnig Dagmar I., Reindl Douglas T., Klein Sanford A (2000), A semi- empirical method for representing domestic refrigerator/freezer compressor calorimeter test data, ASHRAE transaction

20 Klein S A., Reindl D T (1999), Develop Data Base for Determining Optimum Compressor Rating Points for Residential Refrigerator and Freezer Compressors, ASHRAE RP-870

21 Klein S A., Alvarado F.L (2000), EES - Engineering Equation Solver, Middletown

22 Lemmon Eric W., Huber Marcia L., McLinden Mark O (2007), NIST Reference Fluid Thermodynamic and Transport Properties - REFPROP 8.0

23 Li Wenhua (2012), Simplified steady-state modeling for hermetic compressors with focus on extrapolation, International Journal of Refrigeration, 35, 6, pp 1722-1733

24 Li Wenhua (2013), Simplified steady-state modeling for variable speed compressor, Applied Thermal Engineering, 50, 1, pp 318-326

25 Navarro E., et al (2007), A phenomenological model for analyzing reciprocating compressors, International Journal of Refrigeration, 30, 7, pp

26 Popovic Predrag, Shapiro Howard N (1995), A semi-empirical method for modeling a reciprocating compressor in refrigeration systems, ASHRAE Transaction, 101, 2, pp 367-382

27 Systat (2007), TableCurve 3D, Systat Software Inc

28 Xin R., Hatzikazakis P (2000), Reciprocating Compressor Performance Simulation, International compressor engineering conference, pp 1360

29 Yang Liang, et al (2009), Loss-efficiency model of single and variable-speed compressors using neural networks, International Journal of Refrigeration, 32,

MÁY NÉN 2.1 không   tk, o C to, o C

3.12 th so sánh Qo c mô hình mô ph ng v i s a   liu th c nghi m (máy nén S6)  

{CHUONG TRINH CON TINH HE SO CUA ETAS VA LAMDA}

Procedure Tinhtoanetasetav(R$, t_k[1 3], t_o[1 3], Q_o_dot[1 3], N_e[1 3], eta_me, eta_el, DELTAt_qn,DELTAt_ql, N_dd, V_disp, vlt, N_el:lamda[1 3], eta_s[1 3], p_o[1 3], p_k[1 3] ) i:=1

Repeat p_k[i]=PRESSURE(R$,T=T_k[i],x=0) p_o[i]=PRESSURE(R$,T=T_o[i],x=1) t_1[i]=T_o[i]+DELTAT_qn {[C]} p_1[i]=p_o[i] i_1[i]=ENTHALPY(R$,T=t_1[i],P=p_1[i]) s_1[i]=ENTROPY(R$,T=t_1[i],P=p_1[i]) v_1[i]=VOLUME(R$,T=t_1[i],P=p_1[i]) s_2[i]=s_1[i] p_2[i]=p_k[i] t_2[i]=TEMPERATURE(R$,s=s_2[i],P=p_2[i]) i_2[i]=ENTHALPY(R$,s=s_2[i],P=p_2[i]) t_3[i]=T_k[i]-DELTAt_ql i_3[i]=ENTHALPY(R$,T=t_3[i],x=0) i_4[i]=i_3[i] q_o[i]=i_1[i]-i_4[i] m_tt[i]=Q_o_dot[i]/q_o[i] v[i]=m_tt[i]*v_1[i] lamda[i]=v[i]/vlt

DELTA_me=(eta_el*N_el*(1-eta_me))

N_r[i]=N_e[i]*eta_el-DELTA_me l_r[i]=N_r[i]/m_tt[i] l_s[i]=i_2[i]-i_1[i] eta_s[i]=l_s[i]/l_r[i] i:=i+1 until (i>3) end

{inputs 3 diem mau} t_k[1]=lookup(link$,1,1) t_o[1]=lookup(link$,1,2)

N_e[2]=lookup(link$,2,4) t_k[3]=lookup(link$,3,1) t_o[3]=lookup(link$,3,2)

N_e[3]=lookup(link$,3,4) call Tinhtoanetasetav(R$, t_k[1 3], t_o[1 3], Q_o_dot[1 3], N_e[1 3], eta_me, eta_el, DELTAt_qn,DELTAt_ql, N_dd, V_disp, vlt, N_el:lamda[1 3], eta_s[1 3], p_o[1 3], p_k[1 3] )

{Nhap dang ham cua etas} eta_s[1]=a+b*(ln(p_k[1]/p_o[1]))^2+c/((p_k[1]/p_o[1])^1.5) eta_s[2]=a+b*(ln(p_k[2]/p_o[2]))^2+c/((p_k[2]/p_o[2])^1.5) eta_s[3]=a+b*(ln(p_k[3]/p_o[3]))^2+c/((p_k[3]/p_o[3])^1.5) eta_s=a+b*(ln(p_k/p_o))^2+c/((p_k/p_o)^1.5)

{Nhap dang ham cua lamda} lamda[1]=a_v+b_v*ln(t_k[1])+c_v*t_o[1] lamda[2]=a_v+b_v*ln(t_k[2])+c_v*t_o[2] lamda[3]=a_v+b_v*ln(t_k[3])+c_V*t_o[3] lamda=a_v+b_v*ln(t_k)+c_v*t_o

{CHUONG TRINH TINH NANG SUAT LANH VA CONG SUAT DIEN}

V_disp=lookup(link$,1,6){[m^3/h]} eta_me=lookup(link$,1,7) eta_el=lookup(link$,1,8)

{output} vlt=V_disp/3600 {[m^3/s]} p_k=PRESSURE(R$,t=t_k,x=0) p_o=PRESSURE(R$,t=t_o,x=1) t_1=t_o+DELTAT_qn {[C]} p_1=p_o i_1=ENTHALPY(R$,T=t_1,P=p_1) s_1=ENTROPY(R$,T=t_1,P=p_1) v_1=VOLUME(R$,T=t_1,P=p_1)

- 113 - s_2=s_1 p_2=p_k t_2=TEMPERATURE(R$,s=s_2,P=p_2) i_2=ENTHALPY(R$,T=t_2,P=p_2) t_3=T_k i_3=ENTHALPY(R$,T=t_3,x=0) {Chon do qua lanh =0} i_4=i_3 q_o=i_1-i_4 l_s=i_2-i_1 m_tt=lamda*vlt/v_1 {[kg/s]} l_r=l_s/eta_s

DELTA_me=(eta_el*N_el*(1-eta_me))

{CHUONG TRINH MO PHONG BOM NHIET}

{MO HINH MO PHONG MAY NEN}

Procedure Tinhtoanetasetav(R$, t_k[1 3], t_o[1 3], Q_o_dot[1 3], N_e[1 3], eta_me, eta_el, DELTAt_qn,DELTAt_ql, N_dd, V_lt, N_el:lamda[1 3], eta_s[1 3], p_o[1 3], p_k[1 3] ) i:=1

Repeat p_k[i]=PRESSURE(R$,T=T_k[i],x=0) p_o[i]=PRESSURE(R$,T=T_o[i],x=1) t_1[i]=T_o[i]+DELTAT_qn {[C]} p_1[i]=p_o[i] i_1[i]=ENTHALPY(R$,T=t_1[i],P=p_1[i]) s_1[i]=ENTROPY(R$,T=t_1[i],P=p_1[i]) v_1[i]=VOLUME(R$,T=t_1[i],P=p_1[i]) s_2[i]=s_1[i] p_2[i]=p_k[i] t_2[i]=TEMPERATURE(R$,s=s_2[i],P=p_2[i]) i_2[i]=ENTHALPY(R$,s=s_2[i],P=p_2[i]) t_3[i]=T_k[i]-DELTAt_ql i_3[i]=ENTHALPY(R$,T=t_3[i],x=0) i_4[i]=i_3[i] q_o[i]=i_1[i]-i_4[i] m_tt[i]=Q_o_dot[i]/q_o[i] v[i]=m_tt[i]*v_1[i] lamda[i]=v[i]/V_lt

DELTA_me=(eta_el*N_el*(1-eta_me))

N_r[i]=N_e[i]*eta_el-DELTA_me l_r[i]=N_r[i]/m_tt[i] l_s[i]=i_2[i]-i_1[i] eta_s[i]=l_s[i]/l_r[i] i:=i+1 until (i>3) end

{inputs 3 diem mau} t_k[1]=lookup(link$,1,1) t_o[1]=lookup(link$,1,2)

N_e[2]=lookup(link$,2,4) t_k[3]=lookup(link$,3,1) t_o[3]=lookup(link$,3,2)

N_e[3]=lookup(link$,3,4) call Tinhtoanetasetav(R$, t_k[1 3], t_o[1 3], Q_o_dot[1 3], N_e[1 3], eta_me, eta_el, DELTAt_qn,DELTAt_ql, N_dd, V_lt, N_el:lamda[1 3], eta_s[1 3], p_o[1 3], p_k[1 3] )

{Nhap dang ham cua etas} eta_s[1]=a+b*(p_k[1]/p_o[1])/(ln(p_k[1]/p_o[1]))+c*((p_o[1])^3) eta_s[2]=a+b*(p_k[2]/p_o[2])/(ln(p_k[2]/p_o[2]))+c*((p_o[2])^3) eta_s[3]=a+b*(p_k[3]/p_o[3])/(ln(p_k[3]/p_o[3]))+c*((p_o[3])^3) eta_s=a+b*(p_k/p_o)/(ln(p_k/p_o))+c*((p_o)^3)

{Nhap dang ham cua lamda} lamda[1]=a_v+b_v*(t_k[1])+c_v*t_o[1] lamda[2]=a_v+b_v*(t_k[2])+c_v*t_o[2] lamda[3]=a_v+b_v*(t_k[3])+c_V*t_o[3] lamda=a_v+b_v*(t_k)+c_v*t_o

{May nen Bitzer Model : 4CC-9.2Y}

V_disp=lookup(link$,1,6){[m^3/h]} eta_me=lookup(link$,1,7) eta_el=lookup(link$,1,8)

V_lt=V_disp/3600 {[m^3/s]} m_tt=V_lt*lamda/v_1 {[kg/s]} p_o=pressure(R$,t=t_o,x=1) t_1=t_o+DELTAT_qn {[C]} i_1=enthalpy(R$,p=p_o,t=t_1) v_1=VOLUME(R$,T=T_1,P=p_o)