Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shear flow”

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

Superior plastic flow stability of self-patterned carbide – amorphous ceramic nanostructures

Amorphous ceramics and carbides exhibit superb strength but poor plasticity. Here, we synthesized TiC-SiOC nanostructures with TiC-nanocarbides embedded in amorphous ceramic SiOC by co-sputtering followed by high-temperature annealing and/or irradiation. TiC-SiOC nanostructures exhibit high strength and good plastic flow stability even after heavy irradiation, 7 GPa at room temperature and 3.6 GPa at 700 ℃ with a uniform strain of about 10%∼18%. The uniform deformation is accommodated by the shearing of amorphous ceramic and the rotation of nanocarbides. Nanocarbides inhibit the propagation of shear banding in amorphous SiOC, and amorphous-crystal interfaces act as sinks to manage irradiation-induced defects.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Gambusia holbrooki Survive Shear Stress, Pressurization and Avoid Blade Strike in a Simulated Pumped Hydroelectric Scheme

Pumped hydroelectric energy storage (PHES) projects are being considered worldwide as a means of achieving political renewable energy targets in a way that stabilises baseload energy supply from often intermittent renewable energy sources. Unlike a conventional hydroelectric system that only pass water in a downstream direction, a feature of PHES is that it relies on the bi-directional flow of water. In some cases this flow can be across different waterbodies or catchments, posing a risk of inadvertently expanding the range of aquatic animals like fish. The risk of this happening depends on the likelihood of survival of individuals, which remains poorly understood for turbines that are pumping rather than generating. This study quantified the survival of a globally widespread and invasive poeciliid fish, Eastern gambusia (Gambusia holbrooki) when exposed hydraulic stresses characteristic of what would be experienced through a PHES during the pumping phase. A shear flume and hyperbaric chamber were used to expose fish to different strain rates and rapid and sustained pressurisation. A blade strike model was also used to predict survival of fish passing through a Francis dual turbine / pump. The ranges simulated were based on design and operational conditions provided for a PHES scheme being proposed in south-eastern Australia. All gambusia tested survived extremely high (up to 7600 kPa gauge pressure) pressurisation, high levels of shear stress (up to 1853 s -1 ), and the majority (> 93 %) were unlikely to be struck by a turbine blade. Given their tolerance to these extreme simulated stresses, we conclude that gambusia will likely survive passage through a PHES scheme. Therefore, where a new PHES poses the risk of inadvertently expanding the range of gambusia or similar poeciliid species, measures to minimise their spread or mitigate their ecosystem impacts should be considered.

25 ENERGY STORAGE↗
DOE OSTI2021

The influence of surface heat fluxes on the growth of idealized monsoon depressions

This study explores the effect of surface sensible and latent heat fluxes on monsoon depressions using a series of idealized convection-permitting simulations. Each experiment is initialized with a small-amplitude wave that is allowed to grow within an environment representative of the South Asian monsoon. Comparing experiments with and without interactive surface heat fluxes, it is found that these fluxes enhance the growth of the simulated vortices. Without interactive surface fluxes, the strengthening period is short and the vortices fail to reach intensities characteristic of stronger monsoon depressions. Using a large set of experiments in which the vertical and meridional shear are systematically varied, it is found that surface heat fluxes enhance intensity the most when the upper-level shear is weak, the lower-level shear and associated moist static energy (MSE) gradient are sufficiently steep, and the lower-level meridional shear is strong. Here, these experiments reveal two different regimes of convection-coupled monsoon depression growth: one in which convection is driven by MSE advection and one in which it is driven by surface heat fluxes and quasigeostrophic forcing for ascent. Both regimes require sufficiently strong meridional shear to achieve initial growth by barotropic instability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Turbulence model reduction by deep learning

A defining problem of turbulence theory is to produce a predictive model for turbulent fluxes. These have profound implications for virtually all aspects of the turbulence dynamics. In magnetic confinement devices, drift-wave turbulence produces anomalous fluxes via cross-correlations between fluctuations. In this work, we introduce an alternative, data-driven method for parametrizing these fluxes. The method uses deep supervised learning to infer a reduced mean-field model from a set of numerical simulations. We apply the method to a simple drift-wave turbulence system and find a significant new effect which couples the particle flux to the local gradient of vorticity. Notably, here, this effect is much stronger than the oft-invoked shear suppression effect. We also recover the result via a simple calculation. The vorticity gradient effect tends to modulate the density profile. In addition, our method recovers a model for spontaneous zonal flow generation by negative viscosity, stabilized by nonlinear and hyperviscous terms. We highlight the important role of symmetry to implementation of the new method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗
DOE OSTI2021

Extreme wind shear events in US offshore wind energy areas and the role of induced stratification

