Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixing layer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Aerosol Vertical Turbulent Mass Flux Retrievals Through Novel Remote Sensing Algorithm

Abstract Integrated measurements of aerosol, radiation, cloud, and turbulent transport in the planetary boundary layer (PBL) are essential for understanding and modeling climate and air quality. Here, we developed a new technique for the identification of convective turbulent regions and deriving the vertical distribution of aerosol turbulent mass fluxes within PBL. The algorithm uses retrievals from coherent Doppler lidars and a high spectral resolution lidar. The technique was applied to study particle mass fluxes over 2 months (November–December 2020) during the campaign conducted at the DOE Atmospheric Radiation Measurement Southern Great Plains (SGP) site in Lamont, Oklahoma. The algorithm developed here is capable of continuously deriving vertically resolved (curtains) aerosol mass fluxes. Our data analysis shows that at the site, the 30‐min averaged fluxes at 135 m above the surface were mainly positive (upward) at ∼1 μg m −2 s −1 , suggesting that the surface is the primary source of the particle mass supplied to the boundary layer at the SGP site. Analyses of the individual case studies have revealed that not all the derived fluxes can be linked to surface emissions. Both positive and negative values in a range of ±5 μg m −2 s −1 can be caused by convective thermals interacting between the residual layer and the mixed layer and by rotation of the horizontal wind with the height. Large erroneous negative fluxes can also be caused by drizzling/precipitating clouds. We anticipate that the application of the current technique will lead to a more realistic representation of aerosol mass budgets and bidirectional mixing rates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multi‐Scale Temporal Variability of Turbulent Mixing During a Monsoon Intra‐Seasonal Oscillation in the Bay of Bengal: An LES Study

A process study using large‐eddy simulations is carried out to explore the dominant 1‐D processes that affect mixed layer (ML) properties during an event of summer Monsoon Intra‐seasonal Oscillations (MISO) in the Bay of Bengal (BOB). These simulations use realistic air‐sea fluxes and initial conditions that were collected during the summer 2018 MISO‐BOB field experiment to explore the roles of thermal inversion layer (TIL) and Langmuir turbulence (LT) in modulating ML properties. The simulations span an active period with heavy rain and strong winds and a break period with strong solar heat flux and little rain. The mixed layer depth (MLD), sea surface temperature (SST) and sea surface salinity (SSS) are most affected by the presence of near‐inertial oscillations, solar heating and precipitation, all of which occur at different timescales. The subsurface warming induced by the TIL reduces the SST variability at the MISO timescale when compared with the simulation without TIL. Comparison of simulations with and without LT indicates that LT enhances subsurface warming during the active phase and reduces diurnal SST modulation during the break phase. Finally, simulations with 1‐D mixing models show a wide disparity in the evolution of MLD, SST, and SSS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Rayleigh–Taylor Instability With Varying Periods of Zero Acceleration

We present our findings from a numerical investigation of the acceleration-driven Rayleigh–Taylor Instability, modulated by varying periods without an applied acceleration field. It is well known from studies on shock-driven Richtmyer–Meshkov instability that mixing without external forcing grows with a scaling exponent as ≈ t 0.20-0.28 When the Rayleigh–Taylor Instability is subjected to varying periods of “zero” acceleration, the structural changes to the mixing layer remain remarkably small. After the acceleration is re-applied, the mixing layer quickly resumes the profile of development it would have had if there had been no intermission. As a result, this behavior contrasts in particular with the strong sensitivity that is found to other variable acceleration profiles examined previously in the literature.

42 ENGINEERING↗
DOE OSTI2021

Self-consistent, high-order spatial profiles in a model for two-fluid turbulent mixing

A Reynolds-averaged Navier-Stokes model is presented with the property that it admits self-consistent, high-order spatial profiles in simulations of two-fluid turbulent mixing layers. Whereas previous models have been limited by the assumption of a linear mixing profile, the present paper relaxes this assumption and, as a result, is shown to achieve much better agreement with experimental profiles. Similarity analysis is presented to derive constraints on model coefficients to enforce desired self-similar growth rates that are fully consistent with the high-order spatial profiles. Through this similarity analysis, it is shown that care must be taken in model construction, as it is possible to construct certain terms in such a way as to leave growth rates unconstrained. This model, termed the k-Φ-L-a-V model, is then applied in simulations of Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz mixing layers. Overall, these simulations confirm that the expected growth parameters are recovered and high-order spatial profiles are maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Two self-similar Reynolds-stress transport models with anisotropic eddy viscosity

Two Reynolds-averaged Navier-Stokes models with full Reynolds-stress transport (RST) and tensor eddy viscosity are presented. These new models represent RST extensions of the $k−2L−a−\mathscr{C}$ and $k−ϕ−L−a−\mathscr{C}$ models by Morgan. Self-similarity analysis is used to derive constraints on model coefficients required to reproduce expected growth parameters for a variety of canonical flows, including Rayleigh-Taylor (RT) and Kelvin-Helmholtz (KH) mixing layers. Both models are then applied in one-dimensional simulation of RT and KH mixing layers, and the expected self-similar growth rates and anisotropy are obtained. Next, models are applied in two-dimensional simulation of the so-called “tilted rocket rig” inclined RT experiment and in simulation of a shock-accelerated localized patch of turbulence. Here it is found that RST is required to capture the qualitative growth of the shock-accelerated patch, and an anisotropic eddy viscosity provides substantial improvement over a Boussinesq treatment for the tilted rocket rig problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Role of Intense Arctic Storm in Accelerating Summer Sea Ice Melt: An In Situ Observational Study

Abstract Intense storms have been more frequently observed in the Arctic during recent years, in coincidence with extreme sea ice loss events. However, it is still not fully understood how storms drive such events due to deficient observations and modeling discrepancies. Here we address this problem by analyzing in situ observations acquired during an Arctic expedition, which uniquely captured an intense storm in summer 2016. The result shows a pronounced acceleration of sea ice loss during the storm process. Diagnostic analysis indicates a net energy loss at the ice surface, not supporting the accelerated melting. Although the open water surface gained net heat energy, it was insufficient to increase the mixed‐layer temperature to the observed values. Dynamic analysis suggests that storm‐driven increase in ocean mixing and upward Ekman pumping of the Pacific‐origin warm water tremendously increased oceanic heat flux. The thermal advection by the Ekman pumping led to a warmed mixed layer by 0.05°C–0.12°C and, in consequence, an increased basal sea ice melt rate by 0.1–1.7 cm day −1 .

58 GEOSCIENCES↗
DOE OSTI2021

Subtropical Eastern North Pacific SST bias in earth system models

Karthik Balaguru (PNNL), Luke Van Roekel (LLNL), L. Ruby Leung (PNNL) and Milena Veneziani (LLNL) Abstract: This study systematically evaluates the warm SST bias in the Subtropical Eastern North Pacific (SENP), a problem plaguing most CMIP6 models, using the Energy Exascale Earth System Model version 1 (E3SM). Compared to the model at its standard resolution (1o atmosphere, 30-60 km ocean), the SST bias improves in several regions at high resolution (0.25o atmosphere, 18-6 km ocean) except in subtropical eastern boundary regions, such as the SENP. A mixed layer heat budget analysis reveals that errors in the surface radiative fluxes and horizontal heat advection are primarily responsible for the bias. Analysis of HighResMIP models indicates that the lack of improvement in the SENP SST bias with resolution is not native to E3SM alone. A larger CMIP6 inter-model correlation of SST bias with biases in surface winds than radiation highlights the important role of large-scale circulation biases in the SENP SST bias in coupled models.

Balaguru, Karthik↗
DOE OSTI2022

Global decline in ocean memory over the 21st century

