Mobile QR Code QR CODE : Korean Journal of Air-Conditioning and Refrigeration Engineering
Korean Journal of Air-Conditioning and Refrigeration Engineering

Korean Journal of Air-Conditioning and Refrigeration Engineering

ISO Journal TitleKorean J. Air-Cond. Refrig. Eng.
  • Open Access, Monthly
Open Access Monthly
  • ISSN : 1229-6422 (Print)
  • ISSN : 2465-7611 (Online)

  1. 인천대학교 대학원 기계공학과 대학원생 (Graduate Student, Graduate School of Mechanical Engineering, Incheon National University, Incheon, 22012, Korea )
  2. 인천대학교 기계공학과 교수 (Professor, Dept. of Mechanical Engineering, Incheon National University, Incheon, 22012, Korea )



Scroll compressor(스크롤 압축기), Ejector(이젝터), Gas injection(가스 인젝션), Combined heat pump cycle(복합 히트 펌프 사이클)

기호설명

$c_{v}$: 유량 계수
$COP$: 성능계수
$h$: 엔탈피 [kJ/kg]
$k$: 폴리트로픽 지수
$M$: 질량 [kg]
$P$: 압력 [bar]
$Q$: 능력 [kW]
$s$: 엔트로피 [kJ/kg․K]
$t$: 시간 [s]
$T$: 온도 [°C]
$u$: 내부에너지 [kJ/kg]
$W_{c}$: 압축기 일 [kW]
$x$: 건도
$\rho$: 밀도 [kg/m3]
$\sigma$: 유량비
$\theta_{c}$: 크랭크 각도 [°]
$\phi$: 스크롤 랩 각도 [°]

하첨자

$BC$: 기본 히트펌프 사이클
$c$: 냉방, 응축기, 압축실
$e$: 증발기
$f$: 액상
$g$: 기상
$h$: 난방
$inj$: 인젝션

1. 서 론

친환경차인 순수 전기차의 주요 과제 중 하나는 주행거리 확장이다. 특히 냉난방 운전을 할 경우 주행거리는 더욱 급격히 감소하는데, 일반적으로 냉방 시 25%~40%, 난방의 경우에는 50% 정도까지 감소한다. 이 같은 냉난방 시스템 가동 시 주행거리가 급격히 감소하는 문제를 해결하기 위해 성능이 향상된 히트펌프 시스템을 전기차에 적용하는 연구가 주요 관련 기업 등을 중심으로 진행되고 있다.(1) 히트펌프 시스템의 성능을 향상시키기 위한 방안들 가운데 이젝터 적용과 가스 인젝션 기술 등이 있다.

이젝터 적용에 대해서 Huang et al.(2), Eames et al.(3)은 이젝터에 대해 수학적 모델링을 제시하였고, Hisham et al.(4)은 이론적인 수학적 모델과 이젝터의 실험적 데이터를 비교하여 이젝터 설계 및 성능평가를 제시 하였다. 또한 Nehdi et al.(5)는 수치해석 시뮬레이션을 통하여 이젝터의 기하학적인 파라미터가 시스템에 미치는 영향을 연구하였다. Sumeru et al.(6)은 이젝터를 적용한 냉동사이클을 수치 해석적 방법을 적용하여 주위온도의 변화에 따라 성능을 분석하였다. 그 결과 제안한 이젝터 냉동사이클의 COP가 4~14% 정도 증가한다고 보고하였다. 또한, Bilir et al.(7)은 R134a 냉매 사이클에서 이젝터 사용시 최대 22.34%의 COP 개선을 보여주었다. 그외 많은 연구에 따르면 이젝터가 히트펌프 시스템의 성능을 크게 향상 시킬 수 있음을 보여주었다.(8)

가스 인젝션 적용에 대해서는 Beeton and Pham(9)는 가스 인젝션을 스크롤 압축기에 적용한 사례를 통해 가스 인젝션 관련 핵심 기반 기술들을 열거하였고, Ma and Li(10)는 가스 인젝션 압력이 내부열교환기를 장착한 히트펌프 시스템에 미치는 영향을 조사하였으며, Wang et al.(11,12)는 스크롤 압축기에서 가스 인젝션이 내부 압축 과정과 압축기 성능에 주는 영향을 열역학적 모델에 기초하여 수치해석적으로 연구하였다. 또한 전기차의 저온 난방 성능 향상을 위해 가스 인젝션 스크롤이 장착된 히트 펌프에서 난방 성능과 COP를 최대화 할 수 있는 인젝션 포트 위치와 인젝션 압력에 대한 실험적 및 수치해석적 연구(13,14) 및 전기차에서 인젝션 홀의 모양이 난방 성능에 미치는 영향에 대한 실험적 연구도 수행되었다.(15)

가스 인젝션 사이클에서는 여전히 팽창밸브를 통과하는 손실이 발생하므로, 이를 회수하기 위한 방안으로 가스인젝션 사이클에 팽창밸브 대신 이젝터를 결합시키는 복합 사이클도 제안되고 실험적 및 이론적 해석 연구도 수행되었다.(16,17,18)

Pang et al.(16)은 이젝터 출구 유량을 스크롤 압축기 압축 중간부에 주입하여 Fig. 3과 같이 유사 2단 압축을 이루는 R22 복합 사이클을 제안하였고, 실험을 통해 가스 인젝션만 있는 사이클 대비 난방능력과 에너지 효율이 각각 10% 및 4% 향상되는 결과를 얻었다. Xu and Ma(17)도 동일한 구성의 복합 사이클에서 비슷한 실험 결과를 얻었다.

