Quantitative Homogenization of a Cahn--Hilliard System with Source Term in Periodically Perforated Domains

기술 요약: 주기적으로 천공된 도메인 내 소스 항을 포함한 칸-힐리야드(Cahn–Hilliard) 시스템의 정량적 균질화

문제 정의
본 논문은 주기적으로 천공된 도메인 내에서 비보존적(nonconservative)이고 단조적인(monotone) 소스 항을 특징으로 하는 칸-힐리야드 시스템의 질적 및 정량적 균질화를 조사한다. 미시적 모델은 유체가 연결된 기공 공간 $\Omega^\varepsilon_p \subset \Omega \subset \mathbb{R}^d$ ($d \in \{2,3\}$)를 차지하는 다공성 매질 내에서의 상분리를 기술하며, 이 공간은 유계 리프시츠(Lipschitz) 도메인으로부터 주기적인 고체 포함물(solid inclusions) 배열을 제거하여 얻어진다. 시스템은 다음과 같이 질서 매개변수(order parameter) $\phi^\varepsilon$와 화학 퍼텐셜 $\mu^\varepsilon$를 결합한다:
$$\begin{cases} \partial_t \phi^\varepsilon - \Delta \mu^\varepsilon + G(\phi^\varepsilon) = 0, \\ \mu^\varepsilon = -\Delta \phi^\varepsilon + F'(\phi^\varepsilon), \end{cases}$$
여기서 외부 경계 $\partial \Omega$와 기공 벽 $\Gamma^\varepsilon_s$ 모두에 대해 균일한 노이만(Neumann) 경계 조건을 따른다. 퍼텐셜 $F$는 다항식 성장 및 코어시비티(coercivity)를 만족하며(예: 전형적인 이중 우물 퍼텐셜), 소스 항 $G$는 전역적으로 립시츠 연속이며 단조적이고 $G(0)=0$을 만족한다. 주요 과제는 $\varepsilon \to 0$일 때 유효한 거시적 방정식을 유도하고, 미시적 해와 균질화된 해 사이의 오차에 대한 날카로운 수렴 속도를 확립하는 것이다.

방법론
분석에는 천공된 도메인에 적응된 **주기적 언폴딩 방법(periodic unfolding method)**을 사용하며, 진동하는 기하학적 구조를 처리하기 위해 언폴딩 연산자 $T^\varepsilon_*$를 활용한다. 방법론은 두 가지 주요 단계로 진행된다:

1. 질적 균질화:

   • 에너지 추정(Energy Estimates): 미시적 해에 대한 균등한 사전 추정치(a priori estimates)가 유도된다. 핵심적인 구조적 관찰은 단조 소스 항 $G$ ($G(0)=0$인 경우)가 에너지 항등식에서 소산(dissipative) 항에 기여한다는 점이다. 구체적으로, 화학 퍼텐셜 방정식에 $G(\phi^\varepsilon)$를 테스트하면 양의 그래디언트 항이 생성되어, 이는 소산과 경쟁하는 것이 아니라 소산을 강화한다. 이를 통해 지수적 시간 인자를 가진 그론월(Gronwall) 논법 없이도 $\varepsilon$에 독립적인 균등 유계성을 확보할 수 있다.
   • 이-스케일 수렴(Two-Scale Convergence): 언폴딩 연산자와 컴팩트성 논증(Aubin–Lions–Simon 르마)을 사용하여, $(\phi^\varepsilon, \mu^\varepsilon)$가 균질화된 쌍 $(\phi, \mu)$로 수렴함을 증명한다. 극한 시스템은 유효 확산 텐서 $B_{\text{hom}}$를 갖는 거시적 도메인 $\Omega$ 상의 칸-힐리야드 방정식이다.
   • 셀 문제(Cell Problems): 유효 텐서 $B_{\text{hom}}$는 참조 기공 셀 $Y_p$ 상의 스칼라 노이만 셀 문제에 의해 특징지어진다. 보정기(correctors) $\chi_i$는 $H^1_{\text{per},0}(Y_p)$에 속한다.

2. 정량적 분석 (보정기 추정):

   • 스케일 분할 연산자(Scale-Splitting Operator): 정량적 속도를 달성하기 위해, $Q_1$ 스케일 분할 연산자 $Q_\varepsilon$(Griso [29, 30]의 것)를 사용하여 1차 보정 근사 $\Phi^\varepsilon$와 $M^\varepsilon$를 구성한다. 이 연산자는 고정된 거시적 도메인 $\Omega$ 상에서 작용하며, 이후 $\Omega^\varepsilon_p$로 제한된다. 결정적으로, 이 접근 방식은 셀 보정기가 $H^1_{\text{per}}(Y_p)$에 머물 수 있게 하여, 고전적인 형식적 전개(formal expansions)에서 요구되는 제한적인 $W^{1,\infty}$ 가정을 피할 수 있게 한다.
   • 변분 잔차 방법(Variational Residual Method): 오차 분석은 음의 노름(negative-norm) 프레임워크에서 수행된다. 위상 오차 $e^\varepsilon_\phi = \phi^\varepsilon - \Phi^\varepsilon$는 노이만 역 라플라시안 $N_\varepsilon e^\varepsilon_\phi$에 대해 테스트되며, 화학 퍼텐셜 오차는 $e^\varepsilon_\phi$에 대해 테스트된다. 이는 잔차(residual)의 더 높은 정칙성(regularity)을 요구할 필요를 없애준다.
   • 일관성 추정(Consistency Estimates): 증명은 $\partial \Omega$ 근처의 불완전한 셀에 의한 경계층을 고려하는 변분 잔차 추정을 기반으로 한다.

