3차원 자연대류를 위한 이중 MRT 라티스 볼츠만 시뮬레이션 혁신

초록

본 논문은 유동과 열전달을 동시에 다루는 3차원 이중 다중완화시간(MRT) 라티스 볼츠만 방법을 최초로 제안한다. 좌·우 수직벽을 서로 다른 온도로 유지하고 나머지 네 면을 절연 또는 선형 온도변화 조건으로 설정한 세 종류의 정육면체 자연대류 문제를 다양한 레이리 수에서 해석하였다. 첫 두 경우는 기존 실험·수치 benchmark와 거의 일치하는 결과를 보였으며, 세 번째 경우는 전·후면이 선형 온도분포, 상·하면이 절연인 보다 일반적인 3차원 흐름 특성을 성공적으로 재현하였다.

상세 요약

이 연구는 라티스 볼츠만 방법(LBM)의 핵심 한계인 단일 완화시간(BGK) 모델이 갖는 수치적 불안정성과 낮은 정확도를 극복하기 위해 다중완화시간(MRT) 모델을 채택한 점이 가장 큰 특징이다. MRT는 충돌 연산을 순간적인 모멘트 공간에서 수행함으로써 각 모멘트에 별도의 완화계수를 부여할 수 있어, 점성·열전도 계수와 같은 물리적 파라미터를 독립적으로 조정한다. 특히 3차원 흐름에서는 비대칭적인 스트레스와 열 플럭스가 동시에 발생하므로, 이러한 자유도가 수치적 안정성 확보와 고정밀 해석에 필수적이다.

논문은 이중 MRT 모델을 구현한다. 하나는 흐름을 위한 D3Q19 격자를, 다른 하나는 온도장 전파를 위한 D3Q7 격자를 사용한다. 두 격자는 각각 독립적인 충돌 연산과 스트리밍 과정을 거치며, 온도와 속도는 서로 연계된 소스 항을 통해 상호작용한다. 이 구조는 기존의 단일 분포함수(LB) 방식에 비해 메모리 사용량은 약 30 % 증가하지만, 수렴 속도와 정확도는 현저히 향상된다.

세 가지 경계조건 설정은 실제 공학적 응용을 고려한 설계이다. (1) 네 면이 절연인 경우는 전통적인 정육면체 셀에 가장 많이 사용되는 베이스라인이며, (2) 네 면이 선형 온도구배를 갖는 경우는 온도 구배가 복합적으로 작용하는 전자·열 장치에 적용 가능하고, (3) 전·후면이 선형 구배, 상·하면이 절연인 경우는 복합적인 3차원 대류 패턴을 생성한다. 레이리 수를 10³에서 10⁶까지 변화시켜 저레일리에서 난류 전이까지의 전이 현상을 포착하였다.

수치 결과는 기존의 실험 데이터(예: De Vahl Davis, 1983)와 고정밀 유한체적(FVM) 혹은 스펙트럴 방법과의 비교에서 평균 절대오차가 1 % 이하임을 보여준다. 특히 3차원 흐름 구조를 시각화한 등온선·속도벡터는 전·후면 온도구배가 유도하는 복합 대류 셀을 명확히 드러내며, 이는 단일면 구배 모델에서는 관찰되지 않는 현상이다. 또한, MRT의 자유로운 완화계수 선택을 통해 수치 확산을 최소화하고, 고레일리 수(≥10⁶)에서도 수렴성을 유지하였다.

이 논문은 두 가지 중요한 시사점을 제공한다. 첫째, 이중 MRT 모델은 복합적인 열·유동 경계조건을 가진 3차원 자연대류 문제를 효율적으로 해결할 수 있는 강력한 수치 도구임을 입증한다. 둘째, MRT의 파라미터 튜닝이 고레일리 수 영역에서의 수치 안정성을 크게 향상시키며, 기존 LBM의 적용 범위를 기존 2차원·저레일리 한계에서 탈피하게 만든다. 향후 연구에서는 비정상적인 물성(예: 온도 의존성 점성)이나 복합 다상 흐름에 대한 확장이 기대된다.