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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The Spatial and Temporal Variability of the Clear Convective Boundary Layer at the ARM SGP Supersite

The convective boundary layer (CBL), also known as the mixing layer, constitutes the critical lower segment of the atmosphere that significantly influences daily human activities. The growing demand for precise weather forecasts is driven by the requirements of agriculture, transportation, and routine societal functions. Here, to enhance understanding of the CBL, this study investigates the spatiotemporal variability in the CBL and its controlling factors using four-year Doppler lidar, surface flux, and profiling measurements at five ARM Southern Great Plains sites within a 100 km radius. This investigation utilizes data collected exclusively under clear-sky conditions or scattered low-cloud conditions. Results reveal significant spatial differences in CBL evolutions. Daily mixing layer heights (MLHs) vary up to 1 km (30% of the mean) in late afternoon. There is a clear east–west contrast: western sites (C1, E32, E37) exhibit higher summer MLH (1.9–2.1 km) and vertical velocity variances (1.0–1.2 m 2 s −2 ) than eastern sites (1.6–1.8 km), reversing in winter. Temporally, the MLH peaks at 70% of the sunrise–sunset interval, the lagging heat flux (HF) peaks at 50%; and the seasonal MLH maxima lag the HF by approximately one month, influenced by nighttime PBL (planetary boundary layer) properties. The HF and lower tropospheric stability are the main factors of influence for the CBL, but site-specific dependencies highlight the critical roles of local factors, underscoring the need for including them in CBL modeling.

ARM SGP site↗
DOE OSTI2021

Connection Between Mass Flux Transport and Eddy Diffusivity in Convective Atmospheric Boundary Layers

Turbulence parameterizations for convective boundary layer in coarse-scale atmospheric models usually consider a combination of the eddy-diffusive transport and a non-local transport, typically in the form of a mass flux term, such as the widely adopted eddy-diffusivity mass-flux (EDMF) approach. These two types of turbulent transport are generally considered to be independent of each other. Using results from large-eddy simulations, here, we show that a Taylor series expansion of the updraft and downdraft mass-flux transport can be used to approximate the eddy-diffusivity transport in the atmospheric surface layer and the lower part of the mixed layer, connecting both eddy-diffusivity and mass-flux transport theories in convective conditions, which also quantifies departure from the Monin-Obukhov similarity (MOS) in the surface layer. This study provides a theoretical support for a unified EDMF parameterization applied to both the surface layer and mixed layer and highlights important correction required for surface models relying on MOS

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Scalar mixing in a Kelvin-Helmholtz shear layer and implications for Reynolds-averaged Navier-Stokes modeling of mixing layers

Large-eddy simulation of a temporally evolving Kelvin-Helmholtz (KH) mixing layer is performed with the tenth-order compact difference code miranda to examine the steady-state behavior of a passive scalar in a shear-driven mixing layer. It is shown that the integral behavior of scalar variance in a KH mixing layer behaves similarly to the integral behavior of scalar variance in a Rayleigh-Taylor (RT) mixing layer, and mixedness of the simulated KH shear layer tends towards a value of about 0.8. It is further shown that if the k-L-a-V Reynolds-averaged Navier-Stokes (RANS) model [B. E. Morgan et al., Phys. Rev. E 98, 033111 (2018)], calibrated to reproduce steady-state mixing in an RT layer, is applied to simulate a KH mixing layer, the RANS model will significantly overpredict the magnitude of scalar variance in the KH layer. A straightforward addition to the k-L-a-V model is then suggested, and self-similarity analysis is applied to determine constraints on model coefficients. Furthermore, it is shown that with the addition of a buoyancy production term in the model equation for scalar variance, it becomes possible to eliminate the model deficiency and match steady-state mixedness in simulations of both RT and KH mixing layers with a single model calibration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

The DOE E3SM version 2.1: overview and assessment of the impacts of parameterized ocean submesoscales

