The Korean Society Of Automotive Engineers
[ Article ]
Transactions of the Korean Society of Automotive Engineers - Vol. 29, No. 7, pp.667-673
ISSN: 1225-6382 (Print) 2234-0149 (Online)
Print publication date 01 Jul 2021
Received 29 Jan 2021 Revised 01 Apr 2021 Accepted 02 Apr 2021
DOI: https://doi.org/10.7467/KSAE.2021.29.7.667

SCR 환원제 통합모듈의 재질변경에 따른 동결 요소수 해동에 관한 수치해석

박준영1) ; 김만영*, 1) ; 이원희2)
1)전북대학교 항공우주공학과
2)세종공업 DEF 개발팀
A Numerical Study on the Urea Thawing with Different Material in a Reducing Agent-Integrated Module
Junyeong Park1) ; Man Young Kim*, 1) ; Wonhee Lee2)
1)Department of Aerospace Engineering, Jeonbuk National University, Jeonbuk 54896, Korea
2)DEF R&D Team, Sejong Ind. Co., LTD, 120 Heungdeokjungang-ro, Giheung-gu, Yongin-si, Gyeonggi 16950, Korea

Correspondence to: *E-mail: manykim@jbnu.ac.kr

Copyright Ⓒ 2021 KSAE / 188-08
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License(http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium provided the original work is properly cited.

Abstract

The stringent emission regulation for diesel engines results in changes on emissions measurement methods(NEDC→WLTP), as per the Euro-6 standard. To satisfy such a strict regulation standard, Urea-SCR devices are most frequently used to ensure NOx reduction. Under a certain temperature(-11 °C), a urea-water solution, which is a reducing agent, becomes frozen and is present in the form of a solid state inside the urea water tank. This makes it impossible to operate the gear pump and to cause problems in the entire Urea-SCR system, including high emissions during cold starts. Thus, to secure a stable supply within the urea-water solution system, heating modules must be introduced to undertake a thawing phenomenon. In this study, a three-dimensional, unsteady numerical analysis of the process of thawing the frozen urea was performed by using the ANSYS Fluent. The urea flow path was investigated during the thawing phenomenon with respect to the heat transfer characteristics in the urea tank. Finally, the thawing performance of the two different models was compared in terms of thawing periods with the flow path and the frozen urea area inside the pump by changing the material and its property.

Keywords:

Urea-SCR, Phase change, Urea path, Gear pump, Heat transfer, PTC heater

키워드:

요소 첨가 선택적 촉매 반응 시스템, 상변화, 요소수 유로, 기어펌프, 열전달, PTC 히터

1. 서 론

내연기관에서 디젤엔진의 경우, 높은 열효율과 좋은 연비의 장점을 갖지만 가솔린 대비 NOx 및 PM과 같은 유해 물질을 많이 배출하는 것이 단점으로 작용하고 있다.1) 따라서 디젤엔진은 Euro-6과 같은 심화된 배기가스 규제와 배출가스 측정 방식의 변경(NEDC→WLTP)을 통해 배기가스 배출을 엄격하게 제한하고 있다.2,3) 이와 같은 규제에 대응하기 위해 EGR과 같은 엔진내부 저감방법과 LNT, Urea-SCR과 같은 후처리 장치에 대한 연구가 진행되고 있다.4) 특히 Urea-SCR은 상대적으로 높은 NOx 저감효율과 넓은 반응온도의 범위를 가지는 장점으로 인해 대형 및 소형 디젤 차량에 장착된다. 하지만 환원제로 사용되는 요소수는 물과 요소수의 혼합물로서 32.5 wt%의 농도를 가지며, 특정온도 (-11 °C)이하에서 동결되는 문제점을 가지고 있다. 이러한 특성으로 인해 낮은 온도에서 요소수 탱크, 유로, 기어 펌프 내부의 요소수가 동결되거나 얼음알갱이 형태로 남아있는 요소수는 기어펌프가 작동할 수 없게 만들며 인젝터로의 요소수 공급이 불가능하게 한다. 따라서 요소수를 신속하고 안정적으로 공급하기 위해서는 환원제 공급모듈 내 해동시스템의 적용과 열원으로부터 유로 및 기어펌프로의 효율적인 열전달에 대한 연구가 요구된다.

