Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shear 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 OSTI2025

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗
DOE OSTI2026

Suppression of ion temperature gradient modes by Alfvén activity above a drive threshold in DIII-D

A recent study demonstrates that the suppression of ion temperature gradient (ITG) modes can occur during the nonlinear evolution of toroidicity-induced Alfvén eigenmodes (TAEs) in the DIII-D tokamak [Du et al., Phys. Rev. Lett. 135, 265101 (2025)]. In that work, ITG was suppressed by the formation of a narrow, TAE-induced shear flow layer, whose shearing rate exceeded the ITG decorrelation rate. Here, the shear flow arises from an imbalance between Reynolds and Maxwell stress forces, as the TAE departs from the conventional shear Alfvén wave polarization. This follow-up paper systematically identifies the plasma conditions required for robust ITG suppression through a series of comparative experiments. The results show that TAEs routinely suppress ITG turbulence in plasmas with higher local safety factor (q), elevated fast-ion beta, and larger populations of fast ions on passing orbits, i.e., conditions consistent with the substantial TAE drive. Database analysis further reveals the existence of a threshold in TAE drive for ITG suppression. That is, once the drive exceeds a critical value within a favorable q window, the system undergoes the nonlinear bifurcation process, characterized by a sharp increase in TAE saturated amplitudes at a nearly fixed fast ion drive, suppression of ITG turbulence, and the formation of an internal thermal transport barrier at the localized radii.

