The role of subsidence in a weakly unstable marine boundary layer: a case study

The diurnal evolution of a cloud free, marine boundary layer is studied by means of experimental measurements and numerical simulations. Experimental data belong to an investigation of the mixing height over inner Danish waters. The mixed-layer height measured over the sea is generally nearly constant, and does not exhibit the diurnal cycle characteristic of boundary layers over land. A case study, during summer, showing an anomalous development of the mixed layer under unstable and nearly neutral atmospheric conditions, is selected in the campaign. Subsidence is identified as the main physical mechanism causing the sudden decrease in the mixing layer height. This is quantified by comparing radiosounding profiles with data from numerical simulations of a mesoscale model, and a large-eddy simulation model. Subsidence not only affects the mixing layer height, but also the turbulent fluctuations within it. By analyzing wind and scalar spectra, the role of subsidence is further investigated and a more complete interpretation of the experimental results emerges

Direct numerical simulation of fog: The sensitivity of a dissipation phase to environmental conditions

The sensitivity of fog dissipation to the environmental changes in radiation, liquid-water lapse rate, free tropospheric temperature and relative humidity was studied through numerical experiments designed based on the 2007-Paris Fog observations. In particular, we examine how much of the stratocumulus-thinning mechanism can be extended to the near-surface clouds or fog. When the free troposphere is warmed relative to the reference case, fog-top descends and become denser. Reducing the longwave radiative cooling via a more emissive free troposphere favors thickening the physical depth of fog, unlike cloud-thinning in a stratocumulus cloud. Drying the free troposphere allows fog thinning and promotes fog dissipation while sustaining the entrainment rate. The numerical simulation results suggest that the contribution of entrainment drying is more effective than the contribution of entrainment warming yielding the reduction in liquid water path tendency and promoting the onset of fog depletion relative to the reference case studied here. These sensitivity experiments indicate that the fog lifting mechanism can enhance the effect of the inward mixing at the fog top. However, to promote fog dissipation, an inward mixing mechanism only cannot facilitate removing humidity in the fog layer unless a sufficient entrainment rate is simultaneously sustained

Transition Delay in Hypervelocity Boundary Layers By Means of CO₂/Acoustic Instability Interaction

The potential for hypervelocity boundary layer stabilization was investigated using the concept of damping Mack’s second mode disturbances by vibrational relaxation of carbon dioxide (CO₂) within the boundary layer. Experiments were carried out in the Caltech T5 hypervelocity shock tunnel and the Caltech Mach 4 Ludwieg tube. The tests used 5-degree half-angle cones (at zero angle of attack) equipped near the front of the cone with an injector consisting of either discrete holes or a porous section. Gaseous CO₂, argon (Ar) and air were injected into the boundary layer and the effect on boundary layer stability was evaluated by optical visualization, heat flux measurements and numerical simulation. In T5, tests were carried out with CO₂ in the free stream as well as injection. Injection experiments in T5 were inconclusive; however, experiments with mixtures of air/CO₂ in the free stream demonstrated a clear stabilizing effect, limiting the predicted amplification N-factors to be less than 13. During the testing activities in T5, significant improvements were made in experimental technique and data analysis. Testing in the Ludwieg tube enabled optical visualization and the identification of a shear-layer like instability downstream of the injector. Experiments showed and numerical simulation confirmed that injection has a destabilizing influence beyond a critical level of injection mass flow rate. A modified injection geometry was tested in the Ludwieg tube and we demonstrated that it was possible to cancel the shock wave created by injection under carefully selected conditions. However, computations indicate and experiments demonstrate that shear-layer like flow downstream of the porous wall injector is unstable and can transition to turbulence while the injected gas is mixing with the free stream. We conclude that the idea of using vibrational relaxation to delay boundary layer transition is a sound concept but there are significant practical issues to be resolved to minimize the flow disturbance associated with introducing the vibrationally-active gas into the boundary layer

Characteristics of transitional multicomponent gaseous and drop-laden mixing layers from direct numerical simulation: Composition effects

Transitional states are obtained by exercising a model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer at larger Reynolds numbers, Re, than in a previous study. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete dynamic and thermodynamic coupling between phases is included. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass, is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of their PDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices, resulting into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions, and the effects of the initial mass loading and initial free-stream gas temperature are explored. For reference, simulations are also performed for gaseous multicomponent mixing layers for which the effect of Re is investigated in the direct-numerical-simulation–accessible regime. The results encompass examination of the global layer characteristics, flow visualizations, and homogeneous-plane statistics at transition. Comparisons are performed with previous pretransitional MC-liquid simulations and with transitional single-component (SC) liquid-drop-laden mixing layer studies. Contrasting to pretransitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF. Unlike in pretransitional situations, regions of large composition standard deviation no longer necessarily coincide with those of large mean molar mass. The rotational and composition characteristics are all liquid-specific and the variation among liquids is amplified with increasing free-stream gas temperature. The classical energy cascade is found to be of similar strength, but the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation for MC flows. The kinetic energy and dissipation are liquid-specific and the variation among liquids is amplified with increasing free-stream gas temperature. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from, a SGPDF. Eulerian and Lagrangian statistics of gas-phase quantities show that the different observation framework may affect the perception of the flow

