Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shear Layers”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Flame Flashback Investigations in Hydrogen-Enriched Low Swirl Flames using High-Speed Hydroxyl (OH) Planar Laser-Induced Fluorescence

This study reports flashback mechanism observed in hydrogen-enriched flames stabilized in a low swirl burner (LSB) at atmospheric pressure and temperature conditions. The fundamentals of hydrogen-rich stable flames and spatiotemporal investigation of flashback phenomena were observed experimentally using a high-repetition-rate nanosecond (ns)-duration hydroxyl radical planar laser-induced fluorescence (OH-PLIF) diagnostic. Testing was conducted in an optically accessible pre-mixing section of LSB with inlet pre-mixing velocities from 5 to 10 m/s for methane and hydrogen (50 – 90% by mole) blends. Swirlers with two different turning angles, 26 and 33 degrees and three different perforated plate hole diameters, 1.08, 1.12 and 1.16 mm with measured swirl numbers varying from 0.43 to 0.49 were used. The flashback propensity showed dependence on the proximity of the lifted flame to the burner exit which was dependent on pre-mixer velocity, hydrogen content and equivalence ratio (ER). High-speed OH-PLIF images revealed lifted flames were first observed at low ER condition which is changed to M-shaped flame attached to burner rim with the increase in ER. Then, further increase in ER triggers flame flashback into the nozzle. Flashback occurred when burning occurs in the outer shear layer and the leading flame brush ingresses into the nozzle.

Parajuli, Pradeep↗
DOE OSTI2024

High-Speed OH-PLIF Diagnostics of Flame Flashback in Low Swirl Hydrogen-Enriched Flames

This paper reports flashback events observed in hydrogen-enriched flames stabilized in a low swirl burner (LSB) at atmospheric pressure and temperature conditions. The fundamentals of hydrogen-rich stable flames and spatiotemporal investigation of flashback phenomena were observed experimentally using a high-repetition-rate nanosecond (ns)-duration hydroxyl radical planar laser-induced fluorescence (OH-PLIF) diagnostic. Testing was conducted in an optically accessible pre-mixing section of a LSB with inlet pre-mixing velocities ranging from 5 to 10 m/s for methane and hydrogen (50 – 90% by mole) blends. Swirlers with two different turning angles, 26o and 33o, and three different perforated plate hole diameters, 1.08, 1.12 and 1.16 mm, with measured swirl numbers varying from 0.43 to 0.49 were used in this study at atmospheric temperature and pressure inlet conditions. The flashback propensity showed a dependence on the proximity of the lifted flame to the burner exit (termed as lift-off length, L) which was dependent on the pre-mixer velocity (V), hydrogen content (XH2) and equivalence ratio (ϕ) at constant temperature and pressure conditions. High-speed OH-PLIF images revealed lifted flames were first observed at low ϕ/ XH2 condition which then changed to a M-shaped flame attached to burner rim with the increase in ϕ/ XH2. Further increases in ϕ/XH2, depending upon other inlet parameters, triggered flame flashback into the pre-mixing section. Flashback occurred when burning began in the outer shear layer and the leading flame brush propagated into the nozzle. Such spatiotemporally resolved high-speed OH-PLIF imaging provides insights on different stages of flashback in a LSB – flashback initiation, transition of flame propagation from burner exit to the swirler and flashback to flame holding transition. Flame ϕ at flashback showed an expected linearly increasing trend with increasing V and decreasing XH2 and the conclusions drawn aligned well with detailed investigations. For identical inlet conditions, flashback propensity decreased with increasing perforated-plate hole diameter and increasing swirler vane angle.

Parajuli, Pradeep↗
DOE OSTI2026

Validation of OpenPronghorn for Periodic Hill Flow Separation

OpenPronghorn is an open-source, MOOSE-based thermal-hydraulics simulation tool used for advanced reactor analysis. As an open-source code, it offers a transparent framework for validating governing equations, assumptions, and numerical methods against established benchmarks. This study evaluates OpenPronghorn's RANS turbulence model against the ERCOFTAC Case 81 periodic hill benchmark, a standard test case for separated turbulent flow featuring curved-wall separation, recirculation, shear-layer development, and reattachment. Streamwise velocity profiles predicted by OpenPronghorn were compared to reference LES data at multiple x/h locations. Results show that OpenPronghorn captures the overall trend of the velocity profiles, but the largest discrepancies occur in the separated-flow region, where turbulence is highly anisotropic and strongly affected by adverse pressure gradients and wall curvature—conditions that are inherently difficult for standard RANS models to resolve. Future work will test alternative k-e model variants and correction terms to improve prediction accuracy in this region.

