skip to main content


Title: Treatment of the polar coordinate singularity in axisymmetric wave propagation using high-order summation-by-parts operators on a staggered grid

In this paper, we develop a high-order finite difference scheme for axisymmetric wave propagation in a cylindrical conduit filled with a viscous fluid. The scheme is provably stable, and overcomes the difficulty of the polar coordinate singularity in the radial component of the diffusion operator. The finite difference approximation satisfies the principle of summation-by-parts (SBP), which is used to establish stability using the energy method. To treat the coordinate singularity without losing the SBP property of the scheme, a staggered grid is introduced and quadrature rules with weights set to zero at the endpoints are considered. Finally, the accuracy of the scheme is studied both for a model problem with periodic boundary conditions at the ends of the conduit and its practical utility is demonstrated by modeling acoustic-gravity waves in a magmatic conduit.
 [1] ;  [2] ;  [3] ;  [4]
  1. Stanford Univ., CA (United States). Dept. of Geophysics
  2. Stanford Univ., CA (United States). Dept. of Geophysics; Linköping Univ. (Sweden). Dept. of Mathematics. Division of Computational Mathematics
  3. Stanford Univ., CA (United States). Dept. of Geophysics. Inst. for Computational and Mathematical Engineering
  4. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States). Center for Applied Scientific Computing
Publication Date:
Report Number(s):
Journal ID: ISSN 0045-7930
Grant/Contract Number:
Accepted Manuscript
Journal Name:
Computers and Fluids
Additional Journal Information:
Journal Volume: 149; Journal ID: ISSN 0045-7930
Research Org:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Stanford Univ., CA (United States)
Sponsoring Org:
Country of Publication:
United States
97 MATHEMATICS AND COMPUTING; 42 ENGINEERING; polar coordinate singularity; summation-by-parts; high-order finite difference method; staggered grid; axisymmetric wave propagation
OSTI Identifier:
Alternate Identifier(s):
OSTI ID: 1396504