My PhD project involves using a Boussinesq numerical model for predicting gravity waves over mountains, adjusted to include turbulent fluxes. In this model, the independant variables are split into three components: a background state, perturbations due to the orography, and averaged turbulent fluxes.
Prognostic equations will be used to solve for the turbulent fluxes, which will require assumptions (or closure schemes) for unsolvable variables such as third moment terms and terms involving turbulent pressure. These unsolvable terms will be approximated mainly by down-gradient theory and return-to-isentropy.
The model will also be compared with observations. Radiosonde soundings will provide the background flow, and the resulting turbulent fluxes will be compared to fluxes calculated from aircraft data.