복합 사이클에 대한 해석적 연구로는 Fig. 3과 같은 구성을 갖는 사이클에서 Wang et al.(18)은 이젝터에 대한 모델링을 포함한 사이클 해석을 통해, 냉매 R22, R290, 그리고 R32를 사용하는 복합 사이클들의 성능을 비교 하였다. 이 해석적 연구에서는 인젝션 압력을 단순히 기하학적 평균압력($P_{m}=\sqrt{P_{s}P_{d}}$)과 같다고 가정하였고, 이때의 인젝션 유량은 단순히 이젝터 출구부에 설치된 플래시 탱크에서의 건도에 해당하는 만큼의 가스량에 해당한다고 간주하였다. 이럴 경우, 실제 압축기가 수용할 수 있는 인젝션 유량에 대한 부분은 고려되지 못하게 되고, 이로 인한 오차는 인젝션 압력이 낮을수록 더 커진다. 즉, 인젝션 압력(인젝터 후단에 위치한 플래시 탱크 압력)이 낮을수록 플래시 탱크 내부 건도는 증가하여 발생되는 가스량은 많게 되지만, 압축기 측에서는 인젝션 압력과 압축기 내부 압축실 압력과의 차이는 작아지므로 압축기가 수용할 수 있는 인젝션 유량은 점점 감소하게 된다.

따라서 본 연구에서는 이젝터와 가스 인젝션이 결합된 복합 사이클에서 특별히 압축기가 수용할 수 있는 가스 인젝션량에 대한 수치해석적 모델을 제시하고, 이를 복합 사이클 전체의 성능 해석에 적용함으로써 복합 히트 펌프 사이클 성능에 대한 더 정확한 해석 결과를 제시하고자 한다.

2. 사이클 성능 향상 방안

2.1 이젝터(EC)

Fig. 1(a)Fig. 1(b)은 팽창장치로써의 이젝터가 포함된 시스템을 도식적으로 보여준다.

이젝터 사이클(EC)의 경우 압축기 토출유량이 응축기(3)를 통과하여 이젝터의 구동노즐(Primary Nozzle)(t)로 유입되며, 고압의 액상 냉매가 구동노즐을 통해 팽창되어 낮은 압력 영역을 생성한다. 또한, 증발기(6)를 통과한 기상의 냉매는 구동노즐로 인해 팽창된 액상의 냉매와의 압력차로 인해 흡입노즐(Secondary Nozzle)(7)를 통하여 혼합구간(mixing chamber)(m)으로 유입된다. 각각의 노즐로 유입된 액상냉매와 기상냉매는 혼합구간(m)에서 혼합되어 디퓨저(Diffuser)(f) 구간을 지나 플래시 탱크(flash tank)로 유입된다. 이때 플래시 탱크에서 기액분리된 기상냉매는 압축기 흡입부(1)로, 액상냉매는 팽창밸브(5)를 지나 증발기(6)로 흐르는 구조를 가진다. 팽창장치로서 이젝터를 이용하므로 증발기의 압력보다 높은 압력의 냉매가 압축기로 유입되어 용적형 압축기의 경우 유량이 증가되는 효과가 있으며, 압축기의 압력비가 작아져 압축기일이 감소되는 효과가 있다. 또한 팽창과정을 등엔트로피 과정으로 유도함으로써 냉방능력의 증가 효과를 얻을 수 있다.

Fig. 1 Ejector cycle.
../../Resources/sarek/KJACR.2019.31.9.420/fig1.png

2.2 가스 인젝션(VIC)

Fig. 2 Vapor injection cycle.
../../Resources/sarek/KJACR.2019.31.9.420/fig2.png

Fig. 2(a), Fig. 2(b)는 가스인젝션 시스템을 도식적으로 보여준다.

가스인젝션 사이클(VIC)의 경우 압축기 토출 유량은 증발기에서 압축기로 흡입된 유량에 인젝션 유량을 더한 값으로서 응축기를 통과한 후 제 1팽창밸브를 통하여 플래시 탱크(flash tank)로 이동하게 된다. 이때 플래시 탱크의 기상 냉매는 인젝션 라인을 통해 압축기 내부의 인젝션 챔버로 주입된다.(5) 플래시 탱크의 액상냉매는 제2 팽창밸브를 통과하여 증발기로 유입되고,(7) 증발기를 통과한 후 압축기로 흡입된다.

중간압인 플래시 탱크의 기상냉매를 압축기 압축실로 인젝션 해주므로, 응축기를 지나는 냉매가 증가되며, 상대적으로 적은 엔탈피를 가진 인젝션 챔버(5)내의 냉매와 압축기 압축부(1a)의 냉매가 혼합되어 중간 엔탈피 (1b)를 가지게 되며, 이로 인해 압축기 토출부의 온도가 감소된다.

2.3 이젝터와 가스 인젝션이 복합된 사이클(EVIC)

Fig. 3 Ejector+Vapor injection cycle
../../Resources/sarek/KJACR.2019.31.9.420/fig3.png

Fig. 3(a), Fig. 3(b)는 이젝터와 가스인젝션이 결합된 복합시스템을 도식적으로 보여준다.

