Navier-Stokes-Cahn-Hilliard system in a $3$D perforated domain with free slip and source term: Existence and homogenization

논문 개요

이 논문은 3 차원 주기적 다공성 매질 (perforated porous medium) 내에서 이종 비압축성 혼합물의 확산 - 인터페이스 (diffuse-interface) 모델을 연구합니다. 구체적으로, 시간 의존적인 Navier-Stokes-Cahn-Hilliard (NSCH) 시스템을 다루며, 이는 점성, 모세관력, 상 분리 메커니즘, 그리고 비보존적 소스 항 (source term) 을 포함합니다. 연구의 핵심은 미시적 스케일 ($\varepsilon$) 에서의 약해 (weak solution) 존재성을 증명하고, $\varepsilon \to 0$ 극한에서의 주기적 동질화 (periodic homogenization) 를 수행하여 유효 거시적 모델을 도출하는 것입니다.

1. 문제 설정 (Problem Setting)

• 물리적 모델:
  • 영역: 주기적으로 배열된 고체 장애물 ($\Omega_s^\varepsilon$) 로 인해 구멍이 뚫린 유체 영역 ($\Omega_p^\varepsilon \subset \mathbb{R}^3$) 에서 정의됩니다.
  • 방정식:
    1. 운동량 방정식 (Navier-Stokes): 비선형 대류 항, 가변 점성 텐서 ($A^\varepsilon$), Korteweg-type 모세관력 항 ($\lambda_\varepsilon \phi^\varepsilon \nabla \mu^\varepsilon$), 외부 힘 ($g^\varepsilon$) 을 포함합니다.
    2. 연속 방정식: 비압축성 조건 ($\nabla \cdot u^\varepsilon = 0$).
    3. Cahn-Hilliard 방정식: 위상 필드 ($\phi^\varepsilon$) 의 진화를 기술하며, 화학적 퍼텐셜 ($\mu^\varepsilon$), 대류 항 ($u^\varepsilon \cdot \nabla \phi^\varepsilon$), 그리고 비보존적 소스 항 ($G(\phi^\varepsilon)$) 을 포함합니다.
  • 경계 조건:
    • 외곽 경계 ($\partial \Omega$): 고정된 No-slip 조건.
    • 고체 장애물 표면 ($\Gamma_s^\varepsilon$): Free-slip (Navier slip with zero tangential stress) 조건. 즉, 법선 방향 속도는 0 이지만 접선 방향 전단 응력은 0 입니다.
  • 모수: 모세관력 강도 $\lambda_\varepsilon$는 $\varepsilon$에 의존하며, $\varepsilon \to 0$일 때 $\lambda$로 수렴합니다. 외부 힘은 $\sqrt{\lambda_\varepsilon}$에 비례하여 스케일링됩니다.

2. 방법론 (Methodology)

연구는 크게 두 단계로 진행됩니다.

A. 미시적 문제의 존재성과 균일 추정 (Existence and A Priori Estimates)

• 약해 정의: 적절한 함수 공간 (Velocity: $U_\varepsilon$, Order Parameter: $C_\varepsilon$, Chemical Potential: $W_\varepsilon$) 에서의 약해 정의를 제시합니다.
• Bogovski˘ı 및 확장 연산자: 다공성 영역 ($\Omega_p^\varepsilon$) 과 고체 영역 ($\Omega_s^\varepsilon$) 에서 발산 연산자의 역 (Bogovski˘ı 연산자) 을 구성하여 압력을 처리하고, 경계 조건을 만족하는 속도장의 확장을 위해 필요한 도구들을 개발합니다.
• 에너지 추정: 에너지 소산 법칙을 유도하여 미시적 해에 대한 균일한 ( $\varepsilon$ 및 $\lambda_\varepsilon$에 무관한) $L^2$ 및 $H^1$ 추정을 증명합니다.
  • 중요한 점: 소스 항 $G(\phi^\varepsilon)$로 인해 질량 보존이 깨지므로, 농도의 공간 평균이 시간에 따라 변합니다. 이로 인해 기존의 Cahn-Hilliard 이론과 달리 Poincaré-Wirtinger 부등식을 직접 적용할 수 없으며, 평균값에 대한 별도의 추정이 필요합니다.
• 해의 존재성: Galerkin 근사 기법 (Stokes-type 고유함수를 기저로 사용) 을 통해 미시적 NSCH 시스템의 약해 존재성을 증명합니다.

B. 동질화 (Homogenization)

