ALBERT

Description: Thermal convection of a Boussinesq fluid in a layer confined between two infinite horizontal walls is investigated by direct numerical simulation (DNS) and by large-eddy simulation (LES) for zero horizontal mean motion. The lower-surface height varies sinusoidally in one horizontal direction while remaining constant in the other. Several cases are considered with amplitude δ up to 0.15H and wavelength Λ of H to 8H (inclination up to 43°), where H is the mean fluid-layer height. Constant heat flux is prescribed at the lower surface of the initially at rest and isothermal fluid layer. In the LES, the surface is treated as rough surface (z0 = 10-4H) using the Monin-Oboukhov relationships. At the flat top an adiabatic frictionless boundary condition is applied which approximates a strong capping inversion of an atmospheric convective boundary layer. In both horizontal directions, the model domain extends over the same length (either 4H or 8H) with periodic lateral boundary conditions. We compare DNS of moderate turbulence (Reynolds number based on H and on the convective velocity is 100, Prandtl number is 0.7) with LES of the fully developed turbulent state in terms of turbulence statistics and characteristic large-scale-motion structures. The LES results for a flat surface generally agree well with the measurements of Adrian et al. (1986). The gross features of the flow statistics, such as profiles of turbulence variances and fluxes, are found to be not very sensitive to the variations of wavelength, amplitude, domain size and resolution and even the model type (DNS or LES), whereas details of the flow structure are changed considerably. The LES shows more turbulent structures and larger horizontal scales than the DNS. To a weak degree, the orography enforces rolls with axes both perpendicular and parallel to the wave crests and with horizontal wavelengths of about 2H to 4H. The orography has the largest effect for Λ = 4H in the LES and for Λ = 2H in the DNS. The results change little when the size of the computational domain is doubled in both horizontal directions. Most of the motion energy is contained in the large-scale structures and these structures are persistent in time over periods of several convective time units. The motion structure persists considerably longer over wavy terrain than over flat surfaces. © 1992, Cambridge University Press. All rights reserved.

Description: Turbulence in the convective boundary layer (CBL) uniformly heated from below and topped by a layer of uniformly stratified fluid is investigated for zero mean horizontal flow using large-eddy simulations (LES). The Rayleigh number is effectively infinite, the Froude number of the stable layer is 0.09 and the surface roughness height relative to the height of the convective layer is varied between 10−6 and 10−2. The LES uses a finite-difference method to integrate the three-dimensional grid-volume-averaged Navier–Stokes equations for a Boussinesq fluid. Subgrid-scale (SGS) fluxes are determined from algebraically approximated second-order closure (SOC) transport equations for which all essential coefficients are determined from the inertial-range theory. The surface boundary condition uses the Monin–Obukhov relationships. A radiation boundary condition at the top of the computational domain prevents spurious reflections of gravity waves. The simulation uses 160 × 160 × 48 grid cells. In the asymptotic state, the results in terms of vertical mean profiles of turbulence statistics generally agree very well with results available from laboratory and atmospheric field experiments. We found less agreement with respect to horizontal velocity fluctuations, pressure fluctuations and dissipation rates, which previous investigations tend to overestimate. Horizontal spectra exhibit an inertial subrange. The entrainment heat flux at the top of the CBL is carried by cold updraughts and warm downdraughts in the form of wisps at scales comparable with the height of the boundary layer. Plots of instantaneous flow fields show a spoke pattern in the lower quarter of the CBL which feeds large-scale updraughts penetrating into the stable layer aloft. The spoke pattern has also been found in a few previous investigations. Small-scale plumes near the surface and remote from strong updraughts do not merge together but decay while rising through large-scale downdraughts. The structure of updraughts and downdraughts is identified by three-dimensional correlation functions and conditionally averaged fields. The mean circulation extends vertically over the whole boundary layer. We find that updraughts are composed of quasi-steady large-scale plumes together with transient rising thermals which grow in size by lateral entrainment. The skewness of the vertical velocity fluctuations is generally positive but becomes negative in the lowest mesh cells when the dissipation rate exceeds the production rate due to buoyancy near the surface, as is the case for very rough surfaces. The LES results are used to determine the root-mean-square value of the surface friction velocity and the mean temperature difference between the surface and the mixed layer as a function of the roughness height. The results corroborate a simple model of the heat transfer in the surface layer. © 1989, Cambridge University Press. All rights reserved.