Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces

Wang, Siyang; Petersson, N. Anders

doi:10.1137/18m1211465

Title: Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces

Full Record
Other Related Research

Abstract

In this work, we analyze two types of summation-by-parts finite difference operators for approximating the second derivative with variable coefficient. The first type uses ghost points, while the second type does not use any ghost points. A previously unexplored relation between the two types of summation-by-parts operators is investigated. By combining them we develop a new fourth order accurate finite difference discretization with hanging nodes on the mesh refinement interface. We take the model problem as the two-dimensional acoustic wave equation in second order form in terms of acoustic pressure, and we prove energy stability for the proposed method. Compared to previous approaches using ghost points, the proposed method leads to a smaller system of linear equations that needs to be solved for the ghost point values. Another attractive feature of the proposed method is that the explicit time step does not need to be reduced relative to the corresponding periodic problem. Numerical experiments, both for smoothly varying and discontinuous material properties, demonstrate that the proposed method converges to fourth order accuracy. A detailed comparison of the accuracy and the time-step restriction with the simultaneous-approximation-term penalty method is also presented. (An erratum is attached.)

Authors:

Wang, Siyang ^[1]; Petersson, N. Anders ^[2]

Chalmers Univ. of Technology, Gothenburg (Sweden); Univ. of Gothenburg (Sweden)
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States). Center for Applied Scientific Computing

Publication Date:: 2019-10-22

Research Org.:: Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

Sponsoring Org.:: USDOE National Nuclear Security Administration (NNSA)

OSTI Identifier:: 1734602

Report Number(s):: LLNL-JRNL-757334
Journal ID: ISSN 1064-8275; 945085

Grant/Contract Number:: AC52-07NA27344

Resource Type:: Accepted Manuscript

Journal Name:: SIAM Journal on Scientific Computing

Additional Journal Information:: Journal Volume: 41; Journal Issue: 5; Journal ID: ISSN 1064-8275

Publisher:: SIAM

Country of Publication:: United States

Language:: English

Subject:: 97 MATHEMATICS AND COMPUTING; wave equation; finite difference methods; summation-by-parts; ghost points; nonconforming; mesh refinement

Citation Formats


                    Wang, Siyang, and Petersson, N. Anders. Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces.  United States: N. p., 2019. 
Web.  doi:10.1137/18m1211465.

Copy to clipboard


                    Wang, Siyang, & Petersson, N. Anders. Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces.  United States.  https://doi.org/10.1137/18m1211465

Copy to clipboard


                    Wang, Siyang, and Petersson, N. Anders. Tue .  
"Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces".  United States.  https://doi.org/10.1137/18m1211465.  https://www.osti.gov/servlets/purl/1734602.

