skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

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

Authors:
; ; ;
Publication Date:
Sponsoring Org.:
USDOE
OSTI Identifier:
1396504
Grant/Contract Number:
AC52-07NA27344
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Computers and Fluids
Additional Journal Information:
Journal Volume: 149; Journal Issue: C; Related Information: CHORUS Timestamp: 2018-01-08 05:13:54; Journal ID: ISSN 0045-7930
Publisher:
Elsevier
Country of Publication:
United Kingdom
Language:
English

Citation Formats

Prochnow, Bo, O’Reilly, Ossian, Dunham, Eric M., and Petersson, N. Anders. Treatment of the polar coordinate singularity in axisymmetric wave propagation using high-order summation-by-parts operators on a staggered grid. United Kingdom: N. p., 2017. Web. doi:10.1016/j.compfluid.2017.03.015.
Prochnow, Bo, O’Reilly, Ossian, Dunham, Eric M., & Petersson, N. Anders. Treatment of the polar coordinate singularity in axisymmetric wave propagation using high-order summation-by-parts operators on a staggered grid. United Kingdom. doi:10.1016/j.compfluid.2017.03.015.
Prochnow, Bo, O’Reilly, Ossian, Dunham, Eric M., and Petersson, N. Anders. 2017. "Treatment of the polar coordinate singularity in axisymmetric wave propagation using high-order summation-by-parts operators on a staggered grid". United Kingdom. doi:10.1016/j.compfluid.2017.03.015.
@article{osti_1396504,
title = {Treatment of the polar coordinate singularity in axisymmetric wave propagation using high-order summation-by-parts operators on a staggered grid},
author = {Prochnow, Bo and O’Reilly, Ossian and Dunham, Eric M. and Petersson, N. Anders},
abstractNote = {},
doi = {10.1016/j.compfluid.2017.03.015},
journal = {Computers and Fluids},
number = C,
volume = 149,
place = {United Kingdom},
year = 2017,
month = 6
}

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on March 27, 2018
Publisher's Accepted Manuscript

Save / Share:
  • We develop a fourth order accurate finite difference method for solving the three-dimensional elastic wave equation in general heterogeneous anisotropic materials on curvilinear grids. The proposed method is an extension of the method for isotropic materials, previously described in the paper by Sjögreen and Petersson (2012) [11]. The method we proposed discretizes the anisotropic elastic wave equation in second order formulation, using a node centered finite difference method that satisfies the principle of summation by parts. The summation by parts technique results in a provably stable numerical method that is energy conserving. Also, we generalize and evaluate the super-grid far-fieldmore » technique for truncating unbounded domains. Unlike the commonly used perfectly matched layers (PML), the super-grid technique is stable for general anisotropic material, because it is based on a coordinate stretching combined with an artificial dissipation. Moreover, the discretization satisfies an energy estimate, proving that the numerical approximation is stable. We demonstrate by numerical experiments that sufficiently wide super-grid layers result in very small artificial reflections. Applications of the proposed method are demonstrated by three-dimensional simulations of anisotropic wave propagation in crystals.« less
  • Cited by 1
  • For modeling scalar-wave propagation in geophysical problems using finite-difference schemes, optimizing the coefficients of the finite-difference operators can reduce numerical dispersion. Most optimized finite-difference schemes for modeling seismic-wave propagation suppress only spatial but not temporal dispersion errors. We develop a novel optimized finite-difference scheme for numerical scalar-wave modeling to control dispersion errors not only in space but also in time. Our optimized scheme is based on a new stencil that contains a few more grid points than the standard stencil. We design an objective function for minimizing relative errors of phase velocities of waves propagating in all directions within amore » given range of wavenumbers. Dispersion analysis and numerical examples demonstrate that our optimized finite-difference scheme is computationally up to 2.5 times faster than the optimized schemes using the standard stencil to achieve the similar modeling accuracy for a given 2D or 3D problem. Compared with the high-order finite-difference scheme using the same new stencil, our optimized scheme reduces 50 percent of the computational cost to achieve the similar modeling accuracy. This new optimized finite-difference scheme is particularly useful for large-scale 3D scalar-wave modeling and inversion.« less
  • The rotated staggered-grid (RSG) finite-difference method is a powerful tool for elastic-wave modeling in 2D anisotropic media where the symmetry axes of anisotropy are not aligned with the coordinate axes. We develop an improved RSG scheme with fourth-order temporal accuracy to reduce the numerical dispersion associated with prolonged wave propagation or a large temporal step size. The high-order temporal accuracy is achieved by including high-order temporal derivatives, which can be converted to high-order spatial derivatives to reduce computational cost. Dispersion analysis and numerical tests show that our method exhibits very low temporal dispersion even with a large temporal step sizemore » for elastic-wave modeling in complex anisotropic media. Using the same temporal step size, our method is more accurate than the conventional RSG scheme. In conclusion, our improved RSG scheme is therefore suitable for prolonged modeling of elastic-wave propagation in 2D anisotropic media.« less