As the offshore wind industry emerges on the US East Coast, a comprehensive understanding of the wind resource – particularly extreme events – is vital to the industry's success. Such understanding has been hindered by a lack of publicly available wind profile observations in offshore wind energy areas. However, the New York State Energy Research and Development Authority recently funded the deployment of two floating lidars within two current lease areas off the coast of New Jersey. These floating lidars provide publicly available wind speed data from 20 to 200 m height with a 20 m vertical resolution. In this study, we leverage a year of these lidar data to quantify and characterize the frequent occurrence of high-wind-shear and low-level-jet events, both of which will have a considerable impact on turbine operation. In designing a detection algorithm for these events, we find that the typical, non-dimensional power-law-based wind shear exponent is insufficient to identify many of these extreme, high-wind-speed events. Rather, we find that the simple vertical gradient of wind speed better captures the events. Based on this detection method, we find that almost 100 independent events occur throughout the year with mean wind speed at 100 m height and wind speed gradient of 16 m s -1 and 0.05 s -1 , respectively. The events have strong seasonal variability, with the highest number of events in summer and the lowest in winter. A detailed analysis reveals that these events are enabled by an induced stable stratification when warmer air from the south flows over the colder mid-Atlantic waters, leading to a positive air–sea temperature difference.

17 WIND ENERGY↗
DOE OSTI2022

Large-eddy simulation of helical- and straight-bladed vertical-axis wind turbines in boundary layer turbulence

Here, turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines (VAWTs) in boundary layer turbulence are numerically studied using the large-eddy simulation (LES) method combined with the actuator line model. Based on the LES data, systematic statistical analyses are performed to explore the effects of blade geometry on the characteristics of the turbine wake. The time-averaged velocity fields show that the helical-bladed VAWT generates a mean vertical velocity along the center of the turbine wake, which causes a vertical inclination of the turbine wake and alters the vertical gradient of the mean streamwise velocity. Consequently, the intensities of the turbulent fluctuations and Reynolds shear stresses are also affected by the helical-shaped blades when compared with those in the straight-bladed VAWT case. The LES results also show that reversing the twist direction of the helical-bladed VAWT causes the spatial patterns of the turbulent wake flow statistics to be reversed in the vertical direction. Moreover, the mass and kinetic energy transports in the turbine wakes are directly visualized using the transport tube method, and the comparison between the helical- and straight-bladed VAWT cases show significant differences in the downstream evolution of the transport tubes.

17 WIND ENERGY↗
DOE OSTI2026

A Multiscale Simulation Study of the Effects of Nantucket Island on Atmospheric Boundary Layer Flow

Multiscale atmospheric simulations were conducted to investigate how Nantucket Island modifies the marine atmospheric boundary layer (ABL) flow, how the mechanisms responsible are represented in mesoscale versus large-eddy simulation (LES) domains, and potential impacts on wind energy production downstream, in moderate and high wind speed cases. Here, to isolate the effects of the island, each case study was simulated (i) with and without island topography included and (ii) with observed versus spatially uniform sea surface temperature (SST) values, the latter removing the influence of spatial variability in SST, while also enabling the examination of how the island effects vary with near-surface thermodynamic stability over the surrounding ocean. Island effects were found to be strongest during daytime heating, which generates convective rolls that redistribute momentum vertically over the island and downstream. Colder SSTs and weaker winds strengthen these rolls and extend their influence. Flow statistics within a representative wind turbine rotor-swept area downstream show that while vertical momentum redistribution enhances near-surface wind speeds, those aloft decrease, resulting in minimal net change to rotor-equivalent wind speed and slightly reduced vertical shear. Comparisons of mesoscale and LES domains reveal that both capture the roughness-driven near-surface slowdown, but mesoscale models underestimate vertical momentum transport compared to LES, which explicitly resolves convective rolls. Increasing LES resolution beyond the grid spacing required to begin capturing the largest convective roll structures ∼ O(100) m yields minimal additional momentum transport but increases resolved turbulence kinetic energy, relevant for stress loading.

Atmosphere-land interaction↗
DOE OSTI2026

Spherical tokamak physics research in preparation for the operation of NSTX-U

