Large-Eddy Simulations of the Atmospheric Boundary Layer Using a New Subgrid-Scale Model – II. Weakly and moderately stable cases

A series of simulations under weakly to moderately stable boundary layers (SBLs) have been performed using the proposed subgrid-scale (SGS) model implemented into the Terminal Area Simulation System (TASS). The proposed SGS model incorporates some aspects of the two-part eddy viscosity SGS model of Sullivan et al. (1994) and further refinements which include the dependence of SGS mixing length on stratification, two-part separation of the SGS eddy diffusivity of heat, and more realistic empirical forms of Monin–Obukhov similarity functions. The potential temperature profiles from simulations clearly show a three-layer structure: a stable surface layer of strong gradients, a middle layer of small gradients, and an inversion layer on the top. The wind speed profiles show the formation of low level jet (LLJ). However, the sub-layer structures under moderately SBLs differ from those under weakly SBLs. Both the momentum and heat fluxes decrease almost linearly in the lower part of the SBL. The near surface values of the normalized turbulent kinetic energy (TKE/u _* ²) in all simulations are about 4 which is much less than the typical value of 5.5 under the neutral condition. The decay of turbulence first occurs in the area with large values of Richardson number (R _i<0.2). Generally, instantaneous values of the TKE and R _i at the various grid points are negatively correlated, but there is not a unique relationship between the two parameters.     

Large-Eddy Simulations of the Atmospheric Boundary Layer Using a New Subgrid-Scale Model – I. Slightly unstable and neutral cases

Subgrid-scale (SGS) modeling is a long-standing problem and a critical component in the large-eddy simulation (LES) of atmospheric boundary layer. A variety of SGS models with different levels of sophistication have been proposed for different needs, such as Smagorinsky's (1963) eddy viscosity model, Mason and Thomson's (1992) stochastic backscattering model, and Sullivan et al.'s (1994) near surface model. A modified Smagorinsky SGS model has been used in the LES version of Terminal Area Simulation System (TASS-LES). It has successfully simulated the buoyancy-dominated, convective atmospheric boundary layer flows, while simulations of the shear-dominated, slightly unstable, neutral, and stably stratified boundary layer flows are not so good. For the later, we used a simpler version of Sullivan et al.'s subgrid-scale model in which turbulent kinetic energy equation is not included and the model is still the first-order closure. A momentum profile matching approach is adopted in the proposed model. A series of simulations for shear-dominated, slightly unstable and neutral boundary layers are performed using different subgrid-scale models and different grid resolutions. The results are compared with those of Sullivan et al. (1994) and with empirical similarity relations for the surface layer. The simulations with the new SGS model appear to be far more satisfactory than those with the modified Smagorinsky model.