Options
A Nodal Integral Scheme for Acoustic Wavefield Simulation
Journal
Pure and Applied Geophysics
ISSN
00334553
Date Issued
2022-10-01
Author(s)
Kumar, Niteen
Shekar, Bharath
Singh, Suneet
Abstract
The nodal integral method (NIM) has been applied to a variety of both linear and non-linear problems. Here, we outline the development of an NIM scheme for the indefinite Helmholtz equation to simulate acoustic wave propagation. NIM has a simple formulation, and can yield accurate results with coarse grids. In the development of the scheme, the partial differential equation is averaged over a node by a transverse integration process (TIP). While performing the TIP, the product of the wave number and the function are incorporated in pseudo source terms, leading to a simple scheme. The ordinary differential equations resulting from the TIP are solved analytically, constituting the distinguishing step in NIM that leads to accurate solutions. We present the dispersion analysis of the NIM scheme for the very first time to quantify numerical dispersion. We then benchmark the wavefields and time domain seismograms computed from NIM with the analytical solution for a homogeneous model. We find that the NIM scheme yields second-order accurate results for the homogeneous model. Finally, we compute time domain seismograms for the Marmousi model and compare the results with the second-order central finite difference method in the time domain. The second-order scheme presented here is the very first application of NIM to the solution of the acoustic wave equation.
Subjects