Abstract/Summary

The representation of sub-grid scale convection is a weak aspect of weather and climate prediction models and the assumption that no net mass is transported by convection in parameterisations is increasingly unrealistic as models enter the grey zone, partially resolving convection. The solution of conditionally averaged equations of motion (multi-fluid equations) is proposed in order to avoid this assumption. Separate continuity, temperature and momentum equations are solved for inside and outside convective plumes which interact via mass transfer terms, drag and by a common pressure. This is not a convection scheme that can be used with an existing dynamical core -- this requires a whole new model. This paper presents stable numerical methods for solving the multi-fluid equations including large transfer terms between the environment and plume fluids. Without transfer terms the two fluids are not sufficiently coupled and solutions diverge. Two transfer terms are presented which couple the fluids together in order to stabilise the model: diffusion of mass between the fluids (similar to turbulent entrainment) and drag between the fluids. Transfer terms are also proposed to move buoyant air into the plume fluid and vice-versa as would be needed to represent initialisation and termination of sub-grid-scale convection. The transfer terms are limited (clipped in size) and solved implicitly in order to achieve bounded, stable solutions. Results are presented of a well resolved warm bubble with rising air being transferred to the plume fluid. For stability, equations are formulated in advective rather than flux form and solved using bounded finite volume methods. Discretisation choices are made to preserve boundedness and conservation of momentum and energy when mass is transferred between fluids. The formulation of transfer terms in order to represent sub-grid convection is the subject of future work.