The U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) version 2.1 builds on E3SMv2 with several changes, with the most notable being the addition of the Fox-Kemper et al. (2011) mixed-layer eddy parameterization. This parameterization captures the effect of finite-amplitude, mixed-layer eddies as an overturning streamfunction and has the primary function of restratification. Herein, we outline the changes to the mean climate state of E3SM that were introduced by the addition of this parameterization. Overall, the presence of the submesoscale parameterization improves the fidelity of the v2.1 simulation by reducing the ocean surface biases in the North Atlantic present in v2, as illustrated by changes in the climatological sea surface temperature and salinity and the Arctic sea-ice extent. Other impacts include a slight shoaling of the mixed-layer depths in the North Atlantic and a small improvement in the Atlantic Meridional Overturning Circulation (AMOC). We note that the expected shoaling due to the parameterization is regionally dependent in our coupled configuration. In addition, we investigate why the parameterization and its impacts on mixed-layer depth have little impact on the simulated AMOC: despite increased dense-water formation in the Norwegian Sea, only a small fraction of the water formed makes its way south into the North Atlantic basin. Version 2.1 also exhibits small improvements in the atmospheric climatology, with smaller biases in many notable quantities and modes of variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantifying the Role of Ocean Dynamics in SST Variability across GCMs and Observations

Abstract Midlatitude SSTs forced by mesoscale oceanic processes can affect the large-scale atmosphere, pointing to the ocean’s crucial role outside the tropics. Previous studies have shown oceanic mesoscale processes’ effect on global and regional climate variability. This study quantifies the local contribution of ocean dynamics to mixed-layer temperature across the globe by directly estimating the ocean heat flux divergence resolved by state-of-the-art ocean reanalysis, eddy-resolving, and eddy-parameterized versions of two U.S. national climate models and indirectly from air–sea flux satellite-based estimates. Our results show that the eddy-resolving climate simulations resolve mixed-layer temperature variances that are larger and closer to those inferred from observations than both their eddy-parameterized counterparts and ECCO over much of the extratropics. The observations and the eddy-resolving models indicate a more significant role of ocean dynamics in the mixed-layer temperature variability than the surface fluxes over most extratropics compared to their eddy-parameterized versions. A frequency domain analysis shows that the better-resolved ocean mesoscale and thermal gradients enhance the variance over a time scale from 2 months to 30 years. Results show agreement in the ocean’s contribution among satellite-based estimates, ocean reanalysis products, and ocean eddy-resolving simulations. At the same time, differences emerge for ECCO and the eddy-parameterized models, suggesting that surface fluxes account for a larger fraction of the mixed-layer temperature variability in most of the extratropics.

Siqueira, L.↗
DOE OSTI2024

Convective Boundary-Layer Spatial Heterogeneity Experiment Field Campaign Report

Earth system models (ESMs) require accurate heat, mass, and momentum exchanges between components, which requires accurate observations and modeling of the atmospheric boundary layer. Over the land, the daytime convective atmospheric boundary layer (CBL, also called the convective mixing layer) develops and evolves daily, driven by solar surface heating. Doppler lidar (DL) measurements of vertical velocities provide an effective way to document the diurnal cycle of CBL (Chu et al. 2023). With the support of the U.S. Department of Energy (DOE) Atmospheric Research (ASR) program, we studied convective mixing-layer heights (MLH) across multiple DOE Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observatory sites based on DL measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Large-eddy simulation and Reynolds-averaged Navier-Stokes modeling of three Rayleigh-Taylor mixing configurations with gravity reversal

In this work, high-fidelity large-eddy simulation (LES) is performed of Rayleigh-Taylor (RT) mixing in three different configurations involving gravity reversal. In each configuration, LES results are compared with one-dimensional Reynolds-averaged Navier-Stokes (RANS) results, and a deficiency in a commonly used transport equation for the mass-flux velocity, $a_j$, is identified. In the first configuration, a classical two-component RT mixing layer is allowed to develop before it is subjected to rapid acceleration reversal. In the second configuration, a three-component RT mixing layer with an intermediate density layer is allowed to develop before being subjected to rapid acceleration reversal. Finally, in the third configuration, a light layer is interposed between two heavy layers; in this configuration, only one interface is RT-unstable at a time as it undergoes rapid acceleration reversal. In all cases, a commonly used buoyancy production closure in the $a_j$ transport equation is shown to lead to significant over-prediction of mixing layer growth after gravity reversal. An alternative formulation for this closure is then presented which is shown to more accurately capture the stabilization effect of gravity reversal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Experiments and simulations on the Richtmyer-Meshkov instability with a thin intermediate layer

