Conservative-finite-difference-toolbox

Authors

• Mohammad Afzal Shadab ([email protected])
• Marc Hesse ([email protected])

Affiliation

Oden Institute for Computational Engineering and Sciences / Jackson School of Geosciences / University of Texas Institute for Geophysics The University of Texas at Austin

3-line Summary

This work is on developing a second order, staggered grid based finite difference solver in 1/2/3D. The information is given in the comments and test scripts are also enclosed. The solver has been validated with analytical results.

Spreading of a mound inside a porous domain: Conservative vs Non-Conservative schemes

Cartesian coordinates (free surface height, h vs location, x)

test_cartesian.mov

Spherical shell (free surface height, h vs location, Θ) coordinates

test_spherical.mov

Dependences

The codes and live scripts requires the Matlab.

Quick Usage

After cloning the repository, run the Matlab files directly.

Acknowledgement

I would like to acknowledge Professor Marc Hesse for teaching the course in 2021 and providing the basic codes. The course website has most of the theory circumscribing these codes which can be found here.