다중상 흐름과 상변화 열전달을 위한 격자볼츠만 방법

초록

본 리뷰는 격자볼츠만(LB) 방법을 이용한 다중상 흐름·열전달 및 상변화 현상의 최신 연구 동향을 정리한다. 이론적 틀, 열·다중상 모델의 일관성 문제, 스퓨리어스 전류 억제, 수치 안정성 향상 등을 다루며, 고밀도비·고레놀즈수·가변 표면장력 구현을 위한 최신 개선책을 제시한다. 연료전지·배터리, 물방울 충돌, 끓음·증발 등 다양한 응용 사례와 향후 연구 방향도 조명한다.

상세 분석

격자볼츠만(LB) 방법은 미시적 입자 충돌 과정을 확률론적으로 모델링함으로써 연속체 방정식을 유도하는 메소스로, 전통적인 유한차분·유한요소 방식에 비해 격자 구조와 충돌 연산이 병렬화에 유리하고 복잡한 경계조건을 손쉽게 구현할 수 있다는 장점을 가진다. 특히 다중상 흐름과 열전달, 그리고 상변화 현상을 동시에 다루는 경우, 전통적인 CFD는 인터페이스 추적·재구성 비용이 크게 증가하는 반면, LB는 인터페이스를 암시적으로 처리하는 페이즈-필드 혹은 샤프 인터페이스 모델을 자연스럽게 통합한다는 점에서 효율적이다.

하지만 초기 LB 다중상 모델은 열역학적 일관성 부족, 비물리적 스퓨리어스 전류(spurious currents), 제한된 밀도비와 레이놀즈수 등 여러 결함을 안고 있었다. 최근 연구는 자유에너지 기반 모델에 정확한 EOS(Equation of State)를 결합하고, 다중분포 함수(Multi‑distribution)와 이중분포 함수(Dual‑distribution) 체계를 도입해 열·운동량을 별도로 진화시킴으로써 열역학적 일관성을 크게 개선하였다. 또한, 고차 순간(High‑order moments)와 중앙 차분(Central differencing) 스킴을 활용한 충돌 연산은 스퓨리어스 전류를 10⁻⁴ 수준으로 감소시켰으며, 멀티레벨 그리드(Multilevel grid)와 엔트로피 기반 충돌 연산은 수치 안정성을 향상시켜 밀도비 10⁴, 레이놀즈수 10³ 수준까지 시뮬레이션이 가능하도록 만들었다.

표면장력과 접촉각을 조절하는 방법으로는 자유에너지 포텐셜을 직접 수정하거나, 위상장벽(Phase‑field) 파라미터를 동적으로 변환하는 기법이 제안되었다. 이러한 기법들은 물리적 접촉각을 정확히 재현하면서도 인터페이스 두께를 얇게 유지해 계산 효율을 높인다. 열전달 모델에서는 온도 구분 함수와 에너지 분포 함수를 결합한 이중분포 LB가 주류를 이루며, 비등점 근처의 잠열(latent heat) 처리를 위해 가변 열전도도와 잠열 소스 항을 명시적으로 포함한다. 이때, 잠열 항을 충돌 단계에 삽입하거나 스트리밍 단계에 보정함으로써 에너지 보존성을 유지한다.

응용 측면에서, 연료전지의 물-가스 전이 현상, 배터리 전해질 내 기포 발생, 미세 물방울 충돌 시 인터페이스 재구성, 끓음 및 증발 과정에서의 온도·밀도 구배 등 복합 현상을 LB가 효과적으로 재현한다는 사례가 다수 보고되었다. 특히, 고밀도비와 고레놀즈수를 동시에 만족하는 모델은 실제 산업 현장의 복잡한 열·유동 현상을 실시간에 가깝게 예측할 수 있게 하여 설계 최적화와 성능 향상에 기여한다.

향후 과제는 다중물리 현상의 통합 시뮬레이션(전기·화학·기계·열)과, GPU·FPGA 기반 초고속 병렬 구현, 그리고 머신러닝을 이용한 파라미터 자동 튜닝이다. 이러한 방향은 LB 방법을 전통적인 CFD와 차별화된 차세대 멀티스케일 시뮬레이션 플랫폼으로 자리매김하게 할 것이다.