Logo of Science Foundation Ireland  Logo of the Higher Education Authority, Ireland7 Capacities
Ireland's High-Performance Computing Centre | ICHEC
Home | News | Infrastructure | Outreach | Services | Research | Support | Education & Training | Consultancy | About Us | Login


Title:3D rotated and standard staggered finite-difference solutions to Biot's poroelastic wave equations: Stability condition and dispersion analysis
Authors:O'Brien, G.S., 2010
Abstract: A fourth-order in space and second-order in time 3D staggered 􏰀SG􏰁 and rotated-staggered-grid 􏰀RSG􏰁 method for the solution of Biot’s equation are presented. The numerical dispersion and stability conditions are derived using a von Neumann analysis. The exact stability condition is calculated from the roots of a 12th-order polynomial and therefore no nontrivial expression ex- ists. To overcome this, a 1D stability condition is usually general- ized to three dimensions. It is shown that in certain cases, the 1D approximate stability condition is violated by a 3D SG method. The RSG method obeys the approximate 1D stability condition for the material properties and spatiotemporal scales in the exam- ples shown. Both methods have been verified against an analyti- cal solution for an infinite homogeneous porous medium with a misfit error of less than 0.5%. A free surface has been implemented to test the accuracy of this boundary condition. It also serves as a test of the methods to include high material contrasts. The methods have been compared with a quasi-analytical solution. For the specific material properties, spatial grid scaling, and propagation distance used in the test, a maximum error of 3.5% for the SG and 4.1% for the RSG was found. These errors depend on the propagation distance, temporal and spatial scales, and ac- curacy of the quasi-analytical solution. No discernable differ- ence was found between the two methods except for time steps comparable with the stability-criteria threshold time step, the SG was found to be unstable. However, the RSG remained stable for a homogeneous half-space. Time steps, comparable to the stabili- ty criteria, reduce the computational time at the cost of a reduc- tion in accuracy. The methods allow wave propagation to be modeled in a porous medium with a free surface.
ICHEC Project:Seismic source modelling and propagation in complex 3D Earth models
Publication:Geophysics, Vol 75, No. 4, T111-119,
URL: http://dx.doi.org/10.1190/1.3432759
Keywords: porosity, porous materials, seismic waves, seismology
Status: Published

return to publications list