에너지 보존 DPD를 이용한 레이리 베르트 대류 불안정성 연구

초록

본 논문은 에너지 보존 해석적 입자 동역학(Dissipative Particle Dynamics, DPD) 방법을 활용해 이상적인 DPD 유체의 레이리 베르트(Rayleigh‑Bénard) 대류 현상을 최초로 시뮬레이션한다. 3차원 격자에서 온도 구배를 조절해 전도보다 대류 항이 우세하도록 파라미터를 설정하고, 임계 온도 차이에서 두 개의 대류 셀(두 셀 패턴)이 형성되는 것을 확인하였다. 온도 차이를 더욱 증가시키면 냉각벽 근처에 고밀도 층이 급격히 형성되고, 도메인 길이를 두 배로 늘리면 셀 직경은 변하지 않은 채 네 개의 셀이 나타나 패턴의 규모 불변성을 입증한다. 높이를 변화시켰을 때는 보다 균일한 셀 구조가 관찰되며, 이러한 결과는 희박 기체에 대한 기존 연구와 정성적으로 일치한다.

상세 분석

본 연구는 에너지 보존 DPD(E‑DPD) 모델을 레이리 베르트 대류 문제에 적용함으로써 미세유체역학과 열전달 분야에서 새로운 수치 실험 플랫폼을 제시한다. E‑DPD는 전통적인 DPD가 보존하지 못했던 내부 에너지와 온도 변동을 입자 수준에서 직접 다루며, 입자 간의 보존력(보존력, 점성력, 랜덤력) 외에 열전도와 열교환을 담당하는 에너지 교환 항을 포함한다. 이러한 구조는 연속체 열전달 방정식의 미분 형태를 입자 기반으로 재현할 수 있게 하며, 특히 비평형 상태에서의 온도 구배와 흐름 결합을 자연스럽게 포착한다.

시뮬레이션 설정에서는 3차원 직육면체(길이 Lx, 높이 Lz, 깊이 Ly) 내부에 상하 경계에 각각 고정 온도(뜨거운 하부, 차가운 상부)를 부여하고, 측면은 주기적 경계조건을 적용하였다. 파라미터 선택—입자 밀도, 무작위 힘 계수, 점성 계수, 열전도 계수—은 무차원 레이리 수(Ra)와 프란틀 수(Pr)를 목표값(예: Ra≈10⁴–10⁵, Pr≈1)으로 맞추도록 조정되었다. 특히, 전도 항보다 대류 항이 우세하도록 열전도 계수를 낮추고 점성 계수를 적절히 감소시켜, 전통적인 Boussinesq 근사 하에서의 임계 Ra를 초과하도록 설계하였다.

결과 분석에서는 온도 차이 ΔT가 임계값을 초과할 때 흐름이 정지 상태에서 불안정해져 두 개의 대류 셀이 형성되는 것을 확인하였다. 셀의 중심축은 수직 방향으로 정렬되며, 상승 흐름은 뜨거운 하부에서 시작해 상부로 이동하고, 하강 흐름은 반대 방향으로 진행한다. ΔT를 추가로 증가시키면 셀 내부의 온도와 밀도 구배가 급격히 강화되어, 차가운 상부 벽 근처에 고밀도 층이 두드러지게 나타난다. 이는 희박 기체에서 관찰되는 “밀도 급증층” 현상과 유사하며, 입자 기반 모델이 이러한 비선형 현상을 재현할 수 있음을 시사한다.

도메인 길이 Lx를 두 배로 확장했을 때, 기존 두 셀 구조가 가로 방향으로 복제되어 네 개의 셀이 형성되었으며, 각 셀의 직경은 변하지 않았다. 이는 대류 셀의 특성 길이(scale)가 물리적 파라미터(Ra, Pr)와 직접 연관되고, 도메인 크기에 독립적임을 확인하는 중요한 검증이다. 반면, 높이 Lz를 변화시켰을 때는 셀 형태가 보다 균일하고 대칭적으로 변형되었으며, 셀 내부의 온도 구배가 완만해져 전도와 대류의 비율이 조정되는 모습을 보였다. 이러한 결과는 전통적인 선형 안정성 이론과 비선형 전이 현상의 연속성을 입증한다.

마지막으로, 본 연구는 E‑DPD가 희박 기체 흐름을 모델링하는 기존 DSMC(Direct Simulation Monte Carlo) 방법과 비교했을 때, 입자 수와 시간 스텝이 상대적으로 적음에도 불구하고 정성적인 패턴 재현에 성공했음을 강조한다. 다만, 입자 간 상호작용 파라미터의 물리적 해석과 스케일링, 그리고 고온 고압 조건에서의 정확도 향상은 향후 연구 과제로 남는다.