Wall plume and stratification model
This model is presented in details in Caudwell et al 2016, Convection at an isothermal wall in an enclosure and establishment of stratification, J. Fluid Mech or in my Phd Thesis.
It computes the evolution of the flow and the ambient stratification that appear in a closed cavity when a vertical wall is heated at a given constant temperature. It is based on the Morton et al (1965) entrainment theory and Germeles (1975) numerical scheme as per the turbulent part, and on Worster and Leitch (1985) similarity solutions as per the laminar part.
The input parameters are:
- the height of the cavity
- the global Rayleigh number
- the Prandtl number
The model mainly computes :
- the plume volume flux
- the plume momentum flux
- the plume buoyancy flux
- the ambient density field
This code is distributed under CeCILL License (GPL compatible licence).
git clone https://framagit.org/caudwell/wall-plume-stratification-model.git
cd wall-plume-stratification-model python wall_plume_model.py
The main file is called
wall_plume_model.py. It uses the parameters given in its last part, and calls the two attached routines
It possible to run the model as "fully turbulent" (
is_hybrid = False) or as "hybrid" including the laminar modeling (
is_hybrid = True).
It is possible to choose a constant (
is_alpha_var = False) or a variable (
is_alpha_var = True) entrainement coefficient.