주요 결과

• 질적 수렴 (정리 3.7): 미시적 해는 다음의 균질화된 시스템으로 (부분 수열로서) 수렴한다:
  $$\begin{cases} \partial_t \phi - \text{div}(B_{\text{hom}}\nabla \mu) + G(\phi) = 0, \\ \mu = -\text{div}(B_{\text{hom}}\nabla \phi) + F'(\phi), \end{cases}$$
  여기서 $B_{\text{hom}}$는 대칭적이고 균등 타원형(uniformly elliptic)이다.

• 정량적 속도 (정리 4.7): 균질화된 해에 대한 추가적인 $H^2$-정칙성 가정과 잘 준비된 초기 데이터(well-prepared initial data) 하에서, 본 논문은 $\varepsilon^{1/2}$ 차수의 오차 추정치를 확립한다:
  $$\|\phi^\varepsilon - \phi\|_{L^2(0,T; L^2(\Omega^\varepsilon_p))} + \|\nabla \phi^\varepsilon - (\nabla \phi + \sum Q_\varepsilon(\partial_{x_i}\phi)\nabla_y \chi_i(\tfrac{x}{\varepsilon}))\|_{L^2((0,T)\times\Omega^\varepsilon_p)} + \|\mu^\varepsilon - \mu\|_{L^2(0,T; H^1(\Omega^\varepsilon_p)')} \leq C\varepsilon^{1/2}.$$

  • 보정되지 않은 질서 매개변수 오차는 $L^2(0,T; L^2(\Omega^\varepsilon_p))$에서 제어된다.
  • 보정된 그래디언트 오차는 $L^2((0,T)\times\Omega^\varepsilon_p)$에서 제어된다.
  • 화학 퍼텐셜 오차는 듀얼 노름(dual norm) $L^2(0,T; H^1(\Omega^\varepsilon_p)')$에서 제어된다.

• 경계층 효과(Boundary Layer Effects): 유계 도메인에서의 $\varepsilon^{1/2}$ 속도는 $\partial \Omega$ 근처의 불완전한 셀에 의한 경계층 때문이다. 불완전한 셀이 존재하지 않는 평평한 토러스 $\mathbb{T}^d$에서는 속도가 $O(\varepsilon)$로 개선된다 (비고 4.10).

의의 및 기여
본 논문은 천공된 매질 내 상-장(phase-field) 방정식의 균질화 이론에 대해 두 가지 주요 기여를 한다:

1. 질적 확장: 일반적인 단조 비보존적 소스 항을 갖는 균질화된 칸-힐리야드 시스템을 엄밀하게 유도하였으며, 소스 항의 소산적 구조가 일반적인 비보존적 강제력(forcing)에 비해 에너지 추정을 단순화함을 보여주었다.
2. 정량적 개선: 천공된 도메인 내 칸-힐리야드 시스템에 대해 기존에 알려진 $\varepsilon^{1/4}$ 수렴 속도(Schmuck–Kalliadasis [49]에 의해 확립됨)를 천공된 도메인의 2계 타원형 문제에 대한 자연스러운 속도인 $\varepsilon^{1/2}$로 개선하였다.
   • 방법론적 발전: $W^{1,\infty}$ 내의 셀 보정기와 $k \geq 4$인 $C^1(0,T; W^{k,\infty})$ 내의 거시적 해를 요구했던 이전 연구들과 달리, 이 접근 방식은 스케일 분할 연산자 $Q_\varepsilon$를 활용한다. 이를 통해 거시적 해에 요구되는 정칙성을 $H^2$로 낮추고, 셀 보정기를 $H^1_{\text{per}}$에 머물게 하였다.
   • 변분 일관성: 일관성 추정은 전적으로 변분 형식으로 유도되었으며, 셀 문제로부터 상속된 약한 노-플럭스(no-flux) 조건을 활용한다. 이는 기존에 속도 저하($\varepsilon^{1/4}$)를 일으켰던 경계항 기여를 피할 수 있게 하였다.

저자들은 현재의 변분 프레임워크 내에서 경계층으로 인해 $\varepsilon^{1/2}$ 속도가 날카롭다고(sharp) 언급하며, 유계 도메인에서 최적의 $O(\varepsilon)$ 속도를 달라는 것은 Kenig–Lin–Shen [40]과 같은 쌍대성 논법(duality argument)이나 Birman–Suslina [7–9]와 같은 연산자 추정법이 필요할 것이며, 이는 향후 연구 과제로 남겨두었다고 명시하였다. 또한, 현재의 잔차 정칙성으로는 필요한 테스팅을 정당화할 수 없으므로, 이 차수에서 화학 퍼텐셜 $\mu^\varepsilon$의 그래디언트 속도를 얻는 데 대한 주장은 하지 않는다.