42 - ENGINEERING↗
DOE OSTI2020

Variable-density effects in incompressible non-buoyant shear-driven turbulent mixing layers

The asymmetries that arise when a mixing layer involves two miscible fluids of differing densities are investigated using incompressible (low-speed) direct numerical simulations. The simulations are performed in the temporal configuration with very large domain sizes, to allow the mixing layers to reach prolonged states of fully turbulent self-similar growth. Imposing a mean density variation breaks the mean symmetry relative to the classical single-fluid temporal mixing layer problem. Unlike prior variable-density mixing layer simulations in which the streams are composed of the same fluids with dissimilar thermodynamic properties, the density variations are presently due to compositional differences between the fluid streams, leading to different mixing dynamics. Variable-density (non-Boussinesq) effects introduce strong asymmetries in the flow statistics that can be explained by the strongest turbulence increasingly migrating to the lighter fluid side as free-stream density difference increases. Interface thickness growth rates also reduce, with some thickness definitions particularly sensitive to the corresponding changes in alignment between density and streamwise velocity profiles. Additional asymmetries in the sense of statistical distributions of densities at a given position within the mixing layer reveal that fine scales of turbulence are preferentially sustained in lighter fluid, which also is where fastest mixing occurs. These effects influence statistics involving density fluctuations, which have important implications for mixing and more complicated phenomena that are sensitive to the mixing dynamics, such as combustion.

42 ENGINEERING↗
DOE OSTI2021

Environmental sensitivities of shallow-cumulus dilution – Part 2: Vertical wind profile

This second part of a numerical study on shallow-cumulus dilution focuses on the sensitivity of cloud dilution to changes in the vertical wind profile. Insights are obtained through large-eddy simulations of maritime and continental cloud fields. In these simulations, the speed of the initially uniform geostrophic wind and the strength of geostrophic vertical wind shear in the cloud and subcloud layer are varied. Increases in the cloud-layer vertical wind shear (up to 9 ms -1 km -1 ) lead to 40%–50% larger cloud-core dilution rates compared to their respective unsheared counterparts. When the background wind speed, on the other hand, is enhanced by up to 10 ms -1 and subcloud-layer vertical wind shear develops or is initially prescribed, the dilution rate decreases by up to 25%. The sensitivities of the dilution rate are linked to the updraft strength and the properties of the entrained air. Increases in the wind speed or vertical wind shear result in lower vertical velocities across all sets of experiments with stronger reductions in the cloud-layer wind shear simulation (27%–47%). Weaker updrafts are exposed to mixing with the drier surrounding air for a longer time period, allowing more entrainment to occur (i.e., the “core-exposure effect”). However, reduced vertical velocities, in concert with increased cloud-layer turbulence, also assist in widening the humid shell surrounding the cloud cores, leading to entrainment of more humid air (i.e., the “core–shell dilution effect”). In the experiments with cloud-layer vertical wind shear, the core-exposure effect dominates and the cloud-core dilution increases with increasing shear. Conversely, when the wind speed is increased and subcloud-layer vertical wind shear develops or is imposed, the core–shell dilution effect dominates to induce a buffering effect. The sensitivities are generally stronger in the maritime simulations, where weaker sensible heat fluxes lead to narrower, more tilted, and, therefore, more suppressed cumuli when cloud-layer shear is imposed. Moreover, in the experiments with subcloud wind shear, the weaker baseline turbulence in the maritime case allows for a larger turbulence enhancement, resulting in a widening of the transition zones between the cores and their environment, leading to the entrainment of more humid air.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Environmental Controls on MCS Lifetime Rainfall Over Tropical Oceans