The National Spherical Torus Experiment Upgrade (NSTX-U) is preparing to resume operation, representing a crucial step toward realizing compact, cost-effective fusion pilot plants. In advance of this, extensive modeling and data analysis have been conducted to advance the physics basis for low-aspect-ratio, high-performance plasma regimes, focusing on three core objectives: confinement and stability, power and particle handling, and steady-state operation. Significant progress has been made in understanding the electron temperature flattening in high-β plasmas, which is shown to be driven by a complex interplay of magnetohydrodynamic instabilities (e.g. non-resonant infernal modes), fast-ion-driven Alfvén eigenmodes, and electron and ion-scale micro-instabilities, particularly Kinetic Ballooning Modes (KBMs), whose destabilization is strongly dependent on parallel magnetic field fluctuations (δB ∥ ). Furthermore, a new gyrokinetic critical pedestal model was developed, accurately predicting pedestal structure by identifying KBMs as the primary stability limit, offering a critical constraint for future high-confinement scenarios. To address the challenge of high heat flux, novel liquid lithium plasma-facing components were modeled. The analysis confirmed that lithium vapor shielding is a self-regulating mechanism for heat mitigation, while also emphasizing that strong main ion parallel flow is essential to minimize core lithium contamination. Finally, progress toward steady-state operation was anchored by developing the required physics basis and control tools. This includes predictive modeling for reversed magnetic shear sustainment, demonstrating that magnetic island-induced bootstrap current reduction is negligible in STs, and advancing real-time control and disruption avoidance capabilities. The development of high-speed surrogate models (e.g. MMMNet) provides computationally efficient tools vital for non-inductive scenario optimization and integrated, low-disruptivity operations planned for NSTX-U.

NSTX-U↗
DOE OSTI2023

Near-Cloud Atmospheric Ingredients for Deep Convection Initiation

A lack of routine environmental observations located near deepening cumulus congestus clouds limits verification of important theorized and simulated updraft–environment interaction processes occurring during deep convection initiation (CI). We analyze radiosonde profiles collected during several hundred CI events near a mountain range in central Argentina during the CACTI field campaign. Statistical analyses illustrate environmental conditions supporting radar-observed CI outcomes that span a spectrum of convective cell depths, widths, and durations, as well as events lacking precipitating convection. Tested environmental factors include a large variety of sounding-derived measurements of CAPE, CIN, moisture, terrain-relative winds, vertical shear, and lifted parcel properties, with supplemental model reanalysis of background larger-scale vertical motion. CAPE and CIN metrics do not consistently differentiate CI success from failure. Only a few environmental factors contain consistent monotonic relationships among the spectrum of cloud depths achieved during CI: (i) the depth and strength of background ascent, and (ii) the component of low-level flow oriented parallel to the ridgeline. These metrics suggest that the ability of the surrounding flow to lift parcels to their LFC and terrain-modified flow are consistently relevant processes for CI. Low- to midlevel relative humidity strongly discriminated between CI and non-CI events, likely reflecting entrainment-driven dilution processes. However, we could not confidently conclude that relative humidity similarly discriminated robust from marginal CI events. Circumstantial evidence was found linking cell width, an important cloud property governing the probability of CI, to LCL height, boundary layer depth, depth and magnitude of the CIN layer, and ambient wind shear.

Storm environments↗
DOE OSTI2020

The Sensitivity of Joint Inversions of Seismic and Geodynamic Data to Mantle Viscosity

Seismic tomography has revealed the existence of large-scale velocity heterogeneities in the mantle. The interpretation of seismic velocity anomalies in terms of temperature and chemical composition is nonunique. We use geodynamic observations including gravity, plate motions, dynamic topography, and excess ellipticity of the core-mantle boundary combined with seismic observations to investigate the thermo-chemical structure of the mantle through joint inversions. An outstanding issue, however, is the physical connection between mantle density anomalies and the surface geodynamic observations, which requires knowledge of the mantle viscosity structure. Here we perform joint inversions assuming different viscosity profiles and examine the dependence of the results on the viscosity. We first assume that mantle heterogeneity is due to thermal variations, which places a constraint on the relation between seismic velocity and density, and we subsequently relax the constraint to allow for potential nonthermal effects. In all of our joint inversions, a nonthermal origin of density anomalies is required to explain the geodynamic data, though the amount varies with the assumed viscosity structure. A common observation is a high-density chemical signal in the center of the large low-shear-velocity provinces at the base of the mantle resulting in a near neutral or slightly dense overall buoyancy there. Using the derived density models and their corresponding viscosity profiles, we also calculate instantaneous mantle flow fields. The predicted flow fields derived from joint inversions are generally similar but are quite different from flow fields using density models derived from a posteriori scaling of pure seismic tomography models.

58 GEOSCIENCES↗
DOE OSTI2023

New Insights from Legacy Seismic Data regarding Basalt Elevations and Variability on the Hanford Site

