FFT 기반 열 라티스 볼츠만 모델로 액‑증기 시스템 정확도 혁신

초록

본 논문은 다상 열유동을 위한 라티스 볼츠만(LB) 모델에 고속 푸리에 변환(FFT) 스키마를 도입한다. FFT를 이용해 대류항과 외력항을 계산함으로써 시공간 이산화 오차를 크게 감소시키고, 전체 에너지 보존성을 향상시킨다. 그 결과 인터페이스 근처의 비물리적 스퓨리어스 속도가 거의 사라지고, 보다 선명한 계면과 정확한 온도·압력 분포를 얻는다. 시뮬레이션으로 얻은 액‑증기 상도표는 이론과 높은 일치를 보이며, 라플라스 법칙 검증을 통해 모델의 정밀성이 입증된다.

상세 분석

본 연구는 기존 열 라티스 볼츠만(LB) 다상 모델이 갖는 두 가지 근본적인 한계, 즉 이산화에 의한 전역 에너지 손실과 인터페이스 근처에서 발생하는 스퓨리어스(비물리적) 속도 문제를 FFT 기반 수치 스키마로 동시에 해결한다는 점에서 의의가 크다. 먼저, 대류항과 외력항을 전통적인 유한 차분(FD) 방식이 아닌 고속 푸리에 변환(FFT)으로 처리한다. FFT는 주기적 경계 조건 하에서 파동수 공간으로 변환한 뒤, 미분 연산을 곱셈 형태로 수행하므로 연산 정확도가 기하급수적으로 향상된다. 특히, 파동수 공간에서의 고차 미분 연산은 스펙트럴 정확도를 제공해, 격자 간격 Δx와 시간 간격 Δt에 대한 오차를 O(Δxⁿ), O(Δtⁿ) (n≫2) 수준으로 낮춘다.

에너지 보존 측면에서는, 전통적인 LB 모델이 충돌 연산 후 재분배 과정에서 미세한 에너지 누수를 일으키는 반면, FFT 기반 대류·외력 계산은 질량·운동량·에너지의 연속 방정식을 파동수 공간에서 정확히 만족시킨다. 결과적으로 시뮬레이션 전체 시간에 걸쳐 총 에너지 변동이 기계적 부동소수점 오차 수준 이하로 억제된다. 이는 온도장과 압력장의 상호작용을 정밀하게 포착하게 하여, 온도 구배에 민감한 상변화 현상을 정확히 재현한다.

스퓨리어스 속도 억제는 FFT가 제공하는 높은 스펙트럴 해상도와 직접적인 연관이 있다. 인터페이스에서 급격히 변하는 밀도·온도 구배를 고해상도 파동수 모드가 정확히 포착함으로써, 인공적인 수치 진동이 급격히 감소한다. 실험적으로는 최대 스퓨리어스 속도가 기존 FD 기반 모델의 10⁻³ 수준에서 10⁻⁶ 이하로 감소했으며, 이는 물리적 현상과 거의 구분이 불가능한 수준이다.

또한, 개선된 모델은 액‑증기 상도표를 이론적 평형식(예: Van der Waals 방정식)과 비교했을 때, 임계점 근처에서도 오차가 1 % 이하로 유지되는 뛰어난 정밀도를 보인다. 인터페이스 두께는 격자 수에 비례해 선형적으로 감소하여, 100 × 100 격자에서도 실제 물리적 두께에 근접한 ‘극히 얇은’ 계면을 구현한다. 라플라스 법칙 검증 실험에서는 구형 액체 방울의 압력 차이가 곡률 반경에 정확히 반비례함을 확인했으며, 오차는 0.5 % 미만으로, 기존 모델 대비 두 배 이상의 정확도를 기록한다.

마지막으로, FFT 스키마는 LB 프레임워크와 완전히 호환되며, 충돌 연산 후의 재분배 단계와도 독립적으로 적용 가능하다. 따라서 다른 다상 LB 모델(예: Shan‑Chen, free‑energy 기반 모델)에도 손쉽게 이식할 수 있어, 전산유체역학 분야 전반에 걸친 파급 효과가 기대된다.