Du, X. D. [General Atomics, San Diego, CA (United ↗
DOE OSTI2024

Overview of the recent experimental research on the J-TEXT tokamak

The J-TEXT capability is enhanced compared to two years ago with several upgrades of its diagnostics and the increase of electron cyclotron resonance heating (ECRH) power to 1 MW. With the application of electron cyclotron wave (ECW), the ECW assisted plasma startup is achieved; the tearing mode is suppressed; the toroidal injection of 300 kW ECW drives around 24 kA current; fast electrons are generated with toroidal injected ECW and the runaway current conversion efficiency increases with ECRH power. The mode coupling between 2/1 and 3/1 modes are extensively studied. The coupled 2/1 and 3/1 modes usually lead to major disruption. Their coupling can be either suppressed or avoided by external resonant magnetic perturbation fields and hence avoids the major disruption. It is also found that the 2/1 threshold of external field is significantly reduced by a pre-excited 3/1 mode, which can be either a locked island or an external kink mode. The disruption control is studied by developing prediction methods capable of cross tokamak application and by new mitigation methods, such as the biased electrode or electromagnetic pellet injector. The high-density operation and related disruptions are studied from various aspects. Approaching the density limit, the collapse of the edge shear layer is observed and such collapse can be prevented by applying edge biasing, leading to an increased density limit. The density limit is also observed to increase, if the plasma is operated in the poloidal divertor configuration or the plasma purity is increased by increasing the pre-filled gas pressure or ECRH power during the start-up phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development and Verification of an Improved Wake-Added Turbulence Model in FAST.Farm

We introduce a generalized wake-added turbulence (WAT) model in the multiphysics, multiturbine simulation tool FAST.Farm. The WAT model introduces additional small-scale turbulence that represents the breakdown of vortical structures and shear layers in the wake. The article describes the development, implementation, calibration, and verification of the model. The novelties of the model include support for wake asymmetry, buildup of WAT across the wind farm, and secondary effects of wake-induced turbulence (e.g., wake meandering) driven by smaller-scale turbulence structures that arise from wake breakdown. Large-eddy simulations were run to support the calibration of the WAT parameters and verification of the model. Previous studies hypothesized that the lack of WAT modeling was the source of underprediction of fatigue loads, in particular for cases at low turbulence intensities and/or stable atmospheric boundary layers. This study confirms that the newly implemented WAT model enhances the loads predictions in these cases.

17 WIND ENERGY↗
DOE OSTI2025

Visualizing the mass transfer flow in direct-impact accretion

We use a variety of visualization techniques to display the interior and surface flows in a double white dwarf binary undergoing direct-impact mass transfer and evolving dynamically to a merger. The structure of the flow can be interpreted in terms of standard dynamical, cyclostrophic, and geostrophic arguments. We describe and showcase some visualization and analysis techniques of potential interest for astrophysical hydrodynamics. In the context of R Coronae Borealis stars, we find that mixing of accretor material with donor material at the shear layer between the fast accretion belt and the slower rotating accretor body will always result in some dredge-up. We also discuss briefly some potential applications to other types of binaries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Magnetic flux emergence in binary neutron star remnants

Using high-resolution AthenaK simulations of a twisted toroidal flux tube, we study the flux emergence of magnetic structures in the shear layer of a hot massive neutron star typical of a binary neutron star remnant. High-resolution simulations demonstrate that magnetic buoyant instabilities allow for emergence only for extremely large magnetic fields significantly exceeding 10 17 G, and more typical fields around 10 16 G are instead dominated by hydrodynamic effects. Because merger remnants tend to be stable against hydrodynamic convection, our work places strong limitations on the mechanisms by which massive binary neutron star remnants can produce the magnetically-driven outflows needed to power jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Scalar Power Spectra and Scalar Structure Function Evolution in the Richtmyer–Meshkov Instability Upon Reshock

The Richtmyer–Meshkov instability of a twice-shocked gas interface is studied using high-speed planar laser-induced fluorescence in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. Furthermore, this initial condition is accelerated by a shock of nominal strength M = 1.9, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. Three individual experiments are analyzed, the energy spectrum and the structure functions of the light gas mole fraction field are calculated and compared.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The Effect of Pressure on NO x Entitlement and Reaction Timescales in a Premixed Axial Jet-In-Crossflow

This paper investigates the pressure dependency of a lean premixed jet injected into a lean vitiated crossflow with an experimentally verified detailed chemistry computational fluid dynamics (CFD) model and 53 species considered. Experimental data were taken in an axially staged combustor with an optically accessible test section, allowing the use of particle image velocimetry (PIV) and CH* chemiluminescence techniques as well as point measurement of species concentration, temperature, and pressure. The experimental data cases at one, three, and five atmospheres were selected to describe the flame stabilization dependency on pressure and gain the required knowledge for an extrapolation to engine condition. Simulated exit nitrogen oxide levels were validated with experimental emission data, and a global emission trend for the NO reduction at elevated pressure and constant turbine inlet temperature level was defined. The nitrogen oxide benefit at elevated operating pressure was justified with the significantly smaller flame surface area: the analysis of the simulated spanwise and top-view profiles showed a relatively short receded core flame with nitrogen oxide production in the center at high pressure relative to a longer and larger shear layer flame at atmospheric condition that produced NO toward the inner and outer side of the flame. Furthermore, decomposition of the Damköhler number revealed the strong influence of the reaction timescales with higher reaction rates at elevated pressure, along with a moderate influence of the turbulent timescales, showing higher turbulence intensity in the lee-side recirculation zone at lower pressure.

03 NATURAL GAS↗
DOE OSTI2020

Are Supercells Resistant to Entrainment because of Their Rotation?

This research investigates a hypothesis posed by previous authors, which argues that the helical nature of the flow in supercell updrafts makes them more resistant to entrainment than nonsupercellular updrafts because of the suppressed turbulence in purely helical flows. It was further supposed that this entrainment resistance contributes to the steadiness and longevity of supercell updrafts. A series of idealized large-eddy simulations were run to address this idea, wherein the deep-layer shear and hodograph shape were varied, resulting in supercells in the strongly sheared runs, nonsupercells in the weakly sheared runs, and variations in the percentage of streamwise vorticity in updrafts among runs. Fourier energy spectrum analyses show well-developed inertial subranges in all simulations, which suggests that the percentages of streamwise and crosswise vorticity have little effect on turbulence in convective environments. Additional analyses find little evidence of updraft-scale centrifugally stable flow within updrafts, which has also been hypothesized to limit horizontal mass flux across supercell updrafts. Results suggest that supercells do have smaller fractional entrainment rates than nonsupercells, but these differences are consistent with theoretical dependencies of entrainment on updraft width, and with supercells being wider than nonsupercells. Thus, while supercells do experience reduced fractional entrainment rates and entrainment-driven dilution, this advantage is primarily attributable to increased supercell updraft width relative to ordinary convection, and has little to do with updraft helicity and rotation.

Peters, John M.↗
DOE OSTI2021

Spatially and Temporally Resolved Velocimetry for Hypersonic Flows

The development of new hypersonic flight vehicles is limited by the physical understanding that may be obtained from ground test facilities. This has motivated the present development of a temporally and spatially resolved velocimetry measurement for Sandia National Laboratories (SNL) Hypersonic Wind Tunnel (HWT) using Femtosecond Laser Electronic Excitation Tagging (FLEET). First, a multi-line FLEET technique has been created for the first time and tested in a supersonic jet, allowing simultaneous measurements of velocities along multiple profiles in a flow. Secondly, two different approaches have been demonstrated for generating dotted FLEET lines. One employs a slit mask pattern focused into points to yield a dotted line, allowing for two- or three-component velocity measurements free of contamination between components. The other dotted-line approach is based upon an optical wedge array and yields a grid of points rather than a dotted line. Two successful FLEET measurement campaigns have been conducted in SNL’s HWT. The first effort established optimal diagnostic configurations in the hypersonic environment based on earlier benchtop reproductions, including validation of the use of a 267 nm beam to boost the measurement signal-to-noise ratio (SNR) with minimal risk of perturbing the flow and greater simplicity than a comparable resonant technique at 202 nm. The same FLEET system subsequently was reconstituted to demonstrate the ability to make velocimetry measurements of hypersonic turbulence in a realistic flow field. Mean velocity profiles and turbulence intensity profiles of the shear layer in the wake of a hypersonic cone model were measured at several different downstream stations, proving the viability of FLEET as a hypersonic diagnostic.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Control of Large-Amplitude Combustion Oscillations Using Nanosecond Repetitively Pulsed Plasmas

This paper details the use of nanosecond repetitively pulsed discharges to attenuate combustion instabilities in a 14 kW swirl-stabilized methane/air combustor. The combustor exhibits large-amplitude pressure oscillations ranging from 1 to 4% of the mean pressure during which the flame exhibits bulk motion in each instability cycle, upstream and downstream, as revealed by high-speed chemiluminescence. Control is accomplished with an electrode comprising a pin anode at the centerline of the combustor, allowing a nanosecond spark to be generated in a region spanning close to the flame base, through the shear layers of the swirling flow and ending at the metallic combustor wall. The discharges are generated using 20 kV, 9 kHz pulses; and they correspond to about 120 W of mean power. This results in a suppression of the peak amplitude of the pressure oscillations by a factor of two to four, and 5 dB in the rms value. Using phase-averaged visualizations of the flame with and without plasma, we detail the sequence of flame motion in the course of the instability. With the plasma active, this reveals significant interactions between the flame and the plasma during the suppression. Finally, we present a state-space model of the thermoacoustic system, and we demonstrate open-loop control of the instabilities.

Engineering↗
DOE OSTI2023

Hypersonic Fluid–Structure Interaction on a Cone–Slice–Ramp Geometry

In this study, fluid–structure interactions were measured between a representative control surface and the hypersonic flow deflected by it. The control surface is simplified as a spanwise finite ramp placed on a longitudinal slice of a cone. The front surface of the ramp contains a thin panel designed to respond to the unsteady fluid loading arising from the shock-wave/boundary-layer interactions. Experiments were conducted at Mach 5 and Mach 8 with ramps of different angles. High-speed schlieren captured the unsteady flow dynamics and accelerometers behind the thin panel measured its structural response. Panel vibrations were dominated by natural modes that were excited by the broadband aerodynamic fluctuations arising in the flowfield. However, increased structural response was observed in two distinct flow regimes: 1) attached or small separation interactions, where the transitional regime induced the strongest panel fluctuations. This was in agreement with the observation of increased convective undulations or bulges in the separation shock generated by the passage of turbulent spots, and 2) large separated interactions, where shear layer flapping in the laminar regime produced strong panel response at the flapping frequency. In addition, panel heating during the experiment caused a downward shift in its natural mode frequencies.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Particle Image Velocimetry and Total Temperature Characterization of Arc-Heated Mach 4.5 Free Jet

