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

The formation of an radial edge electric field due to finite ion orbit width effects is the possible root cause of the H-mode edge

Abstract Full orbit-following simulations of thermal ions show that finite ion-orbit width effects create charge separation near the last closed flux surface (LCFS) which generates a localized radial electric field. Experimentally, edge electric fields are observed in H-mode plasmas and they are necessary for the edge turbulence suppression via the E × B flow shear mechanism. Confined trapped (and to a lesser extent co-passing) ions near the plasma edge form a positive charge distribution outside the LCFS, while thermal electrons are tied more tightly to field lines owing to their small mass and are poorly confined outside the LCFS, hence charge neutrality is violated outside the LCFS. A large number of reported observations from spherical and conventional tokamaks support the results from the simulations although the simulations were not performed fully self consistently. The results suggest ways to lower the H-mode power threshold and optimize the H-mode plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Saturation of fishbone instability through zonal flows driven by energetic particle transport in tokamak plasmas

Gyrokinetic and kinetic-MHD simulations are performed for the fishbone instability in the DIII-D discharge #178631, chosen for validation of first-principles simulations to predict the energetic particle (EP) transport in an ITER prefusion baseline scenario. Fishbone modes are found to generate zonal flows, which dominate the fishbone saturation. The underlying mechanisms of the two-way fishbone-zonal flows nonlinear interplay are discussed in details. Numerical and analytical analyses identify the fishbone-induced EP redistribution as the dominant generation mechanism for zonal flows. The zonal flows modify the nonlinear dynamics of phase space zonal structures, which reduces the amount of EPs able to resonate with the mode, leading to an early fishbone saturation. Simulation results including zonal flows agree quantitatively with DIII-D experimental measurements of the fishbone saturation amplitude and EP transport, supporting this novel saturation mechanism by self- generated zonal flows. Moreover, the wave-particle mode-locking mechanism is shown to determine quantitatively the fishbone frequency down-chirping, as evident in GTC simulation results in agreement with predictions from analytical theory. Finally, the
fishbone-induced zonal flows are possibly responsible for the formation of an ion-ITB in the DIII-D discharge. Based on the low EP transport and the large zonal flow shearing rates associated with the fishbone instability in gyrokinetic simulations of the ITER scenario, it is conjectured that high performance scenarios could be designed in ITER burning plasmas through fishbone-induced ITBs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗
DOE OSTI2021

On the shear-current effect: toward understanding why theories and simulations have mutually and separately conflicted

ABSTRACT The shear-current effect (SCE) of mean-field dynamo theory refers to the combination of a shear flow and a turbulent coefficient β21 with a favourable negative sign for exponential mean-field growth, rather than positive for diffusion. There have been long-standing disagreements among theoretical calculations and comparisons of theory with numerical experiments as to the sign of kinetic ($\beta ^u_{21}$) and magnetic ($\beta ^b_{21}$) contributions. To resolve these discrepancies, we combine an analytical approach with simulations, and show that unlike $\beta ^b_{21}$, the kinetic SCE $\beta ^u_{21}$ has a strong dependence on the kinetic energy spectral index and can transit from positive to negative values at $\mathcal {O}(10)$ Reynolds numbers if the spectrum is not too steep. Conversely, $\beta ^b_{21}$ is always negative regardless of the spectral index and Reynolds numbers. For very steep energy spectra, the positive $\beta ^u_{21}$ can dominate even at energy equipartition urms ≃ brms, resulting in a positive total β21 even though $\beta ^b_{21}\lt 0$. Our findings bridge the gap between the seemingly contradictory results from the second-order-correlation approximation versus the spectral-τ closure, for which opposite signs for $\beta ^u_{21}$ have been reported, with the same sign for $\beta ^b_{21}\lt 0$. The results also offer an explanation for the simulations that find $\beta ^u_{21}\gt 0$ and an inconclusive overall sign of β21 for $\mathcal {O}(10)$ Reynolds numbers. The transient behaviour of $\beta ^u_{21}$ is demonstrated using the kinematic test-field method. We compute dynamo growth rates for cases with or without rotation, and discuss opportunities for further work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Kinematics and dynamics of disclination lines in three-dimensional nematics