EVIC은 압축기, 응축기, 증발기, 이젝터, 두 개의 플래시 탱크, 두 개의 팽창밸브로 구성된다. 구동노즐(t)로 응축기(3)를 통과한 액상의 냉매가 유입되고 흡입노즐(8)로 제2 플래시 탱크(7)에서 기액분리된 기상의 냉매가 유입된다. 각각의 노즐로 유입된 냉매는 혼합구간(m)에서 혼합되어 디퓨저 구간(f)을 지나 제1 플래시 탱크로 유입된다. 제1 플래시 탱크에서 기액분리된 기상 냉매는 인젝션 챔버(5)로, 액상 냉매는 제1 팽창밸브(6)를 지나 제2 플래시 탱크(7)로 이동한다. 제 2 플래시 탱크(7)에서 분리된 기상 냉매는 흡입노즐(8)로 이동하고, 액상 냉매는 제2 팽창밸브(9)를 지나 증발기(10)로 유입되며, 증발기를 통과한 냉매는 압축기 흡입부(1)로 이동한다. EVIC는 이젝터 효과로 제 1 플래시 탱크의 건도를 보다 높여주어 압축기로 인젝션 할 수 있는 기상 냉매유량을 늘릴 수 있으며, 팽창과정을 등엔트로피 과정으로 유도함으로써 팽창 손실을 감소시키고 냉방 효과의 증가를 얻을 수 있다.

3. 성능 해석

본 논문에서는 앞서 설명한 이젝터 사이클과 가스 인젝션 사이클을 결합한 복합사이클(EVIC)에 대한 성능 해석을 진행하였다. 해석 대상 시스템은 R134a를 냉매로 사용하는 차량용 히트펌프 시스템이며, 압축기는 행정체적 27 cc의 전동식 스크롤 압축기이다.

3.1 스크롤 압축기 성능해석

Fig. 4는 가스 인젝션 홀의 위치가 나타나 있는 스크롤 랩 형상 조합을 보여준다. 스크롤 랩 기본 형상은 폴리노미얼(polynomial) 곡선이며(19), 스크롤 랩 주요 형상 치수는 Table 1에 나타내었다. 인젝션 탱크에서 압축기로 들어오는 인젝션 라인은 먼저 압축기 내부에 설치된 인젝션 챔버로 들어오고 이 인젝션 쳄버로부터 2개의 가스 인젝션 홀이 고정스크롤 경판면에 수직으로 뚫려져 압축실로 연통된다. 인젝션 홀들은 선회스크롤의 선회 운동에 따라 선회스크롤 랩 선단면에 의해 여닫히는 크랭크 각도 구간을 갖게 된다. Fig. 4(a)에서 선회스크롤(O/S) 랩 끝단에 점선으로 그려진 원은 선회 운동에 따라 선회스크롤 전체가 이러한 동일한 원을 그리는 선회운동하는 것을 나타낸다. 스크롤 랩 외곽 끝단이 닫혀서 밀폐 포켓이 형성되는 시점($\theta_{c}=360^{o}$)에서는 닫혀 있던 인젝션 홀이 95o 더 진행된 $\theta_{c}=455^{o}$(Fig. 4(b))에서 완전히 열리게 되고, 이러한 열림은 $\theta_{c}=659^{o}$(Fig. 4(c)) 까지 지속되다가 그 이후로는 닫히기 시작한다. 이러한 열림 면적 변화는 Fig. 6에 나타내었다.

압축기 압축부에서 압축 과정 계산은 먼저 식(1)~(2)의 에너지 방정식과 질량 보존의 법칙으로부터 압축실 가스의 내부 에너지($u_{c}$) 및 질량($M_{c}$)을 구하고, 식(3)에서 밀도($\rho_{c}$)를 구한다. 이렇게 구한 내부 에너지와 밀도로부터 식(4a)~(4c)와 같이 Refprop 9.0을 통하여 가스의 압력과 온도, 그리고 엔탈피를 구할 수 있다.

인젝션 홀에서 압축기 압축실로 인젝션 되는 유동은 압축성 유동으로 간주하여 식(5)과 같은 노즐 유동식을 적용하였다.

Fig. 4 Injection hole opening state (a) $\theta_{c}$ = 360˚, (b) $\theta_{c}$ = 455˚, (c) $\theta_{c}$ = 659˚.
../../Resources/sarek/KJACR.2019.31.9.420/fig4.png

Table 1. Scroll wrap parameter

Notation

Description

Value

Unit

$t_{\max}@\phi_{s}$

Maximum Thickness @ Inner Wrap starting angle

4.35 @ 90˚

[mm]

$t_{\min}@\phi_{e}$

Minimum Thickness @ Wrap end angle

2.72 @ 1110˚

[mm]

$h$

Wrap height

15.91

[mm]

$r_{s}$

Orbiting radius

4.75

[mm]

$V_{s}$

Volume

27

[cc]

$R_{1}$

Inner Circle radius1

5.8458

[mm]

$R_{2}$

Inner Circle radius2

1.1337

[mm]

(1)
$\dfrac{d}{dt}(u_{c}M_{c})=-P\dfrac{d V}{dt}+\sum\dot m_{in }h_{in }-\sum\dot m_{out}h_{c}+\dot m_{in j}h_{in j}$

(2)
$\dfrac{d}{dt}M_{c}=\sum\dot m_{in }-\sum\dot m_{out}+\dot m_{in j}$

(3)
$\rho_{c}=\dfrac{M_{c}}{V_{c}}$

(4a),(4b),(4c)
$P_{c}=f(u_{c},\:\rho_{c}),\:$ $T_{c}=f(u_{c},\:\rho_{c}),\:$ $h_{c}=f(u_{c},\:\rho_{c})$

(5)
$\dot m_{in j}=c_{v}P_{in j}A_{in j}\sqrt{\dfrac{2k}{(k-1)RT_{in j}}}\sqrt{Pr^{\dfrac{2}{k}}-Pr^{\dfrac{k+1}{k}}}$

