makeabhishek / WaveProp_in_MATLAB

Single-file codes in MATLAB for acoustic and elastic wave propagation in time domain using finite-differences(FDTD). As simple formulation and implementation as possible.

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).

DISCRETIZATION DETAILS:

• Finite-Differences in Time Domain (FDTD)
• 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

[e-mail]

[e-mail]