This study describes the derivation of an eddy diffusivity for convective decaying turbulence in the residual layer. The approach is based on the budget equation for turbulent kinetic energy(TKE) where buoyant and shear terms are disregard and turbulence is assumed to be homogeneous and isotropic. In such equation the inertial energy transfer term is closed using a model suggested by Pao on the basis of dimensional considerations. In doing this, we obtain a budget equation for TKE and an analytical solution for the 3-D Energy Density Spectrum(EDS). Furthermore by the isotropy condition, the time decaying vertical velocity variance is calculated and compared with LES data.
The comparison is pretty good until t*=tzi/w*~1. For larger times, the vertical velocity variance calculated from our model decays with the time according to a power law with the exponent -2. Employing an expression suggested by Hanna, described in terms of the spectral wave-number maximum, two vertical decaying eddy diffusivities are derived. One utilizing the decaying vertical velocity variance calculed by LES and the other calculated from our model. Both eddy diffusivities show a good agreement until t*~1. As a consequence of the exponent -2 obtained for the decaying vertical velocity variance calculated from LES, for larger times the LES eddy diffusivity decays more strongly. Finally, the eddy diffusivities are fitted by empirical curves. These algebric fittings can be used in atmospheric diffusion models to parameterize the dispersion in a decaying convective turbulence in the residual layer.
Supplementary URL: