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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effect of Ice Number Concentration on the Evolution of Boundary Layer Clouds During Arctic Marine Cold‐Air Outbreaks

Abstract Marine cold‐air outbreaks (MCAOs) are crucial for Arctic Ocean heat loss, featuring convective cloud rolls that transition into convection cells downstream. Understanding factors controlling this transformation is the key for improving MCAO cloud representation in climate models. This study employs large‐eddy simulations to investigate how cloud ice number concentrations () affect cloud evolution using a case from the Cold‐Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) campaign. The simulations, performed in a Lagrangian framework following an air mass trajectory, are driven by ERA5 reanalysis data. Initially, all simulations produce similar cloud patterns, but higher leads to earlier breakup of cloud rolls. Between 4 and 10 hr, surface precipitation rates are similar across simulations, but precipitation initiates earlier, and the cloud‐base precipitation rates are higher when is higher. The stronger precipitation evaporation leads to increased stability of the boundary layer and reduced intensity of vertical mixing between the surface and cloud layer. An increased sink of cloud layer moisture via precipitation and decreased source through diminished vertical transport result in earlier cloud breakup in higher conditions. Simulations with different sea surface temperatures (SST) indicate that this cloud breakup mechanism remains valid for MCAOs of different strengths, although the cloud organization is more sensitive to SST changes in low environments. This work highlights the importance of accurate representations of ice processes in simulating MCAO clouds and suggests the need for observational constraints of ice nucleating particles and over the mixed‐phase cloud regimes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Stabilization of Miscible Aqueous Phases via Diffusion‐Controlled Multifunctional Nanoparticle‐Ligand Complexation

Liquid‐in‐liquid structuring by harnessing miscible aqueous domains can be achieved by inducing thermodynamically defined phase separation using aqueous solutions of incompatible solutes, where immiscibility is dictated by the concentrations of the solutes. For instance, in aqueous two‐phase systems (ATPSs), the interfacial tension between the two different aqueous solutions is so small that it stabilizes the phase‐separated domains, but not their shape, thus falling short of designing tunable, robust structures. Here, we introduce a diffusion‐controlled strategy that enables liquid‐in‐liquid compartmentalization absent a thermodynamically defined interface or bulk liquid‐liquid phase separation using a barrier formed in situ at the initial contact boundary of an aqueous dispersion of multifunctional nanoparticles and a solution of multifunctional ligands in pure water or a water/alcohol mixture. The nanoparticles and ligands are entropically driven to disperse uniformly but interact at the initial contact boundary more rapidly than bulk mixing, forming a layer that establishes a barrier akin to a semi‐permeable membrane and allows the formation of tubular compartments. A wealth of opportunities is enabled by the proposed technique, where examples in separations, soft electronic devices, and ultralight porous electromagnetic shielding materials are shown.

cryogels↗
DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Molecular simulation using transfer-learned potentials for the disordered nanoscale structure of nitrogen-doped nanoporous carbons

Machine learning (ML)-based molecular dynamics (MD) simulations of the formation of a class of N-doped nanoporous carbons are performed to assess their disordered partially graphitized nanoscale structure. The study is motivated by the effectiveness of so-called nitrogen assembly carbons (NACs) for catalysis applications. Benchmark simulations for pure-C disordered graphitic systems reveal the importance of reliably capturing the vdW component of the potentials in order to accurately describe the tendency for layering of disordered graphene-like sheets. In our modeling, this is achieved by a transfer learning strategy incorporating features of the energetics from the optB88-vdW DFT functional into potentials initially trained with a less expensive functional, thereby providing a superior description of the pure-C systems. Generation from MD simulations of realistic partially graphitized structures is significantly more challenging for N-doped versus for pure C systems. However, such structures are achieved by a tailored MD simulation protocol mimicking the experimental synthesis process and in particular incorporating an annealing and subsequent quenching stages. Simulated PXRD patterns effectively reproduce the features of experimental observations for NACs, including the appearance of a prominent but broad (002) peak at around 25, and the development of another weaker feature associated with in-layer ordering of mixed C-N graphene-like sheets.

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

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

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

Phase stabilization and oxidation of a continuous composition spread multi-principal element (AlFeNiTiVZr) 1-x Cr x alloy

