Chin-Hoh Moeng

A Large-Eddy-Simulation Model for the Study of Planetary Boundary-Layer Turbulence

Abstract A large-eddy-simulation (LFS) model explicitly calculates the large-eddy field and parameterizes the small eddies. The large eddies in the atmospheric boundary layer are believed to be much more important and more flow-dependent than the small eddies. The LES model results are therefore believed to be relatively insensitive to the parameterization scheme for the small eddies. Deardorff first applied this type of numerical model to boundary-layer turbulence. In order to continue his important work, and to take advantage of the fast Fourier transformation algorithm, a new LES model code which uses a mixed pseudospectral finite-difference method was developed. This LES model is described here and tested with a simple vortex flow and with the Wangara day-33 data. This model will be used to systematically investigate fundamental problems in the area of boundary-layer turbulence. It is hoped that three-dimensional simulations will give useful statistical information about turbulence structural and impr...

A comparison of shear- and buoyancy-driven planetary boundary layer flows

Abstract Planetary boundary layer (PBL) flows are known to exhibit fundamental differences depending on the relative combination of wind shear and buoyancy forces. These differences are not unexpected in that shear instabilities occur locally, while buoyancy force sets up vigorous thermals, which result in nonlocal transport of heat and momentum. At the same time, these two forces can act together to modify the flow field. In this study, four large-eddy simulations (LESs) spanning the shear and buoyancy flow regimes were generated; two correspond to the extreme cases of shear and buoyancy-driven PBLs, while the other two represent intermediate PBLs where both forces are important. The extreme cases are used to highlight and quantify the basic differences between shear and convective PBLs in 1) flow structures, 2) overall statistics, and 3) turbulent kinetic energy (TKE) budget distributions. Results from the two intermediate LES cases are used to develop and verify a velocity scaling and a TKE budget mode...

Langmuir turbulence in the ocean

Solutions are analysed from large-eddy simulations of the phase-averaged equations for oceanic currents in the surface planetary boundary layer (PBL), where the averaging is over high-frequency surface gravity waves. These equations have additional terms proportional to the Lagrangian Stokes drift of the waves, including vortex and Coriolis forces and tracer advection. For the wind-driven PBL, the turbulent Langmuir number, La tur = ( U ∗/ U s ) 1/2 , measures the relative influences of directly wind-driven shear (with friction velocity U ∗) and the Stokes drift U s . We focus on equilibrium solutions with steady, aligned wind and waves and a realistic La tur = 0.3. The mean current has an Eulerian volume transport to the right of the wind and against the Stokes drift. The turbulent vertical fluxes of momentum and tracers are enhanced by the presence of the Stokes drift, as are the turbulent kinetic energy and its dissipation and the skewness of vertical velocity. The dominant coherent structure in the turbulence is a Langmuir cell, which has its strongest vorticity aligned longitudinally (with the wind and waves) and intensified near the surface on the scale of the Stokes drift profile. Associated with this are down-wind surface convergence zones connected to interior circulations whose horizontal divergence axis is rotated about 45° to the right of the wind. The horizontal scale of the Langmuir cells expands with depth, and there are also intense motions on a scale finer than the dominant cells very near the surface. In a turbulent PBL, Langmuir cells have irregular patterns with finite correlation scales in space and time, and they undergo occasional mergers in the vicinity of Y-junctions between convergence zones.

Evaluation of Large-Eddy Simulations via Observations of Nocturnal Marine Stratocumulus

