Finite difference method

Abstract: A particle velocity-strain, finite-difference (FD) method with a perfectly matched layer (PML) absorbing boundary condition is developed for the simulation of elastic wave propagation in multidimensional heterogeneous poroelastic media. Instead of the widely used second-order differential equations, a first-order hyperbolic leap-frog system is obtained from Biot’s equations. To achieve a high accuracy, the first-order hyperbolic system is discretized on a staggered grid both in time and space. The perfectly matched layer is used at the computational edge to absorb the outgoing waves. The performance of the PML is investigated by calculating the reflection from the boundary. The numerical method is validated by analytical solutions. This FD algorithm is used to study the interaction of elastic waves with a buried land mine. Three cases are simulated for a mine-like object buried in “sand,” in purely dry “sand” and in “mud.” The results show that the wave responses are significantly different in these cases. The target can be detected by using acoustic measurements after processing.

Abstract: The present investigation is concerned with a fourth order accurate finite difference method and its application to the study of the temporal and spatial stability of the three-dimensional compressible boundary layer flow on a swept wing. This method belongs to the class of compact two-point difference schemes discussed by White (1974) and Keller (1974). The method was apparently first used for solving the two-dimensional boundary layer equations. Attention is given to the governing equations, the solution technique, and the search for eigenvalues. A general purpose subroutine is employed for solving a block tridiagonal system of equations. The computer time can be reduced significantly by exploiting the special structure of two matrices.