3DFD_DVS
The program is designed for computation of seismic wavefields in 3D heterogeneous surface geological structures with planar free surface due to surface and near-surface point doublecouple sources
Authors: | Jozef Kristek, Peter Moczo |
---|---|
General: | Wave propagation |
Code level: | Production Code |
Language: | Fortran 90/95 or later |
Associated groups: | Local Scale |
Supercomputing: | Shared memory (e.g., openMP) |
Grid dimensions: | 3D |
Coordinate system: | Cartesian |
Method: | Finite differences |
Tested operating system: | Linux, Unix (e.g. SunOS, HP Unix), Windows |
Rheology: | Viscoelastic isotropic |
Boundary conditions: | Free surface (planar), Absorbing boundaries, Reflecting boundaries |
Grid type: | Rectangular, regular, Other |
Solution type: | Numerical |
Total hits: | 9094 |
Latest Releases
- 3DFD_DVS 8.0 alpha (2005/11/24 13:17:55.098 GMT+1)
- Show all releases
Project Description
Program 3DFD_DVS is designed for the finite-difference simulation of seismic wave propagation and seismic ground motion in a 3D surface heterogeneous viscoelastic structure with a planar free surface. The computational algorithm is based on the explicit heterogeneous finite-difference scheme solving equations of motion in the heterogeneous viscoelastic medium with material discontinuities. The scheme is 4th-order accurate in space and 2nd-order accurate in time. The displacement-velocity-stress scheme is constructed on a staggered finite-difference grid. The computational region is represented by a volume of a parallelepiped with the top side representing a planar free surface, and bottom, rear, front, left and right sides representing either non-reflecting boundaries or planes of symmetry. Different types of non-reflecting boundaries can be chosen on different sides of the computational region. The discontinuous spatial grid is used to cover the computational region. The upper part of the grid has three times smaller grid spacing than the lower part. Each part itself is a uniform rectangular grid. The rheology of the medium corresponds to the generalized Maxwell body. This makes it possible to account both for spatially varying quality factors of the P and S waves and for arbitrary Q-omega law. The wavefield is excited by a set of point double-couple sources. The core memory optimization is applied in order to significantly reduce requirements of the computer’s core memory.