Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Damping”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Zonal flow screening in negative triangularity tokamaks

Here, this paper presents a comparative study of zonal flow screening in positive and negative triangularity tokamaks. Neoclassical screening sets the strength of zonal flow shear for fixed drive and damping. Orbit calculations show that the radial excursions of trapped particles are reduced in negative triangularity tokamaks, as compared to positive triangularity. Yet surprisingly, the neoclassical dielectric susceptibility actually increases with decreasing triangularity, such that the susceptibility is higher for negative triangularity than for positive triangularity. This is because the reduction in neoclassical polarization by shrinking the banana width is offset by the increase in neoclassical polarization by the enhancement of trapped fraction for negative triangularity. As a result, the zonal flow screening length is actually enhanced for negative triangularity, as compared to positive triangularity. Hence, the zonal flow residual is smaller for negative triangularity than for positive triangularity. Results from gyrokinetic simulations support the analytic calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nonlinear saturation of reverse shear Alfvén eigenmodes induced by kinetic electrons

A strong nonlinear saturation mechanism of the reverse shear Alfvén eigenmode (RSAE) induced by thermal electrons is observed in gyrokinetic $\delta$$f$ particle-in-cell simulation. This nonlinear effect occurs at moderate mode amplitude, $\delta$$B$ ⟂ /$B$ < 10 -3 , and is associated with the electron parallel streaming along the perturbed magnetic field lines, the term v $\|$ $\frac{δΒ_{⟂}}{B}$ ∙ ∇$\delta$$f$, in the electron drift-kinetic equation. In the case of an n = 4 RSAE, where n is the toroidal mode number, this magnetic fluttering nonlinearity leads to strong mode coupling and broadening of the saturation spectrum. The nonlinear components, which are highly damped through the interaction between thermal electrons and the parallel electric field, modifies the RSAE mode structure and suppresses energetic particle drive significantly. The split-weight scheme (Chen and Parker 2007 J. Comput. Phys. 220 839–55) is used to simulate kinetic electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Magnetic island formation and rotation braking induced by low-Z impurity penetration in an EAST plasma

Abstract Recent observations of the successive formations of the 4 / 1 , 3 / 1 , and 2 / 1 magnetic islands as well as the subsequent braking of the 2 / 1 mode during a low- Z impurity penetration process in EAST experiments are well reproduced in our 3 D resistive MHD simulations. The enhanced parallel current perturbation induced by impurity radiation predominately contributes to the tearing mode growth, and the 2 / 1 island rotation is mainly damped by the impurity accumulation as results of the influence from high n modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Analysis of beam ion driven Alfvén eigenmode stability induced by Tungsten contamination in EAST

Alfvén eigenmodes (AE) activity is observed in the EAST high β N and low B T discharge 93910, operation scenario dedicated to explore the ITER baseline scenario. AEs are triggered after the plasma is contaminated by Tungsten that causes an abrupt variation of the thermal plasma and energetic particles (EPs) profiles. The aim of the present study is to analyze the AE stability in the 93910 discharge using the gyro-fluid code FAR3d, identifying the AE stability trends by comparing the plasma before and after the Tungsten contamination. Tungsten contamination causes the destabilization of Toroidal AEs (TAE) and Energetic particle modes (EPMs) in the same frequency range and radial location with respect to the experimental observation and M3D-K/GTAW code results. Next, a set of parametric studies are performed to analyze the effect of the thermal plasma and EP parameters on the AE stability. The analysis indicates a lower EP β threshold for the AEs destabilization if the EP energy increases, an improved AE stability of on-axis NBI configurations due to the stronger continuum damping in the inner plasma region as well as a large enhancement of the EP drive as the thermal ion density increases due to a higher ratio of the EP and Alfven velocities. Consequently, the simulations indicate the increment of the thermal ion density after the Tungsten contamination could be the main cause of the AE/EPM destabilization.

AE↗
DOE OSTI2024

Quasi-mode evolution in a stochastic magnetic field

Abstract We present a multi-scale model of quasi-mode evolution in a stochastic magnetic field. The similarity between a quasi-mode and a ballooning mode enables us to address the challenges arising from the disparate geometries in the theories of ballooning modes in the presence of resonant magnetic perturbations. We obtain useful insights into our understanding of ballooning mode dynamics in a stochastic background. To maintain quasi-neutrality at all scales, the beat between the quasi-mode and the stochastic magnetic field drives microturbulence, which drives the turbulent background that promotes mixing and damps the quasi-mode. As a result of the broad mode structure of the quasi-mode, the turbulent viscosity and the turbulent diffusivity produced by the microturbulence are larger than those in our related study on resistive interchange modes. The stochastic magnetic field can also enhance the effective plasma inertia and reduce the effective drive, thereby slowing the mode growth. A nontrivial correlation between the microturbulence and the magnetic perturbations is shown to develop. This could account for the reduction in the Jensen–Shannon complexity of pedestal turbulence in the Resonant Magnetic Perturbation Edge-Localized Mode suppression phase observed in recent experiments. Directions for future experimental and theoretical studies are suggested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overview of Large Helical Device experiments of basic plasma physics for solving crucial issues in reaching burning plasma conditions

