Numerical Study on the Effects of Wall Wettability on Heat Transfer and Flow Characteristics during Flow Boiling in Microchannels

기술 요약: 마이크로채널 내 유동 비등 중 벽면 젖음성이 열전달 및 유동 특성에 미치는 영향에 관한 수치 연구

문제 정의
고성능 칩의 소형화는 국부 열유속 밀도를 1000 W/cm² 이상으로 급격히 증가시켰으며, 이는 기존의 냉각 방식으로는 해소하기 어려운 수준이다. 마이크로채널 유동 비등은 높은 열전달 계수와 잠열 흡수 능력 덕분에 유망한 해결책이다. 그러나 이러한 시스템의 성능는 기포 역학 및 유동 패턴 진화에 의해 크게 좌우된다. 이 과정에서 결정적이지만 복잡한 요소는 벽면 젖음성이다. 기존 문헌들이 친수성, 소수성 및 혼합 젖음 표면에 대해 탐구해 왔으나, 극단적인 젖음성 조건이 기포 거동, 유동 안정성, 열전달 효율 및 열역학적 비가역성(엔트로피 생성) 사이의 상호작용에 어떻게 영향을 미치는지에 대한 결합된 메커니즘은 여전히 충분히 이해되지 않고 있다. 특히, 다양한 젖음성 영역 하에서 향상된 열전달과 유동 안정성 사이의 절충 관계를 정량화할 필요가 있다.

연구 방법론
본 연구는 그루브(groove)가 있는 마이크로채널 내 유동 비등을 시뮬레이션하기 위해 체적 유체법(Volume of Fluid, VOF)에 기반한 3차원 수치 모델을 사용한다. 시뮬레이션은 표면 장력 효과를 포착하기 위해 연속체 표면 힘(Continuum Surface Force, CSF) 모델을 사용하며, 상변화 처리를 위해 Lee 등의 증발-응축 모델을 활용한다.

• 물리적 모델: 기하학적 구조는 기포 핵 생성 및 성장을 연구하기 위해 설계된 측벽 삼각 그루브를 포함하는 마이크로채널로 구성된다. 작동 유체는 탈이온수이며, 고체 영역은 실리콘이다.
• 젖음성 조건: 두 가지 극단적인 표면 조건을 비교한다: 친수성 벽면(접촉각 $\theta = 30^\circ$)과 소수성 벽면($\theta = 150^\circ$).
• 수치적 접근법: 질량, 운동량 및 에너지 보존 법칙을 따르는 지배 방정식을 유한 부피법을 사용하여 해결하며, 압력-속도 결합을 위해 PISO 알고리즘을 사용한다. 격자 독립성 검증을 통해 약 275,000개의 요소를 가진 격자가 충분한 정확도를 제공함을 확인하였다. 모델은 실험 데이터와 비교 검증되었으며, 칩 베이스 온도 상승치에서 최대 오차는 8% 미만으로 나타났다.
• 분석 지표: 본 연구는 기포 역학, 유동 패턴, 출구 온도 변동, 누셀 수(Nusselt number, Nu), 압력 강하 및 엔트로피 생성(열적 및 마찰적 성분 모두 포함)을 평가한다.

주요 기여 및 결과

1. 기포 역학 및 유동 패턴 진화:

   • 친수성 표면 ($\theta = 30^\circ$): 강한 액체-고체 상호작용으로 인해 기포가 벽면에 부착되어 합쳐지며 액체 막 구조를 형성한다. 이는 환상 유동(annular flow)으로의 전이를 촉진하지만, 기포의 이탈을 방해한다.
   • 소수성 표면 ($\theta = 150^\circ$): 기포가 빠르게 탈착되어 주류로 운반된다. 이는 증기 축적을 가속화하며, 기포 유동(bubbly flow)에서 슬러그 유동(slug flow) 및 환상 유동으로의 조기 전이를 유도한다.

2. 열적 성능:

   • 열전달 효율: 소수성 표면이 우수한 열전달 성능을 보인다. 빠른 기포 탈착은 잠열 제거를 용이하게 하여 낮은 벽면 온도와 높은 누셀 수를 유지한다. 반면, 친수성 표면은 부착된 기포로 인한 열 축적을 겪으며, 이로 인해 벽면 온도가 상승하고 누셀 수가 지속적으로 감소한다.
   • 온도 안정성: 소수성 표면은 평균 열전달을 개선하지만, 친수성 표면에 비해 출구 온도 불안정성이 낮다. 친수성 표면은 지속적인 액체 막에 의한 열 방산 저해로 인해 더 큰 온도 변동을 보인다.

3. 유동 특성 및 불안정성:

   • 압력 강하: 소수성 채널은 현저히 높은 압력 강하와 더 뚜렷한 압력 진동을 경험한다. 이는 높은 증기 체적 분율과 기포의 주기적인 배출이 국부 유동 구조를 변화시키기 때문이다.
   • 속도장: 소수성 채널은 활발한 기포 역학으로 인해 더 강한 국부 유동 교란과 높은 속도 구배를 보이는 반면, 친수성 채널은 더 완만하고 낮은 속도 변화를 보인다.

4. 엔트로피 생성 분석:

   • 열적 vs 마찰적: 본 연구는 열적 엔트로피 생성(온도 구배에 의해 구동)과 마찰적 엔트로피 생성(점성 소산에 의해 구동)을 구분한다.
   • 지배적 메커니즘: 열적 엔트로피 생성이 총 엔트로피 생성의 주요 기여 요인이다. 친수성 표면은 열전달 저하로 인한 큰 온도 구배 때문에 현저히 높은 열적 엔트로피를 생성한다.
   • 절충 관계: 소수성 표면이 강한 속도 구배로 인해 더 많은 마찰적 엔트로피를 생성함에도 불구하고, 상당한 열적 비가역성 감소 덕분에 총 엔트로피 생성은 친수성 표면보다 낮게 유지된다.

의의 및 주장
저자들은 본 연구가 고열유속 조건 하에서 마이크로채널 방열판의 젖음성 설계 및 열-수력 시스템 최적화를 위한 가치 있는 이론적 지침을 제공한다고 주장한다. 연구 결과는 근본적인 절충 관계를 강조한다: 소수성 표면은 비등 열전달을 크게 향상시키고 열적 비가역성을 줄이지만, 더 큰 유동 저항과 압력 변동을 초래한다. 반면, 친수성 표면은 보다 안정적인 압력 강하를 제공하지만, 심각한 열전달 저하와 높은 열역학적 손실을 겪는다.

본 논문은 젖음성 최적화가 상변화 열전달의 향상과 유동 불안정성 증가에 따른 페널티 사이의 균형을 맞추는 작업임을 결론짓는다. 결과는 열 성능과 엔트로피 최소화를 우선시하는 응용 분야의 경우, 관련 압력 강하 및 유동 진동을 수용할 수 있다면 소수성 표면이 유리함을 시사한다. 본 연구는 현재의 수치 조사 범위를 넘어선 특정 미래 응용 분야나 실험적 검증을 제안하지는 않으나, 효과적인 열 관리 설계를 위한 결합된 메커즘 이해의 필요성을 강조한다.