Ocean memory, the persistence of ocean conditions, is a major source of predictability in the climate system beyond weather time scales. We show that ocean memory, as measured by the year-to-year persistence of sea surface temperature anomalies, is projected to steadily decline in the coming decades over much of the globe. This global decline in ocean memory is predominantly driven by shoaling of the upper-ocean mixed layer depth in response to global surface warming, while thermodynamic and dynamic feedbacks can contribute substantially regionally. As the mixed layer depth shoals, stochastic forcing becomes more effective in driving sea surface temperature anomalies, increasing high-frequency noise at the expense of persistent signals. Reduced ocean memory results in shorter lead times of skillful persistence-based predictions of sea surface thermal conditions, which may present previously unknown challenges for predicting climate extremes and managing marine biological resources under climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Global ocean dimethyl sulfide climatology estimated from observations and an artificial neural network

Marine dimethyl sulfide (DMS) is important to climate due to the ability of DMS to alter Earth's radiation budget. Knowledge of the global-scale distribution, seasonal variability, and sea-to-air flux of DMS is needed in order to improve understanding of atmospheric sulfur, aerosol/cloud dynamics, and albedo. Here we examine the use of an artificial neural network (ANN) to extrapolate available DMS measurements to the global ocean and produce a global climatology with monthly temporal resolution. A global database of 82 996 ship-based DMS measurements in surface waters was used along with a suite of environmental parameters consisting of latitude–longitude coordinates, time of day, time of year, solar radiation, mixed layer depth, sea surface temperature, salinity, nitrate, phosphate, and silicate. Linear regressions of DMS against the environmental parameters show that on a global-scale mixed layer depth and solar radiation are the strongest predictors of DMS. These parameters capture ~9 % and ~7 % of the raw DMS data variance, respectively. Multilinear regression can capture more of the raw data variance (~39 %) but strongly underestimates DMS in high-concentration regions. In contrast, the artificial neural network captures ~66 % of the raw data variance in our database. Like prior climatologies our results show a strong seasonal cycle in surface ocean DMS with the highest concentrations and sea-to-air fluxes in the high-latitude summertime oceans. We estimate a lower global sea-to-air DMS flux (20.12±0.43 Tg S yr –1 ) than the prior estimate based on a map interpolation method when the same gas transfer velocity parameterization is used. Our sensitivity test results show that DMS concentration does not change unidirectionally with each of the environmental parameters, which emphasizes the interactions among these parameters. The ANN model suggests that the flux of DMS from the ocean to the atmosphere will increase with global warming. Given that larger DMS fluxes induce greater cloud albedo, this corresponds to a negative climate feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Interface models for three-dimensional Rayleigh–Taylor instability

We derive interface models for three-dimensional Rayleigh–Taylor instability (RTI), making use of a novel asymptotic expansion in the non-locality of the fluid flow. These interface models are derived for the purpose of studying universal features associated with RTI such as the Froude number in single-mode RTI, the predicted quadratic growth of the interface amplitude under multi-mode random perturbations, the optimal (viscous) mixing rates induced by the RTI and the self-similarity of horizontally averaged density profiles and the remarkable stabilization of the mixing layer growth rate which arises for the three-fluid two-interface heavy–light–heavy configuration, in which the addition of a third fluid bulk slows the growth of the mixing layer to a linear rate. Our interface models can capture the formation of small-scale structures induced by severe interface roll-up, reproduce experimental data in a number of different regimes and study the effects of multiple interface interactions even as the interface separation distance becomes exceedingly small. Compared with traditional numerical schemes used to study such phenomena, our models provide a computational speed-up of at least two orders of magnitude.

Mechanics↗
DOE OSTI2022

Investigating the Impacts of Daytime Boundary Layer Clouds on Surface Energy Fluxes and Boundary Layer Structure During CHEESEHEAD19