Recently, experiments on basic plasma physics issues for solving future problems in fusion energy have been performed on a Large Helical Device. There are several problems to be solved in future devices for fusion energy. Emerging issues in burning plasma are: alpha-channeling (ion heating by alpha particles), turbulence and transport in electron dominant heating helium ash exhaust, reduction of the divertor heat load. To solve these problems, understanding the basic plasma physics of (1) wave–particle interaction through (inverse) Landau damping, (2) characteristics of electron-scale (high-k) turbulence, (3) ion mixing and the isotope effect, and (4) turbulence spreading and detachment, is necessary. This overview discusses the experimental studies on these issues and turbulent transport in multi-ion plasma and other issues in the appendix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Cross-scale interaction between microturbulence and meso-scale reversed shear Alfvén eigenmodes in DIII-D plasmas

Abstract This paper reports global nonlinear gyrokinetic simulations that couple meso-scale reversed shear Alfvén eigenmodes (RSAEs) driven by energetic particles (EPs) and ion temperature gradient (ITG) microturbulence driven by thermal plasma, using equilibrium and profiles from DIII-D discharge #159243. In simulations focusing only on the ITG, electrostatic ITG drives a huge thermal ion heat transport, which is reduced by a factor of 10 to a level close to the experimental value in electromagnetic simulation due to finite β effect. In the simulations coupling the RSAE and ITG, ITG can scatter the resonant EP nonlinearly trapped by the RSAE and damp the zonal flows generated by the RSAE. The regulation of the RSAE by the ITG greatly reduces the initial saturation amplitude of the RSAE but increases the RSAE amplitude and associated EP transport to experimental levels in the quasi-steady state. The RSAE effects on the ITG, specifically the stronger zonal flows generated by the RSAE and the RSAE frequency modulation of the ITG-induced thermal ion heat transport, in turn, leads to a reduction of the thermal ion heat transport by more than a factor of 2 . For a stronger background ITG, the regulation of the RSAE by the ITG is stronger, while the RSAE effects on the ITG are weaker. This work highlights the importance of cross-scale coupling in the dynamics of the AE turbulence and EP transport.

Physics↗
DOE OSTI2024

Validation study of RWM stability in DIII-D high- β N plasmas

The n = 1 (n is the toroidal mode number) resistive wall mode (RWM) stability is numerically investigated for two DIII-D high-β N discharges 176440 and 172461, utilizing the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) codes. Systematic validation efforts are attempted, for the first time, for discharges with very slow or vanishing toroidal flow for a large fraction of the plasma volume. While gaining physics insights in accessing stable operation regime at β N exceeding the Troyon no-wall limit in these slow-rotation experiments, the predictive capability of fluid and non-perturbative magnetohydrodynamic-kinetic hybrid models for the RWM is further confirmed. The MARS-F fluid model, with a strong but numerically tunable viscosity mimicking ion Landau damping of parallel sound waves, finds complete stabilization of the n = 1 RWM in the considered DIII-D plasmas under the experimental flow conditions. Similarly, either full stabilization (for discharge 176440) or marginal stability (for discharge 172461) of the mode is computed by the MARS-K hybrid model, which is first-principle based without free model parameters. In particular, all drift kinetic resonances, including those of thermal and energetic particles, are found to synergistically act to marginally stabilize the RWM in discharge 172461. These MARS-F/K modeling results explain the experimentally observed stable operational regime in DIII-D, as far as the RWM stability is concerned. Extensive numerical sensitivity studies, with respect to the plasma toroidal flow speed as well as the radial location of the resistive wall, are also carried out to further support the validation study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Energetic particle physics: Chapter 7 of the special issue: on the path to tokamak burning plasma operation

We review the physics of energetic particles (EPs) in magnetically confined burning fusion plasmas with focus on advances since the last update of the ITER Physics Basis (Fasoli et al 2007 Nucl. Fusion 47 S264). Topics include basic EP physics, EP generation, diagnostics of EPs and instabilities, the interaction of EPs and thermal plasma instabilities, EP-driven instabilities, energetic particle modes (EPMs), and turbulence, linear and nonlinear stability and simulation of EP-driven instabilities and EPMs, 3D effects, scenario optimization strategies based on EP phase-space control, EPs in reduced field scenarios in ITER before DT, and the physics of runaway electrons. We describe the simulation and modeling of EPs in fusion plasmas, including instability drive and damping as well as EP transport, with a range of approaches from first-principles to reduced models, including gyrokinetic simulations, kinetic-MHD models, gyrofluid models, reduced models, and semi-analytical approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Double-null power-sharing dynamics in MAST-U

