Physics – Classical Physics
Scientific paper
2007-06-26
Geophysical Journal International 172 (2008) 252-261
Physics
Classical Physics
accepted and to be published in Geophys. J. Int
Scientific paper
10.1111/j.1365-246X.2007.03620.x
A method is proposed for accurately describing arbitrary-shaped free boundaries in single-grid finite-difference schemes for elastodynamics, in a time-domain velocity-stress framework. The basic idea is as follows: fictitious values of the solution are built in vacuum, and injected into the numerical integration scheme near boundaries. The most original feature of this method is the way in which these fictitious values are calculated. They are based on boundary conditions and compatibility conditions satisfied by the successive spatial derivatives of the solution, up to a given order that depends on the spatial accuracy of the integration scheme adopted. Since the work is mostly done during the preprocessing step, the extra computational cost is negligible. Stress-free conditions can be designed at any arbitrary order without any numerical instability, as numerically checked. Using 10 grid nodes per minimal S-wavelength with a propagation distance of 50 wavelengths yields highly accurate results. With 5 grid nodes per minimal S-wavelength, the solution is less accurate but still acceptable. A subcell resolution of the boundary inside the Cartesian meshing is obtained, and the spurious diffractions induced by staircase descriptions of boundaries are avoided. Contrary to what occurs with the vacuum method, the quality of the numerical solution obtained with this method is almost independent of the angle between the free boundary and the Cartesian meshing.
Gélis Céline
Lombard Bruno
Piraux Joël
Virieux Jean
No associations
LandOfFree
Free and smooth boundaries in 2-D finite-difference schemes for transient elastic waves does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Free and smooth boundaries in 2-D finite-difference schemes for transient elastic waves, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Free and smooth boundaries in 2-D finite-difference schemes for transient elastic waves will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-506111