An exact kinematic law for the motion of disclination lines in nematic liquid crystals as a function of the tensor order parameter Q is derived. Unlike other order parameter fields that become singular at their respective defect cores, the tensor order parameter remains regular. Following earlier experimental and theoretical work, the disclination core is defined to be the line where the uniaxial and biaxial order parameters are equal, or equivalently, where the two largest eigenvalues of Q cross. This allows an exact expression relating the velocity of the line to spatial and temporal derivatives of Q on the line, to be specified by a dynamical model for the evolution of the nematic. By introducing a linear core approximation for Q, analytical results are given for several prototypical configurations, including line interactions and motion, loop annihilation, and the response to external fields and shear flows. Behaviour that follows from topological constraints or defect geometry is highlighted. Finally, the analytic results are shown to be in agreement with three-dimensional numerical calculations based on a singular Maier–Saupe free energy that allows for anisotropic elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Gyrokinetic electrostatic turbulence close to marginality in the Wendelstein 7-X stellarator

The transition from strong (fluidlike) to nearly marginal (Floquet-type) regimes of ion-temperature-gradient (ITG) driven turbulence is studied in the stellarator Wendelstein 7-X by means of numerical simulations. Close to marginality, extended (along magnetic field lines) linearly unstable modes are dominant, even in the presence of kinetic electrons, and provide a drive that results in finite turbulent transport. A total suppression of turbulence above the linear stability threshold of the ITG modes, commonly present in tokamaks and known as the “Dimits shift,” is not observed. We show that this is mostly due to the peculiar radial structure of marginal turbulence, which is more localized than in the fluid case and therefore less likely to be stabilized by shearing flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Torque and Angular-Momentum Transfer in Merging Rotating Bose-Einstein Condensates

When rotating classical fluid drops merge together, angular momentum can be advected from one to another due to the viscous shear flow at the drop interface. It remains elusive what the corresponding mechanism is in inviscid quantum fluids such as Bose-Einstein condensates (BECs). In this work, we report our theoretical study of an initially static BEC merging with a rotating BECin three-dimensional space along the rotational axis. We show that a solitonlike sheet, resembling a corkscrew, spontaneously emerges at the interface. Furthermore, rapid angular momentum transfer at a constant rate universally proportional to the initial angular-momentum density is observed. Strikingly, this transfer does not necessarily involve fluid advection or drifting of the quantized vortices. We reveal that the corkscrew structure can exert a torque that directly creates angular momentum in the static BEC and annihilates angular momentum in the rotating BEC. Uncovering this intriguing angular momentum transport mechanism may benefit our understanding of various coherent matter-wave systems, spanning from atomtronics on chips to dark matter BECs at cosmic scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Structure-mechanical properties correlation in bulk LiPON glass produced by nitridation of metaphosphate melts

The glassy solid electrolyte Lithium phosphorous oxynitride (LiPON) has been widely researched in thin film solid state battery format due to its outstanding stability when cycled against lithium. In addition, recent reports show thin film LiPON having interesting mechanical behaviors, especially its ability to resist micro-scale cracking via densification and shear flow. Here, in the present study, we have produced bulk LiPON glasses with varying nitrogen contents by ammonolysis of LiPO 3 melts. The resulting compositions were determined to be LiPO 3-3z/2 N z , where 0 ≤ z ≤ 0.75, and the z value of 0.75 is among the highest ever reported for this series of LiPON glasses. The short-range order structures of the different resulting compositions were characterized by infrared, Raman, 31 P magic angle spinning nuclear magnetic resonance, and X-ray photoelectron spectroscopies. Instrumented nano-indentation was used to measure mechanical properties. It was observed that similar to previous studies, both trigonally coordinated (N t ) and doubly bonded (N d ) N co-exist in the glasses in about the same amounts for z ≤ 0.36, the limit of N content in most previous studies. For glasses with z > 0.36, it was found that the fraction of the N t increased significantly while the fraction of N d correspondingly decreased. The incorporation of nitrogen increased both the elastic modulus and hardness of the glass by approximately a factor of 1.5 when N/P ratio reaches 0.75. At the same time, an apparent embrittlement of the glass was observed due to nitridation, which was revealed by nanoindentation with an extra sharp nanoindenter tip.