요소수 탱크내 해동현상에 대한 상변화 연구는 Voller 등5)이 제안한 엔탈피 방법이 많이 사용 되고 있다. 또한 권길식 등6)은 앞서 제시한 엔탈피 방법을 이용하여 상변화 과정에 대한 수치해석을 진행하고 이를 온도 수식화 방법과 비교하였다. 한편, Oliveira와 Issa7)는 수정된 PISO 알고리즘을 통해 2차원 공동 내 부력에 의한 유동현상을 관찰하고, 이를 SIMPLEC 알고리즘과 비교하였다. 한편, 이승엽 등8)은 전산해석을 활용하여 냉각수 순환가열 방식 모델에 대한 요소수 해동 및 열전달 특성을 고찰하여 초기 전도에 의한 열전달에 의해 요소수가 확보된 후에는 자연대류에 의한 해동현상이 진행되는 것을 확인하였다. 안태현 등9)은 냉각수 및 전기 가열 방식에 따른 요소수 해동 현상을 수치해석을 통해 비교하였으며, 전기가열 방식이 냉각수 순환 가열 방식보다 더 많은 요소수를 해동시키는 것을 확인하였다. 또한 김태욱 등10)은 히트 스프레더의 형상변경에 따른 요소수 탱크 내 해동현상을 비교하여, 열을 집중하는 형상과 대류현상을 빠르게 유도할 수 있는 형상의 중요성을 제시하였다. 이두기 등11)은 소형 상용차에 사용되는 전기 가열 모듈에 의한 요소수 해동현상을 고찰하여 열원 및 흡입구 위치와 같은 최적화 설계를 위한 기초 데이터를 확보하였다.

본 연구에서는 이전 연구들을 기초로 하여 전산해석을 통한 요소수 탱크, 유로, 기어펌프 내 동결요소수의 해동현상을 고찰하였다. 전산 해석에 사용된 수치해석 모델을 검증하기 위해서 Gau와 Viskanta12)가 수행한 순수 갈륨을 이용한 융해 실험데이터와 비교 분석을 진행하였다. 또한 요소수가 해동되어 이동하는 유로를 따라 포인트를 설정한 뒤, 각 포인트에서 액상분율, 온도, 속도 분포를 분석하여 상 경계면 이동 및 열전달 특성을 고찰하였다. 최종적으로 재질변경에 따른 유로 및 펌프내부의 동결된 요소수가 모두 해동되는 시간을 비교하여 두 가지 모델의 해동성능을 분석하였다.


2. 수치해석

2.1 지배방정식

본 연구에서는 요소수 탱크, 유로 및 기어펌프 내 해동현상에 관해서 Ansys Fluent13)를 이용한 전산 수치해석을 수행하였다. 또한 고정 격자계를 적용한 엔탈피 방법을 사용하여 고상과 액상 그리고 상 경계면의 상태를 0과 1 사이의 값으로 정의하였다. 시간에 따른 열전달 및 상 변화 현상을 계산하기 위하여 사용된 지배 방정식은 다음과 같다.

uit+uixi=0(1) 
ρuit+ujuixj=-pxi+xjμuixj-ρgiβT-Tref+Si(2) 

상 변화가 일어나는 해동현상의 경우 온도구배로 인한 유체의 밀도변화가 일어나 자연대류현상이 발생한다. 이를 고려하기 위해 체적항(Body force)에 Boussinesq 근사를 적용하여 밀도변화를 온도와 열팽창계수(Thermal expansion coefficient)β로 모델링하였다. 또한 속도 생성항인 Si를 도입하여 상 변화로 인한 속도 변화를 고려하였다. Si는 매개변수와 속도성분으로 구성되어 있는데, 매개변수 A는 액상 분율 f의 함수로 표현되며 다음과 같다.

