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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Non-local model for surface tension in fluid-fluid simulations

We propose a non-local model for surface tension obtained in the form of an integral of a molecular-force-like function with support 3.5e added to the Navier-Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non-local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference Pe,in - Pe,out satisfies the Young-Laplace law for droplet radius greater than 3.5e and other- wise deviates from the Young-Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non-local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier-Stokes model subject to the Young-Laplace boundary condition at the fluid-fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non- local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid-size dependent.

Howard, Amanda A.↗
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 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 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 OSTI2021

Suppression of ITG turbulence due to spectral shift during biasing induced H-mode on HBT-EP

Investigations of biasing induced H-mode discharges on HBT-EP show that the edge turbulence is consistent with the ion temperature gradient) mode and have allowed for the controlled observation of the effect of applied flow shear on the turbulence. Measurements of the radial wavenumber spectrum of floating potentials at the edge show that the turbulence intensity decreases with increasing shift in the spectrum average ⟨ k r ⟩ when increasing amounts of bias probe voltage (and increasing amounts of flow shear) is applied. This is in agreement with the spectral shift model [Staebler et al., Phys. Rev. Lett. 110, 055003 2013] for turbulence suppression via sheared flow. A shift in the wavenumber spectrum occurs at applied electrode voltages and currents below the threshold needed for an L–H transition, and a dithering transition is obtained when biasing near the threshold. Suppression of blob-filament turbulence in the scrape-off layer (SOL) precedes the L–H transition, with the SOL turbulence remaining low throughout the dithering phase, despite the modulation of turbulence levels in the nearby edge. This demonstrates that the SOL turbulence “decouples” from the edge turbulence. The spectral shift in the measured radial wavenumber is corroborated by the direct measurement of eddy tilt angle using a novel time delay analysis technique first developed for Doppler reflectometry [Pinzon et al., Plasma Phys. Controlled Fusion 61, 105009 (2019)] but adapted here for floating potential measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonequilibrium design strategies for functional colloidal assemblies

Here, we use a nonequilibrium variational principle to optimize the steady-state, shear-induced interconversion of self-assembled nanoclusters of DNA-coated colloids. Employing this principle within a stochastic optimization algorithm allows us to identify design strategies for functional materials. We find that far-from-equilibrium shear flow can significantly enhance the flux between specific colloidal states by decoupling trade-offs between stability and reactivity required by systems in equilibrium. For isolated nanoclusters, we find nonequilibrium strategies for amplifying transition rates by coupling a given reaction coordinate to the background shear flow. We also find that shear flow can be made to selectively break detailed balance and maximize probability currents by coupling orientational degrees of freedom to conformational transitions. For a microphase consisting of many nanoclusters, we study the flux of colloids hopping between clusters. We find that a shear flow can amplify the flux without a proportional compromise on the microphase structure. This approach provides a general means of uncovering design principles for nanoscale, autonomous, functional materials driven far from equilibrium.

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

Effects beyond ideal MHD on stability of wide and enhanced pedestal regimes in NSTX

Stability of edge-localized modes (ELMs) in spherical tokamaks is explored using the extended MHD model. Linear NIMROD simulations have been performed for three NSTX discharges 132543, 132588, and 141133, to investigate the role of resistivity, diffusivity, and shear flows on the onset of ELMs. The first discharge represents the wide pedestal regime and the later two discharges represent ELM-free enhanced pedestal H-mode. We first present the effect of toroidal rotation shear and find a flow shear destabilizing effect in these NSTX discharges. Simulations are also extended to include the two-fluid and ion gyroviscosity effects. Simulations show that the flow shear can shift the mode spectrum and alter the critical condition of ELM onset. We also uncover that ELM onset prediction in spherical tokamaks requires effects beyond MHD, in particular gyroviscosity and diamagnetic terms could stabilize Peeling-Ballooning modes consistent with the experimental observation of ELM-free regimes in NSTX. The findings give new insight into the nature of the interplay between resistivity, flows, and diamagnetic stabilization in ELM suppression and have potential applications to ELM control schemes in NSTX-U and next-generation spherical tokamaks. This study identifies the essential physical effects that must be included in future predictive and validation simulations.

NSTX↗
DOE OSTI2025

Mechanically Accelerated Depolymerization of Entangled Linear Polymer Melts

Mechanical forces can enhance the chemical depolymerization of synthetic polymers when shear flow accelerates chain scission. To quantify the extent of mechanically-accelerated scission, the effect of simple shear flow (duration and strength) with low Weissenberg and Deborah numbers was investigated by considering the impact of applied work in both simple shear and shear dominated mixed flows. Hydrogenated polyisoprene was chosen as a model linear, entangled system. The conditions (strain amplitude, frequency, and shearing time) necessary to increase chain scission were assessed in the rubbery melt. Shear flow accelerated chain scission at higher temperatures, suggesting an activated process. Isothermal scission versus work curves were superposed by applying shift factors a T,S , whose Arrhenius-like temperature dependence gave an apparent activation energy for chain scission of ~ 110 kJ/mol, which is likely a combination of the activation energy of viscosity and bond energy. This work provides a base for quantifying the impact of shear on depolymerization of polymer melts and highlight the connection between viscous dissipation and scission chemistry.

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