• 주기적 펼침 연산자 (Periodic Unfolding Operator): $T_\varepsilon^*$를 사용하여 미시적 구조를 거시적 변수와 미시적 변수로 분리합니다.
• 2-스케일 수렴 (Two-scale Convergence): 미시적 해의 확장을 통해 거시적 한계장 (macroscopic limit fields) 과 2-스케일 한계장을 도출합니다.
• 모세관력 강도 ($\lambda$) 에 따른 경우 분석:
  1. $\lambda = 0$ (소멸하는 모세관력 regime): 모세관력이 너무 약하여 거시적 대류가 발생하지 않는 경우.
  2. $\lambda \in (0, +\infty)$ (균형 regime): 모세관력과 점성/관성이 균형을 이루는 경우.

3. 주요 결과 (Key Results)

A. 미시적 해의 존재성

• 주기적 다공성 매질에서 가변 점성, 소스 항, 그리고 자유 미끄럼 (free-slip) 경계 조건을 가진 NSCH 시스템의 약해 존재성을 최초로 증명했습니다.
• 질량 비보존 소스 항이 있더라도 에너지 소산 법칙을 만족하며 해가 존재함을 보였습니다.

B. 동질화된 유효 모델 (Homogenized Effective Models)

$\lambda_\varepsilon \to \lambda$의 극한 값에 따라 두 가지 서로 다른 유효 모델이 도출됩니다.

1. Case 1: $\lambda = 0$ (Stokes-Cahn-Hilliard 시스템)

   • 특징: 거시적 대류 항 $(u \cdot \nabla)u$와 위상 필드의 대류 항 $u \cdot \nabla \phi$가 모두 사라집니다.
   • 방정식: Stokes-type 운동량 방정식과 Cahn-Hilliard 방정식이 결합된 형태입니다.
   • 물리적 의미: 모세관력 구동력이 약하여 거시적 흐름을 생성하지 못하므로, 유동은 매우 느린 크리핑 흐름 (creeping flow) 으로 근사됩니다.

2. Case 2: $\lambda \in (0, +\infty)$ (Navier-Stokes-Cahn-Hilliard 시스템)

   • 특징: 거시적 대류 항과 위상 필드 대류 항이 $\sqrt{\lambda}$ 계수와 함께 유지됩니다.
   • 방정식: 완전한 Navier-Stokes-Cahn-Hilliard 시스템으로, 비선형 대류와 인터페이스 운동에 의한 유동이 거시적 스케일에서 발생합니다.
   • 물리적 의미: 모세관력이 점성 및 인터페이스 운동과 경쟁하여 유의미한 거시적 유동을 유도합니다.

C. 유효 텐서 및 에너지 수렴

• 유효 텐서: 점성 ($A^{hom}$), 확산 ($B^{hom}$), 모세관력 ($C^{hom}$) 에 대한 유효 텐서가 세포 문제 (cell problems) 를 통해 정의되며, 이들은 거시적 좌표에 무관하지만 시간에 의존할 수 있습니다.
• 에너지 수렴: 미시적 총 에너지가 동질화된 에너지 함수수로 수렴하며, 이는 미시적 시스템과 유사한 에너지 소산 법칙을 만족함을 증명했습니다.

4. 공헌 및 의의 (Significance)

1. 이론적 확장: 기존 연구들이 주로 고정된 영역이나 No-slip 경계 조건, 질량 보존을 가정했던 것과 달리, 다공성 매질, Free-slip 경계 조건, 그리고 비보존적 소스 항을 모두 포함한 NSCH 시스템에 대한 엄밀한 수학적 분석을 제공했습니다.
2. 경계 조건의 중요성: 고체 장애물 표면에서의 Free-slip 조건이 동질화 과정에서 미시적 유동과 거시적 유동의 상호작용에 어떻게 영향을 미치는지를 규명했습니다.
3. 모세관력의 역할 규명: 모세관력 강도 $\lambda$가 거시적 모델의 구조 (Stokes vs. Navier-Stokes) 를 결정하는 핵심 파라미터임을 보여주었습니다. 이는 다공성 매질 내에서의 상 분리 흐름을 모델링할 때 모세관력의 스케일링이 얼마나 중요한지를 시사합니다.
4. 응용 가능성: 지하수 정화, 석유 회수, 생체 조직 공학 등 다공성 매질에서의 다상 유동 현상을 더 정확하게 모델링하는 데 이론적 기반을 제공합니다.

결론

이 논문은 복잡한 다공성 환경에서의 상 분리 유체 흐름을 기술하는 NSCH 시스템에 대해, 해의 존재성을 rigorously 증명하고, 모세관력의 강도에 따라 두 가지 다른 거시적 거동 (Stokes-type vs. Navier-Stokes-type) 으로 수렴함을 보이는 획기적인 연구입니다. 특히 비보존적 소스 항과 Free-slip 경계 조건을 동시에 고려한 점은 기존 동질화 이론의 한계를 넘어서는 중요한 기여로 평가됩니다.