히터 열전도율이 핵형 풀보일링에 미치는 영향: 다중스케일 오토마타 시뮬레이션

초록

**
본 연구는 구리와 실리콘 옥사이드를 히터 재료로 선택하여, 열확산도와 전도도의 차이가 핵형 풀보일링의 버블 발생 주기, 마이크로레이어 회복, 온도 변동성 등에 미치는 영향을 다중스케일 셀룰러 오토마타 모델로 정량화하였다. 구리는 높은 열확산도로 냉각 효율을 크게 향상시키며, 실리콘 옥사이드는 저열전도성으로 인한 이웃 버블 간 상호작용이 강조되어 온도 변동이 뚜렷하게 나타난다. 시뮬레이션 결과는 기존 반경험식과 일치함을 보여 모델의 신뢰성을 확인한다.

**

상세 분석

**
본 논문은 히터 내부의 열전도와 액체와의 복합 열전달을 동시에 다루는 2‑scale 오토마타 프레임워크를 제시한다. 히터는 1 × 1 cm², 두께 0.25 cm의 정육면체 셀(250 µm³) 16 000개로 이산화되며, 각 셀에서 수정된 ADI 방법을 이용해 Fourier 방정식을 풀어 시간 단계 10⁻³ ~ 10⁻⁴ s를 적용한다. 열손실 메커니즘은 자연대류, 마이크로컨벡션, 마이크로레이어·과열액 기화, 복사 네 가지로 분리되고, 각각의 열유량은 셀의 젖음 상태와 버블 존재 여부에 따라 동적으로 전환된다. 버블 핵생성은 평균 6 µm, 표준편차 3 µm인 가우시안 분포를 따르는 20개의 무작위 캐비티에서 시작되며, 접촉각 φ = 20°가 고정된다. 버블 성장·분리 조건은 표면 초과열, 표면 장력, 열전달 계수 등을 포함한 식(3)‑(5)로 정의되며, 마이크로레이어 두께 δ는 κ/h 로 계산한다.

핵심 결과는 구리(α = 116.6 mm²/s, κ = 401 W/mK)와 실리콘 옥사이드(α = 0.834 mm²/s, κ = 1.38 W/mK)의 열확산도 차이가 버블 성장 시간과 마이크로레이어 회복 시간에 직접적인 영향을 미친다는 점이다. 구리에서는 열이 빠르게 퍼져 버블 성장 주기가 짧아지고, 인접 버블 간 열간섭이 최소화돼 온도 파형이 비교적 평탄하다. 반면 실리콘 옥사이드는 열전도와 확산이 느려 버블 성장·분리 과정이 길어지고, 인접 버블이 만든 국부 과열 영역이 지속돼 온도 변동이 크게 나타난다. 이러한 현상은 시뮬레이션이 제시한 버블 주파수 차이와 마이크로레이어 회복 시간 차이에서 명확히 드러난다. 또한, 모델은 기존의 반경험식(예: Rohsenow, Cooper 등)과 정량적으로 일치함을 보이며, 복잡한 표면 패턴이나 비균일 습윤성에서도 신뢰할 수 있는 예측 능력을 입증한다.

계산 효율성 측면에서, 전통적인 DNS가 요구하는 수십억 셀·시간 단계 대비 현재 오토마타 모델은 2‑3 order magnitude 적은 연산량으로 동일한 물리 현상을 재현한다. 이는 설계 단계에서 다양한 재료·표면 구조를 빠르게 탐색할 수 있게 해, 전자·핵산업 등 고열밀도 냉각 요구 분야에 실용적인 도구가 될 가능성을 시사한다. 다만, 모델은 비임계 상태의 풀보일링에만 적용되며, 버블 탈착 후의 후유동이나 임계 근처 현상은 다루지 않는다. 향후 연구에서는 3‑D 셀 구조 확대, 비정상적인 접촉각 분포, 그리고 고압·고온 조건을 포함한 확장 모델링이 필요하다.

**