R. Martin1, D. Komatitsch1,2, S. D. Gedney3, E. Bruthiaux1,4
CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.1, pp. 17-42, 2010, DOI:10.3970/cmes.2010.056.017
Abstract Unsplit convolutional perfectly matched layers (CPML) for the velocity and stress formulation of the seismic wave equation are classically computed based on a second-order finite-difference time scheme. However it is often of interest to increase the order of the time-stepping scheme in order to increase the accuracy of the algorithm. This is important for instance in the case of very long simulations. We study how to define and implement a new unsplit non-convolutional PML called the Auxiliary Differential Equation PML (ADE-PML), based on a high-order Runge-Kutta time-stepping scheme and optimized at grazing incidence. We demonstrate More >