Si=Aui(3) 
A=Amush1-f2/f3+ϵ(4) 

여기에서 Amush는 상 경계면 영역의 움직임을 나타내는 상수값이며 1×10-5로 설정하였다. 또한 ϵ은 상수값 1×10-3을 사용하여 분모가 0이 되지 않도록 설정하였다. 한편, 에너지 방정식은 다음과 같다.

ρHt+ujHxj=xjkTxj+Sh(5) 

여기에서, Sh는 상변화시 잠열의 흡수 및 방출을 나타내는 생성항이며 다음과 같이 나타낼 수 있다.

Sh=ρΔHt(6) 

식 (5)에서 H는 현열 엔탈피 h와 잠열 엔탈피 ∆H 의 합으로서 H=h+∆H 로 표현되고, 각 엔탈피의 표현식은 다음과 같다.

h=href+TrefTcpdT(7) 
ΔH=fL(8) 

여기에서 L은 상 변화가 일어날 때 필요한 에너지인 융해잠열을 나타내고, f는 액상 분율을 의미하며 온도 범위에 따라 다음과 같은 값을 가진다.

TTs        f=0(9a) 
TTi        f=1(9b) 
Ts<T<Ti        f=T-TsTl-Ts(9c) 

여기에서 f=0는 고체 상태, 0<f<1는 머시영역(Mush region) 즉 상 변화 과정인 상태이며 f=1은 액체 상태를 의미한다.

2.2 수치 해석 모델

본 연구에서 고려한 요소수 펌프 및 히터 통합 모듈을 모델링하여 Fig. 1에 나타내었다. 요소수 통합 모듈은 크게 상부 요소수 탱크 영역과 하부 펌프 모듈로 나눌수 있는데 상부에는 히터 모듈, 하우징, 펌프 모듈로 구성되어 있으며 히터모듈에는 열원인 PTC thermistor가 장착되어있다. 그리고 하부 펌프모듈은 요소수 공급을 위한 기어펌프와 요소수 유로가 포함된다. 동결된 요소수는 탱크와 Inlet pipe, Urea line, Gear pump 내 존재하며 요소수 탱크의 용량은 약 3.5 L이다. 또한 Inle pipe와 Urea line 그리고 Gear pump 내 존재하는 동결요소수의 체적은 각각 0.5 ml, 1.67 ml이다.

Fig. 1

Modeling of the integrated urea supply module and pump module

환원제로 사용되는 요소수는 32.5 wt%의 농도에서 약 –11.7 °C의 어는점을 가지는 공용혼합물이며 물성치는 Table 1에 나타내었다.

Properties of the 32.5 % AdBlueTM

초기온도는 요소수를 완전히 동결하고자 –20 °C로 설정하고 요소수를 저장하는 탱크의 벽면은 단열로 가정하였다. 또한 동결요소수의 상 변화가 일어날 때 온도에 따른 밀도와 비열의 변화는 선형적으로 변한다고 가정하였다. 한편, 열원인 PTC thermistor의 재질은 Al6061이며 160 W의 가열조건으로 설정하였다. 그리고 Gear bed plate, Pump gear, Pump body는 각각 SUS304, Steel, PPS+30 % 재질의 물성치를 적용하였다. Table 2는 각 Case에 따른 Inlet pipe, Gear housing 의 재질을 도시하였고 Table 3에서는 각 재질에 따른 주요 물성치를 나타내었다.

Material of each case

Property of each material

수치해석을 위한 격자는 품질이 우수한 Polyhedral mesh를 사용하여 최종적으로 약 110만개의 격자와 0.812의 품질을 가지는 격자를 구성하였다.


3. 결과 및 고찰

3.1 수치해석 검증

