Generation of inflow turbulence in large-eddy simulations of nonneutral atmospheric boundary layers with the cell perturbation method

Realistic multiscale simulations involve coupling of mesoscale and large-eddy simulation (LES) models, thus requiring efficient generation of turbulence in nested LES domains. Herein, we extend our previous work on the cell perturbation (CP) method to nonneutral atmospheric boundary layers (ABLs). A modified Richardson number scaling is proposed to determine the amplitude of the potential temperature perturbations in stableABLs, withRi_m ∼ -1.0 overall providing optimum turbulence transition to a fully developed state (fetch reduced by a factor of 4-5, compared to the unperturbed cases). In the absence of perturbations, turbulence onset is triggered by a Kelvin- Helmholtz instability, typically occurring in the vicinity of the low-level jet maximum. It is found that a turbulent length scale l=q^1/2/N can be used tomore accurately estimate the optimumRi_m, where q is the turbulence kinetic energy, and N is the Brunt-Väisälä frequency. In convective ABLs, a perturbation amplitude based on mixed layer temperature variance scaling is proposed:σθ. For that criterion to be optimum, the ratioU_ci/w*,whereU_ci is thewind speed at the top of the capping inversion, and w* is the convective velocity scale, needs to be incorporated: σθ(U_ci/w*). This allows us to account for the competing roles of the surface thermal instability and the mean flow advection. For U_ci/w*∼10, the development fetch is reduced by a factor of ∼6, while whenU_ci/w*≤3, the use of the CPmethod does not provide a significant advantage in the ability to generate turbulence, provided a smoothmesoscale inflow.