Experiments and simulations of the Richtmyer-Meshkov instability (RMI) in two- and three-layer configurations are presented. The two-layer case utilizes a light-over-heavy configuration and consists of air as the light gas and sulfur hexafluoride (SF 6 ) as the heavy gas. The three-layer case utilizes a light-intermediate-heavy configuration, with helium (He) as the light gas, air as the intermediate gas, and SF 6 as the heavy gas. Statistically significant differences in the mixing layer width of the lower interface are not observed between the two cases. This differs from the experiments of Schalles et al., where a small, though statistically significant, difference in mixing layer growth was observed between the two- and three-layer cases with a nominally two-dimensional, single mode perturbation. Notably, the perturbations on the lower interface in the present work do not grow large enough to significantly interact with the upper interface during the duration of the experiments. In conclusion, this suggests that the differences in mixing layer growth observed by Schalles et al. may be due to interactions of the perturbations on one interface with the other interface rather than being inherent to the three-layer problem.

42 ENGINEERING↗
DOE OSTI2024

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Demons in the North Atlantic: Variability of Deep Ocean Ventilation

Abstract Translation of atmospheric forcing variability into the ocean interior via ocean ventilation is an important aspect of transient climate change. On a seasonal timescale in the subtropics, this translation is mediated by a so‐called “Demon” that prevents access to all except late‐winter mixed‐layer water. Here, we use an eddy‐permitting numerical circulation model to investigate a similar process operating on longer (interannual) timescales in the subpolar North Atlantic. We find that variations in atmospheric forcing are mediated in their translation to the ocean interior, with year‐to‐year changes in the late‐winter mixed layer depth being the critical factor. The signature of persistent strong atmospheric forcing driving deep mixed layers is preferentially ventilated to the interior when the forcing is ceased. Susceptibility to this effect depends on the location and density of subduction—with the rate at which newly ventilated water escapes its region of subduction being the crucial factor.

MacGilchrist, G. A.↗
DOE OSTI2026

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗
DOE OSTI2024

Assessment of the literature about Be-W mixed material layer formation in the fusion reactor environment

All plasma facing surfaces in a fusion reactor, whether initially pure or an alloy, will rapidly evolve into a mixed material due to plasma-induced erosion, migration and redeposition. Beryllium (Be) erosion from the main chamber, and its transport and deposition on to a tungsten (W) divertor results in the growth of mixed Be-W layers, which can evolve to form beryllides. These Be-W mixed materials exhibit generally less desirable properties than pure tungsten or pure beryllium, such as lower melting points. In order to better understand the parameter space for growth of these alloys, this paper reviews the literature on Be-W mixed material formation experiments—in magnetically confined fusion reactors, in linear plasma test stands, and during thin-film deposition—and on computational modeling of Be-W interactions, as well as briefly assesses the Be-W growth kinetics. We conclude that the following kinetic steps drive the material mixing: adsorption of the implanted/deposited ion on the metal surface; diffusion of the implanted/deposited ion from surface into the bulk, which is accelerated by defects; and loss of deposited material through erosion. Adsorption dominates (or prevents) material mixing in thin-film deposition experiments, whereas diffusion drives material mixing in plasma exposures due to the energetic ion implantation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Assessment of the literature about Be-W mixed material layer formation in the fusion reactor environment

Abstract All plasma facing surfaces in a fusion reactor, whether initially pure or an alloy, will rapidly evolve into a mixed material due to plasma-induced erosion, migration and redeposition. Beryllium (Be) erosion from the main chamber, and its transport and deposition on to a tungsten (W) divertor results in the growth of mixed Be-W layers, which can evolve to form beryllides. These Be-W mixed materials exhibit generally less desirable properties than pure tungsten or pure beryllium, such as lower melting points. In order to better understand the parameter space for growth of these alloys, this paper reviews the literature on Be-W mixed material formation experiments—in magnetically confined fusion reactors, in linear plasma test stands, and during thin-film deposition—and on computational modeling of Be-W interactions, as well as briefly assesses the Be-W growth kinetics. We conclude that the following kinetic steps drive the material mixing: adsorption of the implanted/deposited ion on the metal surface; diffusion of the implanted/deposited ion from surface into the bulk, which is accelerated by defects; and loss of deposited material through erosion. Adsorption dominates (or prevents) material mixing in thin-film deposition experiments, whereas diffusion drives material mixing in plasma exposures due to the energetic ion implantation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗
DOE OSTI2022

A Hierarchy of Global Ocean Models Coupled to CESM1