본 연구에서는 Gau와 Viskanta12)가 수행한 순수 갈륨 융해 실험을 분석한 후, 전산해석을 통한 결과와 실험값의 비교를 통해 본 모델의 수치해석 방법을 검증하였다. 순수 갈륨의 융해현상을 해석하기 위해 88.9 mm × 63.5 mm의 사각평면에 충분한 수의 격자를 생성하고 초기 갈륨의 온도를 301.45 K로 설정 하였다.좌측 벽면은 열원으로서 311.15 K의 온도 값을 가지고 나머지 모든 벽면은 단열로 가정하여 2차원 비정상상태 해석을 수행하였다. 수치해석 및 실험에 사용된 순수갈륨의 물성치는 Table 4에 나타내었다.

Properties of the gallium

가열 시간이 경과함에 따라 고체에서 액체로의 상 변화가 발생하게 되고 그에 따른 상 경계면의 이동에 대한 실험값과 수치해석 값을 Fig. 2에 도시하였다. 열원인 좌측벽면으로부터 갈륨의 융해현상이 발생하게 되고 상 경계면이 시간이 경과하면서 오른쪽으로 이동하는 현상을 볼 수 있다. 또한 Gau와 Viskanta12)가 수행한 실험데이터와 전산해석 결과 값이 비슷한 것으로 보아 본 연구에서 사용된 수치해석 방법이 타당한 것을 알 수 있다.

Fig. 2

Comparison of liquid fractions for melting of the pure gallium with schematic of the calculated domain

3.2 요소수 탱크 내 해동현상

Case1, 2에서 탱크 내 동결 요소수의 해동현상을 파악하기 위해 우레아 탱크의 중심을 기준으로 하여 시간에 따른 액상 분율의 변화를 Fig. 3에 나타내었다. 초기에는 요소수가 히터 모듈 주변부터 녹기 시작하고 시간이 경과하면서 점차 하부영역보다 상부 영역이 더 많이 해동되는 것을 볼 수 있다. 이는 히터 모듈 부근의 요소수가 가열되어 상승하는 유동과 동결 요소수 부근에서 냉각된 요소수가 하강하는 유동으로 인한 자연대류 현상 때문인 것으로 사료된다. 자연대류 현상으로 탱크 상부에는 온도가 상대적으로 높은 요소수가 위치하게 되고, 하부에는 온도가 낮은 요소수가 분포하게 된다. 따라서 시간이 경과할수록 탱크 상부영역이 하부 영역보다 더 많은 해동이 일어나게 된다.

Fig. 3

Liquid fraction distributions at center plane with time

한편, 요소수 탱크 내 열전달 및 상 변화 현상을 자세히 고찰하기 위해 Fig. 4에 시간에 따른 각 포인트의 온도 및 액상분율 분포를 도시하였다. Fig. 4(a)는 우레아 탱크 중심을 기준으로 Heat sink의 바닥면으로부터 각각 1/4, 3/4의 높이에 대해 y+방향으로 3 mm 간격으로 총 8개를 설정하였다. Fig. 4(b)는 h1영역 Point h1-1, h1-2, h1-3, h1-4에 대한 온도 및 액상분율 분포를 도시하였다. 그래프를 보면 가열 초기, 히터 싱크 에 가까울수록 (h1→ h4)온도 그래프의 기울기가 더 급격한 것을 볼 수 있다. 이는 히트싱크와 동결 요소수가 접촉되어 전도 열전달이 발생하고, 히트 싱크에 가까울수록 많은 에너지가 더 빠르게 전달되기 때문이다. –12∼–11 °C 부근에서 액상 분율이 1에 도달하는 것으로 보아 상 변화가 일어나는 것을 알 수 있고, 온도가 일정하게 유지되는 현상은 요소수가 용융점에 도달하여 잠열의 흡수가 일어나기 때문이다.

Fig. 4

Liquid fraction distributions at center plane with time