Data from the first research flight (RF01) of the second Dynamics and Chemistry of Marine Stratocumulus (DYCOMS-II) field study are used to evaluate the fidelity with which large-eddy simulations (LESs) can represent the turbulent structure of stratocumulus-topped boundary layers. The initial data and forcings for this case placed it in an interesting part of parameter space, near the boundary where cloud-top mixing is thought to render the cloud layer unstable on the one hand, or tending toward a decoupled structure on the other hand. The basis of this evaluation consists of sixteen 4-h simulations from 10 modeling centers over grids whose vertical spacing wa s5ma t thecloud-top interface and whose horizontal spacing was 35 m. Extensive sensitivity studies of both the configuration of the case and the numerical setup also enhanced the analysis. Overall it was found that (i) if efforts are made to reduce spurious mixing at cloud top, either by refining the vertical grid or limiting the effects of the subgrid model in this region, then the observed turbulent and thermodynamic structure of the layer can be reproduced with some fidelity; (ii) the base, or native configuration of most simulations greatly overestimated mixing at cloud top, tending toward a decoupled layer in which cloud liquid water path and turbulent intensities were grossly underestimated; (iii) the sensitivity of the simulations to the representation of mixing at cloud top is, to a certain extent, amplified by particulars of this case. Overall the results suggest that the use of LESs to map out the behavior of the stratocumulus-topped boundary layer in this interesting region of parameter space requires a more compelling representation of processes at cloud top. In the absence of significant leaps in the understanding of subgrid-scale (SGS) physics, such a representation can only be achieved by a significant refinement in resolution—a refinement that, while conceivable given existing resources, is probably still beyond the reach of most centers.

A subgrid-scale model for large-eddy simulation of planetary boundary-layer flows

A long-standing problem in large-eddy simulations (LES) of the planetary boundary layer (PBL) is that the mean wind and temperature profiles differ from the Monin-Obukhov similarity forms in the surface layer. This shortcoming of LES has been attributed to poor grid resolution and inadequate sub-grid-scale (SGS) modeling. We study this deficiency in PBL LES solutions calculated over a range of shear and buoyancy forcing conditions. The discrepancy from similarity forms becomes larger with increasing shear and smaller buoyancy forcing, and persists even with substantial horizontal grid refinement. With strong buoyancy forcing, however, the error is negligible.In order to achieve better agreement between LES and similarity forms in the surface layer, a two-part SGS eddy-viscosity model is proposed. The model preserves the usual SGS turbulent kinetic energy formulation for the SGS eddy viscosity, but it explicitly includes a contribution from the mean flow and a reduction of the contributions from the turbulent fluctuations near the surface. Solutions with the new model yield increased fluctuation amplitudes near the surface and better correspondence with similarity forms out to a distance of 0.1–0.2 times the PBL depth, i.e., a typical surface-layer depth. These results are also found to be independent of grid anisotropy. The new model is simple to implement and computationally inexpensive.

Eddy Diffusivity and Countergradient Transport in the Convective Atmospheric Boundary Layer

Abstract To describe the heat and scalar fluxes in the convective boundary layer, we propose expressions for eddy diffusivities and countergradient terms. The latter expressions can be used in a modified flux-gradient approach, which takes account for nonlocal convective vertical exchange. The results for heat are based on a derivation similar to that of Deardorff by utilizing the turbulent heat-flux equation, but the closure assumptions applied to the heat-flux budgets are different. As a result, the physical interpretation for the countergradient term differs; our countergradient term results from the third-moment transport effect, while Deardorffs results from the buoyancy production term. On the basis of our analysis, we are able to calculate an eddy diffusivity for heat, using large-eddy simulation results. The results are presented in the form of a similarity profile, using the convective velocity scale w* and the inversion height zi. It is shown that the latter profile is well behaved and that it ...

Evaluation of Turbulent Transport and Dissipation Closures in Second-Order Modeling

Abstract We show that the turbulence statistics from our (96)3 large-eddy-simulation (LES) studies of a convective boundary layer are in excellent agreement with those from the Deardorff–Willis laboratory convection tank. Using these LES data, we evaluate contemporary parameterizations for turbulent transport and dissipation in second-order closure models of the convective boundary layer. The gradient-diffusion parameterization for turbulent transport fares poorly, due in large part to the direct influence of buoyancy. This leads to poor predictions of the vertical profiles of some turbulence statistics. We also find that the characteristic length scales for the mechanical and thermal dissipation rates typically used in second-order closure models are a factor of 2–3 too small; this leads to underpredictions of turbulence kinetic energy levels. Finally, we find that the flux and variance budgets for conservative scalars are substantially different in top-down and bottom-up diffusion. In order to reproduce...

