다공성 매질과 유체가 결합된 시스템의 난류 대류·상변화 통합 수치 프레임워크

초록

**
본 논문은 Darcy‑Brinkman 방정식과 수정된 2차 위상장 모델을 기반으로, 유체와 다공성 매질이 동시에 존재하는 영역에서 난류 대류와 고체‑액체 상변화를 일관되게 시뮬레이션할 수 있는 단일 도메인 수치 해석 체계를 제시한다. 열·용질 전달, 열분산, 토러스 효과 등을 포함한 완전한 비정규화 방정식을 도입하고, 연산 효율성을 위한 연산자 분할·투영 방법을 구현하였다. 다양한 검증 사례(채널 난류·다공성 기판, 이중층 열대류, 1‑D 스테판, 2‑D 물 동결, 해수 동결에 의한 무스 아이스 성장)에서 기존 결과와 일치함을 보이며, 복합 현상 연구에 활용 가능한 범용 도구임을 입증한다.

**

상세 분석

**
이 연구는 기존의 두 영역(순수 유체와 Darcy 흐름) 사이에 인위적인 경계조건을 부과하는 방식과 달리, Darcy‑Brinkman 방정식을 전체 영역에 적용함으로써 속도와 압력의 연속성을 자연스럽게 보장한다. 특히, 다공성 매질 내부의 점성항을 유지하면서도 Darcy 저항을 가변 퍼미어빌리티 K(ϕ) 함수로 표현해, ϕ가 0→1으로 변할 때 매끄러운 전이 구역을 구현한다. 이는 무스 레이어와 같은 고체‑액체 혼합 영역을 모델링할 때 필수적인 특성이다.

위상장 방정식은 전통적인 Cahn‑Hilliard 형태 대신 2차형식(대류‑확산 + 상변화 반응항)으로 단순화했으며, 이때 이중우물 잠재력을 생략함으로써 시간 단계 제한을 크게 완화한다. 반응계수 G와 액상 온도 Tϕ (액상 온도 T₀ − λS) 를 통해 용질 농도에 의한 어는점 변화를 직접 연결한다. 이는 해수 동결과 같은 용질‑온도 상호작용이 중요한 현상을 정확히 포착한다.

열전달 모델에서는 고체·액체의 비열·전도도를 ϕ에 대한 선형 가중 평균으로 정의하고, 고속 흐름에서 발생하는 열분산을 ε_disp(ϕ) |U| 형태의 추가 전도항으로 보정한다. 이는 다공성 매질 내부의 미세 흐름이 전체 열전달에 미치는 영향을 정량화한다는 점에서 의미가 크다. 용질 전달은 Archie 법칙을 변형해 토러스 τ = ϕ − γ (γ=1) 로 설정, ϕ가 낮을수록 확산이 크게 억제되는 현실적인 거동을 재현한다.

수치 해법은 시간 적분에 2차 정확도 Runge‑Kutta와 연산자 분할을 적용해, 비선형 대류‑확산·위상장·모멘텀 방정식을 각각 독립적으로 풀면서도 전체 연계성을 유지한다. 압력–속도 결합은 전통적인 투영(프레셔‑코렉션) 방법을 사용해 연속 방정식을 만족시킨다. 격자 독립성 검증과 함께, 채널 난류·다공성 기판 사례에서 평균 속도 프로파일과 전단 응력 분포가 실험·DNS 결과와 일치함을 보였다.

검증 사례는 점진적으로 복잡성을 높이며 프레임워크의 범용성을 확인한다. 1‑D 스테판 문제에서는 고정된 경계와 이동 경계 모두에서 라인타임 정확도가 기존 고정밀 해법에 필적한다. 2‑D 물 동결 시뮬레이션에서는 인터페이스 형성·전파 속도가 실험 영상과 일치하고, 무스 아이스 성장 모델에서는 해수 농도와 온도 프로파일이 실험 데이터와 거의 동일하게 재현된다.

이러한 결과는 (1) 단일 도메인 접근법이 복합 경계조건 없이도 다중 물리 현상을 일관되게 처리할 수 있음을, (2) 수정된 위상장 모델이 높은 고체‑액체 열전도도 대비를 효율적으로 다룰 수 있음을, (3) 열·용질·상변화가 동시에 일어나는 자연·공학 시스템(예: 화산, 금속 주조, 해양 빙하) 연구에 바로 적용 가능함을 시사한다. 다만, 현재 구현은 Boussinesq 근사와 등온·등농도 가정에 의존하므로, 강한 비선형 밀도 변동이나 비등온·비등농도 상황에서는 추가 모델링이 필요할 것이다.

**