이후, 상 변화에 사용된 에너지가 현열현상을 통해 온도가 급격히 상승한다. 그리고 일정구간 변동(Fluctuation)이 발생하는 것을 볼 수 있는데, 이는 상 경계면을 기준으로 전도 및 대류 열전달이 복합적으로 발생하고 동시에 자연대류로 인한 유동이 형성되어 불안정한 상태가 형성되기 때문이다. 또한 가열시간이 지속되면서 온도가 점차 일정해지는 것은 해동된 요소수가 열원으로부터 얻은 에너지를 동결된 요소수로 전달하여 열평형이 이루어지기 때문인 것으로 사료된다. Fig. 4(c)는 h2영역에서 Point h2-1, h2-2, h2-3, h2-4의 온도 및 액상분율 분포를 도시하였다. 가열초기, h1영역의 온도분포와는 다르게 h2-1을 제외한 나머지 Point에서 서서히 온도가 상승되는 것을 볼 수 있다. 열원 근처에 있는 Point h2-1은 가열 초기 전도에 의한 열전달이 주로 일어나기 때문에 point h1-1과 온도분포의 차이가 없다. 하지만 point h2-2, h2-3, h2-4는 히트 싱크를 기준으로 h1영역보다 하부에 위치하고 있어 주위에 상대적으로 낮은 온도의 액상 요소수가 분포하게 된다. 따라서 대류 열전달과 자연대류 현상이 상부에 비해 비교적 약하기 때문에 온도가 서서히 상승하고, 상변화는 느리게 진행되는 것을 확인 할 수 있다. 또한 자연대류의 구심력(Driving force)이 약하기 때문에 이전 Fig. 4(a)의 h2영역에서 발생한 변동은 본 그래프에서 약하게 발생하는 것으로 사료된다.

3.3 요소수 유로 및 기어펌프 내 해동현상

요소수가 이동하는 유로 내 해동현상을 파악하기 위하여 Fig. 5와 같이 Inlet pipe와 Gear housing의 중심 Plane에서 시간에 따른 각 Case의 액상분율 분포를 도시하였다. Case1의 경우 가열시간이 90 s 일 때, Inlet pipe 내 요소수의 해동이 이루어지고 Urea line과 Gear pump 내 해동은 일어나지 않는다. 그리고 가열시간이 지속되면서 Gear pump 중심을 기준으로 좌측에서 우측으로 Urea line과 Gear pump 내 요소수 해동이 진행되는 것을 볼 수 있다. 한편, Case2에서는 90 s 일 때, Case1보다 Inlet pipe 내 요소수 해동이 조금 더 이루어진 것을 볼 수 있고, 180 s일 때는 Inlet pipe 뿐만 아니라 Urea line과 Gear pump 내 전체 영역에 걸쳐 해동이 진행되고 있는 것을 알 수 있다. 그리고 300 s일 때는 Gear pump 내 중심 두 영역을 제외하고 모든 영역에서 해동이 이루어 졌으며, 그 이후에는 모든 영역의 요소수가 해동된 것으로 보인다. 이와 같이 Case1과 Case2의 해동현상은 각 Case에 사용된 재질의 전도 특성에 따라 다른 것으로 판단된다. 열전도율이 높은 재질을 사용한 Case2는 Inlet, Outlet pipe, Gear housing에서 빠르게 열이 전달되어 전체영역에서 동시에 해동이 일어난다. 반면, 낮은 전도율의 재질을 사용한 Case1은 열원과 가까운 영역(Inlet, Outlet pipe)에서 점차적으로 열이 전달되며, Case2보다 늦은 해동속도를 보인다. Fig. 6은 시간에 따른 하부영역의 온도변화를 나타낸 것으로서, 열원으로부터 열이 상부에서 하부로 전달되면서 온도가 상승되는 모습을 볼 수 있다. 또한 상대적으로 높은 전도율의 재질을 사용한 Case2가 Case1보다 더 빨리 온도가 상승하는 것을 알 수 있다. 한편, 요소수 유로 내 요소수 해동속도를 비교하기 위해서 Fig. 7과 같이 Inlet pipe, Urea line, Gear pump 내 총 16개의 Point를 설정하고, Fig. 8에는 시간에 따른 각 포인트의 액상분율 변화를 그래프로 도시하였다. 여기서 유로 내 동결요소수의 해동속도는 요소수 공급시점 및 펌프작동시점을 결정하므로 해동성능에 있어 중요한 결과중 하나이다. Inlet pipe의 중간영역인 Point u-4와 Gear pump 출구 영역인 Point u-16을 보면, Case1은 124 s, 553 s, Case2는 117 s, 310 s에 액상분율이 1에 도달하는 것을 볼 수 있다. 그리고 Point u-5 이후부터 Case1과 Case2의 액상분율이 1에 도달하기까지 걸리는 시간 차이가 이전 Inlet pipe 영역에서보다 더 큰 것을 볼 수 있다. 이는 Case1과 Case2의 Inlet pipe에 사용된 ADC12와 AL6061의 전도율 차이보다 Gear housing에 사용된 SUS304와 AL6061의 전도율 차이가 훨씬 크기 때문이다. 결과적으로 두 가지 Case 모두 Point u-16에서 해동이 완료되면 Inlet pipe, Urea line, Gear pump 내 모든 동결요소수가 해동되므로 펌프작동을 통한 요소수 공급이 가능할 것이라 사료된다.

