Adjoint Methods for Guiding Adaptive Mesh Refinement in Wave Propagation Problems
by Brisa N. Davis and R. J. LeVeque

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.

This paper was not accepted for publication.
Some of the material from the preprint has been incorporated into the new manuscript Adjoint Methods for Guiding Adaptive Mesh Refinement in Tsunami Modeling that focuses on tsunami modeling. Other material from this preprint will be used in a forthcoming manuscript on the adjoint method for AMR that will also discuss adjoint error estimation in more detail.

Preprint: http://arxiv.org/abs/1511.03645

Code: https://github.com/clawpack/adjoint

bibtex entry:

@misc{DavisLeVeque:adjoint1,
  author="B. N. Davis and R. J. LeVeque",
  title="Adjoint Methods for Guiding Adaptive Mesh Refinement in Wave Propagation Problems",
  howpublished="\url{http://arxiv.org/abs/1511.03645}",
  year="2015"
}