다공성 매질에서 주기적 불혼합 흐름을 위한 유효 거시 스톡스칸하밀턴 방정식 유도

초록

본 연구는 두 불혼합 비압축성 유체의 상변화 현상을 다공성 매질 내에서 기술하기 위해, 다중 스케일 방법과 드리프트를 포함한 변분 원리를 적용한다. 주기적 흐름과 충분히 큰 페크레 수를 가정하고, Ginzburg‑Landau/칸‑하밀턴 방정식에 대한 새로운 분할 전략을 이용해 유효 거시 상장 모델을 엄밀히 유도한다. 결과적으로 전통적인 대류‑확산 문제에서 나타나는 테일러‑아리스 분산 효과를 포함한 확산‑분산 관계를 얻으며, 다공성 매질에서 인터페이스를 거시적으로 추적할 수 있는 실용적인 수학·수치 프레임워크를 제공한다.

상세 분석

이 논문은 다공성 매질 내에서 두 비압축성 유체가 서로 섞이지 않는(immiscible) 상황을 기술하기 위해, 전통적인 상장 모델인 스톡스‑칸‑하밀턴(또는 Cahn‑Hilliard) 방정식을 기반으로 새로운 거시 모델을 도출한다. 핵심은 ‘다중 스케일 방법 with drift’를 적용하여 미세 구조(주기적 셀)와 거시 흐름을 동시에 고려한다는 점이다. 먼저, 강하게 천공된(strongly perforated) 영역을 작은 매개변수 ε으로 스케일링하고, 물리량을 ε‑전개(series expansion)한다. 여기서 드리프트 항은 평균 흐름 속도와 연계된 페크레 수(Pe)와 직접적인 관계를 가지며, Pe≫1인 경우 대류가 확산을 압도하는 상황을 모델링한다.

저자들은 기존 Ginzburg‑Landau/Cahn‑Hilliard 방정식의 비선형성 때문에 직접적인 평균화가 어려운 점을 인식하고, ‘분할 전략(splitting strategy)’을 도입한다. 이는 화학 퍼텐셜 μ와 순수 상장 변수 φ를 각각 독립적인 문제로 분리한 뒤, 각각에 대해 다중 스케일 전개를 수행하는 방법이다. 이렇게 하면 미세 셀 문제는 선형(또는 준선형) 형태가 되며, 셀 함수(정기 함수)와 보조 함수(cell problem)들을 정의할 수 있다. 특히, 셀 문제는 주기적 경계조건과 무공극(Neumann) 조건을 만족하도록 설정되어, 유효 확산 텐서와 유효 대류 텐서를 명시적으로 계산한다.

유도 과정에서 얻어지는 주요 결과는 두 가지이다. 첫째, 거시 상장 방정식은 기존의 Cahn‑Hilliard 형태를 유지하지만, 확산 계수 D가 단순 상수 대신 ‘확산‑분산 텐서 D_eff’를 갖는다. 이 텐서는 미세 구조와 평균 흐름 속도에 의존하며, 테일러‑아리스(Taylor‑Aris) 분산과 유사한 형태를 보인다. 둘째, 압력-속도 관계는 스톡스 방정식의 형태를 유지하면서, 유효 점성 μ_eff가 미세 구조에 의해 수정된다. 따라서 최종 거시 모델은 스톡스‑칸‑하밀턴 연동 방정식으로, 두 상의 인터페이스는 φ의 등위면(level set)으로 자연스럽게 표현된다.

수학적 엄밀성 측면에서, 저자들은 에너지 법칙과 엔트로피 생산성을 보존하는 변분 구조를 유지함을 증명한다. 또한, 존재와 유일성, 그리고 정규성(regularity) 결과를 기존의 Cahn‑Hilliard 이론에 맞추어 확장한다. 수치적 검증을 위해, 저자들은 2‑차원 주기적 셀 구조와 간단한 채널 흐름을 시뮬레이션하고, 유효 텐서가 직접 계산된 값과 비교하여 높은 일치도를 보인다. 이는 제안된 모델이 실제 다공성 매질(예: 토양, 연료전지 전극)에서의 다상 흐름을 정확히 포착할 수 있음을 시사한다.

전반적으로, 이 연구는 다공성 매질에서의 다상 흐름을 거시적으로 기술하기 위한 체계적인 수학적 프레임워크를 제공한다. 특히, 고페크레 수 regime에서 대류‑확산 상호작용을 정확히 반영한 확산‑분산 관계를 도출함으로써, 기존의 경험적 모델을 대체하거나 보완할 수 있는 강력한 도구를 제시한다.