A very simple 1D CFD tool that solves the compressible euler equations
It is using a Multi-Stage "MUSTA" solution scheme, see [1] for details.
The equations are formulated such we simply solve:
dq/dt + df(q)/dx = 0
where q and f(q) are the vectors:
q = [rho, rho*v, E]
f(q) = [rho*v, rho*v^2, (E+p)v]
It was developed as a test/learning/demonstration project to see how small a CFD program in python can be, also demonstrating usage of a thermodynamic library for equation of state calculations. In this form here, the thermodynamic library usage is removed and a simple ideal gas equation of state is used.
TODO:
- Add multi-component support
- Extend to two-phase, e.g. like in [2]
- Find more insanely cool ways to visualize simulations
[1] Titarev, V. A., and E. F. Toro. "MUSTA schemes for multi-dimensional hyperbolic systems: analysis and improvements." International journal for numerical methods in fluids 49.2 (2005): 117-148.
[2] Munkejord, Svend Tollak, Steinar Evje, and Tore Flåtten. "A MUSTA scheme for a nonconservative two-fluid model." SIAM Journal on Scientific Computing 31.4 (2009): 2587-2622.
Usage instructions:
python run.py
python view.py