[论文解读] 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

一句话总结

本文提出使用时变伴随方程来指导双曲守恒律中的自适应网格加密（AMR），仅针对影响特定关注区域的区域进行加密。通过计算前向解与伴随解之间的内积，该方法实现了精确且高效的加密——在保持波传播问题（如声学和海啸）精度的同时，显著降低了计算成本。

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中现有的Riemann求解器和Richardson误差估计器，以保持精度与效率。
• 在海啸应用中，使用线性化浅水方程进行伴随建模，并对海岸线相互作用采用壁边界条件。

实验结果

研究问题

• RQ1 是否可以利用伴随解识别波传播问题中仅影响特定目标区域的区域？
• RQ2 与基于表面高程或梯度的传统加密方法相比，伴随引导的AMR在计算成本和精度方面有何差异？
• RQ3 伴随方法是否能在复杂波动力学（如海啸传播）中减少模拟时间，同时保持精度？
• RQ4 仅对影响相关区域进行加密，对目标位置解的保真度有何影响？
• RQ5 如何实现基于伴随的加密的自动化，并将其集成到Clawpack和GeoClaw等现有AMR软件中？

主要发现

• 伴随方法通过仅对影响目标区域的区域进行加密，显著降低了计算成本，同时未牺牲解的精度。
• 在海啸算例中，基于伴随的加密在容差为0.09时达到与表面标志法相当的准确结果，但计算时间显著减少。
• 该方法成功捕捉了表面标志法在更高容差（0.14）下遗漏的后期到达波，如边缘波和反射海啸。
• 伴随方法在精度与效率之间实现了更优平衡，优于粗略和精细的表面标志策略。
• 在AMRClaw和GeoClaw中成功演示了基于伴随的加密集成，涵盖一维和二维声学问题以及全尺度海啸模型。
• 作者指出，未来改进可包括使用Richardson外推法进行误差界估计，以及采用自适应检查点技术以优化伴随解的加密。