[論文レビュー] A phase-field approach for detecting cavities via a Kohn-Vogelius type functional

Andrea Aspri|arXiv (Cornell University)|Feb 24, 2022

Numerical methods in inverse problems参考文献 65被引用数 3

ひとこと要約

本稿では、境界上の応力と変位の測定値から線形弾性体内の空洞を再構築するための新しいフェーズフィールド法を提案する。この手法は、周囲長ペナルティによる正則化を施したKohn-Vogelius型汎関数に基づく。周囲長のModica-Mortola緩和と、小さな弾性テンソルを用いた仮想材料アプローチを組み合わせることで、形状最適化問題がH¹(Ω)の凸部分集合上の滑らかで微分可能な最適化問題に変換され、反復的プライマルデュアルアルゴリズムにより、ノイズのあるデータに対しても正確な空洞形状に収束する、安定的かつ効率的な数値再構築が可能になる。

ABSTRACT

We deal with the geometrical inverse problem of the shape reconstruction of cavities in a bounded linear isotropic medium by means of boundary data. The problem is addressed from the point of view of optimal control: the goal is to minimize in the class of Lipschitz domains a Kohn-Vogelius type functional with a perimeter regularization term which penalizes the perimeter of the cavity to be reconstructed. To solve numerically the optimization problem, we use a phase-field approach, approximating the perimeter functional with a Modica-Mortola relaxation and modeling the cavity as an inclusion with a very small elastic tensor. We provide a detailed analysis showing the robustness of the algorithm through some numerical experiments.

研究の動機と目的

• 境界上の応力と変位の測定値から線形弾性体内の空洞を同定する、不適切に定式化された逆問題に対処すること。
• 弱い安定性推定値がある中でも、安定性と収束性を保証する、空洞再構築のための頑健な数値アルゴリズムを開発すること。
• フェーズフィールド法を、これまでこの文脈で応用されていなかったKohn-Vogelius型汎関数へと拡張すること。
• フェーズフィールド緩和と周囲長正則化を用いて、形状再構築問題の数値的安定性と微分可能性を確保する定式化を提供すること。

提案手法

• 2つの境界値問題の解のL²エネルギー差として定義されるKohn-Vogelius汎関数を用いて、逆空洞検出問題を制約付き最適化問題として定式化する。
• 不適切な問題を安定化させ、複雑な空洞境界を罰するため、周囲長正則化項を導入する。
• 空洞を表すスカラーのフェーズ変数v ∈ [0,1]を導入することでフェーズフィールドアプローチを適用し、C1 = δC0（δ > 0は小さい）として、空洞を低剛性インクルージョンとしてモデル化する。
• 周囲長の近似にModica-Mortola汎関数を用い、エネルギー項ε|∇v|² + (1/ε)v(1−v)により、空洞境界の周りに幅εの拡散界面を形成する。
• Kohn-Vogeliusエネルギーとフェーズフィールドによる周囲長近似を組み合わせた緩和汎関数Jδ,ε(v)を導出し、H¹(Ω)上でFrechét微分可能であることを保証する。
• 緩和汎関数の一次最適性条件に基づく、反復的プライマルデュアルアクティブセット法を実装し、数値的再構築を実現する。

実験結果

リサーチクエスチョン

• RQ1 フェーズフィールド法を、線形弾性における空洞再構築に適応できるか？
• RQ2 Modica-Mortola周囲長緩和と仮想材料モデリングの組み合わせが、形状再構築のための安定的かつ収束性のある数値スキームをもたらすか？
• RQ3 標準的な不一致汎関数と比較して、ノイズのある境界測定値下での本手法の性能はいかがなものか？
• RQ4 非凸形状に対して、フェーズフィールド緩和が真の空洞の幾何的特徴をどの程度保持できるか？

主な発見

• 周囲長正則化付きKohn-Vogelius汎関数のフェーズフィールド緩和は、H¹(Ω)上でFrechét微分可能な汎関数をもたらし、効率的な数値最適化を可能にする。
• 数値実験の結果、δとεが十分に小さいとき、緩和汎関数Jδ,εの極小値は元のJregの極小値を高い精度で近似する。
• 凸空洞（例：円、楕円、長方形）に対しては、境界測定値に最大6.5％のノイズが加わっても、正確な再構築が達成される。
• 複数の空洞（例：2つの円形インクルージョン）の再構築において、5％のノイズ下でn = 35,318回の反復で収束が観察された。
• 非凸領域の再構築では、Modica-Mortola周囲長近似のため、凸化が生じる傾向にあり、非凸幾何形状に対しては限界があることが示された。
• 不一致汎関数（Jmisfitδ,ε）との比較では、複数インクルージョンの場合はKohn-Vogeliusに基づく手法がわずかに優れた結果を示したが、境界にアーチファクトが生じる傾向にあった。