Fig. 5는 가스 인젝션이 포함된 압축기의 성능해석 계산 절차를 보여준다. 압축기 형상 치수와 운전조건 으로부터 식(1)~식(4)에 따라 가스 인젝션이 반영된 압축실 압력을 계산하고, 이로부터 선회스크롤에 작용하는 가스력을 구한 뒤, 선회 스크롤 및 크랭크 축 등과 같은 운동 부재들에 작용하는 각종 반력과 마찰 손실 등을 구하여 최종적으로 압축기 소요동력, 냉매 토출량, 가스 인젝션 유량 등을 얻는다.(19)

Fig. 6에는 크랭크 각에 따른 인젝션 홀 개폐도와 인젝션 질량 유량을 보여준다. 인젝션 압력에 따라 인젝션 유량은 증가하는데, 인젝션 홀이 열려 있는 동안 인젝션 유량이 항상 양(+)의 값을 갖는 것은 아니다. 이는 압축실 압력($P_{c}$)이 인젝션 압력($P_{in j}$)보다 높아지는 영역에서는, 압축실 내부 가스가 인젝션 라인으로 역류되어 인젝션 유량이 음(-)의 값을 갖게 된다. 인젝션 압력이 높을수록 주입(+) 유량도 증가하고, 역류는 지연되어 더 많은 인젝션이 이루어진다. 인젝션 홀 직경이 커질수록 더 많은 인젝션 유량을 얻을 수 있지만 선회 스크롤 랩 두께를 넘어가게 되면 주위의 다른 압축실과의 연통이 일어나므로 홀 직경은 랩 두께 이하로 제한된다.

Fig. 5 Flow chart for scroll compressor performance analysis.
../../Resources/sarek/KJACR.2019.31.9.420/fig5.png

Fig. 6 Injection mass flow profiles.
../../Resources/sarek/KJACR.2019.31.9.420/fig6.png

3.2 EVIC 성능해석

Fig. 7 Ejector geometry.
../../Resources/sarek/KJACR.2019.31.9.420/fig7.png

이젝터에 대한 성능해석은 다양한 해석들이 존재한다. Xu et al.(18)는 혼합구간을 등압과정으로 가정 하였고, 식을 단순화하여 이젝터 해석을 진행하였다. 또한 Liu and Groll.(20)은 혼합구간 입출구의 압력을 고려한 정적과정 으로 이젝터 해석을 진행하였다. 본 연구에서는 혼합구간을 정적과정으로 가정하여 이젝터 해석을 진행하였다.

Fig. 7은 이젝터 각부의 명칭을 보여준다. 이젝터 성능해석을 위해서는 이젝터 각 부위에서 질량 보존 법칙과 에너지 보존 법칙, 그리고 운동량 방정식을 적용한다. 각 지점들의 주요 열역학적 값들은 Refprop 9.0을 이용하여 구하였다.

Fig. 3(a), Fig. 3(b)와 같은 가스인젝션과 이젝터가 복합된 EVIC에서 각점의 질량 보존 관계식은 식(6)~식(12)와 같다.

(6)
$\dot m_{1}=\dot m_{10}=\dot m_{9}$

(7)
$\dot m_{4}=\dot m_{3}+\dot m_{8}$

(8)
$\dot m_{5}=x_{f}\dot m_{f}$

(9)
$\dot m_{6}=(1-x_{f})\dot m_{f}$

(10)
$\dot m_{7}=\dot m_{6}$

(11)
$\dot m_{8}=x_{7}\dot m_{7}$

(12)
$\dot m_{9}=(1-x_{7})\dot m_{7}$

이젝터에서 흡입노즐(secondary nozzle)을 통하여 유입되는 냉매유량($\dot m_{8}$)과 구동노즐(primary nozzle)을 통하여 유입되는 냉매유량($\dot m_{3}$)의 비를 $\sigma$라고 하면 식(13)과 같이 정의된다.

(13)
$\sigma =\dfrac{\dot m_{8}}{\dot m_{3}}$

(14)
$\dot m_{3}\left(h_{3}+\dfrac{V_{3}^{2}}{2}\right)=\dot m_{t}\left(h_{t}+\dfrac{V_{t}^{2}}{2}\right)$

식(14)에서 $V_{3}$는 $V_{t}$에 비해 매우 작은 값을 가지므로, 무시 할 수 있다. 또한 $\dot m_{3}=\dot m_{t}$이므로 식(15)와 같이 정리 할 수 있다.

(15)
$V_{t}=\sqrt{2\left(h_{3}-h_{t}\right)}$

여기서 구동노즐 출구 엔탈피 $h_{t}$는 구동노즐 입구(3)에서 출구(t)까지를 등엔트로피 과정($s_{3a}=s_{3}$)으로 간주하여 구한 $h_{t,\:is}$에 효율을 적용하여 구할 수 있다.

(16)
$h_{t,\:is}=f(P_{t},\:h_{3})$

(17)
$h_{t}=h_{3}-\eta_{t}(h_{3}-h_{t,\:is})$

(18)
$\dot m_{t}=\rho_{t}A_{t}V_{t}$

흡입 노즐에서의 질량유량($\dot m_{8}$)은 식(19)에서 얻을 수 있다. 식(20)은 제2 플래시 탱크에서 이젝터 내부(b) 까지의 압력강하를 보여주며 이때, 식(21)와 같이 유입되는 유량의 유동은 등엔탈피 과정으로 가정하였다. 식(22)에서 이젝터 혼합구간으로 들어가는 흡입노즐에서 유입된 유량의 속도($V_{mt}$)를 얻을 수 있다.

(19)
$\dot m_{8}=\rho_{8}A_{8}V_{8}$