Linear and Nonlinear Kelvin–Helmholtz Instability and Magnetohydrodynamic Wave Emission in Sheared Astrophysical Plasma Flows

The evolution of the Kelvin–Helmholtz instability (KHI) and magnetohydrodynamic (MHD) wave emission is investigated at shear-flow boundaries of magnetized plasmas. While MHD wave emission has been suggested to be only possible during the nonlinear stages, we find that there is also significant wave emission during the KHI's linear stages. These emitted MHD waves may have stronger impacts than KHI surface waves since they can act to transport energy away from the local region of the shear flow. The removal of energy from the shear-flow region, instead of just the local redistribution considered in previous studies, and its propagation away from the interface could have major implications for the evolution of astrophysical objects characterized by fast plasma flow shears.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Gouge formation and dilation impacts to flow during fracture shearing

We report flow in fractures is of critical importance to fluid flow in low-permeability shales, particularly when they restrict vertical migration from reservoirs. The Mt. Simon Sandstone is a target for geologic CO 2 sequestration and the overlying Eau Claire Formation is the sealing formation capping this reservoir. The lower Eau Claire, a primary barrier to upward migration of CO 2 from the Mt. Simon, is lithologically heterogenous with a mixture of sandstone and shale. We use a series of flow-through experiments on fractured samples from the Eau Claire to document the morphological changes taking place in the fracture and surrounding rock as it is sheared in a stepwise fashion. The samples were imaged using computed tomography after each shearing episode and flow resistance was calculated based on the pressure differential across each core. We found that gouge accumulations in the fracture plane led to specific instances of self-sealing behavior despite the fact that most samples increased in average fracture aperture and transmissivity. Gouge limited the increase or, in some cases, decreased flow through these samples. Fractures adjacent to shale-rich zones and concordant with bedding were more likely to show evidence of self-sealing, while increased sandstone content and fracture orientation perpendicular to bedding planes was more likely to result in fracture dilation and fracture permeability increases. Reactivation and shearing of vertical fractures in lithologically heterogenous caprocks is more likely to pose risks to long term CO 2 storage scenarios.

58 GEOSCIENCES↗
DOE OSTI2024

Stability and transport of gyrokinetic critical pedestals

Abstract A gyrokinetic threshold model for pedestal width–height scaling prediction is applied to multiple devices. A shaping and aspect ratio scan is performed on National Spherical Torus Experiment (NSTX) equilibria, finding Δ ped = 0.92 A 1.04 κ − 1.24 0.38 δ β θ , ped 1.05 for the wide-pedestal branch with pedestal width Δ ped , aspect ratio A , elongation κ , triangularity δ , and normalized pedestal height β θ , ped . The width–transport scaling is found to vary significantly if the pedestal height is varied either with a fixed density or fixed temperature, showing how fueling and heating sources affect the pedestal density and temperature profiles for the kinetic-ballooning-mode (KBM) limited profiles. For an NSTX equilibrium, at fixed density, the wide branch is Δ ped = 0.028 ( q e / Γ e − 1.7 ) 1.5 ∼ η e 1.5 and at fixed temperature Δ ped = 0.31 ( q e / Γ e − 4.7 ) 0.85 ∼ η e 0.85 , where q e and Γ e are turbulent electron heat and particle fluxes and η e = ∇ ln ⁡ T e / ∇ ln ⁡ n e for an electron temperature T e and density n e . Pedestals close to the KBM limit are shown to have modified turbulent transport coefficients compared to the strongly driven KBMs. The role of flow shear is studied as a width–height scaling constraint and pedestal saturation mechanism for a standard and lithiated wide pedestal discharge. Finally, the stability, transport, and flow shear constraints are combined and examined for an NSTX experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2023

Statistical Analysis of the Limiting Dynamics of Two dimensional Boussinesq Turbulent Systems

In the current study, we investigate rotational and stratified turbulent systems within the framework of the Boussinesq equations. The effects of rotation and stratification give rise to a variety of dynamical regimes in atmospheric and oceanic turbulence. To develop a better understanding of the fundamental flow properties through a more straightforward mathematical framework, we commence our research with a simplified two-dimensional Boussinesq model that incorporates rotation and stratification effects. We explore two distinct limiting dynamics, considering pure rotation and pure stratification separately. We examine these intriguing dynamics through numerical investigations, applying the exact solution theory of Boussinesq equations and examining small-scale perturbations to the exact solution under rotation and stratification. We conduct an analysis of statistical quantities associated with the turbulent state variables. The results indicate that, under weak rotation at the equilibrium statistical state, the flow field exhibits vortex flow. However, strong rotation transforms the flow into vertical shear flow. We observe both downscale and upscale energy transfers with a decay rate of k −3 in the absence of external forcing and dissipation. Furthermore, we delve into the study of internal gravity wave mode interactions in the simplified two-dimensional Boussinesq system, which is a crucial aspect of geophysical turbulence due to gravity. Our findings reveal that the time series of the mode coefficients exhibit wave-like interactions between wave modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Different responses of the Rayleigh–Taylor type and resistive drift wave instabilities to the velocity shear

