A finite element solution to turbulent diffusion in a convective boundary layer

A second‐order closure turbulence model is used to simulate the plume behaviour of a passive contaminant dispersed in a convective boundary layer. A time‐splitting finite element scheme is used to solve the set of partial differential equations. It is shown that the second‐order closure model compares favourably with recent findings from laboratories, wind‐tunnel experiments and large‐eddy simulations. We also compare the second‐order closure model with the commonly used K‐diffusion model for the same meteorological conditions. Case studies also show the effects of model parameters and turbulence variables on the plume behaviour.