(20)
$P_{b}=P_{8}-\dfrac{\rho_{8}V_{8}^{2}}{2}$

(21)
$h_{b}=h_{8}$

(22)
$\dot m_{s}=\dot m_{b}=\rho_{b}A_{mt}V_{mt}$

혼합구간 입구(t)에서 출구(m)까지 질량보존 방정식과 모멘텀 방정식 그리고 에너지 방정식은 각각 식(23)~ 식(25)와 같다.

(23)
$\rho_{t}A_{t}V_{t}+\rho_{b}A_{b}V_{b}=\rho_{m}A_{m}V_{m}$

(24)
$P_{t}A_{t}+\eta_{m}\rho_{t}A_{t}V_{t}^{2}+P_{b}(A_{m}-A_{t})+\eta_{m}\rho_{b}A_{b}V_{b}V_{mt}=P_{m}A_{m}+\rho_{m}A_{m}V_{m}^{2}$

(25)
$\dot m_{t}\left(h_{t}+\dfrac{V_{t}^{2}}{2}\right)+\dot m_{b}\left(h_{b}+\dfrac{V_{b}^{2}}{2}\right)=\dot m_{m}\left(h_{m}+\dfrac{V_{m}^{2}}{2}\right)$

혼합구간에서 혼합된 작동유체가 디퓨저를 지나는 동안 유체의 운동에너지는 정압으로 변환되어 압력이 상승하게 된다. 디퓨저에 의한 압력상승은 식(26)을 통하여 구할 수 있으며,(21) 식(27)을 통하여 디퓨저 출구(f)의 엔탈피를 구할 수 있다. 디퓨저 출구(f)의 엔탈피와 압력을 이용하여 건도를 구할 수 있으며, 디퓨저에서의 건도를 통해 인젝션유량 또한 구할 수 있다.

(26)
$\dfrac{P_{d}-P_{m}}{\dfrac{1}{2}\rho_{m}V_{m}^{2}}=\eta_{d}\rho_{m}\left[1-\left(\dfrac{A_{m}}{A_{d}}\right)^{2}\right]\left[\dfrac{x_{m}^{2}}{\rho_{g,\:m}}+\dfrac{\left(1-x_{m}\right)^{2}}{\rho_{f,\:m}}\right]$

(27)
$\dot m_{3}h_{3}+\dot m_{8}h_{8}=(\dot m_{3}+\dot m_{8})h_{d}$

한편 식(6)~식(12)의 질량 보존 관계식들로부터 식(28)과 같은 관계식을 구할 수 있다.

(28)
$x_{7}(1-x_{f})=\dfrac{\sigma}{1+\sigma}$

각종 식들에서 사용된 이젝터 효율들은 각각 $\eta_{t}=0.9$, $\eta_{b}=0.9$, $\eta_{m}=0.85$, $\eta_{f}=0.85$이 사용되었으며,(22) 이상의 관계식들로부터 냉매유량비($\sigma$), 가스 인젝션 압력($P_{in j}$), 가스 인젝션 유량($\dot m_{in j}$) 등을 다음의 Fig. 8과 같은 계산 과정을 통해 구할 수 있다.

Fig. 8 Flow chart for EVIC performance analysis.
../../Resources/sarek/KJACR.2019.31.9.420/fig8.png

앞서 진행한 계산에 의하여, EVIC 사이클은 제 2 플래시탱크 압력($P_{8}$)을 제어함으로써 가스 인젝션 압력 ($P_{in j}$)을 제어할 수 있다. 가스 인젝션 압력($P_{in j}$)은 실제 압축기가 수용할 수 있는 인젝션 유량에 대한 부분을 고려하여 선정하여야 한다. Fig. 9는 인젝션 압력에 따른 인젝션 유량 변화($\dot m_{"in "j}$)를 보여준다. 인젝션 압력이 높을수록 압축기에서 받을 수 있는 인젝션 유량은 증가한다. 하지만 플래시 탱크에서 발생시킬 수 있는 가스 유량($\dot m_{flash}$)은 플래시 탱크의 압력(즉, 인젝션 압력)과 건도($x_{v}$)에 의해 결정된다. 따라서 인젝션 홀을 통해 압축실 내로 주입될 수 있는 유량의 최대치는 플래시 탱크에서의 가스 발생량에 의해 제한받게 되므로 두 곡선이 만나는 지점에서 최적 인젝션 압력($P_{in j}$)이 정해진다.

참고문헌(18)에서는 복합 사이클에서 인젝션 압력을 저압과 고압의 기하학적 평균치 $P_{in j}=\sqrt{P_{s}P_{d}}$로 가정했는데, 이는 사이클에서 제공할 수 있는 가스량과 압축기에서 수용할 수 있는 가스량의 균형을 고려하지 않은 가정이다. Fig. 9에서 가스 인젝션 압력을 기하학적 평균치 $P_{in j}=\sqrt{P_{s}P_{d}}$로 잡을 경우, 사이클에서의 가스 공급량이 부족하게 되어 실제로는 가스 인젝션이 제대로 일어날 수 없음을 보여준다. 본 해석에서는 복합사이클(EVIC)에서의 인젝션 압력이 기하학적 평균치 인젝션 압력보다는 낮게 형성된다. 또한 복합사이클(EVIC) 에서의 인젝션 압력(E점)이 인젝션만 있는 사이클(VIC)에서의 인젝션 압력(B점)보다 더 높게 형성되는데, 이는 복합사이클에서 인젝션 시점(Fig. 3의 f 점)의 건도($x_{f}$)가 순수 인젝션 사이클에서의 인젝션 시점(Fig. 2의 4점)의 건도($x_{4}$)보다 높기 때문이다.