Studies of land-atmosphere interactions under a clear sky and low cumulus cloud conditions are common from long-term observatories like at the southern great plains. How well the relationships and responses of surface radiative and turbulent heat fluxes determined from these investigations hold for more heterogeneous surfaces in other climate regimes, however, is uncertain. In this study, detailed observations of the surface energy budget and daytime boundary layer properties are analyzed using measurements from the Chequamegon Heterogenous Ecosystem Energy-Balance Study Enabled by a High-Density Extensive Array of Detectors 2019 (CHEESEHEAD19) field campaign, July-October 2019, across a heterogeneous forested landscape of northern Wisconsin. A cloud regime framework is employed to classify consecutive periods of clear skies from lower atmosphere stratiform and cumulus clouds. A seasonal transition from low cumulus to low stratiform periods occurred, together with a diurnal pattern in cloudy or clear sky period dominance. Radiative forcing was highly dependent on sky conditions, leading to changes in the redistribution efficiency of radiative energy by the surface turbulent heat fluxes. During CHEESEHEAD19, small Bowen ratios dominated with daytime latent heat fluxes three times as large as sensible heat fluxes for all sky conditions studied; the forested region, therefore, falls within an energy-limited regime. The depth of the daytime mixed layer depended upon the sky condition and thermodynamic setting; deeper mixed layers occurred during periods of low cumulus and not clear skies. Profiles of vertical velocity were found to have enhanced variance under low cumulus compared to clear sky periods, suggesting potential for cloud feedbacks on boundary layer structure and surface energy fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Time-resolved particle image velocimetry measurements of the turbulent Richtmyer–Meshkov instability

Experiments are presented on the Richtmyer–Meshkov instability (RMI) with a three-dimensional, multi-mode initial perturbation. The experiments use a vertical shock tube, where a stably stratified interface is formed between air and sulphur hexafluoride (SF 6 ) via counterflow. A perturbation is imposed at the interface by vertical oscillation of the gas column, forming Faraday waves. The interface is accelerated by a Mach 1.17 (in air) shock wave, and the development of the mixing region between the gases is investigated using particle image velocimetry. Following shock acceleration, a reflected shock wave from the bottom of the shock tube interacts with the mixing layer a second time (reshock). The experiment is initialized with both high and low amplitude perturbations to examine the effect of the perturbation amplitude on measured quantities. The instability growth exponent (θ) is determined from the kinetic energy field using the width of the mixing layer and the decay of kinetic energy, which are found to be in agreement when the flow is most strongly excited. A growth exponent of θ ≈ 0.5 is found for all cases except the high-amplitude reshocked regime (where θ ≈ 0.33). High-amplitude experiments exhibit the transitional outer Reynolds number (Re ≡ hh˙/ν >10 4 ) required for mixing transition following the incident shock, and both experiments are elevated well above this threshold following reshock. However, neither set of experiments meet the more stringent requirements proposed by Zhou et al. (Phys. Rev.E, vol. 67, issue 5, 2003) which include the time dependent aspect of the RMI, an observation which is also made when examining the spectra.

42 ENGINEERING↗
DOE OSTI2020

The αs and θs in Rayleigh–Taylor and Richtmyer–Meshkov instabilities

The alphas (α b , α s ) and thetas (θ b , θ s ) are the essential parameters describing the asymptotic growth behavior of the mixing layer driven by Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities at the interface between the fluids. Extensive simulations and experiments have been conducted in the determination of these parameters. The results depend on numerical algorithms, resolutions, and the accuracy of experimental measurements. In this article, we present a unique relationship between the alphas and thetas by applying a simple buoyancy-drag mixing model. Our first main result is that, for the same fluids, these parameters are not independent because they are fundamentally connected by physics principles and they all are functions of the drag coefficients of the fluids. As a consequence, the alphas and thetas are intrinsically coupled. All are functions of any single one of these. Additionally, our second main result is that the RT bubble growth rate (α b ) is determined by the bubble competition process – a merger model – at the edge of the mixing layer. Together, these two main results determine all the alphas and thetas by strictly theoretical reasoning. Applications to both hydro and laser driven experiments at OMEGA are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A Scaling Analysis of Post-Detonation Mixing With Detailed Chemical Kinetics

