Fast high-order solvers for the Lippmann--Schwinger equation in piecewise-smooth heterogeneous media

技術要約：区分的に滑らかな不均質媒体におけるリップマン・シュワルツ方程式のための高速高次ソルバー

問題提起
本論文は、透過性のある区分的に滑らかな不均質媒体による時間調和波散乱から生じる、二次元リップマン・シュワルツ方程式の数値解法を扱う。主な課題は、既存の文献ではしばしば相反する関係にある以下の2つの目的を同時に達成することである：

1. 高次の精度： 屈折率（コントラスト）が区分的に滑らかであり、材料界面において跳躍不連続性を伴う場合でも、解を正確に解像すること。
2. 周波数に対して頑健な反復収束： 得られる密な線形系を、波数が増大しても有界なままとなる既約部分空間（Krylov部分空間、GMRES）反復回数で解くこと。

既存の高次手法は通常、一様なデカルト格子に依存しており、界面が格子に適合しない場合に精度が低下する。一方、界面適合型手法（適合メッシュを用いたNyströmスキームなど）は高次の精度を提供するが、周波数に対して頑健な前処理手法を欠いているため、高周波において反復回数が急速に増大する場合が多い。

手法
著者らは、転送作用素を用いて、高次な界面適合型離散化と構造格子前処理を結合させたソルバーを提案している。

1. 非構造離散化（ソルバー）：

   • ドメイン $\Omega$（コントラストの支持集合）は、材料界面に適合した適合三角形メッシュを用いて離散化される。
   • 高次のNyström法が、Vioreanu–Rokhlin (VR) 求積ノードを用いて採用されている。
   • 体積ポテンシャルは、体積密度補間法 (VDIM) を用いて評価される。VDIMは、各三角形上で密度の多項式補間を構成することにより、弱特異カーネルを正則化し、不連続を横断する高次精度の積分を可能にする。
   • 得られる線形系は密であるが、高速多重極展開 (FMM) または H行列圧縮を用いた体積ポテンシャル演算子の行列フリー適用が可能であり、反復あたりの計算量は準線形な複雑さを実現する。

2. 前処理戦略：

   • 収束を加速するために、一様なデカルト格子用に設計された「疎化前処理（sparsifying preconditioner）」（具体的にはYing [46] による手法）を用いて、非構造系を前処理する。
   • 転送作用素： 2つの作用素が、非構造求積ノードと一様なデカルト格子を橋渡しする：
     • デカルトから求積への転送 ($T^Q_C$): 格子値を求積ノードへ写像する双線形補間（prolongation）。
     • 求積からデカルトへの転送 ($T^C_Q$): 本論文では、この制限作用素（restriction operator）の2つの構成を提案し、分析している：
       • アプローチ I: 区分的線形関数への局所的な $L^2$ 射影。
       • アプローチ II: 補間作用素（prolongation operator）の正規化された転置（幾何学的マルチグリッドにおける標準的な手法）。
   • 前処理 $P$ は、$P = T^Q_C S^{(C)} T^C_Q + (I - T^Q_C T^C_Q)$ として構成される。ここで、$S^{(C)}$ はデカルト格子上の疎化前処理である。補正項 $(I - T^Q_C T^C_Q)$ は、埋め込み項のランク不足を補償するために不可欠である。

3. 理論的解析：

   • 論文では、前処理 $P$ の可逆性に関する十分条件を確立している。
   • アプローチ I については、メッシュ比 $r_h = h/h_C$（非構造メッシュサイズとデカルトステップサイズの比）が十分に小さいという条件下で可逆性が証明されている。
   • アプローチ II については、特定の内積におけるスペクトル条件（特定の非正の実軸を避ける条件）を満たす場合、メッシュサイズ $h$ に関わらず、前処理 $S^{(C)}$ が連続的な散乱問題のエネルギーバランスに関連している限り、可逆性が証明されている。

主な貢献

1. 新規ソルバー： 界面において高次の精度を維持しながら、高速な畳み込みベースのソルバーを利用できる、2Dリップマン・シュワルツ方程式のための初の高速高次Nyströmソルバー。
2. ハイブリッド前処理： 転送作用素を介して適合メッシュ離散化と構造格子前処理を結合させる新しい戦略を提供し、周波数に対して頑健なソルバー（疎化前処理など）を非構造メッシュ上で利用可能にした。
3. 可逆性の証明： メッシュ細分化に関する明示的な条件および基礎となるデカルトソルバーの特性の下で、前処理の可逆性に関する厳密な理論的保証を与えた。
4. 数値検証： 提案された前処理系が、前処理なしの場合と比較して大幅に少ないGMRES反復回数を必要とすることを実証する広範な実験。反復回数は、高周波散乱や高コントラストの区分的に滑らかな屈折率に対しても、メッシュサイズおよび波数に対してほぼ独立していることが示されている。

結果と意義
本論文は、高次の幾何学的柔軟性と周波数に対して頑健な反復収束との間の長年の対立を解消することを主張している。数値結果は、提案手法が以下を実現することを示している：

• 跳躍不連続を伴う区分的に滑らかなコントラストに対して、高次の収束率を達成する。
• 波数が増大してもGMRESの反復回数が有界に保たれる「周波数に対して頑健な」挙動を示す。これは、非構造高次手法では通常失われる特性である。
• メッシュ細分化に対する安定性を維持する。

著者らは、本研究を、鋭い材料界面やサブドメイン内で連続的に変化する屈折率を持つ複雑な人工構造物（ナノフォトニクスや集積光学など）の高周波散乱問題を解くための重要な一歩として位置づけている。本手法は、適合メッシュの精度と構造格子のソルバー効率の両方の利点を、どちらかを犠牲にすることなく活用している。