Homogenization of an optimal control problem for nonlocal semilinear elasticity with soft inclusions

技術要約：ソフトインクルージョンを伴う非局所半線形弾性における最適制御問題の均質化

1. 問題定式化
本論文は、ソフトな周期的なインクルージョン（包含物）を含む高コントラスト弾性媒体上で定義された最適制御問題（OCP）の漸近解析を調査している。物理的領域 $\Omega \subset \mathbb{R}^N$ は、連結した硬いマトリックス相（$\Omega^2_\varepsilon$）とソフトなインクルージョン相（$\Omega^1_\varepsilon$）に分解される。系は、非局所項を持つ半線形弾性方程式によって支配されており、これは変形振幅のグローバルな尺度（例：溶媒の再分配や平均場フィードバックによるもの）に変形への機械的応答が依存する材料をモデル化している。

変位 $u_{\varepsilon\delta}$ に関する支配的な状態方程式は以下の通りである：
$$-\nabla \cdot \sigma_{\varepsilon\delta}(u_{\varepsilon\delta}) + \alpha \|u_{\varepsilon\delta}\|_{L^p(\Omega)}^p u_{\varepsilon\delta} = f(x) + \mathbb{1}_{\Omega^1_\varepsilon} \theta_\varepsilon \quad \text{in } \Omega,$$
ここで、境界条件としてディリクレおよびノイマン条件が課される。ここで：

• $\varepsilon > 0$ は周期スケールを表す。
• $\delta > 0$ はコントラストパラメータであり、ソフト相における剛性は硬い相に対して $\delta^2$ のオーダーでスケールする。
• $\theta_\varepsilon \in L^2(\Omega^1_\varepsilon)^N$ は、ソフトインクルージョン内でのみ作用する分布制御である。
• 非局所項 $\alpha \|u\|_{L^p}^p u$ （ただし $p \ge 2$）は、変位のグローバルな $L^2$ 振幅を通じて、すべての空間点を結合させる。

目的は、変位の $L^2$ トラッキング誤差を目標プロファイル $u_d$ とバランスさせ、制御エネルギーに対する正則化項を組み合わせた、簡約されたコスト汎関数 $j_{\varepsilon\delta}(\theta_\varepsilon)$ を最小化することである：
$$j_{\varepsilon\delta}(\theta_\varepsilon) = \frac{1}{2} \int_\Omega |u_{\varepsilon\delta}(\theta_\varepsilon) - u_d|^2 dx + \frac{\gamma}{2} \int_{\Omega^1_\varepsilon} |\theta_\varepsilon|^2 dx.$$

解析は、スケーリング・レジーム $\kappa = \lim \delta/\varepsilon \in (0, +\infty]$ の下での $(\varepsilon, \delta) \to (0, 0)$ の極限に焦点を当てている。

2. 手法
著者らは、極限系を導出するために、**周期展開演算子（periodic unfolding operator）と$\Gamma$-収束（$\Gamma$-convergence）**の議論を組み合わせて用いている。

• 展開演算子： 振動する領域と係数を扱うために、周期展開演算子 $T_\varepsilon$ が利用されている。この手法は、積分を $\Omega$ 上の積分から $\Omega \times Y$ 上の積分（ここで $Y$ は単位周期セル）へと変換することで、マクロ的およびミクロ的スケールの分離を可能にする。
• エネルギー分解と先験的評価： 高コントラストの性質（ソフト相における剛性の退化）のため、標準的な $H^1$ 有界性は成立しない。著者らは、変位場をマクロ的部分とミクロ的な補正項の部分へと分解する特定の分解法を導入している（補題 3.2–3.3）。これにより、コンパクトネスに必要な一様な評価の導出が可能となる。
• 非局所項による強収束： 非局所半線形項 $\alpha \|u\|_{L^p}^p u$ において極限をとることは、決定的な方法論的課題である。この項に対しては弱収束では不十分である。著者らは、全エネルギーの収束（$\Gamma$-収束経由）およびエネルギー汎関数の厳密凸性を利用することで、展開された状態変数の強 $L^2$ 収束を確立している。
• 随伴系解析： 最適制御を特徴付けるために、著者らは随伴状態方程式を導出する。彼らは一意な随伴状態の存在を証明し、その漸近挙動を確立している。これは、極限の最適性条件を導出するために不可欠である。

