Navier-Stokes with Mimetic Methods

2D Lock exchange test case using mimetic methods

Single, extensively commented file that solves NS equations assuming small density and temperature variations (Boussinesq). I have opted for mimetic methods (via MOLE) and explicit forward time schemes. More elaborated and precise time discretization schemes contribute very little to this particular scenario. Numerical diffusion could be avoided by computing the transient terms using a symplectic method such as Leapfrog.

For sharp Kelvin-Helmholtz billows, use a dx == 0.0625, keep in mind the CFL condition imposed on dt for a two-dimensional advection problem. Picture was obtained using m = 1600, n = 320, and dt = 0.1

The most time-consuming part (as expected) is the computation of the pressure field, Laplacian matrix is highly sparse but not positive definite.