search

daniel-koehn / DENISE-Black-Edition

2D time-domain isotropic (visco)elastic FD modeling and full waveform inversion (FWI) code for P/SV-waves
GNU General Public License v2.0
121 stars 66 forks source link

Topography #27

Closed gaoyajian closed 3 years ago

gaoyajian commented 3 years ago

Hey, dear Daniel,

I am a new user for your codes, I am wondering can I set a relief topography for the simulations?

Many thanks!

daniel-koehn commented 3 years ago

Hi gaoyalian,

The implementation of a relief topography by using vertically deformed meshes, like the one described in this paper:

Puente et al. (2014) Mimetic seismic wave modeling including topography on deformed staggered grids

is on my ToDo list, but not implemented, yet. You can force a free-surface boundary condition in the DENISE code by using the Improved Vaccumm Formulation (IVF):

Zeng et al. (2012) An improved vacuum formulation for 2D finite-difference modeling of Rayleigh waves including surface topography and internal discontinuities

and spatially discretize the free surface topography sufficiently.

Best regards,

Daniel