Abstract Mesoscale convective systems (MCSs) contribute a majority of rainfall over tropical oceans. However, our understanding of the environmental controls on tropical oceanic MCS precipitation remains incomplete. Using 20‐year of satellite observations, reanalysis data, and MCS tracking, we found that MCSs initiating in a mesoscale environment with enhanced lower‐free‐tropospheric moisture, warmer middle troposphere, stronger low‐level ascent, and stronger deep‐layer (surface‐400 hPa) wind shear tend to produce more precipitation during their lifetimes. While most of these environmental factors are correlated with one another, the deep‐layer shear is not. A rapid pickup in MCS lifetime rainfall is found when the lower‐free‐tropospheric specific humidity exceeds 10 g kg −1 . This nonlinearity is mostly dominated by the nonlinear increase in MCS area. On the other hand, both MCS area and rain rate increase quasi‐linearly with the deep‐layer shear. The increase in rain rate is related to the enhancement of heavy precipitating convective activity with deep‐layer shear.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Polarimetric Signatures in Landfalling Tropical Cyclones

Polarimetric radar observations from the NEXRAD WSR-88D operational radar network in the contiguous United States, routinely available since 2013, are used to reveal three prominent microphysical signatures in landfalling tropical cyclones: 1) hydrometeor size sorting within the eyewall convection, 2) vertical displacement of the melting layer within the inner core, and 3) dendritic growth layers within stratiform regions of the inner core. Size sorting signatures within eyewall convection are observed with greater frequency and prominence in more intense hurricanes, and are observed predominantly within the deep-layer environmental wind shear vector-relative quadrants that harbor the greatest frequency of deep convection (i.e., downshear and left-of-shear). Melting-layer displacements are shown that exceed 1 km in altitude compared to melting-layer altitudes in outer rainbands and are complemented by analyses of archived dropsonde data. Dendritic growth and attendant snow aggregation signatures in the inner core are found to occur more often when echo-top altitudes are low (≤10 km MSL), nearer the -15°C isotherm commonly associated with dendritic growth. These signatures, uniquely observed by polarimetric radar, provide greater insight into the physical structure and thermodynamic characteristics of tropical cyclones, which are important for improving rainfall estimation and the representation of tropical cyclones in numerical models.

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

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Indications of a Decrease in the Depth of Deep Convective Cores with Increasing Aerosol Concentration during the CACTI Campaign

Abstract An aerosol indirect effect on deep convective cores (DCCs), by which increasing aerosol concentration increases cloud-top height via enhanced latent heating and updraft velocity, has been proposed in many studies. However, the magnitude of this effect remains uncertain due to aerosol measurement limitations, modulation of the effect by meteorological conditions, and difficulties untangling meteorological and aerosol effects on DCCs. The Cloud, Aerosol, and Complex Terrain Interactions (CACTI) campaign in 2018–19 produced concentrated aerosol and cloud observations in a location with frequent DCCs, providing an opportunity to examine the proposed aerosol indirect effect on DCC depth in a rigorous and robust manner. For periods throughout the campaign with well-mixed boundary layers, we analyze relationships that exist between aerosol variables (condensation nuclei concentration > 10 nm, 0.4% cloud condensation nuclei concentration, 55–1000-nm aerosol concentration, and aerosol optical depth) and meteorological variables [level of neutral buoyancy (LNB), convective available potential energy, midlevel relative humidity, and deep-layer vertical wind shear] with the maximum radar-echo-top height and cloud-top temperature (CTT) of DCCs. Meteorological variables such as LNB and deep-layer shear are strongly correlated with DCC depth. LNB is also highly correlated with three of the aerosol variables. After accounting for meteorological correlations, increasing values of the aerosol variables [with the exception of one formulation of aerosol optical depth (AOD)] are generally correlated at a statistically significant level with a warmer CTT of DCCs. Therefore, for the study region and period considered, increasing aerosol concentration is mostly associated with a decrease in DCC depth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

New insights into creep characteristics of calcium silicate hydrates at molecular level

The fundamental mechanisms under concrete creep are far from being fully understood, especially at the molecular level. Hereby, a calcium-silicate-hydrate (C-S-H) molecular model is developed to explain, for the first time, the creep characteristics at various stress states, temperature levels and water contents, which are not accessible experimentally. Rather tensile and compressive loadings, C-S-H only creeps under shear loadings originating from the sliding of the calcium silicate layers over each other as the interlayer component (water and ions) acts as a lubricator. A heterogeneous creep characteristic is observed. Elevated temperature reduces the interlayer lubricator viscosity and weakens the interfacial adhesion between the layers and the interlayer lubricator, which accelerates C-S-H creep. The removal of interlayer water enhances the creep resistance, resulting from the reduced interlayer space and enhanced interfacial adhesion. The atomic-level mechanisms explain the inter-CSH-particle behaviours at the microscale, which bridges the gap between atomistic simulation and microcosmic phenomenon.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nanoscale phase transition in impact indented Mo/Cu thin films

