Plume Simulation in Python/NumPy

This short video presents the flow simulation using the Paper by Jos Stam and her algorithm/method presented in the SIGGRAPH 1999.

In this problem, a source of momentum (a force source) is added in a small region in the lower central part of the domain, which generates an upward plume represented by the velocity vector field. The contour map represents the vorticity field (curl of velocity), which in 2D has only one non-zero component, the normal to the image plane.

For CFD engineers, this methodology seems strange at first sight. However, it intends to generate a reliable and convincing solution of the Navier-Stokes Equation. This approach is adequate for games, movies FX, Visual Effects, and Computer Graphics.
I actually do not know if it is possible to calculate engineering quantities with precision such as CFD classic methods.

Nonetheless, I think really fun to create code that solves a flow with a simple algorithm. Very instructive, and didactic to my classes.

I have implemented in Google COLAB based on the video of @Machine Learning and Simulation:

I advise this channel as a good reference for Python and Simulation.

[QuippersUnited]

*sees Navier-Stokes in description*
*blacks out momentarily*
Wha-what just happened?