25 ENERGY STORAGE↗
DOE OSTI2020

Theory of the tertiary instability and the Dimits shift within a scalar model

The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry--Horton system has been proposed by St-Onge [J. Plasma Phys. {\bf 83}, 905830504 (2017)] as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$, where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U''$. Then, depending on $U''$, the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $U'$. Our analytic predictions are in agreement with direct numerical simulations of the modified Terry--Horton system.

fusion plasmas↗
DOE OSTI2023

Optimization and Physics Informed AI/ML on the Cooling Tower Power Output and the One-Dimensional Model

Graduate Research at Missouri S&T Thesis: Theoretical Study of Magnetic Particles in a Shear Flow Subjected to a Uniform Magnetic Field. • Published several journal articles. • Presented at the Indiana University Purdue University Indianapolis graduate seminar. • Presented research findings at American Physical Society: Division of Fluid Dynamics and poster sessions at Missouri S&T for the Chancellor's Distinguished Fellowship.

Sobecki, Christopher A.↗
DOE OSTI2024

Portable Soft X-Ray Diagnostics for Transformative Fusion-Energy Concepts (Final Scientific/Technical Report)

A time-gated extreme-ultraviolet spectrometer was fielded on a transformative fusion-energy device, FuZE-Q at Zap Energy, Inc. (ZEI), a sheared-flow-stabilized z-pinch. Plasma emission spectra were collected in the 5–40 nm wavelength (30–250 eV) range. This was the first time a direct-vacuum-coupled, high-voltage, deposition-sensitive diagnostic was fielded on this type of dense-plasma fusion device. To overcome the special challenges of such a coupling, a unique, highly effective implementation was developed, with fast (10-ns risetime) pulsed high-voltage electronics and a multi-stage differential pumping system. This enabled the UNR spectrometer to collect three independently timed spectra per FuZE-Q shot while also protecting sensitive internal components. Analysis of line emission identified oxygen (N-, C-, B-, Be-, Li-, and He-like O), peaking in intensity shortly after maximum current (>500 kA). Spectral databases and collisional-radiative modeling have been used to identify impurities and estimate edge-plasma temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Global Non-axisymmetric Hall Instabilities in a Rotating Plasma

Non-axisymmetric, flow-driven instabilities in the incompressible Hall-MHD model are studied in a differentially rotating cylindrical plasma. It is found that, in the Hall-MHD regime, both whistler waves and ion-cyclotron waves can extract energy from the flow shear, resulting in two distinct branches of global instability. The non-axisymmetric whistler modes grow significantly faster than non-axisymmetric, ideal MHD modes. A discussion of the global whistler instability mechanism is presented in the large-ion-skin-depth, “electron-MHD” limit. When the magnetic field is azimuthal, a subset of the whistler modes having zero axial wave number are uncovered to be destabilized by the “corotation amplifier” mechanism. It is observed that the effect of the Hall term on the non-axisymmetric modes can be appreciable when d i is on the order of a few percent of the width of the cylindrical annulus. Distinct global modes emerge in the strong Hall-MHD regime at significantly stronger magnetic fields than those required for unstable global MHD modes, as the Hall effect weakens the stabilizing “field-line bending” by decoupling ion motion from the magnetic field. These global non-axisymmetric modes may play an important role in weakly ionized accretion disks.

Alfven waves↗
DOE OSTI2026

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Mega-Gauss Plasma Jet Creation Using a Ring of Laser Beams

Using 20 OMEGA laser beams at the Laboratory for Laser Energetics, University of Rochester, to irradiate a flat plastic target in a hollow ring configuration, we created supersonic cylindrical stable plasma jets with self-generated megagauss magnetic fields extending out to >4 mm. These well-collimated magnetized jets possess a number of distinct and novel properties that will allow us to study the dynamics, physical processes, and scaling properties of astrophysical jets with a dynamic range exceeding those of previous laboratory settings. The dimensionless parameters of these laboratory jets fall in the same regime as those of young stellar object jets. These jets will also provide new versatile laser-based platforms to study magnetized shocks, shear flows, and other plasma processes under controllable conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