Multi-Principal Element Alloys (MPEA) have been hypothesized to maintain a stable single alloy phase and resist corrosion at high temperatures, which are desirable in various functional applications, including material of turbine blades. However, studies into oxidation resistant alloys have been constrained to individual MPEA systems or a limited selection of compositions within those systems. Here, High-Throughput Experimental (HTE) techniques were utilized to rapidly assess phase stability and resistance to oxidation of an (AlFeNiTiVZr) 1-x Cr x continuous composition spread (CCS) thin-film alloy system. High-throughput synchrotron diffraction and Raman spectroscopy were used to identify phases, especially high symmetry structures and oxides, as a function of annealing time and composition. The formation of two distinct phase regions, R 1 and R 2 , across the CCS was observed. One of these phase regions, R 2 , exhibited a stable high symmetry BCC-like phase at Cr concentrations greater than 25.4 at%, and grew a passivating mixed metal surface oxide; whereas, the R 1 phase region, decomposed on oxidation to multiple oxides and reacted with the substrate to form silicides. The passivating mixed metal oxide layer appears to prevent destabilization of the (AlFeNiTiVZr) 1-x Cr x MPEA BCC solid solution. Furthermore, at longer annealing times the passivating phase was seen only above 29.2 at% Cr. High Cr-doped (AlFeNiTiVZr) 1-x Cr x appear to be a good candidate for high temperature, oxidation-resistant structural alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Self-healing Growth of LaNiO 3 on a Mixed-Terminated Perovskite Surface

Developing atomic-scale synthesis control is a prerequisite for understanding and engineering the exotic physics inherent to transition-metal oxide heterostructures. Thus, far, however, the number of materials systems explored has been extremely limited, particularly with regard to the crystalline substrate, which is routinely SrTiO 3 . Here, we investigate the growth of a rare-earth nickelate-LaNiO 3 -on (LaAlO 3 )(Sr 2 AlTaO 6 ) (LSAT) (001) by oxide molecular beam epitaxy (MBE). Whereas the LSAT substrates are smooth, they do not exhibit the single surface termination usually assumed necessary for control over the interface structure. Performing both nonresonant and resonant anomalous in situ synchrotron surface X-ray scattering during MBE growth, we show that reproducible heterostructures can be achieved regardless of both the mixed surface termination and the layer-by-layer deposition sequence. The rearrangement of the layers occurs dynamically during growth, resulting in the fabrication of high-quality LaNiO 3 /LSAT heterostructures with a sharp and consistent interfacial structure. This is due to the thermodynamics of the deposition window as well as the nature of the chemical species at interfaces-here, the flexible charge state of nickel at the oxide surface. This has important implications regarding the use of a wider variety of substrates for fundamental studies on complex oxide synthesis.

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

Dependence of Cloud Macrophysical Properties and Phase Distributions on Environmental Conditions Over the North Atlantic and Southern Ocean: Results From COMBLE and MARCUS

The accurate representation of Cold Air Outbreaks (CAOs) and affiliated mixed-phase boundary layer (BL) clouds in models is challenging. How BL cloud properties evolve during CAOs and their dependence on meteorological conditions is not well understood but is important for the simulation of Earth's energy budgets. Here the properties of polar BL clouds over the North Atlantic (NA) and Southern Ocean (SO) are compared using observations from the Measurements of Aerosol Radiation and CloUds over the SO (MARCUS) and CAOs in the Marine BL Experiment (COMBLE) conducted over the NA. MARCUS observations show a stronger BL inversion than COMBLE, with a higher mean EIS (estimated inversion strength)/LTS (lower tropospheric stability) of -0.03 K/13 K compared to COMBLE’s -3.2 K/9.3 K. 39% of CAOs observed during COMBLE were intense with M > 5 K, while MARCUS only had 1.3%. 78%/72% of clouds sampled in CAOs during COMBLE/MARCUS had cloud top heights <4 km. The mean BL cloud top height was over 400 m higher, and the BL was over 500 m deeper for M of 10 K compared to 0 K for both regions. MARCUS observed a 27% moister BL structure than COMBLE when M > 5 K due to stronger BL inversion trapping more moisture within the BL. Under the same LTS, EIS, and M conditions, MARCUS observed a 12% drier BL structure, and clouds were 46% more turbulent than COMBLE. During CAOs, 54% of single-layer BL clouds sampled during MARCUS had liquid-dominated bases compared to 39% during COMBLE.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Generating a stratocumulus-like cloud top in a convection-cloud chamber

Stratocumulus-topped boundary layers play a crucial role in influencing daily weather and earth energy balance. Entrainment at the stratocumulus cloud top affects the cloud’s lifetime, precipitation, and radiative properties, but our understanding remains limited due to the lack of resolution in both field observations and numerical simulations. A recently proposed convection-cloud chamber with detailed control of sidewall temperatures can provide a unique opportunity to explore this mechanism in a laboratory setting. In this work, we use numerical simulations to demonstrate that this design can produce a cloud top that mimics the entrainment interfacial layer in a stratocumulus cloud. Our results show that a steady-state cloud can be formed by cooling the lower portions of the sidewalls and warming the bottom surface, while a temperature inversion at the cloud top can be generated by keeping the upper sidewalls and top surface warmer than the bottom. The turbulent kinetic energy profile and budget are similar to those found in a convective boundary layer, and inhomogeneous mixing near the cloud top can be observed. These findings significantly enhance the scientific value of constructing the tall convection-cloud chamber.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

