A quasi-incompressible Cahn-Hilliard-Darcy model for two immiscible fluids in porous media

기술 요약: 다공성 매질 내 두 가지 비혼합 유체를 위한 준-비압축성 Cahn–Hilliard–Darcy 모델

문제 정의
본 연구는 다공성 매질 내 두 가지 비혼합 유체의 이상 유동(two-phase flows)에 관한 수학적 모델링을 다룬다. 고전적인 무스카트(Muskat) 문제는 유체 계면을 무한히 얇은 불연속면으로 효과적으로 묘사하지만, 매끄러운 프로파일의 무제한적 성장이나 겉보기에 안정적인 영역에서도 발생하는 불안정성과 같은 유한 시간 특이점(finite-time singularities)에 취약하다. 이러한 한계를 극복하기 위해, 급격한 불연속을 좁은 전이층으로 대체하여 역학을 정규화하는 확산 계면(diffuse-interface, phase-field) 정식화가 채택된다.

기존 문헌은 밀도가 일치하거나 일치하지 않는 경우를 포함한 Cahn–Hilliard–Navier–Stokes (CHNS) 시스템과 Cahn–Hilliard–Hele–Shaw/Darcy (CHHS/CHD) 모델을 광범위하게 다루어 왔다. 그러나 특이 로그형 플로리–허긴스(Flory–Huggins) 퍼텐셜과 비균일한 유체 특성(구체적으로 서로 다른 밀도)을 통합하는 준-비압축성(quasi-incompressible) Cahn–Hilliard–Darcy (qCHD) 모델에 대한 완전한 이론적 프레임워크는 여전히 미해결 과제로 남아 있다. 대부분의 기존 Darcy 유형 모델은 엄격한 비압축성($\nabla \cdot u = 0$)을 가정하거나, 상 경계 근처에서의 물리적 실재성을 결여한 다항식 자유 에너지 퍼텐셜을 사용한다.

방법론
저자들은 qCHD 모델을 유도하고 다음의 방법론적 단계를 통해 약해(weak solution)의 전역적 존재성을 확립한다:

1. 변분 유도 (Variational Derivation):

   • 모델은 온사거의 극치 원리(Onsager's extremum principle)(Rayleigh의 최소 에너지 소산 원리와 동등)를 사용하여 유도된다.
   • 시스템은 부피 분율 차이 $\phi$를 질서 매개변수(order parameter)로 사용하며, 밀도 차이로 인해 비용적(non-solenoidal, $\nabla \cdot u \neq 0$)인 질량 평균 혼합 속도 $u$를 사용한다.
   • 총 에너지 범함수 $E(\phi)$는 계면 에너지를 나타내는 구배 항, Flory–Huggins 균질 자유 에너지 $F(\phi)$, 그리고 중력 퍼텐셜 항을 포함한다.
   • 연속체 제약을 완화하기 위해 라그랑주 승수 $p$를 도입하고 Rayleighian을 최소화함으로써, 저자들은 속도, 압력, 화학적 퍼텐셜 및 중력을 지배하는 결합된 방정식계를 유도한다.

2. 모델 정식화 (Model Formulation):

   • 도출된 qCHD 시스템은 다음으로 구성된다:
     • 속도, 압력, 화학적 퍼텐셜 및 중력의 관계를 나타내는 수정된 Darcy 방정식.
     • 혼합물 밀도에 대한 연속 방정식.
     • Atwood 수 $\alpha$에 의해 구동되는 질서 매개변수 $\phi$에 대한 Cahn–Hilliard 유형 방정식 ($\mu_p = \mu + \alpha p$).
   • 자유 에너지 밀도 $F(\phi)$는 다음과 같은 특이 로그 함수이다: $F(\phi) = \frac{\theta_0}{2}[(1+\phi)\ln(1+\phi) + (1-\phi)\ln(1-\phi)] - \frac{\theta_c}{2}\phi^2 + C$.

3. 존재성을 위한 해석적 프레임워크 (Analytical Framework for Existence):

   • 시간 이산화 (Time-Discretization): 자유 에너지 퍼텐셜의 볼록-오목 분할($F = F_v + F_n$)을 사용하여 시간 진행 스킴을 구축한다.
   • 고정점 논법 (Fixed Point Argument): 관련 연산자의 엄격한 단조성(strict monotonicity)과 응집성(coercivity)에 의존하여 시간 이산 시스템의 해의 존재성을 Browder-Minty 렘마를 통해 증명한다.
   • 사전 추정 (A Priori Estimates): 에너지 소산 법칙과 Darcy 방정식으로부터 유도된 압력에 대한 특정 추정을 활용하여 근사해에 대한 균등 추정치를 도출한다.
   • 컴팩트성 (Compactness): Aubin-Lions-Simon 렘마를 사용하여 질서 매개변수의 강한 수렴을 확립함으로써, 특히 특이 로그 퍼텐셜을 포함한 비선형 항에서 극한으로의 통과를 가능하게 한다.

주요 기여 및 결과

• qCHD 모델의 유도: 본 논문은 국부 질량 보존과 에너지 소산 법칙을 만족하는 준-비압축성 Cahn–Hilliard–Darcy 모델의 체계적인 유도를 제시한다. 이 모델은 Atwood 수를 통해 밀도 차이를 고려하며 질량 평균 속도를 사용한다.
• 약해의 전역적 존재성: 주요 수학적 결과는 2D 및 3D 영역 모두에서 q-CHD 시스템의 약해의 전역적 존재성을 증명하는 것이다. 이는 물리적으로 타당한 Flory–Huggins 로그 퍼텐셜에 대해 달성되었으며, 이를 통해 질서 매개변수가 물리적 범위 $(-1, 1)$ 내에 엄격하게 유지됨을 보장한다.
• 급격 계면 극한 (Sharp Interface Limit): 형식적 일치 점근 전개(formal matched asymptotic expansion)를 통해, 저자들은 qCHD 모델이 급격 계면 극한에서 수정된 Muskat 문제로 수렴함을 입증하여 모델이 고전적인 계면 역학과 일치함을 확인한다.
• 에너지 안정성: 도출된 시스템은 엄격한 에너지 소산 법칙을 만족하며, 이는 점성 소산 및 화학적 확산에 의해 총 에너지가 시간이 지남에 따라 감소함을 보장한다.

의의
본 논문은 특이 로그 퍼텐셜과 비균일한 유체 특성(불일치 밀도)을 동시에 통합하는 준-비압축성 Cahn–Hilliard–Darcy 시스템을 제시하는 첫 번째 연구라고 주장한다.

이러한 결과의 의의는 그러한 퍼텐셜을 가진 준-비압축성 Cahn–Hilliard–Navier–Stokes 방정식에 대한 전역적 약해 존재성이 현재 알려져 있지 않다는 점에 있다. Darcy 유사 모델에 대한 적절성을 확립함으로써, 본 연구는 다공성 매질 내 이상 유동의 수학적 이론에서 중요한 공백을 메우며, 불포화 토양 흐름에서 지질학적 CO2 저장 및 연료 전지 수 관리에 이르는 응용 분야를 위한 견고한 프레임워크를 제공한다. 로그 퍼텐셜의 사용은 모델이 순수 상(pure phases) 근처에서도 물리적으로 유효하게 유지되도록 하여, 다항식 근사에서 발생할 수 있는 비물리적 행동을 방지한다.