Maintaining an effective double-null (DN) configuration is expected to be challenging in reactor-scale tokamak devices. As divertor power-sharing is closely linked to the magnetic topology, even minor variations can lead to fast power-sharing fluctuations which exacerbate the already daunting exhaust challenge. While the static aspects of DN power-sharing have been extensively studied across various devices, this paper presents the first detailed investigation of its dynamic behaviour. We employ dedicated H-mode experiments in MAST-U, in Super-X divertor configuration, featuring perturbation frequencies up to 200 Hz. Our results clearly show no significant dynamic damping of the power-sharing within this frequency range: the divertor responds equally to both fast and slow perturbations. Moreover, the dynamic response also aligns with quasi-static results from slow ramps, implying that static power-sharing models remain valid even for fast fluctuations. Occasionally, some deviations from the otherwise mainly linear behaviour are observed, alongside notable scatter and asymmetries between upwards and downwards trajectories. These observations are likely linked to changes in core conditions, though the underlying mechanisms remain unclear and require further study.

divertor↗
DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: I. Stability

Ion cyclotron emission (ICE) driven by fast ions has the potential to be a reactor relevant diagnostic of confined and escaping fast ions; moreover, controlled experiments in magnetic fusion devices can shed light on similar magnetospheric instabilities below and above the ion cyclotron frequency f ci . Dedicated DIII-D experiments studied instabilities below f ci and at cyclotron harmonics l f ci by injecting thirteen different neutral beam populations into plasmas with different thermal compositions of H, D, and 3He at five different values of toroidal field. High-frequency magnetic loops diagnose mode properties. Each beam type is individually injected for an approximate slowing-down time (~100 ms), resulting in a database of 2529 conditions that are analyzed in frequency bands corresponding to ICE harmonics and < f ci global Alfvén eigenmodes. The sub-cyclotron modes are more unstable at low field and in plasmas with large H concentration. For D injection, first harmonic ICE at f cD is more unstable in hydrogen plasmas but, for H injection, the converse is true (ICE at f cH is more unstable in D), suggesting that thermal same-species ions damp fundamental ICE. ICE stability also depends upon the pitch-angle anisotropy gradient and upon the ‘bump-on-tail’ gradient ∂f/∂v. Although the frequency of unstable ICE in these plasmas always occurs near the central cyclotron frequency, the precise value also depends on the Alfvén speed and beam direction, consistent with the idea that instability occurs for modes that simultaneously satisfy ω ≃ 2πlf ci , the magnetosonic dispersion relation, and the Doppler-shifted resonance condition.

electromagnetic ion cyclotron instabilities↗
DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: III. Novel sub-cyclotron modes

Dedicated DIII-D experiments studied instabilities above and below the ion cyclotron frequency by injecting thirteen different neutral beam populations into L-mode plasmas with different thermal compositions of H, D, and 3 He at five different values of toroidal field. Two companion papers describe the results of these experiments in detail. This paper focuses on data and analysis of two instabilities that occur at ${\sim}0.6$ of the deuterium and hydrogen cyclotron frequencies $f_\mathrm {cD}$ and $f_\mathrm {cH}$ in mixed hydrogen and deuterium plasmas. The modes below $f_\mathrm {cD}$ are often driven unstable by deuterium beams but are rarely destabilized by hydrogen beams. The modes below $f_\mathrm {CD}$ are reported here for the first time in a fusion plasma. These $\mathrm{sub-}f_\mathrm {cH}$ modes are only driven by deuterium beams but, unlike the $\mathrm{sub-}f_\mathrm {cD}$ modes, they are only unstable in plasmas with deuterium concentrations that are ${\lesssim}20\%$; they also have a smaller amplitude than the sub-deuterium modes. Addition of fully-ionized 3 He to the plasma has little effect on the sub-cyclotron deuterium modes but is stabilizing for the sub-cyclotron hydrogen modes. The modes are excited in the core and have toroidal mode numbers around 20. These trends are compared with calculations of the growth rate in a simplified model. The theory successfully predicts global Alfvén eigenmode instability with frequencies and mode numbers close to experimental values, the stronger amplitude of the sub-deuterium modes, and the species-mix dependence of the sub-hydrogen modes. Increased damping associated with the ion-ion hybrid, cutoff, and crossover frequencies plausibly accounts for the stabilization of the $\mathrm{sub-}f_\mathrm {cH}$ modes when 3 He is puffed. However, the relative stability of $\mathrm{sub-}f_\mathrm {cH}$ modes during hydrogen injection remains unexplained.

electromagnetic ion cyclotron instability↗
DOE OSTI2026

Auxiliary heating and current drive physics for the ST-E1 fusion power plant