The study reveals a size effect in deformation mechanisms via unique structural transitions in sputter deposited fcc/bcc nanolayers under high strain rate impact nanoindentation as layer thickness is reduced from around 100-200 nm to around 5 nm. Here, the indented 5 nm Cu/5 nm Mo multilayer films characterized via cross-sectional scanning/transmission electron microscopy have revealed a phase transition in the 5 nm thick Mo layers from BCC to a coherent metastable FCC structure and fine deformation twins in the 5 nm Cu layers inside a shear band that is approximately 500 nm long and inclined to the layer interfaces. In contrast, thicker layers, 80 nm Cu/180 nm Mo, show different deformation behavior with uneven layer thickness reduction, shearing of Cu layers, and bending of the columnar Mo grains beneath the indentation without any sign of shear band and localized phase transformation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microturbulence Suppression by Alfvén Eigenmodes in the DIII-D Tokamak

Mitigation and suppression of low-𝑘 turbulence are observed during the nonlinear evolution of toroidicity-induced Alfvén eigenmodes (TAEs) in DIII-D experiments. Turbulence mitigation begins when the dominant TAE starts to depart from the typical shear Alfvén wave polarization. TAE activity then evolves into a new state, characterized by more discrete coherent modes, increased total amplitude, and more localized radial structures. During this transition, a narrow shear flow layer forms, driven by an enhanced Reynolds stress force, with a shearing rate that exceeds the local turbulence decorrelation rate, leading to a full turbulence suppression. Furthermore, these observations indicate the imbalance between Reynolds and Maxwell stress forces during the nonlinear evolution of the TAE and their important roles in shear flow generation and turbulence reduction.

Alfvén waves↗
DOE OSTI2024

Influences of the South American Low-Level Jet on the Convective Environment in Central Argentina Using a Convection-Permitting Simulation

This study documents the spatial and temporal distribution of the South American low-level jet (SALLJ) and quantifies its impact on the convective environment using a 6.5-month convection-permitting simulation during the Remote Sensing of Electrification, Lightning, And Mesoscale/Microscale Processes with Adaptive Ground Observations and Clouds, Aerosols, and Complex Terrain Interactions (RELAMPAGO-CACTI) campaigns. Overall, the simulation reproduces the observed SALLJ characteristics in central Argentina near the Sierras de Córdoba (SDC), a focal point for terrain-focused upscale growth. SALLJs most frequently occur in the summer with maxima to the northwest and east of the SDC and minima over the higher terrain. The shallower SALLJs (<1750 m) have a strong overnight skew, while the elevated jets are more equally spread throughout the day. SALLJ periods often have higher amounts of low-level moisture and instability compared to non-SALLJ periods, with these impacts increasing over time when the SALLJ is present and decreasing afterward. The SALLJ may enhance low-level wind shear magnitudes (particularly when accounting for the jet height); however, enhancement is somewhat limited due to the presence of speed shear in most situations. SALLJ periods are associated with low-level directional shear favorable for organized convection and an orientation of cloud-layer wind shear parallel to the terrain, which could favor upscale growth. A case study is shown in which the SALLJ influenced both the magnitude and direction of wind shear concurrent with convective upscale growth near the SDC. This study highlights the complex relationship between the SALLJ and its impacts during periods of widespread convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Coupling between Alfven wave and Kelvin-Helmholtz waves in the low latitude boundary layer

