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

技術要約：周期的に穿孔された領域におけるソース項を伴うカーン・ヒリアード系の定量的ホモジナイゼーション

問題設定
本論文は、非保存的で単調なソース項を持つ、周期的に穿孔された領域におけるカーン・ヒリアード系の定性的および定量的なホモジナイゼーションを調査している。微視的モデルは、有界なリプシッツ領域 $\Omega \subset \mathbb{R}^d$ ($d \in \{2,3\}$) 内の連結な空隙空間 $\Omega^\varepsilon_p$ における相分離を記述しており、これは周期的な固相包含物の配列を $\Omega$ から除去することによって得られる。系は、秩序パラメータ $\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$ の両方に対して均一なノイマン境界条件が課されている。ポテンシャル $F$ は多項式成長および強凸性を満たし（例えば、古典的なダブルウェル・ポテンシャル）、ソース項 $G$ は全域リプシッツ連続であり、単調かつ $G(0)=0$ を満たす。主な課題は、スケール $\varepsilon \to 0$ の極限として有効なマクロ的方程式を導出し、微視的解とホモジナイズされた解の間の誤差に関する鋭い収束率を確立することにある。

手法
解析には、振動する幾何学的構造を扱うために、アンフォールディング演算子 $T^\varepsilon_*$ を用いて穿孔領域に適応させた**周期的なアンフォールディング法（periodic unfolding method）**を用いている。手法は主に以下の2つの段階で進められる：

1. 定性的ホモジナイゼーション:

   • エネルギー見積もり: 微視的な解に対する一様な先験的評価が導出される。重要な構造的観察として、単調なソース項 $G$（$G(0)=0$）はエネルギー恒等式において散逸項に寄与することである。具体的には、化学ポテンシャルの式を $G(\phi^\varepsilon)$ でテストすると、正の勾配項が得られ、これが散逸と競合するのではなく散逸を強化する。これにより、指数関数的な時間因子を伴うグロンウォールの議論を必要とせずに、$\varepsilon$ に依存しない一様な有界性が得られる。
   • 二スケール収束: アンフォールディング演算子とコンパクト性議論（Aubin–Lions–Simon補題）を用いることで、$(\phi^\varepsilon, \mu^\varepsilon)$ がホモジナイズされたペア $(\phi, \mu)$ に収束することを証明する。極限系は、有効拡散テンソル $B_{\text{hom}}$ を持つ、マクロ領域 $\Omega$ 上のカーン・ヒリアード方程式となる。
   • セル問題: 有効テンソル $B_{\text{hom}}$ は、参照空隙セル $Y_p$ 上のスカラ・ノイマン・セル問題によって特徴付けられる。補正項 $\chi_i$ は $H^1_{\text{per},0}(Y_p)$ に属する。

2. 定量的解析（補正器見積もり）:

   • スケール分離演算子: 定量的なレートを得るために、$Q_1$ スケール分離演算子 $Q_\varepsilon$（Griso [29, 30] による）を用いて、一次補正近似 $\Phi^\varepsilon$ および $M^\varepsilon$ を構築する。この演算子は固定されたマクロ領域 $\Omega$ 上で作用し、その後 $\Omega^\varepsilon_p$ に制限される。決定的なのは、このアプローチにより、セル補正項が $H^1_{\text{per}}(Y_p)$ に留まることが可能となり、古典的な形式展開で必要とされる制限的な $W^{1,\infty}$ の仮定を回避できる点である。
   • 変分残差法: 誤差解析は負のノルムの枠組みで行われる。秩序パラメータの誤差 $e^\varepsilon_\phi = \phi^\varepsilon - \Phi^\varepsilon$ はノイマン逆ラプラシアン $N_\varepsilon e^\varepsilon_\phi$ に対してテストされ、化学ポテンシャルの誤差は $e^\varepsilon_\phi$ に対してテストされる。これにより、残差の高い正則性を必要とせずに済む。
   • 整合性見積もり: 証明は、$\partial \Omega$ 付近の不完全なセルの境界層を考慮した変分残差見積もり（補題 4.6）の導出に依拠している。

主要な結果

• 定性的収束 (定理 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}}$ は対称かつ一様に楕円型である。

• 定量的レート (定理 4.7): ホモジナイズされた解に対する追加の $H^2$ 正則性の仮定と、適切に準備された初期データの下で、本論文は $\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)$ で制御される。
  • 化学ポテンシャルの誤差は、双対ノルム $L^2(0,T; H^1(\Omega^\varepsilon_p)')$ で制御される。

• 境界層効果: 有界領域における $\varepsilon^{1/2}$ のレートは、$\partial \Omega$ 付近の不完全なセルの境界層に起因する。境界層が存在しない平坦なトーラス $\mathbb{T}^d$ 上では、レートは $O(\varepsilon)$ に改善される（Remark 4.10）。

意義および貢献
本論文は、穿孔媒体における相場方程式のホモジナイゼーション理論に対して、主に2つの貢献を行っている：

1. 定性的拡張: 一般的な単調非保存ソース項を持つホモジナイズされたカーン・ヒリアード系を厳密に導出し、ソース項の散逸構造が一般的な非保存強制力と比較してエネルギー見積もりを簡素化することを実証した。
2. 定量的改善: 穿孔領域におけるカーン・ヒリアード系の既知の収束レートを、Schmuck–Kalliadasis [49] による $\varepsilon^{1/4}$ から、二次楕円問題における自然なレートである $\varepsilon^{1/2}$ へと向上させた。
   • 手法的な進展: セル補正項に $W^{1,\infty}$ を要求し、マクロ解に $k \geq 4$ なる $C^1(0,T; W^{k,\infty})$ を要求する従来の手法とは異なり、本アプローチはスケール分離演算子 $Q_\varepsilon$ を利用している。これにより、マクロ解に必要な正則性を $H^2$ へと下げ、セル補正項を $H^1_{\text{per}}$ に留めることが可能となった。
   • 変分的一貫性: 整合性見積もりは、セル問題から継承された弱な無フラックス条件を用いて、完全に変分形式で導出されている。これにより、以前のレート低下の原因となっていた境界項の寄与を回避している。

著者らは、現在の変分フレームワーク内では境界層のために $\varepsilon^{1/2}$ のレートが鋭いものであるが、有界領域において最適な $O(\varepsilon)$ のレートを達成するには、Kenig–Lin–Shen [40] のような双対性議論や、Birman–Suslina [7–9] のような演算子見積もりアプローチが必要になるであろうと述べている（これらは今後の課題としている）。また、現在の残差の正則性では必要なテストが正当化されないため、本論文ではこの次数における化学ポテンシャル $\mu^\varepsilon$ の勾配レートについては主張していない。