Structure of the Entrainment Zone Capping the Convective Atmospheric Boundary Layer

The authors use large-eddy simulation (LES) to investigate entrainment and structure of the inversion layer of a clear convectively driven planetary boundary layer (PBL) over a range of bulk Richardson numbers, Ri. The LES code uses a nested grid technique to achieve fine resolution in all three directions in the inversion layer. Extensive flow visualization is used to examine the structure of the inversion layer and to illustrate the temporal and spatial interaction of a thermal plume and the overlying inversion. It is found that coherent structures in the convective PBL, that is, thermal plumes, are primary instigators of entrainment in the Ri range 13.6 # Ri # 43.8. At Ri 5 13.6, strong horizontal and downward velocities are generated near the inversion layer because of the plume‐interface interaction. This leads to folding of the interface and hence entrainment of warm inversion air at the plume’s edge. At Ri 5 34.5, the inversion’s strong stability prevents folding of the interface but strong horizontal and downward motions near the plume’s edge pull down pockets of warm air below the nominal inversion height. These pockets of warm air are then scoured off by turbulent motions and entrained into the PBL. The structure of the inversion interface from LES is in good visual agreement with lidar measurements in the PBL obtained during the Lidars in Flat Terrain field experiment. A quadrant analysis of the buoyancy flux shows that net entrainment flux (or average minimum buoyancy flux wu min) is identified with quadrant IV w2u1 , 0 motions, that is, warm air moving downward. Plumes generate both large negative quadrant II w1u2 , 0 and positive quadrant III w2u2 . 0 buoyancy fluxes that tend to cancel. The maximum vertical gradient in potential temperature at every (x, y) grid point is used to define a local PBL height, zi(x, y). A statistical analysis of zi shows that skewness of zi depends on the inversion strength. Spectra of zi exhibit a sensitivity to grid resolution. The normalized entrainment rate we/w * , where we and w * are entrainment and convective velocities, varies as ARi21 with A 0.2 in the range 13.6 # Ri # 43.8 and is in good agreement with convection tank measurements. For a clear convective PBL, the authors found that the finite thickness of the inversion layer needs to be considered in an entrainment rate parameterization derived from a jump condition.

Large‐eddy simulation of the diurnal cycle of shallow cumulus convection over land

SUMMARY Large-eddy simulations of the development of shallow cumulus convection over land are presented. Many characteristics of the cumulus layer previously found in simulations of quasi-steady convection over the sea are found to be reproduced in this more strongly forced, unsteady case. Furthermore, the results are shown to be encouragingly robust, with similar results obtained with eight independent models, and also across a range of numerical resolutions. The datasets produced are already being used in the development and evaluation of parametrizations used in numerical weather-prediction and climate models.

Spectral Analysis of Large-Eddy Simulations of the Convective Boundary Layer

Abstract To examine the fidelity of the simulated turbulent flow, we analyzed the spectra of 96×&96×96 large-eddy-simulation results. We derived expressions for the inertial-range spectra of a filtered wind field, compared our computed spectra with the theoretical predictions, and drew two main conclusions. First, the wave cutoff filter is more appropriate for our large-eddy-simulation model than the Gaussian filter. Second, only a certain combination of the subgrid-scale parameters for the dissipation rate and eddy viscosity provides good inertial-range spectra. We offer an explanation for the unreasonably large sugrid temperature and moisture variances reported in Deardorffs 1974 lame-eddy-simulation results, and show that moment statistics up to third order are not sensitive to moderate changes in the subgrid-scale parameters.