WaveProp_in_MATLAB icon indicating copy to clipboard operation
WaveProp_in_MATLAB copied to clipboard

Single-file implementations of 2D and 3D acoustic and elastic wave propagation in time domain using finite-differences(FDTD). Simple formulation and implementation

SIMPLE FDTD WAVE PROPAGATION IN MATLAB

Single-file vectorized implementations of wave propagation in MATLAB. We solve second-order wave equation in displacement formulation in time domain(FDTD). We don't account for derivatives of elastic parameters to keep it simple. Medium is surrounded by simple absorbing sponge boundaries with exponential decay (Cerjan, 1985).

Wavefield examples

DISCRETIZATION DETAILS:

  • Finite-Differences in Time Domain (FDTD)
  • Regular collocated grid
  • Explicit time stepping
  • O(2,2)
  • Conventional stencils derived from Taylor series:
    • in space [1: -2 :1]/dx2 and [1: -1: -1:1]/4dxdz
    • in time [1: -2 :1]/dt2

HOW TO?

Open folder with a scenario of interest e.g. acoustic_2D_FDTD_wave_propagation
Open .m file in MATLAB
Press "Run" button 

[email protected]

[email protected]