[論文レビュー] Adjoint Methods for Guiding Adaptive Mesh Refinement in Wave Propagation Problems

Brisa N. Davis, Randall J. LeVeque|arXiv (Cornell University)|Nov 11, 2015

Seismic Imaging and Inversion Techniques被引用数 5

ひとこと要約

本稿では、波動伝播問題（音響および津波を含む）において、特定の領域に影響を与える領域にのみメッシュを適応的に細分化するため、時間に依存する随伴方程式を用いたアプローチを提案する。前向き解と随伴解の内積を計算することで、影響を及ぼす領域を正確かつ効率的に特定可能となり、計算コストを著しく削減しながらも解の精度を維持する。

ABSTRACT

One difficulty in developing numerical methods for hyperbolic systems of conservation laws is the fact that solutions often contain regions where much higher resolution is required than elsewhere in the domain, particularly since the solution may contain discontinuities or other localized features. The Clawpack software deals with this issue by using block-structured adaptive mesh refinement to selectively refine around propagating waves. For problems where only a target area of the total solution is of interest, a method that allows identifying and refining the grid only in regions that influence this target area would significantly reduce the computational cost of finding a solution. In this work, we show that solving the time-dependent adjoint equation and using a suitable inner product with the forward solution allows more precise refinement of the relevant waves. We present acoustics examples in one and two dimensions and a tsunami propagation example. To perform these simulations, the use of the adjoint method has been integrated into the adaptive mesh refinement strategy of the open source Clawpack and GeoClaw software. We also present results that show that the accuracy of the solution is maintained and the computational time required is significantly reduced through the integration of the adjoint method into AMR.

研究の動機と目的

• 特定の領域に影響を与える領域にのみメッシュを細分化することで、波動伝播シミュレーションにおける計算コストを削減すること。
• 目標領域に影響を与えるグリッドセルを自動的に同定・細分化する手法を開発し、関係のない波動特徴の無駄な細分化を回避すること。
• ハイパボリックPDEに適したAMRClawおよびGeoClawフレームワークに、随伴に基づく細分化基準を統合すること。
• 従来のAMR戦略と比較して、著しくシミュレーション時間を短縮しながらも解の精度を維持すること。
• 津波のような地球物理的波動の効率的かつ的確なモデリングを可能にし、海岸付近の領域にのみ関心が向けられる状況を想定すること。

提案手法

• ハイパボリック系の時間に依存する随伴方程式を解き、目標関数に影響を与える領域を特定する。
• 前向き解と随伴解の内積を計算し、目標領域への影響度を定量化する。
• 内積ノルムに許容誤差を設け、前向き問題における細分化対象セルを特定する。
• AMRClawおよびGeoClawのブロック構造化AMRフレームワークに、随伴に基づく細分化基準を統合する。
• AMRClawに既存のリーマンソルバーやリチャードソン誤差推定器を活用し、精度と効率を維持する。
• 津波応用では、線形化された浅水域方程式を随伴モデリングに用い、海岸線での相互作用には壁境界条件を適用する。

実験結果

リサーチクエスチョン

• RQ1 随伴解を用いることで、波動伝播問題において特定の目標領域に影響を与える領域のみを同定できるか？
• RQ2 従来の表面高さ基準または勾配基準の細分化と比較して、随伴指導致しAMRは計算コストと精度の面でどのように異なるか？
• RQ3 複雑な波動ダイナミクス（例：津波伝搬）において、随伴法によりシミュレーション時間を短縮しつつも精度を維持できるか？
• RQ4 影響を及ぼす領域にのみ細分化することの、目標位置における解の忠実度に与える影響は何か？
• RQ5 随伴基準に基づく細分化を、ClawpackおよびGeoClawのような既存のAMRソフトウェアに自動化して統合するにはどうすればよいか？

主な発見

• 随伴法により、目標領域に影響を与える領域にのみメッシュを細分化することで、著しい計算コストの削減が達成され、解の精度を損なわない。
• 津波の例では、随伴に基づく細分化は、許容誤差0.09で表面フラグ付きの結果と同等の精度を達成したが、著しく短い計算時間で実現した。
• 高い許容誤差（0.14）で表面フラグ付きの手法では見逃された後から到着する波（縁波や反射津波）を、随伴法により正確に捉えた。
• 随伴法は、粗いおよび細かい表面フラグ付き戦略を組み合わせた精度と効率の両立に優れている。
• AMRClawおよびGeoClawへの随伴基準に基づく細分化の統合は、1次元および2次元の音響問題および本格的な津波モデルにおいて成功裏に実証された。
• 著者らは、将来的にはリチャードソンの外挿法による誤差境界推定や、随伴解の細分化を最適化するためのアダプティブチェックポイント技術の導入が可能になると指摘している。