Particle image velocimetry (PIV) measurements were conducted in the Arc-Heated Combustion Tunnel (ACT-II) on a high-enthalpy Mach 4.5 nozzle flow. The primary motivation for the current work is to establish the feasibility of PIV as a high-quality flow diagnostic technique in hypersonic arc-heated flows and produce a high-resolution hypersonic nozzle flow characterization, including nozzle free jet shear layer velocity profiles. Due to the uncertain nature of the gas composition in the flowfield of interest, PIV measurements of a supersonic wedge flow were used to couple the flow velocities to the oblique shock relations. This enabled a PIV-derived total temperature calculation that does not require a priori assumptions of the exact gas composition. A PIV expansion module was developed for this study, and proper seed sizing and material were selected for the present flow. An analysis of the seed particle response from the PIV measurements provided an initial validation of the seed performance while revealing potential issues with consistent seed sizing. These measurements and derived conditions represent a means of acquiring high-resolution velocimetry and temperature datasets that consider uncertainties in the experimental generation of hypersonic high-enthalpy flows.

Engineering↗
DOE OSTI2025

100 kHz OH PLIF of High-Pressure H2/NH3 and H2/CH4 Flames

The reaction zones in a multi-element, multi-stage, micromix injector are investigated using planar laser-induced fluorescence (PLIF) imaging of the hydroxyl (OH) radical at a 100 kHz repetition rate. The combustor is operated at a pressure of ∼10 bar. Multiple fuel blends of hydrogen, methane, and ammonia are considered to study the effect of changing fuel composition on the flame structure. Compact regions of high OH signal are observed in each flame for fuel mixtures with high hydrogen content, indicating rapid reactant consumption in the shear layers. Instantaneous images exhibit wrinkled flame fronts and the formation of unreacted pockets due to highly turbulent burning. The addition of methane is observed to reduce the reactivity and flame speed, resulting in more diffuse regions of OH between each flame. The shape of the reactant core is altered, and the length scale of the flame surface wrinkling is increased. Ammonia addition has a more severe effect on the distribution of OH and the local flame structure than methane addition. A time series of images shows the presence of OH vortical structures and significant turbulent fluctuations in the flames near the wall, suggesting the presence of a recirculation zone that could drive the observed disparity in flame lengths.

