Juan Pedro Mellado

The evaporatively driven cloud-top mixing layer

Direct numerical simulations of the turbulent temporally evolving cloud-top mixing layer are used to investigate the role of evaporative cooling by isobaric mixing locally at the stratocumulus top. It is shown that the system develops a horizontal layered structure whose evolution is determined by molecular transport. A relatively thin inversion with a constant thickness h =κ/ω e is formed on top and travels upwards at a mean velocity ω e ≃0.1(κ|b s |χ 2 c ) 1/3 , where κ is the mixture-fraction diffusivity, b s < 0 is the buoyancy anomaly at saturation conditions χ s and χ c is the cross-over mixture fraction defining the interval of buoyancy reversing mixtures. A turbulent convection layer develops below and continuously broadens into the cloud (the lower saturated fluid). This turbulent layer approaches a self-preserving state that is characterized by the convection scales constructed from a constant reference buoyancy flux B, = |b s |ω e /χ s. Right underneath the inversion base, a transition or buffer zone is defined based on a strong local conversion of vertical to horizontal motion that leads to a cellular pattern and sheet-like plumes, as observed in cloud measurements and reported in other free-convection problems. The fluctuating saturation surface (instantaneous cloud top) is contained inside this intermediate region. Results show that the inversion is not broken due to the turbulent convection generated by the evaporative cooling, and the upward mean entrainment velocity ω e is negligibly small compared to the convection velocity scale ω* of the turbulent layer and the corresponding growth rate into the cloud.

Global Intermittency and Collapsing Turbulence in the Stratified Planetary Boundary Layer

Direct numerical simulation of the turbulent Ekman layer over a smooth wall is used to investigate bulk properties of a planetary boundary layer under stable stratification. Our simplified configuration depends on two non-dimensional parameters: a Richardson number characterizing the stratification and a Reynolds number characterizing the turbulence scale separation. This simplified configuration is sufficient to reproduce global intermittency, a turbulence collapse, and the decoupling of the surface from the outer region of the boundary layer. Global intermittency appears even in the absence of local perturbations at the surface; the only requirement is that large-scale structures several times wider than the boundary-layer height have enough space to develop. Analysis of the mean velocity, turbulence kinetic energy, and external intermittency is used to investigate the large-scale structures and corresponding differences between stably stratified Ekman flow and channel flow. Both configurations show a similar transition to the turbulence collapse, overshoot of turbulence kinetic energy, and spectral properties. Differences in the outer region resulting from the rotation of the system lead, however, to the generation of enstrophy in the non-turbulent patches of the Ekman flow. The coefficient of the stability correction function from Monin–Obukhov similarity theory is estimated as

Gradient trajectory analysis of a scalar field with external intermittency

\beta \approx 5.7

Large-eddy simulation of Rayleigh-Taylor turbulence with compressible miscible fluids

β≈5.7 in agreement with atmospheric observations, theoretical considerations, and results from stably stratified channel flows. Our results demonstrate the applicability of this set-up to atmospheric problems despite the intermediate Reynolds number achieved in our simulations.

Direct Numerical Simulations of a Smoke Cloud-Top Mixing Layer as a Model for Stratocumuli

Large-eddy simulation of combustion problems involves highly nonlinear terms that, when filtered, result in a contribution from subgrid fluctuations of scalars, Z, to the dynamics of the filtered value. This subgrid contribution requires modeling. Reconstruction models try to recover as much information as possible from the resolved field Z, based on a deconvolution procedure to obtain an intermediate field ZM. The approximate reconstruction using moments (ARM) method combines approximate reconstruction, a purely mathematical procedure, with additional physics-based information required to match specific scalar moments, in the simplest case, the Reynolds-averaged value of the subgrid variance. Here, results from the analysis of the ARM model in the case of a spatially evolving turbulent plane jet are presented. A priori and a posteriori evaluations using data from direct numerical simulation are carried out. The nonlinearities considered are representative of reacting flows: power functions, the dependence of the density on the mixture fraction (relevant for conserved scalar approaches) and the Arrhenius nonlinearity (very localized in Z space). Comparisons are made against the more popular beta probability density function (PDF) approach in the a priori analysis, trying to define ranges of validity for each approach. The results show that the ARM model is able to capture the subgrid part of the variance accurately over a wide range of filter sizes and performs well for the different nonlinearities, giving uniformly better predictions than the beta PDF for the polynomial case. In the case of the density and Arrhenius nonlinearities, the relative performance of the ARM and traditional PDF approaches depends on the size of the subgrid variance with respect to a characteristic scale of each function. Furthermore, the sources of error associated with the ARM method are considered and analytical bounds on that error are obtained.