Abstract We develop a hierarchy of simplified ocean models for coupled ocean, atmosphere, and sea ice climate simulations using the Community Earth System Model version 1 (CESM1). The hierarchy has four members: a slab ocean model, a mixed‐layer model (MLM) with entrainment and detrainment, an Ekman MLM, and an ocean general circulation model (OGCM). Flux corrections of heat and salt are applied to the simplified models ensuring that all hierarchy members have the same climatology. We diagnose the needed flux corrections from auxiliary simulations in which we restore the temperature and salinity to the daily climatology obtained from a target CESM1 simulation. The resulting three‐dimensional corrections contain the interannual variability fluxes that maintain the correct vertical gradients of temperature and salinity in the tropics. We find that the inclusion of mixed‐layer entrainment and Ekman flow produces sea surface temperature and surface air temperature fields whose means and variances are progressively more similar to those produced by the target CESM1 simulation. We illustrate the application of the hierarchy to the problem of understanding the response of the climate system to the loss of Arctic sea ice. We find that the shifts in the positions of the mid‐latitude westerly jet and of the Inter‐tropical Convergence Zone (ITCZ) in response to sea‐ice loss depend critically on upper ocean processes. Specifically, heat uptake associated with the mixed‐layer entrainment influences the shift in the westerly jet and ITCZ. Moreover, the shift of ITCZ is sensitive to the form of Ekman flow parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A hierarchy of global ocean models coupled to CESM1

We develop a hierarchy of simplified ocean models for coupled ocean, atmosphere, and sea ice climate simulations using the Community Earth System Model version 1 (CESM1). The hierarchy has four members: a slab ocean model, a mixed-layer model (MLM) with entrainment and detrainment, an Ekman MLM, and an ocean general circulation model (OGCM). Flux corrections of heat and salt are applied to the simplified models ensuring that all hierarchy members have the same climatology. We diagnose the needed flux corrections from auxiliary simulations in which we restore the temperature and salinity to the daily climatology obtained from a target CESM1 simulation. The resulting three-dimensional corrections contain the interannual variability fluxes that maintain the correct vertical gradients of temperature and salinity in the tropics. We find that the inclusion of mixed-layer entrainment and Ekman flow produces sea surface temperature and surface air temperature fields whose means and variances are progressively more similar to those produced by the target CESM1 simulation. We illustrate the application of the hierarchy to the problem of understanding the response of the climate system to the loss of Arctic sea ice. We find that the shifts in the positions of the mid-latitude westerly jet and of the Inter-tropical Convergence Zone (ITCZ) in response to sea-ice loss depend critically on upper ocean processes. Specifically, heat uptake associated with the mixed-layer entrainment influences the shift in the westerly jet and ITCZ. Moreover, the shift of ITCZ is sensitive to the form of Ekman flow parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Flow regime and Reynolds number variation effects on the mixing behavior of parallel flows

The hydraulic single-phase mixing of three parallel rectangular channels is experimentally investigated at various Reynolds numbers (Re) and flow regime combinations. Particle Image Velocimetry results for seven mixing cases are presented and discussed with varying Re combinations ranging from 1,824 to 20,844. While all cases result in the same Re ratio of ~0.69 between the inner and outer flows, two cases represent multi-regime mixing with the inner-outer regime pair of laminar-transitional and transitional-turbulent, while the other 5 cases are all characteristic of turbulent mixing with varying levels of turbulence. The outer channels initially share characteristics with a backward facing step. The center channel is found to initially behave like a slot jet, but then sees a significant increase in velocity decay. This inner flow velocity decay increased dramatically in the laminar-transitional mixing case, whose centerline velocity decay was ~6 times larger than the decay in the turbulent mixing cases. Second order statistics revealed a consistent mixing layer thickness of ~0.1 hydraulic diameters for all the cases but showed more intense shearing in the multi-regime mixing cases. The combined point and thereby the mixing layer length is determined using centerline velocity decay profiles, which show a much more aggressive mixing in multi-regime flows. Multi-regime mixing demonstrated superior characteristics relative to turbulent mixing due to a more dramatic velocity decay in the inner flow and a shorter mixing length. The contributions of this work include communicating the benefits of multi-regime mixing and providing detailed characterization efforts that can serve future efforts for validating computational models. Here this research also lays the groundwork for future studies aimed at achieving high levels of mixing without a severe penalty in pressure drop.

42 ENGINEERING↗
DOE OSTI2023

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