This is a much better approach to CFD / Navier-Stokes and will help you understand the various phenomenon along the way. https://lorenabarba.com/blog/cfd-python-12-steps-to-navier-s...
> So that means, while I know what it does, I don't really know how, since all the work is in that mysterious function.
that I spent the time to work it out myself. (answer: It arises from discretising the Laplacian -- 6 is the number of direct neighbours in 3D)
Lots of points in there like
> "Air is an example of a compressible fluid; you squish it and it gets smaller. Water is an example of an incompressible fluid; you squish it and it pushes back, and doesn't get any smaller" (this only really depends on the Mach number, Ma>~0.3 and you are in compressible territory for any fluid. Incompressibility usually means we assume the divergence of velocity of zero)
> "Incompressible fluids are simpler to simulate because their density and pressure is always constant." (This is only true if you choose to adopt a grad P = 0 approximation)
are incorrect from a physics perspective.
If you look at what real incompressible Navier-Stokes solvers do [1], it's mathematically totally different from what this post shows. In fact, the part that this post omits (handle the pressure gradient term by first approximate time stepping the velocity term by ignoring the pressure term and then correct by solving a Poisson equation for the pressure residual, and then correct the velocity) is the most expensive step in incompressible solvers by far.
[1] https://en.wikipedia.org/wiki/Projection_method_(fluid_dynam...
Iiuc water might compress ~50% at the right place in the Earth’s mantle, maybe just not looking much like liquid.
Has there been anything like this published in the past 20 years but for compressible fluids? I have wanted to make a simple atmospheric model for years but have been unable to because of the complexity and probably my lack of complete understanding.