Hodge, Alexander J. [Purdue University]↗
DOE OSTI2024

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Large Eddy Simulation of Rotationally Induced Ingress and Egress around an Axial Seal between Rotor and Stator Disks

In gas turbines, the hot gas exiting the combustor can have temperatures as high as 2000 °C, and some of this hot gas enter into the space between the stator and rotor disks (wheelspace). Since the entering hot gas could damage the disks, its ingestion must be minimized. This is carried out by rim seals and by introducing a cooler flow from the compressor (sealing flow) into the wheelspace. Ingress and egress into rim seals are driven by the stator vanes, the rotor and its rotation, and the rotor blades. This study focuses on the ingress and egress driven by the rotor and its rotation. This is carried out by performing wall-resolved large eddy simulation (LES) around an axial seal in a rotor–stator configuration without vanes and blades. Results obtained show the mechanisms by which the rotor and its rotation induce ingress, egress, and flow trajectories. Kelvin–Helmholtz instability was found to create a wavy shear layer and displacement thickness that produces alternating regions of high and low pressures around the rotor side of the seal. Vortex shedding on the backward-facing side of the seal and its impingement on the rotor side of the seal also produces alternating regions of high and low pressures. The locations of the alternating regions of high and low pressures were found to be statistically stationary and to cause ingress to start on the rotor side of the seal. Vortex shedding and recirculating flow in the seal clearance also cause ingress by entrainment. With the effects of the rotor and its rotation on ingress and egress isolated, this study enables the effects of stator vanes and rotor blades to be assessed.

42 ENGINEERING↗
DOE OSTI2026

Stability and Characteristics of Lower-hybrid Drift Waves: Dependence on Electron Beta and Cross-field Relative Drift

Lower-hybrid drift waves (LHDWs) are frequently observed microinstabilities in both space and laboratory plasmas. Despite decades of study, the relationship between electrostatic (ES-LHDW) and electromagnetic (EM-LHDW) variants and the plasma parameters controlling their stability remains unclear. Here, we systematically examine LHDW behavior by solving the local linear dispersion relation over a wide range of plasma and field conditions. Our results demonstrate that ES-LHDWs and EM-LHDWs are not distinct modes but rather two different regimes of the same drift wave whose character evolves smoothly with electron beta (β e ) and the cross-field electron drift velocity relative to ions, normalized to the ion sound speed (u 0x /C s ). The nature of the waves changes from electrostatic to electromagnetic when β e increases. Growth rates increase with u 0x /C s but decrease with β e , while the most unstable wavelength remains nearly universal, with kρ e ∼ 0.8 (k is the magnitude of the wave vector and ρ e is the electron gyroradius). We further present quasi-linear estimates of nonlinear saturation properties, including energy partition among electric fields, magnetic fields, and particle kinetic responses. We show that ES-LHDWs reach higher electric-field saturation amplitudes, whereas EM-LHDWs generate strong magnetic perturbations and parallel electric fields that may enable efficient particle heating. Comparisons with the classical model reveal that retaining electromagnetic effects is essential for accurate predictions of frequency, growth rate, and the propagation angle. These findings provide a unified framework for understanding LHDWs across diverse collisionless plasma environments, including current sheets of magnetic reconnection, shear layers, collisionless shocks, and boundary regions.

Solar coronal waves↗
DOE OSTI2024

On the Geometry of the Near-core Magnetic Field in Massive Stars

It is well-known that the cores of massive stars sustain a stellar dynamo with a complex magnetic field configuration. However, the same cannot be said for the field's strength and geometry at the convective–radiative boundary, which are crucial when performing asteroseismic inference. In this Letter, we present 3D magnetohydrodynamic (MHD) simulations of a 7 M ⊙ mid-main-sequence star, with particular attention given to the convective–radiative boundary in the near-core region. Our simulations reveal that the toroidal magnetic field is significantly stronger than the poloidal field in this region, contrary to recent assumptions. Moreover, the rotational shear layer, also important for asteroseismic inference, is specifically confined within the extent of the Brunt–Väisälä frequency peak. These results, which are based on the inferred properties of HD 43317, have widespread implications for asteroseismic studies of rotation, mixing, and magnetism in stars. While we expect our results to be broadly applicable across stars with similar Brunt–Väisälä frequency profiles and stellar masses, we also expect the MHD parameters (e.g., Re m ) and the initial stellar rotation rate to impact the geometry of the field and differential rotation at the convective–radiative interface.

79 ASTRONOMY AND ASTROPHYSICS↗