Higher-level simulations of turbulent flows

The fundamentals of large eddy simulation are considered and the approaches to it are compared. Subgrid scale models and the development of models for the Reynolds-averaged equations are discussed as well as the use of full simulation in testing these models. Numerical methods used in simulating large eddies, the simulation of homogeneous flows, and results from full and large scale eddy simulations of such flows are examined. Free shear flows are considered with emphasis on the mixing layer and wake simulation. Wall-bounded flow (channel flow) and recent work on the boundary layer are also discussed. Applications of large eddy simulation and full simulation in meteorological and environmental contexts are included along with a look at the direction in which work is proceeding and what can be expected from higher-level simulation in the future

On compressibility assumptions in aeroacoustic integrals: a numerical study with subsonic mixing layers

Two assumptions commonly made in predictions based on Lighthill’s formalism are investigated: a constant density in the quadrupole expression, and the evaluation of the source quantity from incompressible simulations. Numerical predictions of the acoustic field are conducted in the case of a subsonic spatially evolving two-dimensional mixing layer at Re = 400. Published results of the direct noise computation (DNC) of the flow are use as reference and input for hybrid approaches before the assumptions on density are progressively introduced. Divergence free velocity fields are obtained from an incompressible simulation of the same flow case, exhibiting the same hydrodynamic field as the DNC. Fair comparisons of the hybrid predictions with the reference acoustic field valid both assumptions in the source region for the tested values of the Mach number. However, in the observer region, the inclusion of flow effects in the Lighthill source term is not preserved, which is illustrated through a comparison with the Kirchhoff wave-extrapolation formalism, and with the use of a convected Green function in the integration process

Transition Delay in Hypervelocity Boundary Layers By Means of CO₂/Acoustic Instability Interaction

The potential for hypervelocity boundary layer stabilization was investigated using the concept of damping Mack’s second mode disturbances by vibrational relaxation of carbon dioxide (CO₂) within the boundary layer. Experiments were carried out in the Caltech T5 hypervelocity shock tunnel and the Caltech Mach 4 Ludwieg tube. The tests used 5-degree half-angle cones (at zero angle of attack) equipped near the front of the cone with an injector consisting of either discrete holes or a porous section. Gaseous CO₂, argon (Ar) and air were injected into the boundary layer and the effect on boundary layer stability was evaluated by optical visualization, heat flux measurements and numerical simulation. In T5, tests were carried out with CO₂ in the free stream as well as injection. Injection experiments in T5 were inconclusive; however, experiments with mixtures of air/CO₂ in the free stream demonstrated a clear stabilizing effect, limiting the predicted amplification N-factors to be less than 13. During the testing activities in T5, significant improvements were made in experimental technique and data analysis. Testing in the Ludwieg tube enabled optical visualization and the identification of a shear-layer like instability downstream of the injector. Experiments showed and numerical simulation confirmed that injection has a destabilizing influence beyond a critical level of injection mass flow rate. A modified injection geometry was tested in the Ludwieg tube and we demonstrated that it was possible to cancel the shock wave created by injection under carefully selected conditions. However, computations indicate and experiments demonstrate that shear-layer like flow downstream of the porous wall injector is unstable and can transition to turbulence while the injected gas is mixing with the free stream. We conclude that the idea of using vibrational relaxation to delay boundary layer transition is a sound concept but there are significant practical issues to be resolved to minimize the flow disturbance associated with introducing the vibrationally-active gas into the boundary layer

Direct numerical simulation of incompressible spatially developing turbulent mixing layers

Turbulent mixing layers are a canonical free shear flow in which two parallel fluid streams of different velocities mix at their interface. Understanding spatial development of a turbulent mixing layer is essential for various engineering applications. However, multiple factors affect physics of this flow, making it difficult to reproduce results in experiments and simulations. The current study investigates sensitivity of direct numerical simulation (DNS) of such a flow to computational parameters. In particular, effects of the computational domain dimensions, grid refinement, thickness of the splitter plate, and the laminar boundary layer characteristics at the splitter plate trailing edge are considered. Flow conditions used in DNS are close to those from the experiments by Bell & Mehta (1990), where untripped boundary layers co-flowing on both sides of a splitter plate mix downstream of the plate. No artificial perturbations are used in simulations to trigger the flow transition to turbulence. DNS were conducted using the spectral-element method implemented in the open-source code Nek5000. Mean flow statistics are presented for the spatially developing self-similar flow, including high-order velocity moments. Such statistics will be used for validation of high-order Reynolds-Averaged Navier-Stokes (RANS) closure models