[論文レビュー] Benchmark problems for transcranial ultrasound simulation: Intercomparison of compressional wave models

Jean‐François Aubry, Oscar Bates+20 more|arXiv (Cornell University)|Feb 9, 2022

Ultrasound Imaging and Elastography被引用数 4

ひとこと要約

本論文は、頭蓋骨ホモジニアスを用いた圧縮波の伝播を用いて、11種類の異なる計算モデルを評価する標準化されたベンチマークスイートを提示する。モデル間の比較により、焦点圧力において10%未満、焦点位置において1 mm未塔の良好な一致が得られ、焦点超音波治療や位相補正の臨床的応用における数値シミュレーションの信頼性が裏付けられる。

ABSTRACT

Computational models of acoustic wave propagation are frequently used in transcranial ultrasound therapy, for example, to calculate the intracranial pressure field or to calculate phase delays to correct for skull distortions. To allow intercomparison between the different modeling tools and techniques used by the community, an international working group was convened to formulate a set of numerical benchmarks. Here, these benchmarks are presented, along with intercomparison results. Nine different benchmarks of increasing geometric complexity are defined. These include a single-layer planar bone immersed in water, a multi-layer bone, and a whole skull. Two transducer configurations are considered (a focused bowl and a plane piston), giving a total of 18 permutations of the benchmarks. Eleven different modeling tools are used to compute the benchmark results. The models span a wide range of numerical techniques, including the finite-difference time-domain method, angular-spectrum method, pseudospectral method, boundary-element method, and spectral-element method. Good agreement is found between the models, particularly for the position, size, and magnitude of the acoustic focus within the skull. When comparing results for each model with every other model in a cross comparison, the median values for each benchmark for the difference in focal pressure and position are less than 10\% and 1 mm, respectively. The benchmark definitions, model results, and intercomparison codes are freely available to facilitate further comparisons.

研究の動機と目的

• 頭蓋内超音波伝搬の計算モデルを評価するための標準化されたベンチマークフレームワークを確立すること。
• 頭蓋内超音波シミュレーションに用いられる異なる数値的手法の整合性と正確性を評価すること。
• 多様なモデリングツールや技術におけるシミュレーション結果のばらつきの原因を特定すること。
• 焦点超音波治療および頭内圧力予測の分野における計算モデリングのベストプラクティスを支援すること。
• 再現性と将来の新しいソルバーの検証を可能にするために、自由に利用可能なベンチマーク、データ、コードを提供すること。

提案手法

• 幾何学的複雑度が増す9つのベンチマーク問題を定義し、水中に浸漬された平面およびマルチレイヤー頭蓋骨モデル、および完全な頭蓋骨を含む。
• 焦点をもつボウル型および平面ピストン型の2種類のトランスデューサー設定を用い、合計18通りのベンチマークの組み合わせが得られる。
• 11種類の異なるモデリングツールが使用され、有限差分時域法、アングルスペクトル法、擬スペクトル法、境界要素法、スペクトル要素法など、多様な数値手法をカバーする。
• 各モデルが音響圧力場を計算し、焦点位置、焦点サイズ、圧力大きさなどの指標を用いて比較される。
• すべてのモデルペア間でのクロス比較が実施され、すべてのベンチマークにおいてモデル間の一致度が定量化される。
• すべてのベンチマーク定義、シミュレーション結果、および相互比較コードが公開され、再現性と拡張性が確保される。

実験結果

リサーチクエスチョン

• RQ1 異なる数値モデルが、幾何学的複雑度が異なる頭蓋骨構造を通じて超音波を伝搬させる際、どのように性能を発揮するか？
• RQ2 頭蓋骨通過後の脳内における焦点圧力場を予測する際、多様なモデリングツール間でどの程度の一致が達成できるか？
• RQ3 数値手法や実装選択の違いが、焦点位置や圧力大きさなどの重要なシミュレーション指標にどのように影響するか？
• RQ4 モデルの差異は、幾何学的複雑度やトランスデューサー設定に応じてどの程度変動するか？
• RQ5 標準化されたベンチマークフレームワークは、頭蓋内超音波治療における計算予測の不確実性を低減し、信頼性を向上させることができるか？

主な発見

• すべてのベンチマークにおいて、任意のモデル同士を比較した場合、焦点圧力の中央値差が10%未塔、焦点位置の中央値差が1 mm未塔であることが示された。
• 特に、頭蓋骨通過後の音響焦点の位置、サイズ、大きさにおいて良好な一致が確認された。
• 近場波挙動の違いに起因し、トランスデューサーと頭蓋骨の間で、特に全領域指標において大きな差が観察された。
• 有限差分法、スペクトル法、境界要素法を含む幅広い数値技術において、強い一貫性が示された。
• ベンチマーク定義、シミュレーション結果、および相互比較コードが公開されており、新しいモデリングツールの再現と将来の検証が可能である。
• 本研究は、弾性波モデルや患者特異的CT由来の材料パrameterを用いたシミュレーションを含む、将来の相互比較の基盤を提供する。