[論文レビュー] A fullwave model of the nonlinear wave equation with multiple relaxations and relaxing perfectly matched layers for high-order numerical finite-difference solutions

Gianmarco Pinton|arXiv (Cornell University)|Jun 22, 2021

Seismic Imaging and Inversion Techniques参考文献 13被引用数 9

ひとこと要約

本稿では、非線形波動方程式に対して空間的に変化する複数の緩和メカニズムを用いて周波数依存の減衰および分散をモデル化する高次有限差分ソルバを提示する。畳み込み緩和を用いて緩和的パーソンマッチドレイヤー（PML）を統合することで、数値分散および減衰を最小限に抑え、長距離超音波シミュレーションを高精度に実現した。2次元および3次元の実験的および理論的ベンチマークで検証され、スペックル輝度の誤差が1 dB未満、7次高調波の精度が保証された。

ABSTRACT

Large-scale simulations of the wave equation in electromagnetism, seismology, and acoustics, can be solved efficiently by finite difference methods. The accuracy of these numerical solutions usually depends on the minimization of dissipation and dispersion error using optimized schemes. However, these optimized schemes are often difficult to implement for more complex formulations of the wave equation. Here we present a model of the quadratically nonlinear wave equation with multiple relaxations as a function of space that can model perfectly matched layers and is designed to be easily used with optimized finite-difference stencils on a staggered grid. Solutions are demonstrated in 2D and 3D for ultrasound imaging applications where the long propagation lengths, heterogeneous media in the human body, large dynamic range, and requirements for computational speed, demonstrate the utility of the proposed methods.

研究の動機と目的

• 数値分散および減衰を最小限に抑えつつ、不均質なヒト組織における長距離・高ダイナミクスレンジ超音波波動伝播をシミュレートする課題に対処する。
• 任意の減衰則および非線形性をサポートする、空間的に変化する緩和メカニズムを用いた柔軟で高次有限差分式を構築する。
• 畳み込み緩和を介して緩和的パーソンマッチドレイヤー（PML）を波動方程式に統合し、反射アーチファクトおよびインピーダンス不整合を最小限に抑える。
• 非線形性、減衰、およびサブリゾリューション散乱体を含む、現実的な物理現象を高精度に再現可能な超音波画像化システムのシミュレーションを可能にする。
• バーガース方程式のRusanovスキームを用いた数値ベンチマークと、ブタの肝臓からの実験データを用いてモデルを検証する。

提案手法

• 圧力-速度変数を用いた線形波動方程式を、空間微分の複素座標変換を用いて、N個の緩和メカニズムにより減衰および分散をモデル化する。
• 軟部組織などの不均質媒体における周波数依存減衰および分散を表現するため、空間的に変化する緩和パラメータ（κ, dν, αν）を実装する。
• 因果性を保証し、解析的減衰則（例：0.5 dB/MHz/cm）に一致させるために、クライマース・クロニッガー関係から導出された畳み込み緩和カーネルを用い、パラメータは最小二乗法によるフィッティングで決定する。
• 境界層に緩和フレームワークを拡張して緩和的PMLを構築し、内部の分散特性と一致させることで、分散が一定であり、反射を最小限に抑える。
• ステaggerドグリッド上に高次最適化有限差分ステンシルを適用し、分散および減衰誤差を最小限に抑え、粗いグリッド上でも効率的かつ安定したシミュレーションを実現する。
• 1次元Rusanov解を用いてバーガース方程式のベンチマークを実施し、1波長あたり15点の周波数および時間領域応答を比較する。

実験結果

リサーチクエスチョン

• RQ1 空間的に変化する複数の緩和メカニズムを有する非線形波動方程式のフルウェーブモデルは、軟部組織における周波数依存減衰および分散を正確にシミュレートできるか？
• RQ2 畳み込み緩和メカニズムに基づく緩和的PMLは、不均質媒体において数値反射をどれほど効果的に低減し、安定性を維持できるか？
• RQ3 ステガードグリッド上に適用された本手法の高次有限差分スキームは、長距離超音波シミュレーションにおいて、スペクトル成分をどれほど正確に保持し、分散および減衰誤差を最小限に抑えることができるか？
• RQ4 本モデルは、特にスペックル輝度および信号対ゴミ比（LOC）に関して、実験的Bモード画像データをどれほど正確に再現できるか？
• RQ5 15点/波長の粗いグリッドを用いた場合、シミュレーションで保持される最大の高調波成分は何か？これはナイキスト周波数限界と比較してどうなるか？

主な発見

• 3つの緩和メカニズムを用いて1–10 MHzの周波数範囲で平均二乗誤差を最小化した結果、解析的減衰および分散則と優れた一致を示した。
• 緩和的PMLは、垂直入射時において反射係数を-60 dB未満に低下させ、グリッドサイズおよび境界層厚さにほとんど依存しなかった。これは、波の効果的吸収を確認するものであった。
• フルウェーブソルバは非線形波の7次高調波を正確に捉えており、ナイキスト周波数限界（7.5f₀）に近い周波数領域でのスペクトルサポートを示した。
• ブタの肝臓の実験的Bモード画像ではスペックル輝度が-57.8 dBであったが、散乱体インピーダンスを背景の4.0%に設定したシミュレーションでは-57.5 dBとよく一致した。
• 信号対ゴミ比（LOC）は実験で0.32 ± 0.093、シミュレーションで0.31 ± 0.037と測定され、チャネル信号ダイナミクスの高精度な再現性を確認した。
• 可視ヒューマンの腹部を用いた3次元波動伝播のシミュレーションでは、20 mmを超える深さで現実的なリバーブレーションリングダウンを再現し、時間領域および周波数領域の両方で実験データと一致した。