The Kelvin-Helmholtz (KH) instability of magnetohydrodynamic surface waves at the low latitude boundary layer is examined using both an eigenfrequency analysis and a time-dependent wave simulation. The analysis includes the effects of sheared flow and Alfven velocity gradient. When the magnetosheath flows are perpendicular to the ambient magnetic field direction, unstable KH waves that propagate obliquely to the sheared flow direction occur at the sheared flow surface when the Alfv\'en Mach number is higher than an instability threshold. Including a shear transition layer between the magnetosphere and magnetosheath leads to secondary KH waves (driven by the sheared flow) that are coupled to the resonant surface Alfven wave. There are remarkable differences between the primary and the secondary KH waves including wave frequency, the growth rate, and the ratio between transverse and the compressional component. The secondary KH wave energy is concentrated near the shear Alfven wave frequency at the magnetosheath with a lower frequency than the primary KH waves. Although the growth rate of the secondary KH waves is lower than the primary KH waves, the threshold condition is lower, so it is expected that these types of waves will dominate at lower Mach number. Because the transverse component of the secondary KH waves is stronger than the primary KH waves, more efficient wave energy transfer from the boundary layer to the inner magnetosphere is also predicted.

Alfven wave↗
DOE OSTI2022

Coupling Between Alfvén Wave and Kelvin–Helmholtz Waves in the Low Latitude Boundary Layer

The Kelvin–Helmholtz (KH) instability of magnetohydrodynamic surface waves at the low latitude boundary layer is examined using both an eigenfrequency analysis and a time-dependent wave simulation. The analysis includes the effects of sheared flow and Alfvén velocity gradient. When the magnetosheath flows are perpendicular to the ambient magnetic field direction, unstable KH waves that propagate obliquely to the sheared flow direction occur at the sheared flow surface when the Alfvén Mach number is higher than an instability threshold. Including a shear transition layer between the magnetosphere and magnetosheath leads to secondary KH waves (driven by the sheared flow) that are coupled to the resonant surface Alfvén wave. There are remarkable differences between the primary and the secondary KH waves, including wave frequency, the growth rate, and the ratio between the transverse and compressional components. The secondary KH wave energy is concentrated near the shear Alfvén wave frequency at the magnetosheath with a lower frequency than the primary KH waves. Although the growth rate of the secondary KH waves is lower than the primary KH waves, the threshold condition is lower, so it is expected that these types of waves will dominate at a lower Mach number. Because the transverse component of the secondary KH waves is stronger than that of the primary KH waves, more efficient wave energy transfer from the boundary layer to the inner magnetosphere is also predicted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Geodetic Evidence for Distributed Shear Below the Brittle Crust of the Walker Lane, Western United States

Abstract Models of active deformation of the Earth's crust are predominantly represented with dislocations having a downdip continuation into the lower crust, where the fault slips continuously. This model predicts surface strain accumulation concentrated near the fault during the interseismic period. In an alternative model, faults do not extend beneath the elastic portion of the crust and are accompanied by a wide zone of distributed shear underneath, predicting a more constant strain rate lacking concentrations at the faults. We use high‐precision GPS data collected across the northern and central Walker Lane, USA— a region of complex faulting near the western edge of the Basin and Range Province to evaluate which model is appropriate. Despite the existence of dense continuous and semi‐continuous geodetic networks that have been surveyed for ∼20 years, the horizontal velocities reveal no evidence of localized strain accumulation across the fault surface expressions. Instead, deformation within the Walker Lane is uniformly linear, suggesting that the surface deformation reflects distributed shear within the ductile crust rather than focused deformation at faults. This suggests no downdip extension of the faults below the seismogenic layer. The shear zone is 172 ± 6 km wide in the northernmost Walker Lane narrowing to 116 ± 4 km in the central Walker Lane. The total velocity budget across the shear zone is 7.2 ± 0.1 mm/yr in the north, increasing to 10.1 ± 0.1 mm/yr in the central Walker Lane. We conclude that assuming the presence of lower crustal dislocations when estimating geodetic faults slip rates may be inappropriate.

Geochemistry & Geophysics↗
DOE OSTI2022

Reduction of transport due to magnetic shear in gyrokinetic simulations of the scrape-off layer

We report the effect of varying magnetic shear on scrape-off layer turbulence and profiles is studied via electromagnetic gyrokinetic simulations of a helical scrape-off layer model. We develop a model helical geometry with magnetic shear and a corresponding field-aligned coordinate system, which is used for simulations with the Gkeyll code. We find that perpendicular transport is reduced in cases with stronger shear, resulting in higher peak particle and heat fluxes to the endplates. Electromagnetic effects slightly increase transport in strong shear cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