Migration of groundwater contaminants in the Gable Gap area of the Hanford Site in southeastern Washington State is strongly influenced by the distribution and permeability of basalts that lie beneath an unconfined aquifer. Locally, folding and faulting of the Columbia River Basalt associated with the Yakima fold and thrust belt followed by erosion due to the Lake Missoula floods resulted in a complex basalt surface that represents either an impermeable lower boundary to the unconfined aquifer system or localized regions of increased permeability that potentially promote communication between the unconfined aquifer system and deeper, confined aquifer systems. Paleo-channels carved into the basalt by floodwaters are thought to provide preferential flow paths for groundwater contaminants. In 2011, a seismic landstreamer campaign was carried out to image the basalt surface and produced pre-stack depth migrated p-wave reflection images. The reflection images identified two large troughs that may represent paleo-channels and several areas of possible faulting. Here, the streamer data are re-analyzed using refraction travel-time and Rayleigh wave dispersion analyses to obtain images of compressional and shear wave velocities within the suprabasalt sediment sections and the upper basalt surface. The combined interpretation of reflection and seismic velocity images shows complexity in the basalt velocity and elevation, which varies by 50 m or more within the study area. These results, along with other ongoing geophysical investigations, will be used to inform the site geologic model and potentially guide placement of future boreholes needed to quantify vertical flow between the confined and unconfined aquifers.

58 GEOSCIENCES↗
DOE OSTI2024

Self-Focused Pulse Propagation Is Mediated by Spatiotemporal Optical Vortices

We show that the dynamics of high-intensity laser pulses undergoing self-focused propagation in a nonlinear medium can be understood in terms of the topological constraints imposed by the formation and evolution of spatiotemporal optical vortices (STOVs). STOVs are born from point phase defects on the sides of the pulse nucleated by spatiotemporal phase shear. These defects grow into closed loops of spatiotemporal vorticity that initially exclude the pulse propagation axis, but then reconnect to form a pair of toroidal vortex rings that wrap around it. STOVs constrain the intrapulse flow of electromagnetic energy, controlling the focusing-defocusing cycles and pulse splitting inherent to nonlinear pulse propagation. Here, we illustrate this in two widely studied but very different regimes, relativistic self-focusing in plasma and nonrelativistic self-focusing in gas, demonstrating that STOVs mediate nonlinear propagation irrespective of the detailed physics.

42 ENGINEERING↗
DOE OSTI2020

Investigation of the Elliptic Flow Fluctuations of the Identified Particles Using the a Multi-Phase Transport Model

A Multi-Phase Transport (AMPT) model is used to study the elliptic flow fluctuations of identified particles using participant and spectator event planes. The elliptic flow measured using the first order spectator event plane is expected to give the elliptic flow relative to the true reaction plane which suppresses the flow fluctuations. However, the elliptic flow measured using the second-order participant plane is expected to capture the elliptic flow fluctuations. Our study shows that the first order spectator event plane could be used to study the elliptic flow fluctuations of the identified particles in the AMPT model. The elliptic flow fluctuations magnitude shows weak particle species dependence and transverse momentum dependence. Such observation will have important implications for understanding the source of the elliptic flow fluctuations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Linking structural and rheological memory in disordered soft materials

Linking the macroscopic flow properties and nanoscopic structure is a fundamental challenge to understanding, predicting, and designing disordered soft materials. Under small stresses, these materials are soft solids, while larger loads can lead to yielding and the acquisition of plastic strain, which adds complexity to the task. In this work, we connect the transient structure and rheological memory of a colloidal gel under cyclic shearing across a range of amplitudes via a generalized memory function using rheo-X-ray photon correlation spectroscopy (rheo-XPCS). Our rheo-XPCS data show that the nanometer scale aggregate-level structure recorrelates whenever the change in recoverable strain over some interval is zero. The macroscopic recoverable strain is therefore a measure of the nano-scale structural memory. We further show that yielding in disordered colloidal materials is strongly heterogeneous and that memories of prior deformation can exist even after the material has been subjected to flow.

Kamani, Krutarth M. [Univ. of Illinois at Urbana-C↗
DOE OSTI2026

Reverse segregation and self-organization in inclined chute flows of bidisperse granular mixtures

In the usual segregation scenario for stable inclined chute flows of bidisperse mixtures of fine and coarse spherical particles, coarse particles rise toward the free surface, forming a coarse-rich region atop the flowing pile. Beyond a threshold coarse-to-fine diameter ratio of approximately 4, conversely, the weight of the coarse particles exceeds the segregation driving forces, causing individual coarse particles to sink within the pile and producing a reversed segregation state. However, an understanding of the collective evolution of the pile structure is still lacking when the particle diameter ratio exceeds 4 and the coarse-particle mass fraction is appreciable. To explore this broadly bidisperse limit, we perform discrete element method simulations considering mean particle diameter ratios of up to 8 and coarse-particle mass fractions spanning 0.1 to 0.9. The steady-state flow profiles reveal several intriguing behaviors that depend on the diameter ratio and mass fraction. These include a previously identified transition from usual to reverse segregation and a newfound tendency to self-organize into alternating coarse- and fine-rich particle layers stacked along the shear gradient direction, with layer thickness dictated by the coarse-particle diameter. A fuller understanding of segregation at this scale could pave the way for enhanced mixing or demixing techniques at the commercial scale.

granular flow↗