Fig. 9 Injection mass flow rates vs. Injection pressure.
../../Resources/sarek/KJACR.2019.31.9.420/fig9.png

4. 성능해석 결과 및 고찰

4.1 이젝터 형상의 영향

이젝터의 형상은 이젝터 성능의 중요한 변수가 된다. 이젝터의 성능은 혼합 구간의 면적($A_{m}$)과 구동노즐의 면적($A_{t}$)의 비 $A_{m}/A_{t}$를 설계 변수로 하여 최적 성능을 보이는 설계점을 얻을 수 있다.(5) 그런데 순수 이젝터 사이클과 복합 사이클에서의 이젝터 최적 설계점은 다르게 나타난다.

Fig. 10은 이젝터 사이클(EC)과 복합 사이클(EVIC)에서 면적비($A_{m}/A_{t}$)에 따른 COP 변화를 보여준다. 이젝터 사이클(EC)에서는 $A_{m}/A_{t}=7.3\sim 7.9$에서, 복합사이클(EVIC)에서는 $A_{m}/A_{t}=1.7\sim 2.0$ 부근에서 최고 성능을 보인다. 이러한 차이는 기본적으로 각 사이클에서 이젝터 혼합구간 압력 $P_{m}$ 수준이 많은 차이가 나기 때문 이다(Fig. 11 참조). 혼합 냉매가 이젝터 혼합 구간을 지나 디퓨저 출구에 이르면서 압력 상승이 일어나 이젝터 출구 압력이 올라가게 되는데, 디퓨저에서의 압력 상승효과는 면적비 $A_{m}/A_{t}$가 커지면서 점점 작아지게 된다. 이는 혼합구간에서의 단면적 증가가 디퓨저 출구 및 입구 면적비($A_{f}/A_{m}$)의 감소로 이어지기 때문이다. 이젝터 혼합구간 압력 $P_{m}$ 변화에 디퓨저 압력 상승을 더한 값이 이젝터 출구 압력 $P_{f}$을 결정하며, 이젝터 출구 압력이 높을수록 높은 성능을 나타낸다.

Table 2에는 이젝터 사이클 및 복합 사이클에서 각각 최적 성능을 내는 면적비 영역에서 설계한 이젝터 주요 치수를 보여준다. 복합사이클에서는 이젝터 출구에서의 인젝션 유량이 압축기에서 받아들일 수 있는 인젝션 유량과 일치해야 하는데, 이는 제2 플래시 탱크 압력 $P_{8}$을 조절하여 이루어 줄 수 있다.

Fig. 10 Effects of area ratio of $A_{m}/A_{t}$ on COP.
../../Resources/sarek/KJACR.2019.31.9.420/fig10.png

Fig. 11 Effects of area ratio of $A_{m}/A_{t}$ on $P_{m}$ and $P_{f}$.
../../Resources/sarek/KJACR.2019.31.9.420/fig11.png

Table 2. Ejector Geometry

Notation

Description

EC

EVIC

Unit

$D_{t}$

Diameter of Primary Nozzle

2

2

[mm]

$D_{b}$

Diameter of Secondary Nozzle

10

10

[mm]

$D_{m}$

Diameter of Mixing Chamber

5.6

2.8

[mm]

$D_{d}$

Diameter of Diffuser

20

10

[mm]

$A_{m}/A_{t}$

Area ratio

7.84

2

-

4.2 기본 운전 조건 해석 결과

Table 3에 나타낸 기본 운전조건에서 기본 히트펌프 사이클(BC), 이젝터 사이클(EC), 인젝션 사이클(VIC), 그리고 복합사이클(EVIC) 각각에 대해 성능해석을 진행하였고, 그 결과를 Table 4에 나타내었다. 난방 능력의 경우, 기본 히트펌프 사이클(BC) 대비, 이젝터 사이클(EC)과 가스 인젝션 사이클(VIC), 그리고 복합사이클 (EVIC)은 각각 7.72%, 34.3%, 40.8%가 증가하였고, 난방성적계수 COPh는 각각 4.18%, 10.3%, 12.6% 증가 하였다.

Table 3. Operating Conditions

Notation

Description

Unit

Value

BC

EC

VIC

EVIC

$P_{s}$

Evaporator pressure/Evaporator Temperature

[bar]/[℃]

2/-10

2/-10

2/-10

2/-10

$P_{d}$

Condenser pressure/Condenser Temperature

[bar]/[℃]

18/62.9

18/62.9

18/62.9

18/62.9

$T_{sh}$

Super heating Temp

[℃]

5

0

5

5

$T_{sc}$

Sub cooling Temp

[℃]

5

5

5

5

이젝터 사이클(EC)에서 압축기 입구 압력이 2.42 bar로 높아진 것은 이젝터 내부에서 압력이 상승된 냉매가 압축기 입구로 흡입되기 때문이다. 압축기 입력의 경우, 이젝터 사이클(EC)에서 압축기 흡입 압력이 높아져서 압축비가 작아졌음에도 불구하고 압축기 입력이 증가한 것은 흡입 압력에 높아짐에 따른 흡입밀도 증가로 인해 기본 사이클 대비 압축기 토출량($\dot m_{c}$)이 늘어났기 때문이다.

EVIC에서 VIC보다 인젝션 압력이 높아진 것은 Fig. 9에서 설명한 것처럼 사이클 내 인젝션 지점에서의 건도가 EVIC에서 더 많기 때문에 이로 인해 증가된 인젝션 유량을 수용하기 위한 것이다.

