SPECFEM / specfem2d

SPECFEM2D simulates forward and adjoint seismic wave propagation in two-dimensional acoustic, (an)elastic, poroelastic or coupled acoustic-(an)elastic-poroelastic media, with Convolution PML absorbing conditions.
https://specfem.org
GNU General Public License v3.0
203 stars 147 forks source link

Tunnel Simulation #1229

Closed zuoxuan-de closed 5 months ago

zuoxuan-de commented 5 months ago

I have designed a model of the palm face of a tunnel with a width of 10 m and a length of 50 m. The simulation area is 115 m long and 60 m wide, in which the thickness of the PML absorbing boundary layer is set to 5 m. The mesh profile is based on a grid of unstructured quadrilateral cells. The earthquake source is located at the center of the palm face. In this model, the longitudinal wave velocity of the normal surrounding rock is 5000 m/s, the transverse wave velocity is 2500 m/s, and the density is 2500 kg/m3; the Ricker sub-wave with a main frequency of 120 Hz is used as the excitation source, and the direction is perpendicular to the palm surface. The time step is taken as 0.2e-4. 10000 moments are taken. In the results, the wave field is zero until 6400 moments, what is the reason and how to solve it. There is also the fact that the wave field propagates very quickly, without any specific details of the process

Translated with DeepL.com (free version) forward_image000006400 forward_image000007600 forward_image000007800 forward_image000008000