This report documents a scaling study of hemispherical pentaerythritol tetranitrate explosives detonated in ambient air. The focus is on the mixing layer and evolution of the chemical species. The HyBurn code developed at the University of Florida is used to simulate center initiated 12, 38, 120, and 1200 mm diameter charges. The results show blast structures commonly found in the literature. Azimuthal averaging is used to radially average the flow variables to get mean values. The charges exhibit self-similar behavior at early times in a primarily gas dynamic regime, but discrepancies between the scales form at later times due to the finite rate chemistry. The turbulent mixing layer is calculated and used to define a region to conditionally average by mixture fraction. The results of the conditional averages imply a transition from a kinetic-limited regime to a mixing-limited regime with increasing scale. The scaling of the total heat release rate confirms an early time gas dynamic regime, and a later time regime where mixing and chemistry becomes impactful.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor mixing induced by power-law accelerations in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing turbulent Rayleigh–Taylor mixing driven by a temporal power-law acceleration are derived in the small Atwood number (Boussinesq) limit. The solutions generalize those previously derived for constant acceleration Rayleigh–Taylor mixing for models based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ [O. Schilling, “Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit,” Phys. Fluids 33, 085129 (2021)]. The turbulent fields are expressed in terms of the model coefficients and power-law exponent, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. Mixing layer growth parameters and other physical observables are obtained explicitly as functions of the model coefficients and parameterized by the exponent of the power-law acceleration. Values for physical observables in the constant acceleration case are used to calibrate the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to a canonical (i.e., constant acceleration) Rayleigh–Taylor turbulent flow. The calibrated four-equation model is then used to numerically reconstruct the mean and turbulent fields, and turbulent equation budgets across the mixing layer for several values of the power-law exponent. Finally, the reference solutions derived here can be used to understand the model predictions for strongly accelerated or decelerated Rayleigh–Taylor mixing in the large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2021

Morning boundary layer conditions for shallow to deep convective cloud evolution during the dry season in the central Amazon

Observations of the boundary layer (BL) processes are analyzed statistically for dry seasons of 2 years and in detail, as case studies, for 4 shallow convective days (ShCu) and 4 shallow-to-deep convective days (ShDeep) using a suite of ground-based measurements from the Observation and Modeling of the Green Ocean Amazon (GoAmazon 2014/5) Experiment. The BL stages in ShDeep days, from the nighttime to the cloudy mixing layer stage, are then described in comparison with ShCu days. Atmospheric thermodynamics and dynamics, environmental profiles, and surface turbulent fluxes were employed to compare these two distinct situations for each stage of the BL evolution. Particular attention is given to the morning transition stage, in which the BL changes from stable to unstable conditions in the early morning hours. Results show that the decrease in time duration of the morning transition on ShDeep days is associated with high humidity and well-established vertical wind shear patterns. Higher humidity since nighttime not only contributes to lowering the cloud base during the rapid growth of the BL but also contributes to the balance between radiative cooling and turbulent mixing during nighttime, resulting in higher sensible heat flux in the early morning. The sensible heat flux promotes rapid growth of the well-mixed layer, thus favoring the deeper BL starting from around 08:00 LST (UTC-4 h). Under these conditions, the time duration of morning transition is used to promote convection, having an important effect on the convective BL strength and leading to the formation of shallow cumulus clouds and their subsequent evolution into deep convective clouds. Statistical analysis was used to validate the conceptual model obtained from the case studies. Despite the case-to-case variability, the statistical analyses of the processes in the BL show that the described processes are very representative of cloud evolution during the dry season.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

Bo, Tong (ORCID:0000000260300561)↗
DOE OSTI2024

Simulations of three-layer Richtmyer–Meshkov mixing in a shock tube

Here, the Richtmyer–Meshkov instability causes perturbations to grow after a shock traverses a fluid density interface. This increases the mixing rate between fluid from either side of the interface. We use the Flash Eulerian hydrodynamic code to investigate alterations when a thin third layer of intermediate density is placed along the interface, effectively creating two adjacent unstable interfaces. This is a common occurrence in engineering applications where a thin barrier initially separates two materials. We find that the width of the mixing layer is similar or slightly reduced; however, the total mass of mixed material can actually increase. The mixing layer becomes more compact and efficient. However, the normalized mixed mass decreases, meaning that finger entrainment becomes more important than in the simple two-layer case. The effect of adding the central layer appears to decrease when the Atwood number is decreased. The Flash results are also benchmarked against two-layer experimental data from a shock tube at the University of Arizona.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