Fig. 5

Liquid fraction distributions at yz and xy plane with time

Fig. 6

Temperature distributions at yz and xy plane with time

Fig. 7

Data points at urea path

Fig. 8

Liquid fraction distributions at each point with time


4. 결 론

본 연구에서는 전산수치해석을 통해 두 가지 재질을 적용한 SCR 통합 모듈 내 동결요소수의 해동 및 열전달 특성을 고찰하였다. 요소수 탱크 내부에서는 가열초기, 히트싱크 주변으로 전도에 의한 열전달이 주로 이루어졌지만, 요소수가 고상에서 액상으로 상변화가 일어나면서 대류에 의한 열전달이 지배적으로 이루어지는 것을 알 수 있었다. 한편, 재질변경에 따른 하부 요소수 유로의 해동완료 시점을 알기위해 요소수 유로 내 각 Point에서 액상분율변화를 관찰하였다. 유로 내 모든 동결 요소수가 해동되기까지 Case1은 553 s, Case2는 310 s가 소요되어 전도율이 높은 재질을 사용한 Case2가 더 빨리 해동되는 것을 확인하였다. 이를 통해 Inlet pipe와 Gear hosing의 재질이 유로 내 요소수 해동속도에 큰 영향을 미치는 것을 알 수 있었다.

Acknowledgments

본 연구는 중소기업청의 재원으로 한 WC300 R&D 사업의 지원(S2307213) 그리고 2018년 정부(미래창조과학부)의 재원으로 한국연구재단의 지원(No. NRF-2018R1D1A1B07048355)을 받아 수행된 연구이며 이에 감사드립니다.

