A staggeredgrid finitedifference scheme optimized in the time–space domain for modeling scalarwave propagation in geophysical problems
Abstract
For modeling scalarwave propagation in geophysical problems using finitedifference schemes, optimizing the coefficients of the finitedifference operators can reduce numerical dispersion. Most optimized finitedifference schemes for modeling seismicwave propagation suppress only spatial but not temporal dispersion errors. We develop a novel optimized finitedifference scheme for numerical scalarwave 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 a given range of wavenumbers. Dispersion analysis and numerical examples demonstrate that our optimized finitedifference 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 highorder finitedifference 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 finitedifference scheme is particularly useful for largescale 3D scalarwave modeling and inversion.
 Authors:
 Formerly Los Alamos National Laboratory, Geophysics Group, Los Alamos, NM 87545 (United States)
 Los Alamos National Laboratory, Geophysics Group, Los Alamos, NM 87545 (United States)
 Publication Date:
 OSTI Identifier:
 22382140
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 276; Other Information: Copyright (c) 2014 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ACCURACY; COMPARATIVE EVALUATIONS; COMPUTERIZED SIMULATION; CONTROL; ERRORS; FINITE DIFFERENCE METHOD; OPTIMIZATION; PHASE VELOCITY; SCALARS; SEISMIC WAVES; WAVE PROPAGATION
Citation Formats
Tan, Sirui, Email: siruitan@hotmail.com, and Huang, Lianjie, Email: ljh@lanl.gov. A staggeredgrid finitedifference scheme optimized in the time–space domain for modeling scalarwave propagation in geophysical problems. United States: N. p., 2014.
Web. doi:10.1016/J.JCP.2014.07.044.
Tan, Sirui, Email: siruitan@hotmail.com, & Huang, Lianjie, Email: ljh@lanl.gov. A staggeredgrid finitedifference scheme optimized in the time–space domain for modeling scalarwave propagation in geophysical problems. United States. doi:10.1016/J.JCP.2014.07.044.
Tan, Sirui, Email: siruitan@hotmail.com, and Huang, Lianjie, Email: ljh@lanl.gov. Sat .
"A staggeredgrid finitedifference scheme optimized in the time–space domain for modeling scalarwave propagation in geophysical problems". United States.
doi:10.1016/J.JCP.2014.07.044.
@article{osti_22382140,
title = {A staggeredgrid finitedifference scheme optimized in the time–space domain for modeling scalarwave propagation in geophysical problems},
author = {Tan, Sirui, Email: siruitan@hotmail.com and Huang, Lianjie, Email: ljh@lanl.gov},
abstractNote = {For modeling scalarwave propagation in geophysical problems using finitedifference schemes, optimizing the coefficients of the finitedifference operators can reduce numerical dispersion. Most optimized finitedifference schemes for modeling seismicwave propagation suppress only spatial but not temporal dispersion errors. We develop a novel optimized finitedifference scheme for numerical scalarwave 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 a given range of wavenumbers. Dispersion analysis and numerical examples demonstrate that our optimized finitedifference 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 highorder finitedifference 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 finitedifference scheme is particularly useful for largescale 3D scalarwave modeling and inversion.},
doi = {10.1016/J.JCP.2014.07.044},
journal = {Journal of Computational Physics},
number = ,
volume = 276,
place = {United States},
year = {Sat Nov 01 00:00:00 EDT 2014},
month = {Sat Nov 01 00:00:00 EDT 2014}
}

This document is intended to serve as a users guide for the timedomain atmospheric acoustic propagation suite (TDAAPS) program developed as part of the Department of Defense HighPerformance Modernization Office (HPCMP) Common HighPerformance Computing Scalable Software Initiative (CHSSI). TDAAPS performs staggeredgrid finitedifference modeling of the acoustic velocitypressure system with the incorporation of spatially inhomogeneous winds. Wherever practical the control structure of the codes are written in C++ using an object oriented design. Sections of code where a large number of calculations are required are written in C or F77 in order to enable better compiler optimization of these sections. Themore »

An improved rotated staggeredgrid finitedifference method with fourthorder temporal accuracy for elasticwave modeling in anisotropic media
The rotated staggeredgrid (RSG) finitedifference method is a powerful tool for elasticwave 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 fourthorder temporal accuracy to reduce the numerical dispersion associated with prolonged wave propagation or a large temporal step size. The highorder temporal accuracy is achieved by including highorder temporal derivatives, which can be converted to highorder 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 » 
An optimized finitedifference scheme for wave propagation problems
Two fullydiscrete finitedifference schemes for wave propagation problems are presented, a maximumorder scheme and an optimized (or spectrallike) scheme. Both combine a sevenpoint spatial operator and an explicit sixstage timemarch method. The maximumorder operator is fifthorder in space and is sixthorder in time for a linear problem with periodic boundary conditions. The phase and amplitude errors of the schemes obtained using Fourier analysis are given and compared with a secondorder and a fourthorder method. Numerical experiments are presented which demonstrate the usefulness of the schemes for a range of problems. For some problems, the optimized scheme leads to a reductionmore » 
Finite difference scheme for parabolic problems on composite grids with refinement in time and space
Finite difference schemes for transient convectiondiffusion problems on grids with local refinement in time and space are constructed and studied. The construction utilizes a modified upwind approximation and linear interpolation at the slave nodes. The proposed schemes are implicit of backward Euler type and unconditionally stable. Error analysis is presented in the maximum norm, and convergence estimates are derived for smooth solutions. Optimal approximation results for ratios between the spatial and time discretization parameters away from the CFL condition are shown. Finally, numerical examples illustrating the theory are given.