A Least-Squares Weak Galerkin Method for the Biharmonic Cauchy Problem

技術要約：双調和コーシー問題に対する最小二乗弱ガラーキン法

問題設定
本論文は、板の曲げ、弾性学、および画像処理などの応用分野で発生する4階の偏微分方程式である双調和コーシー問題を扱う。標準的な境界値問題とは異 달리、コーシー問題は、境界のアクセシブルな部分 $\Gamma_1$ にのみ境界条件（関数の値と法線微分）が与えられ、アクセシブルではない部分 $\Gamma_2$ には情報が存在しない領域 $\Omega$ で定義される。

この問題は、アダマールの意味において**極めて不良設定（ill-posed）**である。一意性の継続原理に基づく特定の仮定の下では一意性を確立できる可能性があるものの、解はコーシーデータに対して連続的に依存しない。その結果、境界測定における微小な摂動が、再構成された解に大きな誤差をもたらす可能性がある。この不安定性は、4階の作用素の性質と、部分的なデータから欠落した境界情報を復元する必要性によって増幅される。

手法
著者らは、**最小二乗弱ガラーキン（LS-WG）**有限要素法を提案している。その手法は以下のステップで進行する。

1. 再定式化: 4階の方程式 $\Delta^2 u = f$ を、補助変数 $w = \Delta u$ を導入することで、2つの2階連立方程式として再定式化する：
   $$\begin{aligned} \Delta w &= f, \\ \Delta u - w &= 0, \end{aligned}$$
   ここで、$\Gamma_1$ 上のコーシーデータに従うものとする。

2. 弱ガラーキン・フレームワーク: 本手法は、関数を内部値、境界値、および法線微分の近似を表すトリプレット $v = \{v_0, v_b, v_n\}$ として定義する弱有限要素空間を利用する。これらの成分は独立した未知数として扱われ、要素境界を越えて連続である必要はなく、また内部値のトレースと一致する必要もない。

3. 離散的弱ラプラシアン: 古典的なラプラシアンは、多項式テスト関数に対して積分を行う弱形式を通じて定義される離散的弱ラプラシアン演算子（$\Delta_{w,r,T}$）によって置き換えられる。これにより、不連続有限要素空間の使用が可能になる。

4. 最小二乗定式化: 連立系は最小二乗アプローチを用いて離散化される。この手法は、支配方程式の残差と、弱関数の内部成分と境界成分の間の整合性を最小化する。これにより、対称正定値線形系が得られる。

   • 安定化: 弱関数における内部成分（$v_0$）と境界成分（$v_b, v_n$）の間の整合性を弱く強制するために、安定化項が導入される。
   • 境界条件: $\Gamma_1$ 上のコーシーデータは、境界成分への $L^2$ 射影を通じて強制される。

5. メッシュの柔軟性: この定式化は、一般の多角形メッシュを自然に受け入れることができ、一般のメッシュ上で構築が困難な（アルギリス要素のような）大域的な $C^1$ 適合要素を必要としない。

主な貢献
論文では、以下の具体的な貢献を主張している。

• 双調和コーシー問題に対する初の LS-WG: 双調和コーシー問題のために特別に設計された最初の最小二乗弱ガラーキン有限要素法を提示している。
• 対称正定値系: 提案された定式化は、対称正定値線形系をもたらす。これにより、混合有限要素法で通常必要とされる離散的 inf-sup 条件を回避し、サドルポイント構造を避けることができる。
• $C^1$ 連続性の排除: 弱演算子を用いることで、不連続空間を利用し、$C^1$ 適合な有限要素空間の構築を回避する。
• 理論的解析: 連続的な問題が一意の解を持つという仮定の下で、著者らは離散 LS-WG 解の一意性を確立している。さらに、離散エネルギーノルムにおける最適次数の誤差評価を導出している。
• 効率的なソルバー: 得られる系の対称正定値な性質により、効率的な反復ソルバー（マルチグリッド・プリコンディショナを用いた共役勾配法など）を利用できる。

数値実験
著者らは、3種類のメッシュ（一様正方形格子、三角形格子、非凸多角形格子）を用いて、単位正方形領域における数値実験を通じて手法を検証している。

• 収束性: 数値実験は理論的な収束率を裏付けている。滑らかな解に対して、本手法は様々な多項式次数（$k=2, 3, 4$）において、離散エネルギーノルムおよび $L^2$ ノルムで最適な収束次数を達成する。
• 堅牢性: 結果は、本手法が異なるメッシュタイプに対して効果的であることを示しており、正方形、三角形、および非凸多角形メッシュの間で性能の差は小さい。
• 課題: 自由境界（アクセシブルではない部分）付近で解の変動が大きい場合、数値解の近似がより困難になり、問題固有の不良設定性に起因する大きな条件数により、より細かい格子レベルでマシン精度に達することが論文内で指摘されている。

意義
本論文は、LS-WG 法を、4階の不良設定問題に対する堅牢な代替案として位置づけている。最小二乗定式化の安定性と弱ガラーキン・フレームワークの柔軟性を組み合わせることで、本手法は、双調和コーシー問題を解く際の2つの主要な障壁、すなわち $C^1$ 適合要素の構築の難しさと、混合定式化における inf-sup 条件を満たす複雑さを克服している。一般の多角形メッシュを扱える能力は、従来の適合的手法が苦戦する複雑な幾何学的形状において、本手法を特に有用なものにしている。既存の弱ガラーキン法がほぼ排他的に良設定な境界値問題のために開発されてきたことを踏まえると、本研究は文献上の空白を埋めるものである。