VIC 대비 EVIC의 압축기 토출량은 5.5% 증가하였으나, 난방능력은 4.8% 증가하였다. 이는 EVIC에서 더 많은 인젝션 유량으로 인해 압축기 출구 엔탈피가 VIC 경우보다 낮아져서 응축기 입출구 엔탈피 차가 줄어들었기 때문이다. 압축기 입력은 EVIC에서 VIC 대비 2.95% 증가하여 결국 COPh는 1.62% 증가하게 된다.

Table 4. Performance comparison of various cycles

Notation

Description

Unit

Value

BC

EC

VIC

EVIC

$P_{c}$

Compressor inlet pressure

[bar]

2

2.42

2

2

$P_{in j}$

vapor Injection pressure

[bar]

-

-

5.25

5.64

$\dot m_{c}$

Compressor out mass flow

[kg/s]

0.01204

0.01484

0.01778

0.01876

$\dot m_{in j}$

Vapor injection mass flow

[kg/s]

-

-

0.00582

0.00680

$W_{c}$

Compressor input work

[kW]

1.029

1.064

1.254

1.291

$Q_{h}$

Heating capacity

[kW]

2.304

2.482

3.096

3.245

$COP_{h}$

Heating Coefficient of performance

-

2.24

2.33

2.47

2.51

4.3 증발 온도 변화

Fig. 12 Capacity increase in EC, VIC, and EVIC for various conditions.
../../Resources/sarek/KJACR.2019.31.9.420/fig12.png

Fig. 12는 응축기 온도를 $T_{c}$ = 62.89℃($P_{d}=18$ bar) 에 고정하고, 난방운전의 경우 증발기 온도를 $T_{e}$= -25℃~ 0℃($P_{s}$ = 1 bar~2.93 bar)로, 냉방운전의 경우 $T_{e}$ = 5℃~15℃($P_{s}$ = 3.5 bar~4.88 bar)로 변경하였을 때, 기본사이클 (BC) 대비 난방 및 냉방 능력의 변화를 보여주며, Fig. 13은 COP의 변화를 보여준다. 여기서 $\triangle Q$는 $Q-Q_{BC}$를 나타내며, $\triangle COP$는 $COP-COP_{BC}$를 나타낸다. 증발 온도가 높아질수록 EVIC, VIC, 그리고 EC 모두에서 난방 및 냉방 능력의 향상분과 성적계수 향상분이 감소하는 경향을 보인다. EVIC에서는 기본 사이클 대비 난방 능력은 52.8~35.1%P, 냉방능력은 53.4~39.9%P 향상되는 성능을 나타내며, VIC 대비해서는 난방능력은 4~6P% 정도, 냉방능력은 15~11P% 정도 높은 성능을 보인다. COP의 경우, 기본 사이클 대비 EVIC에서 난방에서는 0.33~ 0.22, 냉방에서는 0.52~0.44 높은 값을 보이며, VIC 대비해서는 난방에서는 0.044~0.023, 냉방에서는 0.22~0.21 높게 얻어진다. EVIC와 VIC의 비교에서 난방보다는 냉방시 성능에 더 큰 차이를 보이는 것은 난방 시에는 응축기를 통과하는 질량유량의 차이가 성능 차이의 주요 원인이라면, 냉방에서는 증발기를 통과하는 유량은 차이가 거의 없으나 증발기 입구 엔탈피가 EVIC에서 2번에 걸친 감압 과정을 통해 더 낮아지는 것이 주 원인이 되기 때문이다. 증발기를 통과하는 엔탈피 차이가 만드는 냉방능력의 차이가 질량유량의 차이가 만드는 난방능력의 차이보다 더 크기 때문에 EVIC와 VIC의 사이의 냉방 능력 차이가 난방능력 차이보다 상대적으로 더 커지게 된다.

Fig. 13 COP increase in EC, VIC, and EVIC for various conditions.
../../Resources/sarek/KJACR.2019.31.9.420/fig13.png

5. 결 론

차량용 R134a 히트 펌프 사이클에서 성능 향상을 위해 적용되는 이젝터, 가스 인젝션, 그리고 이들이 복합된 사이클에 대한 성능 해석을 수행하여 다음과 같은 결론을 얻었다.

(1) 가스 인젝션 사이클(VIC)과 복합 사이클(EVIC)에서 인젝션 적정 압력은 압축기가 수용할 수 있는 인젝션 유량과 플래시 탱크에서 생성하는 가스량이 같아지는 지점에서 인젝션 압력이 결정되며,

(2) 난방 기본 운전 조건($P_{s}/P_{d}$ = 2/18 bar)에서 복합 사이클의 인젝션 최대 압력은 $P_{in j}$ = 5.64 bar이고, 인젝션 사이클에서의 인젝션 최대 압력은 $P_{in j}$ = 5.25 bar로서 복합 사이클에서의 인젝션 압력이 더 높게 나타난다. 이는 이젝터 효과로 인해 복합 사이클에서 제 1 플래시 탱크의 건도가 증가하기 때문이다.

(3) 이젝터의 성능은 이젝터 혼합구간 면적($A_{m}$)과 구동노즐의 면적($A_{t}$)의 비에 따라 달라지는데, 이젝터만 있는 사이클(EC)에서는 면적비 7.3~7.9에서 그리고 복합 사이클(EVIC)에서는 면적비가 1.7~2.0에서 최고의 성능을 나타낸다.

(4) 난방 기본 운전 조건에서 이젝터 사이클(EC)과 가스 인젝션 사이클(VIC), 복합 사이클(EVIC)의 난방능력은 기본 히트 펌프 사이클 대비 각각 7.72%, 34.3%, 40.8% 증가하였으며, COPh는 각각 4.18%, 10.3%, 12.6% 증가하였다.