The effects of velocity shear on the unstable modes driven by the effective gravity (Rayleigh–Taylor and interchange) and resistive drift wave instabilities for inhomogeneous equilibrium fluid/plasma density are analyzed for the localized eigenmode problems. It is shown that the effect of the velocity shear drastically depends on the type of instability. Whereas the velocity shear can significantly suppress both Rayleigh–Taylor and interchange instabilities, it has only a weak impact on the growth rate of the resistive drift wave. This is directly related to the physical nature of these instabilities. For the Rayleigh–Taylor and interchange instabilities, the shear flow tilts the eddies of the stream functions, while for the resistive drift wave instability the shear flow simply shifts the eddies in the radial direction with no tilting. However, for a large velocity shear, the eigenmode solutions for resistive drift waves cease to exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Development of Integrated Education, Training and Research Program with Interdisciplinary Applications

Plasma models have been formulated including realistic spatial profiles of both flow and radio frequency induced ponderomotive force. With these inclusions the picture of stability of various plasma and fluid instabilities, as expected, changed drastically with ground-breaking consequences. The inhomogeneous parallel flow and the radio frequency waves can actually shown to stabilize turbulence. This is different from the prevalent notion that both parallel flow shear and radio frequency waves are responsible for the excitation (destabilization) of plasma turbulence. This has several ground-breaking consequences:- (1) the stabilization by parallel flow clearly goes against the conventional notion of the origin of ionospheric oscillation which invokes parallel flow destabilization as the origin, (2) the stabilization by parallel flow opens us a new avenue for improved mode formation in fusion devices - which mostly rely on the perpendicular flow shear stabilization for improved mode formation but the perpendicular flow is damped in a tokamak - so the improved mode formed by the parallel flow can sustain longer and has more prospect for ignition, (3) the complete stabilization of the ITG mode (and consequent suppression of transport) not only explain many unknown phenomena in the space physics but it also raises a prospect for transport barrier formation by the RF waves but not by the RF induced flow (as most works suggest) which is never observed in a tokamak of that magnitude to create a barrier. These are indeed ground-breaking consequences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of Integrated Education, Training and Research Program with Interdisciplinary Applications II

Plasma models have been formulated including realistic spatial profiles of both flow and radio frequency induced ponderomotive force. With these inclusions the picture of stability of various plasma and fluid instabilities, as expected, changed drastically with ground-breaking consequences. The inhomogeneous parallel flow and the radio frequency waves can actually shown to stabilize turbulence. This is different from the prevalent notion that both parallel flow shear and radio frequency waves are responsible for the excitation (destabilization) of plasma turbulence. This has several ground-breaking consequences:- (1) the stabilization by parallel flow clearly goes against the conventional notion of the origin of ionospheric oscillation which invokes parallel flow destabilization as the origin, (2) the stabilization by parallel flow opens us a new avenue for improved mode formation in fusion devices - which mostly rely on the perpendicular flow shear stabilization for improved mode formation but the perpendicular flow is damped in a tokamak - so the improved mode formed by the parallel flow can sustain longer and has more prospect for ignition, (3) the complete stabilization of the ITG mode (and consequent suppression of transport) not only explain many unknown phenomena in the space physics but it also raises a prospect for transport barrier formation by the RF waves but not by the RF induced flow (as most works suggest) which is never observed in a tokamak of that magnitude to create a barrier. These are indeed ground-breaking consequences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2022

Control of internal transport barriers in magnetically confined tokamak burning plasmas

The initiation, termination, and control of internal transport barriers associated with E × B flow shear near local minima of magnetic shear are examined for burning plasmas to determine if the positive feedback loops between profiles, instability, transport, and flow shear operate in regimes with fusion self-heating. A five-field transport model for the evolution of profiles of density, ion and electron temperature, ion and electron fluctuations, and radial electric field is utilized to examine the efficacy of controls associated with external inputs of heat and particles, including neutral beam injection, RF, pellets, and gas puffing. The response of the plasma to these inputs is studied in the presence of self-heating. The latter is affected by the external inputs and their modification of profiles and is, therefore, not an external control. Provided sufficient external power is applied, internal transport barriers can be created and controlled, both in ion and electron channels. Barrier control is sensitive to the locations of power deposition and pellet ablation, as well as temporal sequencing of external inputs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