References

  • J. Kim, “Diesel Vehicle NOx Reduction Strategy,” Auto Journal, KSAE, Vol.37, No.9, pp.48-52, 2017.
  • A. Dimaratos, D. Tsokolis, G. Fontaras, S. Tsiakmakis, B. Ciuffo and Z. Samaras, “Comparative Evaluation of the Effect of Various Technologies on Light-duty Vehicle CO2 Emission over NEDC and WLTP,” Transport Research Procedia, Vol.14, pp.3169-3178, 2016. [https://doi.org/10.1016/j.trpro.2016.05.257]
  • I. G. Kim and T. W. Lee, “Environmental Policy to Reduce Emissions of In-use Vehicles,” Auto Journal, KSAE, Vol.40, No.8, pp.20-24, 2018.
  • M. Y. Kim, “Performance Prediction of SCR-DeNOx System for Reduction of Diesel Engine NOx Emission,” Transactions of KSAE, Vol.11, No.3, pp.71-76, 2003.
  • V. R. Voller, M. Cross and N. C. Markatos, “An Enthalpy Method for Convection/Diffusion Phase Change,” International Journal for Numerical Methods in Engineering, Vol.24, No.1, pp.271-284, 1987. [https://doi.org/10.1002/nme.1620240119]
  • K. S. Kwon, B. J. Son, M. S. Yeo and S. J. Yoon, “Numerical Analysis by the Enthalpy Method in the Phase Change Process,” Korean Journal of Air-Conditioning and Refrigeration Engineering, Vol.6, No.2, pp.67-77, 1994.
  • P. J. Oliveira and R. I. Issa, “An Improved PISO Algorithm for the Computation of Buoyancy-Driven Flows,” Numerical Heat Transfer B, Vol.40, No.6, pp.473-493, 2001. [https://doi.org/10.1080/104077901753306601]
  • S. Y. Lee, M. Y. Kim, C. H. Lee and Y. B. Park, “Numerical Investigation of the Urea Melting and Heat Transfer Characteristics with Three Different Types of Coolant Heaters,” Transactions of KSAE, Vol.20, No.4, pp.125-132, 2012. [https://doi.org/10.7467/KSAE.2012.20.4.125]
  • T. H. An, S. Y. Lee, C. H. Lee, Y. B. Park and M. Y. Kim, “Numerical Investigation on the Urea Melting Characteristics with Coolant and Electric Heaters,” KSAE13-H0014, pp.87-93, 2013.
  • T. Kim, B. Choi and S. Jung, “Numerical Analysis of the Melting Characteristics of a Frozen Urea-water Solution by Heat Spreaders with a Positive Temperature Coefficient Heater,” Applied Thermal Engineering, Vol.119, pp.275-282, 2017. [https://doi.org/10.1016/j.applthermaleng.2017.03.069]
  • D. K. Lee, M. Y. Kim, K. C. Oh and T. S. No, “Numerical Study on the Urea Thawing and Heat Transfer Characteristics with Electric Heating Module for Light Commercial Vehicle,” Transaction of the Korean Society of Mechanical Engineers(B), Vol.43, No.3, pp.161-168, 2019. [https://doi.org/10.3795/KSME-B.2019.43.3.161]
  • C. Gau and R. Viskanta, “Melting and Solidification of a Pure Metal on a Vertical Wall,” Journal of Heat Transfer, Vol.108, pp.174-181, 1986. [https://doi.org/10.1115/1.3246884]
  • Ansys Inc., Fluent, Release 12.0 User’s Guide, 2009.

Fig. 1

Fig. 1
Modeling of the integrated urea supply module and pump module

Fig. 2

Fig. 2
Comparison of liquid fractions for melting of the pure gallium with schematic of the calculated domain

Fig. 3

Fig. 3
Liquid fraction distributions at center plane with time

Fig. 4

Fig. 4
Liquid fraction distributions at center plane with time

Fig. 5

Fig. 5
Liquid fraction distributions at yz and xy plane with time

Fig. 6

Fig. 6
Temperature distributions at yz and xy plane with time

Fig. 7

Fig. 7
Data points at urea path

Fig. 8

Fig. 8
Liquid fraction distributions at each point with time

Table 1

Properties of the 32.5 % AdBlueTM

Density (solid), ρs [kg/m3] 1,010
Density (liquid), ρl [kg/m3] 1,090
Dynamic viscosity, μl [kg/ms] 1.4
Reference density, ρref [kg/m3] 1095.6
Reference temperature, Tref [K] 283.16
Freezing(and Melting) point [K] 261.46
Latent heat, L [kJ/kg] 152.86
Thermal expansion coefficient, β [1/K] 4.5×10-4
Thermal conductivity, k [W/mK] 0.57

Table 2

Material of each case

Gear housing Inlet pipe
Case1 SUS304 ADC12
Case2 Al6061 Al6061

Table 3

Property of each material

Specific heat, cp [J/kgK] Thermal conductivity, k [W/mK]
Al6061 896 180
ADC12 962.78 96
SUS304 500 16.2

Table 4

Properties of the gallium

Density (liquid), ρl [kg/m3] 6,093
Dynamic viscosity, μl [kg/ms] 1.8×10-3
Reference density, ρref [kg/m3] 6,095
Reference temperature, T ref [K] 302.93
Melting point [K] 302.93
Latent heat, L [kJ/kg] 801.60
Thermal expansion coefficient, β [1/K] 1.2×10-4
Thermal conductivity, k [W/mK] 32