(5) 응축 온도를 62.9℃(18 bar)에 고정하고, 증발 온도를 -25℃~15℃로 변화시켰을 때, 증발 온도가 높아질수록 EC, VIC, 그리고 EVIC 모두에서 난방 및 냉방 능력의 향상 정도와 성적계수 향상 정도가 감소하는 경향을 보인다. 기본 사이클 대비 EVIC에서 난방능력은 52.8~35.1%P, 냉방능력은 53.4~39.9%P 향상된 값을 보이며, VIC 대비해서는 난방능력은 4~6P% 정도, 냉방능력은 15~11P% 정도 높게 나타난다. COP의 경우, 기본 사이클 대비 EVIC에서 난방에서는 0.33~0.22, 냉방에서는 0.52~0.44 높은 값을 보이며, VIC 대비해서는 난방에서는 0.044~0.023, 냉방에서는 0.22~0.21 높게 얻어진다.

References

1 
2017.12., Automotive reportGoogle Search
2 
Huang B. J., Chang J. M., Wang C. P., Petrenko V. A., 1998, A 1-D analysis of ejector performance, International Journal of Refrigeration, Vol. 22, pp. 354-364DOI
3 
Eames I. W., Aphornratana S., Haider H., 1995, A theoretical and experimental study of a small-scale steam jet refrigerator, International Journal of Refrigeration, Vol. 18, No. 6, pp. 378-386DOI
4 
Hisham E. D., Hisham E., Imad A., Ghada A. N., 2002, Evaluation of steam jet ejectors, Chemical Engineering and Processing, Vol. 41, pp. 551-561DOI
5 
Nehdi E., Kairouani L., Bouzaina M., 2006, Performance analysis of the vapour compression cycle using ejector as an expander, International Journal of Energy Research, Vol. 31, pp. 364-375DOI
6 
Sumeru K., Nasution H., Ani F. N., 2012, A review on two-phase ejector as an expansion device in vapor compression refrigeration cycle, Renewable Sustainable Energy Reviews, Vol. 16, pp. 4927-4937DOI
7 
Bilir N., Ersoy H. K., 2009, Performance improvement of the vapour compression refrigeration cycle by a two-phase constant area ejector, International Journal of Energy Research, Vol. 33, No. 5, pp. 469-480DOI
8 
Sarkar J., 2012, Ejector enhanced vapor compression refrigeration and heat pump systems-A review, Renewable Sustainable Energy Reviews, Vol. 16, pp. 6647-6659DOI
9 
Beeton W. L., Pham H. M., 2003, Vapor-injected scroll compressor, ASHRAE Journal, Vol. 45, pp. 22-27Google Search
10 
Ma G. Li X., 2007, Exergetic optimization of a key design parameter in heat pump systems with economizer coupled with scroll compressor, Energy Conversion and Management, Vol. 48, pp. 1150-1159DOI
11 
Wang B., Shi W., Li X., Yan Q., 2008, Numerical research on the scroll compressor with refrigeration injection, Applied Thermal Engineering, Vol. 28, No. 5, pp. 440-449DOI
12 
Wang B., Shi W., Li X., 2009, Numerical analysis on the effects of refrigerant injection on the scroll compressor, Applied Thermal Engineering, Vol. 29, No. 1, pp. 37-46DOI
13 
Kwon C., Kim M. S., Choi Y., Kim M. S., 2017, Performance evaluation of a vapor injection heat pump system for electric vehicles, Int. J. Refrig., Vol. 74, pp. 138-150Google Search
14 
Choi Y. U., Kim M. S., Kim G. T., Kim M., and Kim M. S., 2017, Performance analysis of vaporinjection heat pump system for electric vehicle in cold startup condition, Int. J. Refrig, Vol. 80, pp. 24-36DOI
15 
Qin F., Zhang G. Y., Xue Q. F., Zou H. M., Tian C. Q., 2017, Experimental investigation and theoretical analysis of heat pump systems with two different injection portholes compressors for electric vehicles, Applied Energy, Vol. 185, pp. 2085-2093DOI
16 
Pang Z. Z., Ma G. Y., and Yu L. H., 2007, Experimental Study on Quasi Two-stage Compression Heat Pump System Coupled with Ejector., 22nd International Congress of Refrigeration, Beijing.Google Search
17 
Xu S. X. and Ma G. Y., 2011, Research on air-source heat pump coupled with economized vapor injection scroll compressor and ejector, Int. J. Refrig, Vol. 34, No. , pp. 1587-1595DOI
18 
Wang X., Yu J., Xing M., 2015, Performance analysis of a new ejector enhanced vapor injection heat pump cycle, Energy Conversion and Management, Vol. 100, pp. 242-248DOI
19 
Lim J. T., Choi J. W., Kim H. J., 2019, Performance analysis of gas injection in an electric scroll compressor, Korean Journal of Air-Conditioning and Refrigeration Engineering, Vol. 30, No. 3, pp. 107-115Google Search
20 
Fang L., Eckhard A. G., 2013, Study of ejector efficiencies in refrigeration cycles, Applied Thermal Engineering, Vol. 52, pp. 360-370DOI
21 
Owen I., Abdul-Ghani A., Amini A. M., 1992, Diffusing a homogenized two-phase flow, Int. J. Multiphase Flow, Vol. 18, pp. 531-540DOI
22 
Ma G. Y., Zhao H. X., 2008, Experimental study of a heat pump system with flash-tank coupled with scroll compressor, Energy and Buildings, Vol. 40, pp. 697-701DOI