3. 主要な貢献と結果

• 均質化された状態系の導出：
  本論文は、スケーリング $\delta/\varepsilon \to \kappa$ に関連する極限（均質化）状態系を導出している。極限系は、以下の二スケール構造を持つ：

  • マクロな変位 $u_0$。
  • ソフト相に関連するミクロな補正項 $\tilde{W}$、および硬い相に関連する $\hat{U}$。
  • 極限エネルギー汎関数 $T_\kappa$ は、パラメータ $\kappa$ に明示的に依存する。
  • $\kappa = +\infty$ の特定の場合、系は完全にスケール分離された均質化システムへとデカップル（分離）される。
  • 極限問題は、均質化された剛性テンソル $A^{hom}$ と、結合された極限場 $u_0 + \frac{1}{\kappa}\tilde{W}$ に作用する非局所項を含む変分不等式によって特徴付けられる。

• 最適制御問題の均質化：
  著者らは、ミクロな最適制御の列 $\Theta_\varepsilon$ が、極限問題に対する最適制御 $\hat{\Theta}$ に収束することを証明している。

  • 存在と一意性： ミクロおよび極限の両方の問題における最適制御の存在が確立されている。極限問題については特定の条件下で一意性が保証されるが、非線形性に起因する簡約コスト汎関数の非凸性により、ミクロ問題における大域的一意性は一般には保証されない。
  • 収束： 簡約コスト汎関数の $\Gamma$-収束および状態変数と随伴変数の強収束を用いることで、ミクロな最適制御の列の任意の極限点は、極限 OCP の最適制御になることが証明されている。
  • 極限最適性条件： 極限の最適制御は、極限の随伴状態を介して明示的に特徴付けられる：$\hat{\Theta} = -\frac{1}{\gamma} (v_0 + \frac{1}{\kappa}\hat{Z})$。ここで、$v_0$ と $\hat{Z}$ は極限の随伴状態の成分である。

• 高コントラストと非局所性の処理：
  本研究は、高コントラスト弾性（ソフト相においてエネルギーが退化する）と非局所半線形性が同時に存在する困難に対処している。著者らは、非局所項が強収束を必要とすることを実証しており、これはエネルギー収束の議論を通じて達成される。これは、弱収束でしばしば十分である標準的な局所半線形均質化の解析とは異なる点である。

4. 重要性と新規性
本論文は、主に非局所的な半線形構造と高コントラストの設定において重要性を主張している。

• 非局所演算子の構造： 局所項 $S(u)$ を用いる従来の半線形方程式を用いた OCP の研究とは異なり、本モデルはヒルベルト空間演算子 $u \mapsto \alpha \|u\|_{L^p}^p u$ を採用している。これは、すべての空間点をグローバルに結合させる。著者らは、この構造が弱収束ではなく強 $L^2$ 収束（ノルムの収束）を必要とすることを指摘しており、これは一様な $H^1$ 評価が通常成立しない高コントラスト設定においては非自明な要求である。
• 高コントラスト弾性： 本解析は、ソフトインクルージョンを伴う弾性における最適制御問題へと均質化理論を拡張している。これは、弾性エネルギーが退化する領域である。本論文は、パラメータ $\kappa$ が有効モデルにどのように影響するかを示すことで、このレジームにおける極限系の導出のための厳密な枠組みを提供している。
• 方法論的な厳密さ： 本研究は、制御から状態への写像とコスト汎関数の収束を扱うために、間接法（最適性システム）と直接的な $\Gamma$-収束の議論を組み合わせ、最適制御の極限が確かに極限 OCP に対して最適であることを保証している。

本論文は、結果を、状態および制御の存在証明、先験的評価、二スケール極限系の導出、そして最終的な OCP の収束へと整理することで、このような複雑なメタマテリアルにおける形状／応答最適化のための完全な数学的枠組みを提供して締めくくられている。