The source codes available here are based on the publication:
The codes solve the Reynolds-averaged Navier-Stokes equaitons for a fully developed turbulent channel flow with varying properties, such as density and viscosity. The code demonstrates how to modify existing turbulence models to properly account these thermophysical property variations. Five models are used for demonstration:
- an algebraic eddy viscosity model (Cess, 1958),
- the Spalart and Allmaras model (1994),
- k-epsilon model based on Myong and Kasagi (1993),
- Menter's SST k-omega model (Menter, 1995)
- V2F model (medic and Durbin, 2012).
They are available as
- matlab (matlab/main.m) and
- python source in form of a jupyter notebook: main.ipynb (might take some time to load) or using the nbviewer
- matlab
- Jupyter Notebook, python3.5
Either run the matlab file main.m, or execute the jupyter notebook. The codes run cases for which DNS data is available to compare with. The DNS data is given in the directory DNS_data, which are based on this publication.