This is a very minimal, acoustic, multidimensional FDTD simulation in under 25 lines of code. The FDTD calculations are actually performed using very few commands (the number of dimensions of your simulation +1). The source is a ‘hard source’ Gaussian pulse and there are no boundary conditions, or anything else to complicate the script as yet. This hopefully makes each process is very obvious. Pressure and Velocity field calculations are vectorized for efficiency and compact code. The output is animated over time.

This script is supplied ready to run. Just hit F5. The code is supplied in 2 dimensions, but, this can be changed to 1D or 3D with a couple of simple copy/paste and edits.

Enjoy! (email: nick@ihr.mrc.ac.uk)