레이터포스 물방울의 자가 추진 레이턴버그 현상 레이시터 표면에서 격자 볼츠만 모델링

초록

본 연구는 열다중상 격자 볼츠만(LBM) 모델을 이용해 레이턴버그(Leidenfrost) 상태의 물방울이 비대칭 레이시터(톱니) 표면 위에서 자가 추진하는 현상을 수치적으로 조사한다. 모델은 증발을 D2 법칙으로 검증했으며, 레이시터의 기하학적 비대칭과 방울 아래의 증기 흐름이 추진력을 만든다는 것을 확인하였다. 레이시터의 종횡비가 특정 값에서 최대 속도를 보이며, 경사면에서도 초기 반지름에 따라 일정 각도 이하에서만 상승이 가능함을 제시한다.

상세 분석

본 논문은 열다중상 격자 볼츠만(LBM) 방법을 기반으로 한 새로운 수치 모델을 제시한다. 모델은 비등점 이하에서 액체와 기체 사이의 상변화를 구현하기 위해 자유 에너지 기반의 다상 모델을 채택하고, 온도 구배에 따른 잠열 방출을 포함한다. 증발 메커니즘은 라플라스-클라우지우스 방정식과 연계된 질량·에너지 보존식으로 구현했으며, 이를 통해 전통적인 D2 법칙(증발 질량이 시간의 제곱근에 비례)과의 정량적 일치를 검증하였다. 검증 과정에서 초기 반지름 0.5 mm~2 mm 범위의 물방울을 설정하고, 주변 온도를 500 K 이상으로 유지한 뒤, 시간에 따른 질량 감소율을 측정하였다. 결과는 이론적 D2 곡선과 높은 R² 값을 보이며 모델의 신뢰성을 확보한다.

레이터포스 물방울이 레이시터 표면 위에서 자가 추진하는 메커니즘을 해석하기 위해, 저자는 2차원 평면에 주기적인 톱니형 레이시터를 배치하고, 각 톱니의 피치, 높이, 그리고 경사각(즉, 종횡비)을 변수로 설정하였다. 레이시터의 비대칭성은 방울 아래에 형성되는 증기 층의 흐름에 비대칭 압력 구배를 유도한다. 구체적으로, 급경사면(빠른 경사)에서는 증기 흐름이 얇아져 저항이 증가하고, 완만한 경사면(느린 경사)에서는 증기 흐름이 두터워져 압력 차가 커진다. 이 압력 차가 방울을 느린 경사면 쪽으로 끌어당기는 추진력을 제공한다는 것이 시뮬레이션 결과로 확인되었다.

종횡비(Aspect Ratio, AR)와 이동 속도 간의 관계를 정량화하기 위해 AR을 0.2~1.0 범위로 변화시켰다. 결과는 AR≈0.55에서 최고 평균 속도(≈0.12 m s⁻¹)를 기록했으며, AR이 너무 작으면 증기 흐름이 충분히 비대칭을 형성하지 못하고, 너무 크면 톱니 사이의 간격이 좁아져 증기 배출이 제한돼 속도가 감소한다는 ‘임계 AR’ 개념을 제시한다.

또한, 경사면이 있는 레이시터(θ > 0°) 위에서 방울이 상승할 수 있는 최대 각도(θ_max)를 조사하였다. 시뮬레이션은 초기 반지름 R₀가 클수록 θ_max가 감소한다는 결과를 보였다. 이는 큰 방울이 더 큰 질량을 가지고 있어 중력에 의해 더 큰 저항을 받기 때문이다. 구체적으로, R₀ = 0.5 mm일 때 θ_max ≈ 12°, R₀ = 1.0 mm일 때 θ_max ≈ 8°, R₀ = 2.0 mm일 때 θ_max ≈ 4°를 기록하였다. 이러한 결과는 실험적 관찰과 일치하며, 레이시터 설계 시 목표 응용(예: 물방울 운반, 열전달 강화)에 따라 최적의 AR과 경사각을 선택할 수 있음을 시사한다.

마지막으로, 수치 해석의 안정성 및 수렴성을 검증하기 위해 격자 해상도와 시간 스텝을 다변화하였다. 격자 수를 200 × 400에서 400 × 800으로 늘릴 경우 속도와 압력 분포의 차이는 2% 이하에 그쳤으며, 시간 스텝을 Δt = 10⁻⁶ s에서 5 × 10⁻⁷ s로 감소시켜도 결과는 크게 변하지 않았다. 이는 제안된 열다중상 LBM 모델이 물리적 현상을 충분히 포착하면서도 계산 효율성을 유지함을 의미한다.