Mixing Driven by Radiative and Evaporative Cooling at the Stratocumulus Top

The passive scalar field of a temporally evolving shear layer is investigated using gradient trajectories as a means to analyse the scalar probability density function and the conditional scalar dissipation rate in the presence of external intermittency. These results are of significance for turbulent combustion, where improved predictions of the statistics of the conditional dissipation rate are needed in several models. First, the variation of the conventional first and second moments of the conditional dissipation rate across the layer is quantitatively documented in detail. A strong dependence of the conditional dissipation rate on the lateral position and on the conditioning value of the scalar is observed. The dependence on the transverse distance to the centre-plane partially explains the double-hump profile usually reported when this dependence is ignored. The variation with the scalar observed in the ratio between the second and first moments would invalidate certain assumptions commonly done in turbulent combustion. It is also seen that conditioning on the scalar does not reduce the fluctuation of the dissipation rate with respect to unconditional values. Next, the role of external intermittency in these results is investigated. For that purpose, the flow is partitioned into different zones based on different types of gradient trajectories passing through each point, thereby introducing non-local information in comparison with the standard turbulent/non-turbulent separation based on the conventional intermittency function. In addition to the homogeneous outer regions, three zones are identified: a turbulent zone, a turbulence interface and quasi-laminar diffusion layers. The relative contribution from each of these zones to the conventional intermittency factor is reported. The statistics are then conditioned on each of these zones, and the spatial variation of the scalar distribution and of the conditional scalar dissipation rate is explained in terms of the observed zonal statistics. For the Reynolds numbers of the present simulation, between 1500 and 3000 based on the vorticity thickness and the velocity difference, and a Schmidt number equal to 1, it results that the major contribution to both statistics is due to the turbulence interfaces. At the same time, the turbulent zone shows a distinct behaviour, being approximately homogeneous but anisotropic.

Wind shear and buoyancy reversal at the stratocumulus top

The entrainment zone (EZ) of a dry, shear-free convective boundary layer growing into a linearly stratified fluid is studied by means of direct numerical simulation. The scale separation between the boundary layer thickness and the Kolmogorov length scale is shown to be sufficient to observe Reynolds number similarity in the statistics of interest during the equilibrium entrainment regime. Contrary to previous considerations, the vertical structure of the entrainment zone is found to be better described by the superposition of two sublayers: 1) an upper EZ sublayer that is dominated by overshooting thermals and is characterized by a penetration depth that scales with the ratio of the convective velocity and the buoyancy frequency of the free troposphereand2)alowerEZsublayerthatisdominatedbytroughsofmixedfluidandischaracterizedbythe integral length scale of the mixed layer. Correspondingly, different buoyancy scales are identified. The consequences of this multiplicity of scales on the entrainment rate parameters are evaluated directly, without resorting to any bulk model, through an exact relation among the mean entrainment rate, the local buoyancy increment, and both the turbulent and the finite-thickness contributions to the entrainment ratio A measured at the height of minimum buoyancy flux. The smaller turbulent contribution to A that is usually observed for relatively thick EZs is found to be compensated by the smaller local buoyancy increment instead of by the finite-thickness contribution. The two-layer structure of the entrainment zone is found to affect the exponent of the power-law relation between the normalized mean entrainment rate and the convective Richardson number such that the exponent deviates from 21 for typical atmospheric conditions, although it asymptotically approaches 21 for higher Richardson numbers.

Scaling laws for the heterogeneously heated free convective boundary layer

Turbulence developed from Rayleigh-Taylor instability between two compressible miscible fluids in an unbounded domain is addressed in this paper. It is demonstrated that the turbulent Mach number in the turbulent core has an upper bound, independent of the density ratio under a broad range of initial mean configurations. The initial thermodynamic state of the system determines the amount of potential energy per unit mass involved in the turbulent mixing stage, and thus the characteristic level of turbulent fluctuations that is achievable is linked to the characteristic speed of sound such that the turbulent Mach number is limited. For the particular case of an ideal gas, this bound on the turbulent Mach number is found to be between 0.25 and 0.6, depending on the particular initial thermodynamic state. Hence, intrinsic compressibility effects (those owing to large Mach number) are likely to be limited in the turbulent stage of a pure Rayleigh-Taylor problem. This result is confirmed by large-eddy simulati...