Copy to clipboard


                    
@article{osti_1734602,

  title        = {Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces},

  author       = {Wang, Siyang and Petersson, N. Anders},

  abstractNote = {In this work, we analyze two types of summation-by-parts finite difference operators for approximating the second derivative with variable coefficient. The first type uses ghost points, while the second type does not use any ghost points. A previously unexplored relation between the two types of summation-by-parts operators is investigated. By combining them we develop a new fourth order accurate finite difference discretization with hanging nodes on the mesh refinement interface. We take the model problem as the two-dimensional acoustic wave equation in second order form in terms of acoustic pressure, and we prove energy stability for the proposed method. Compared to previous approaches using ghost points, the proposed method leads to a smaller system of linear equations that needs to be solved for the ghost point values. Another attractive feature of the proposed method is that the explicit time step does not need to be reduced relative to the corresponding periodic problem. Numerical experiments, both for smoothly varying and discontinuous material properties, demonstrate that the proposed method converges to fourth order accuracy. A detailed comparison of the accuracy and the time-step restriction with the simultaneous-approximation-term penalty method is also presented. (An erratum is attached.)},

  doi          = {10.1137/18m1211465},

  journal      = {SIAM Journal on Scientific Computing},

  number       = 5,

  volume       = 41,

  place        = {United States},

  year         = {2019},

  month        = {10}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1137/18m1211465

Other availability

Search WorldCat to find libraries that may hold this journal

Save / Share:

Export Metadata

Save to My Library

Similar Records in DOE PAGES and OSTI.GOV collections:

Elastic Wave Propagation in Curvilinear Coordinates with Mesh Refinement Interfaces by a Fourth Order Finite Difference Method

Journal Article Zhang, Lu ; Wang, Siyang ; Petersson, N. Anders - SIAM Journal on Scientific Computing

In this work, we develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinementmore »« less
https://doi.org/10.1137/20m1339702

Full Text Available
Super-Grid Modeling of the Elastic Wave Equation in Semi-Bounded Domains

Journal Article Petersson, N. Anders ; Sjögreen, Björn - Communications in Computational Physics

Abstract We develop a super-grid modeling technique for solving the elastic wave equation in semi-bounded two- and three-dimensional spatial domains. In this method, waves are slowed down and dissipated in sponge layers near the far-field boundaries. Mathematically, this is equivalent to a coordinate mapping that transforms a very large physical domain to a significantly smaller computational domain, where the elastic wave equation is solved numerically on a regular grid. To damp out waves that become poorly resolved because of the coordinate mapping, a high order artificial dissipation operator is added in layers near the boundaries of the computational domain. Wemore »« less
https://doi.org/10.4208/cicp.290113.220514a

Full Text Available
A Fourth Order Difference Scheme for the Maxwell Equations on Yee Grid

Journal Article Fathy, Aly E ; Wilson, Joshua L - Journal of Hyperbolic Differential Equations

The Maxwell equations are solved by a long-stencil fourth order finite difference method over a Yee grid, in which different physical variables are located at staggered mesh points. A careful treatment of the numerical values near the boundary is introduced, which in turn leads to a 'symmetric image' formula at the 'ghost' grid points. Such a symmetric formula assures the stability of the boundary extrapolation. In turn, the fourth order discrete curl operator for the electric and magnetic vectors gives a complete set of eigenvalues in the purely imaginary axis. To advance the dynamic equations, the four-stage Runge-Kutta method ismore »« less
https://doi.org/10.1142/S0219891608001623
Dynamic earthquake rupture simulations on nonplanar faults embedded in 3D geometrically complex, heterogeneous elastic solids

Journal Article Duru, Kenneth ; Dunham, Eric M. ; Institute for Computational and Mathematical Engineering, Stanford University, Stanford, CA - Journal of Computational Physics

Dynamic propagation of shear ruptures on a frictional interface in an elastic solid is a useful idealization of natural earthquakes. The conditions relating discontinuities in particle velocities across fault zones and tractions acting on the fault are often expressed as nonlinear friction laws. The corresponding initial boundary value problems are both numerically and computationally challenging. In addition, seismic waves generated by earthquake ruptures must be propagated for many wavelengths away from the fault. Therefore, reliable and efficient numerical simulations require both provably stable and high order accurate numerical methods. We present a high order accurate finite difference method for: a)more »« less
https://doi.org/10.1016/J.JCP.2015.10.021
Stable grid refinement and singular source discretization for seismic wave simulations

Journal Article Petersson, N A ; Sjogreen, B - Communications in Computational Physics, vol. 8, no. 5, May 31, 2010, pp. 1074-1110

An energy conserving discretization of the elastic wave equation in second order formulation is developed for a composite grid, consisting of a set of structured rectangular component grids with hanging nodes on the grid refinement interface. Previously developed summation-by-parts properties are generalized to devise a stable second order accurate coupling of the solution across mesh refinement interfaces. The discretization of singular source terms of point force and point moment tensor type are also studied. Based on enforcing discrete moment conditions that mimic properties of the Dirac distribution and its gradient, previous single grid formulas are generalized to work in themore »« less
Full Text Available

Similar Records