Global weak solutions to the Cahn-Hilliard equation with degenerate mobility and singular diffusion

技術要約：退化型移動度および特異拡散を伴うカーン・ヒリアード方程式の全域的弱解

問題設定
本論文は、ポリマーブレンドにおける相分離をモデル化する、フローリー・ハギンズ・ド・ジェンヌ自由エネルギーに従うカーン・ヒリアード系（Cahn-Hilliard system）の初期境界値問題を扱っている。系は、有界で滑らかな凸領域 $\Omega \subset \mathbb{R}^3$ 上の時間区間 $(0, T)$ で定義される。支配方程式は以下の通りである：
$$\partial_t u = \text{div}(b(u)\nabla\mu), \quad \mu = -a(u)\Delta u - \frac{a'(u)}{2}|\nabla u|^2 + f(u) - \theta_0 u,$$
ここで、斉次ノイマン境界条件および初期データ $u_0$ を伴う。

調査対象となっている特定の物理的レジームは以下の通りである：

1. 退化型移動度（Degenerate Mobility）: $b(u) = 1 - u^2$ であり、純相 $u = \pm 1$ において消失する。
2. 特異拡散（Singular Diffusion）: $a(u) = \frac{2}{1-u^2}$ であり、$u = \pm 1$ において特異となる。
3. 特異ポテンシャル（Singular Potential）: 自由エネルギー密度は、$u = \pm 1$ で特異となるフローリー・ハギンズ・エントロピー $F(u)$ を含む。

主な数学的課題は、退化型移動度と特異拡散が同時に存在することにある。移動度が非退化であるか、あるいは拡散が一定である場合には存在結果が存在するが、本研究のようなフローリー・ハギンズ・ド・ジェンヌ構造を持つ特定のケースにおいて、これら両方の退化と特異性が3次元で共存する場合、高度な解析的困難が生じる。標準的なエントロピーに基づくアプローチは、通常 $L^2(0, T; H^2(\Omega))$ 推定値を導出するが、化学ポテンシャル $\mu$ が標準的なルベーグ空間でウェルディファインド（well-defined）であることを保証せずに特異項を扱うには不十分である。

手法
著者らは、正則化アプローチと、新しい変数変換および新しい楕円推定法を組み合わせて用いている。

1. 近似スキーム: 退化型移動度 $b(u)$ と特異係数 $a(u)$ を、厳密に正の列 $b_\varepsilon(u) = 1 - u^2 + \varepsilon$ および $a_\varepsilon(u) = 2/b_\varepsilon(u)$ によって近似する。決定的なのは、この近似が係数間の基本的な代数的関係を保持している点である：$a_\varepsilon(u)b_\varepsilon(u) = 2$ および $b_\varepsilon(u)f'(u) = \theta$。ポテンシャル $\Psi$ は元の特異な形式のまま維持される。
2. 近似解の存在: 非退化型移動度および正則な拡散（Giorginiら [13] による結果）を用いることで、正則化された系の全域的弱解の存在を確立する。これらの解は、ほとんど至る所（a.e.）で $|u_\varepsilon| < 1$ を満たす。
3. 変数変換: 特異拡散を扱うために、$\phi(u) = \arcsin(u)$ という変換を導入する。この選択は、恒等式 $\frac{a(u)}{2}|\nabla u|^2 = |\nabla \phi(u)|^2$ によって動機付けられている。正則化されたバージョンは $\phi_\varepsilon(u) = \arcsin(u/\sqrt{1+\varepsilon})$ である。
4. 主要な楕楽推定（Key Elliptic Estimate）: $\phi_\varepsilon(u_\varepsilon)$ および $u_\varepsilon$ に関する $\phi_\varepsilon(u_\varepsilon)$ と $u_\varepsilon$ の $L^4(0, T; H^2(\Omega))$ における一様有界性の導出が核心である。従来の $L^2$ 時間可積分性を得るエントロピー推定とは異なり、著者らはより強い $L^4$ 推定値を導出している。これは以下によって達成される：
   • 方程式を $-\Delta u_\varepsilon$ でテストする。
   • $\phi_\varepsilon$ の特定の構造に適合するように設計された補助関数 $h_\varepsilon$ を用いた、ダル・パッソ、ガルケ、およびグルン（Dal Passo, Garcke, and Grun）による積分部分変数化公式を利用する。
   • 境界項の非負性を保証するために、領域 $\Omega$ の凸性を活用する。
   • これらをエネルギー不等式と組み合わせ、特異項を制御する。
5. 極限への移行: コンパクト性議論（オービン・リオンスの補題）を用いて、極限関数 $u$ を抽出する。重要なステップは、$\varepsilon \int f(u_\varepsilon) \text{div}\eta$ が $\varepsilon \to 0$ のときに消滅することを証明することである。これは、$\|f(u_\varepsilon)\|_{L^2(0,T; L^1(\Omega))}$ が高々 $\varepsilon^{-1/2}$ のオーダーで増大するという特定の評価を確立することで、前方の係数 $\varepsilon$ によって相殺されることを通じて達成される。

主要な貢献と結果

• 定理 1.1（全域的弱解の存在）: 本論文は、有限の初期エネルギーを仮定した場合、3次元の有界で滑らかな凸領域における、退化型移動度と特異拡散を伴うカーン・ヒリアード系の全域的弱解の存在を証明する。

  • 正則性: 解 $u$ は $-1 \le u \le 1$ (a.e.) を満たし、$L^\infty(0, T; H^1(\Omega)) \cap L^4(0, T; H^2(\Omega)) \cap C([0, T]; L^2(\Omega))$ に属する。
  • 変換された変数の強化された正則性: 変換された変数 $\phi(u) = \arcsin(u)$ もまた $L^4(0, T; H^2(\Omega))$ に属する。この $L^4$ 時間正則性は、標準的なエントロピー法から得られる $L^2$ 時間正則性に比べて大幅な改善である。
  • 弱定式化: 解は、$b(u)=0$ となる特異性を扱うために分布の意味で定義されたフラックスベクトル $J \in L^2(0, T; L^2(\Omega))$ を含む特定の積分恒等式を満たす。この定式化は、フラックスの定義における明示的な $\mu$ への依存を回避している。
  • エネルギー不等式: 解はエネルギー不等式 $E(u(t)) + \int_0^t \int_\Omega |\tilde{J}|^2 dx ds \le E(u_0)$ を満たす。ここで $\tilde{J}$ はフラックスに関連している。

• 定理 1.3（等価な定式化）: 著者らは、自身の弱解の概念が、非局所系に関する Cancès および Matthes [6] によって導入された定式化と等価であることを示す。彼らは重み付き化学ポテンシャル $q = -2\Delta \phi(u) + \sqrt{b(u)}\Psi'(u)$ を定義し、解が $q$ を含む積分恒等式を満たすことを示す。これにより、局所的な特異モデルを重み付きワッサーシュタイン空間における勾配流構造へと接続している。

意義と主張
古典的なエントロピーベースのアプローチに頼ることなく、解 $u$ と変換された変数 $\phi(u)$ の両方に対して $L^4(0, T; H^2(\Omega))$ 推定値を導出したことが、本論文の主要な新規性であると著者らは主張している。このより強い可積分性は、エネルギー有界性と、移動度および拡散係数の間の特定の代数的関係を十分に活用することによって得られる。

著者らは、領域の凸性が、特異拡散モデルの先行する解析的研究と一致する重要な仮定であることを述べており、これを非凸領域へ拡張することは未解決の問題として残されている。さらに、拡散係数が一定である場合でさえ、一意性と長期的ダイナミクスの特性化は依然として課題であるが、本研究は、3次元における退化型移動度と特異拡散を伴う物理的に関連性の高いケースに対する基礎的な存在論を確立したものである。

開発された手法、特に新しい楕円推定は、より広範な移動度関数（例：$b(u) = (1-u^2)^m$ ただし $m > 1$）や、潜在的に結合されたナビエ・ストークス／カーン・ヒリアード系にも適応可能であることが示唆されているが、これらの拡張は今後の課題として保留されている。