This work describes the physics basis for the proposed auxiliary heating and current drive system on the ST-E1 fusion power plant. The ST-E1 flattop plasma considered here is fully non-inductive with a bootstrap fraction of 0.9 and the remaining current driven by EC waves. Using the recently published physics-based optimization method for EC launchers (Lopez et al 2025 Plasma Phys. Control. Fusion 67 055012), we show that the target flattop ECCD can be achieved with a net efficiency of 52 kA MW −1 using fundamental O-mode (O1) with frequency range 160–200 GHz launched from the low-field side top half of the vacuum vessel (LFS top-launch). From considering two candidate rampup scenarios, we conclude that LFS top-launch O1 ECCD can be equally effective during the early stages of plasma operation, although poloidal steering might be needed. X-mode waves injected from the LFS midplane are also shown to be effective for rampup even when T e < 1 keV. We also present modeling results for the pre-conceptual design of an ICRH system proposed for ST-E1. Using TORIC, we find that an ICRH system aiming for 42–48 MHz and toroidal mode number n φ ~ 10 robustly achieves dominant ion damping via Helium-3 minority heating transitioning to second-harmonic Tritium heating. We then show that such waves can be efficiently generated by a 5-strap traveling-wave antenna (TWA) using the Petra-M code. The TWA has a 40–45 MHz passband within which ~60% of the power entering the TWA is coupled to the plasma with the remaining ~40% of the power being transmitted through the TWA and possibly recirculated; the power reflected back into the transmission lines is negligible. This passband structure persists even when the evanescent distance is increased by a factor of two, or when the magnetic-field angle is increased by 30°, demonstrating inherent load resilience that will be crucial for effective ICRH on ST-E1.

electron cyclotron↗
DOE OSTI2022

Dynamic inflow and unsteady aerodynamics models for modal and stability analyses in OpenFAST

This work presents two aerodynamic models that are used to perform modal and stability analyses of wind turbines. Including aerodynamic effects in stability analyses is important to capture aerodynamic inertia, stiffness, and damping. The two models consist of a dynamic inflow model based on Øye’s model and an unsteady aerodynamics model based on a simplified Beddoes-Leishman dynamic stall model. Both models are written in a continuous-time state-space form, which is the form needed to perform linearizations within OpenFAST. In this article, we present the new models, perform some simple verifications against HAWC2 results, and then perform stability analyses. We investigate the contribution of the aerodynamic effects on the turbine modes, particularly the flapwise modes, which are the most affected. Because of the development presented in this work, Campbell diagrams accounting for dynamic stall and dynamic wake effects were generated for the first time using OpenFAST.

17 WIND ENERGY↗
DOE OSTI2024

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth-order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A generalized wind turbine cross section as a reduced-order model to gain insights in blade aeroelastic challenges

In this work, we present an approach to study the aeroelastic stability of a wind turbine by focusing on the dynamics of a blade cross section. We present a methodology to obtain a reduced-order model of the blade dynamics in the form of generalized cross-sectional quantities that approximates the aerodynamic and structural properties of the full blade. The motivation for the work is to gain a physical understanding of the influence of aerodynamic models such as dynamic wake and dynamic stall on the frequency and damping of the structure using a reduced-order model with low computational cost. The model may be coupled to two-dimensional computational fluid dynamics softwares or engineering unsteady airfoil aerodynamics models accounting for dynamic wake and dynamic stall. In the latter case, we can obtain monolithic state-space forms of the aeroelastic system of equations, which simplifies the determination of the modal parameters and therefore the study of stability. The work investigates wind turbines in operation or at standstill, where vortex-induced vibrations and stall-induced vibrations, respectively, might be an issue. The implementation is made available as part of the open-source Python package WELIB and as part of the open-source unsteady aerodynamic driver of OpenFAST.

17 WIND ENERGY↗
DOE OSTI2022

Beam dynamics challenges in linear colliders based on laser-plasma accelerators

In this paper we discuss design considerations and beam dynamics challenges associated with laser-driven plasma-based accelerators as applied to multi-TeV-scale linear colliders. Plasma accelerators provide ultra-high gradients and ultra-short bunches, offering the potential for compact linacs and reduced power requirements. We show that stable, efficient acceleration with beam quality preservation is possible in the nonlinear bubble regime of laser-plasma accelerators using beam shaping. Ion motion, naturally occuring for dense beams (i.e., low emittance and high energy) severely damps transverse beam instabilities. Coulomb scattering by the background ions is considered and it is shown that the strong focusing in the plasma strongly suppresses scattering-induced emittance growth. Finally, Betatron radiation emission from the transverse motion of the beam in the plasma will result in beam power loss and energy spread growth; however for sub-100 nm emittances, the beam power loss and energy spread growth will be sub-percent for multi-TeV-class plasma linacs.

47 OTHER INSTRUMENTATION↗