C 13 ( n , 2 n γ ) C 12 γ -ray production in the 14–16 MeV incident neutron energy range

We report γ-ray emission from 12 C and 13 C samples irradiated with deuterium-tritium fusion neutrons was experimentally measured at the Omega Laser Facility and at the Ohio University Edwards Accelerator Laboratory. The intent of these measurements was to determine the feasibility of using 13 C-based plastic ablators with embedded 12 C layers for “dark mix” diagnosis of inertial confinement fusion implosions. Spectrally resolved measurements at Ohio University identified significant 4.44-MeV γ-ray emission from the 13 C(n, 2nγ) 12 C-L1 reaction channel. The recorded 4.44-MeV 13 C signal was compared against emission from an identically irradiated 12 C target with known 12 C(n, n'γ) 12 C-L1 cross section, which resulted in an average 13 C(n, 2nγ) 12 C-L1 cross section of 117 ± 17 mb over the incident neutron energy distribution range from 14.4 to 15.8 MeV. Integrated 13C γ -ray signals above 2.9 MeV recorded with the Gas Cherenkov Detector at Omega exceeded MCNP6.1 predictions by a factor of 3. The additional signal was attributed to 4.44-MeV γs resulting in an inferred 13 C(n, 2nγ) 12 C-L1 cross section of 95 ± 11 mb at 14.1-MeV average incident neutron energy. As a result, the 13 C-based “dark mix” diagnostic concept was deemed infeasible.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Detection of multi-modal Doppler spectra – Part 1: Establishing characteristic signals in radar moment data

Vertically pointing millimeter-wavelength radars provide a wealth of information about cloud and precipitation particle properties. Doppler spectral data can inform on how particles of varying vertical velocities contribute to the total backscattered power observed. It is more computationally cost effective to process moment data instead of spectra data, but doing so leaves valuable information on the cutting room floor. To confidently identify a multi-modal spectra event, in which two or more modes are present within a layer, Doppler spectral data are essential. This means long-term identification of layers featuring multi-modal spectra can be cost prohibitive. To address this, we explore three multi-modal spectra cases from winter precipitation events to determine characteristic signatures of these layers in the moment data averaged over short time periods (∼ 145 s) and explore how these layers differ from the rest of the vertical profiles. We find that the mean spectrum width and the standard deviation of mean Doppler velocity can be used to determine whether or not a layer is multi-modal. In particular, multi-modal layers in mixed-phase and ice clouds feature larger mean spectrum width (exceeding 0.17 m s −1 ) and smaller standard deviation of the mean Doppler velocity (below 0.1 m s −1 ). In Part 1 of this study, the identification criteria and methods are described. In Part 2 (Wugofski and Kumjian, 2025), we perform a verification of the method for three years of vertically pointing radar data, and explore the meteorological conditions associated with identified multi-modal spectral events.

Wugofski, Sarah [Pennsylvania State Univ., Univers↗
DOE OSTI2021

Retrieved Number concentration of Cloud Condensation Nuclei (RNCCN) Value-Added Product (VAP)

The objective of Retrieved Number concentration of Cloud Condensation Nuclei (RNCCN) Value-Added Product (VAP) is to provide the vertical distribution of cloud condensation nuclei (CCN) number concentrations to better represent aerosol indirect effects in climate models. The RNCCN VAP allows to collect a large number of independent samples of vertical CCN measurements at various supersaturation values from the ground. This helps to avoid the high cost of aircraft and increases the temporal distribution of retrieved CCN in the planetary boundary layer. The VAP method is based on Ghan and Collins (2004) and Ghan et al. (2006) proposed parameterizations that retrieve the vertical profiles of CCN directly from the ground level CCN, humidification factor, and in-situ lidar measurements. Further, the VAP accuracy is improved by using the quality controlled data that ensures that the boundary layer is well mixed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Dynamics of multilayer Rayleigh–Taylor instability at moderately high Atwood numbers

This paper investigates the multilayer Rayleigh–Taylor instability (RTI) using statistically stationary experiments conducted in a gas tunnel. Employing diagnostics such as particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF), we make simultaneous velocity–density measurements to study how dynamics and mixing are linked in this variable density flow. Experiments are conducted in a newly built, blow-down three-layer gas tunnel facility. Mixing between three gas streams is studied, where the top and bottom streams are comprised of air, and the middle stream is an air–helium mixture. Shear is minimized between these streams by matching their inlet velocities. The four experimental conditions investigated here consist of two different density ratios (Atwood numbers 0.3 and 0.6), each investigated at two instability development times (or equivalently, two streamwise locations), and all experiments are with the same middle stream thickness of 3 cm. The growth of the middle layer is measured using laser-based planar Mie scattering visualization. The mixing width is found to grow linearly with time at late times. Various quantitative measures of molecular mixing indicate a very high degree of molecular mixing at late times in the multilayer RTI flow. The vertical turbulent mass flux $a_y$ is calculated. In addition to mostly negative values of $a_y$ , typical of buoyancy-dominated flows due to negative correlation between velocity and density fluctuations, positive regions are also observed in profiles of $a_y$ due to entrainment and erosion at the lower edge of the mixing region. Global energy budgets are calculated for the multilayer RTI flow at late times and it is found that the majority of potential energy released has been dissipated due to viscous effects, and a large value of mixing efficiency ( $\sim$ 60 %) is observed.

Mechanics↗
DOE OSTI2025

Large Eddy Simulation of the Diurnal Cycle of Shallow Convection in the Central Amazon

Climate models often face challenges in accurately simulating the daily precipitation cycle over tropical land areas, particularly in the Amazon. One contributing factor may be the incomplete representation of the diurnal evolution of shallow cumulus (ShCu) clouds. This study aimed to enhance the understanding of the diurnal cycles of ShCu clouds—from formation to maturation and dissipation—over the Central Amazon (CAMZ). Using observational data from the Green Ocean Amazon 2014 (GoAmazon) campaign and large eddy simulation (LES) modeling, we analyzed the diurnal cycles of six selected pure ShCu cases and their composite behavior. Our results revealed a well-defined cycle, with cloud formation occurring between 10 and 11 local time (LT), maturity from 13 to 15 LT, and dissipation by 17–18 LT. The vertical extent of the liquid water mixing ratio and the intensity of the updraft mass flux were closely associated with increases in turbulent kinetic energy (TKE), enhanced buoyancy flux within the cloud layer, and reduced large-scale subsidence. We further analyzed the diurnal cycles of the convective available potential energy (CAPE), the convective inhibition (CIN), the Bowen ratio (BR), and the vertically integrated TKE in the mixed layer (ITKE-ML), exploring their relationships with the cloud base mass flux (Mb) and cloud depth across the six ShCu cases. ITKE-ML and Mb exhibited similar diurnal trends, peaking at approximately 14–15 LT. However, no consistent relationships were found between CAPE (or BR) and Mb. Similarly, comparisons of the cloud depth with CAPE, BR, ITKE-ML, CIN, and Mb revealed no clear relationships. Smaller ShCu clouds were sometimes linked to higher CAPE and lower CIN. It is important to emphasize that these findings are preliminary and based on a limited sample of ShCu cases. Further research involving an expanded dataset and more detailed analyses of the TKE budget and synoptic conditions is necessary. Such efforts would yield a more comprehensive understanding of the factors influencing ShCu clouds’ vertical development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Study of Intermittent Turbulence in Stable Arctic Boundary Layers

Abstract Polar boundary layers are difficult to model due to the existence of intermittent turbulence within stable layers. Here we present a case study evaluation of coherent structures in a stable boundary layer observed during a series of flights with an uncrewed aircraft system (DataHawk2) on 19 October 2016 at Oliktok Point, Alaska as part of the ERASMUS (Evaluation of Routine Atmospheric Sounding Measurements using Unmanned Systems) field campaign. During a sequence of five flights over a nine-hour period, 57 profiles of atmospheric properties (0–400 m a.g.l) were collected. Turbulence was identified using derived Richardson Number, temperature structure function parameter, and turbulence kinetic energy dissipation. Throughout all flights on this strongly stable day, intermittent turbulence was observed. These turbulent layers showed well-mixed potential temperature profiles embedded within otherwise stable potential temperature profiles; often resulting in a characteristic staircase pattern. Turbulent layers ranged from 1 to 30 m deep, with most individual layers being 1–2 m deep. Vertical propagation velocities of layers in the lower atmosphere were on the order of a few cm s −1 , typical of non-convective environments. In different regions of the profile, turbulence was driven by a different balance of buoyancy and shear forces, with turbulence in the near surface environment driven by strong shear forces overcoming strong resistance to buoyancy, while turbulence in elevated layers characterized by weaker shear forces overcoming weaker resistance to buoyancy. We discuss the potential of such datasets for improving subgrid parameterizations of small-scale turbulence embedded within stable boundary layers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