Search NASA⌕ Search

SEARCH · Search NASA

Results for “eigenmode”

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 163 records · Page 9

Development of high-performance long-pulse discharge in KSTAR

High-performance long-pulse plasma operation is essential for producing economically viable fusion energy in tokamak devices. To achieve such discharges in KSTAR, firstly, the rapid increase in the temperature of plasma-facing components was mitigated. The temperature increase of the poloidal limiter, especially, was associated with beam-driven fast ion orbit loss and the discrepancy of the equilibrium reconstructed with heated magnetic probes of signal drift. The fast ions lost to the poloidal limiter were reduced by optimizing the plasma shape and the composition of neutral beam injection (NBI). This nonlinear signal drift was successfully reduced by a new thermal shielding protector on the magnetic probes. Secondly, a lower loop voltage approach was implemented to reduce a poloidal flux consumption rate. A plasma current of 400 kA and a line-averaged electron density of ~2.0 × 10 19 m –3 were chosen by considering the L–H power threshold, fast ion orbit loss, and beam shine-through power loss for low loop voltage in KSTAR. In addition, the application of electron cyclotron heating also helped maintain the plasma with low loop voltage (~25 mV) by enhancing the NBI-driven current and achieving a high poloidal beta (β P ) state. KSTAR has achieved a long pulse (~90 s) operation with the high performance of β P $\leqslant$ 2.7, thermal energy confinement enhancement factor (H 98y2 ) ~ 1.1, and fraction of non-inductive current (f NI ) $\leqslant$ 0.96. Still, gradual degradation of the plasma performance has been observed over time in the discharges. In one of the long-pulse discharges, β P reduced by ~18% over the time of ~8τ R (current relaxation time, τ R ~5 s) and ~1067τ E,th (thermal energy confinement time, τ E,th ~45 ms). The degradation may be closely associated with weak, yet growing, and persistent toroidal Alfvén eigenmodes and their effect on fast ion confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transport and losses of energetic particles in tokamaks in the presence of Alfvén activity using the new full orbit TAPaS code coupled to FAR3d

Recent developments and tools integrated into the TAPaS code are presented, enabling realistic scenario simulations of particle dynamics within experimental tokamak magnetic equilibria. In particular, the enhanced capabilities of TAPaS enable seamless coupling with external simulations, provided the metric and equilibrium magnetic field of the external code are known. Coupling TAPaS with the gyro-fluid code FAR3d, the transport and losses of energetic particles in the presence Alfvén eigenmodes (AEs) in DIII-D plasma discharge #159243 were investigated. Detailed analyses of prompt losses with and without collisions were performed. Then, further analysis was performed in the presence of electromagnetic perturbations resulting from AEs activity. The results indicate that, for the energies and the initial conditions considered here, the presence of AEs enhances the particle losses.

Alfven eigenmodes↗

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↗

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↗

Shape fluctuations and optical transition of He 2 * excimer tracers in superfluid He 4

Metastable He$_{2}^{*}$ excimer molecules have been utilized as tracer particles of the normal component in superfluid $^4$He (He II) which can be imaged via laser-induced fluorescence. These excimer molecules form tiny bubbles in He II and can bind to quantized vortices at sufficiently low temperatures, thereby allowing for direct visualization of vortex dynamics in an inviscid superfluid. However, the $a^{3}\Sigma^+_{u}$${\rightarrow}$$c^{3}\Sigma^+_{g}$ optical absorption line, which is responsible for the fluorescence imaging of the He$_{2}^{*}$ molecules, is controlled by fluctuations on the bubble shape, and its exact line profile is not known at low temperatures. In this paper, we present a bubble model for evaluating the surface fluctuation eigenmodes of the excimers in He II. The line profile of the $a^{3}\Sigma^+_{u}{\rightarrow}c^{3}\Sigma^+_{g}$ transition is calculated at different temperatures by considering both the zero-point and thermal fluctuations on the bubble shape. We show that, as the temperature drops from 2~K to 20 mK, the peak absorption strength is enhanced by a factor of about five, accompanying a blueshift of the peak location by about 2 nm. A double-peak line profile due to the rotational levels of the molecular core can be resolved. This bubble model also allows us to evaluate the stiffness of the He$_{2}^{*}$ bubbles and hence their diffusion constant in He II due to scattering off thermal phonons. Our results will aid the design of future experiments on imaging quantized vortices in He II using He$_{2}^{*}$ tracers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Energetic particle transport in 3D equilibrium. Final report

In this project, we have further developed the state-of-the-art fusion simulation code, gyrokinetic toroidal code (GTC), for studying the energetic particle (EP) transport in the 3D equilibrium. We have applied the new capabilities to simulate the EP transport by the Alfven eigenmodes in stellarators and by 3D magnetic islands in tokamaks. As part of V&V efforts for 3D capability, the 3D capabilities of GTC have been applied to study thermal transport in stellarators and tokamak with 3D magnetic fields by microturbulence, which could cause additional EP loss. Effects magnetic islands on neoclassical bootstrap current have also been studied. The computational capability developed and physics learned in this project are in support of fusion SciDAC Center for Integrated Simulation of Energetic Particles in Burning Plasmas (ISEP). During this project, the PI has supervised 6 postdoctoral researchers and 4 PhD students at the University of California, Irvine. A complete list of 49 publications and 23 invited talks partially supported by this grant since 7/2015 appears at the end of this report. Five key publications closely related to this project are highlighted in this report.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stability optimization of energetic particle driven modes in nuclear fusion devices: the FAR3d gyro-fluid code

The development of reduced models provide efficient methods that can be used to perform short term experimental data analysis or narrow down the parametric range of more sophisticated numerical approaches. Reduced models are derived by simplifying the physics description with the goal of retaining only the essential ingredients required to reproduce the phenomena under study. This is the role of the gyro-fluid code FAR3d, dedicated to analyze the linear and nonlinear stability of Alfvén Eigenmodes (AE), Energetic Particle Modes (EPM) and magnetic-hydrodynamic modes as pressure gradient driven mode (PGDM) and current driven modes (CDM) in nuclear fusion devices. Such analysis is valuable for improving the plasma heating efficiency and confinement; this can enhance the overall device performance. The present review is dedicated to a description of the most important contributions of the FAR3d code in the field of energetic particles (EP) and AE/EPM stability. FAR3d is used to model and characterize the AE/EPM activity measured in fusion devices as LHD, JET, DIII-D, EAST, TJ-II and Heliotron J. In addition, the computational efficiency of FAR3d facilitates performing massive parametric studies leading to the identification of optimization trends with respect to the AE/EPM stability. This can aid in identifying operational regimes where AE/EPM activity is avoided or minimized. This technique is applied to the analysis of optimized configurations with respect to the thermal plasma parameters, magnetic field configuration, external actuators and the effect of multiple EP populations. In addition, the AE/EPM saturation phase is analyzed, taking into account both steady-state phases and bursting activity observed in LHD and DIII-D devices. The nonlinear calculations provide: the induced EP transport, the generation of zonal structures as well as the energy transfer towards the thermal plasma and between different toroidal/helical families. Finally, FAR3d is used to forecast the AE/EPM stability in operational scenarios of future devices as ITER, CFETR, JT60SA and CFQS as well as possible approaches to optimization with respect to variations in the most important plasma parameters.

Alfv én Eigenmodes↗

Tunable Phonon Polariton Hybridization in a Van der Waals Hetero‐Bicrystal

Abstract Phonon polaritons, the hybrid quasiparticles resulting from the coupling of photons and lattice vibrations, have gained significant attention in the field of layered van der Waals heterostructures. Particular interest has been paid to hetero‐bicrystals composed of molybdenum oxide (MoO 3 ) and hexagonal boron nitride (hBN), which feature polariton dispersion tailorable via avoided polariton mode crossings. In this work, the polariton eigenmodes in MoO 3 ‐hBN hetero‐bicrystals self‐assembled on ultrasmooth gold are systematically studied using synchrotron infrared nanospectroscopy. It is experimentally demonstrated that the spectral gap in bicrystal dispersion and corresponding regimes of negative refraction can be tuned by material layer thickness, and these results are quantitatively matched with a simple analytic model. Polaritonic cavity modes and polariton propagation along “forbidden” directions are also investigated in microscale bicrystals, which arise from the finite in‐plane dimension of the synthesized MoO 3 micro‐ribbons. The findings shed light on the unique dispersion properties of polaritons in van der Waals heterostructures and pave the way for applications leveraging deeply sub‐wavelength mid‐infrared light‐matter interactions.

36 MATERIALS SCIENCE↗

Analysis of the Spontaneous Emission Limited Linewidth of an Integrated III–V/SiN Laser

Abstract This article describes a calculation of the spontaneous emission limited linewidth of a semiconductor laser consisting of hybrid or heterogeneously integrated, silicon and III–V intracavity components. Central to the approach are a) description of the multi‐element laser cavity in terms of composite laser/free‐space eigenmodes, b) use of multimode laser theory to treat mode competition and multiwave mixing, and c) incorporation of quantum‐optical contributions to account for spontaneous emission effects. Application of the model is illustrated for the case of linewidth narrowing in an InAs quantum‐dot laser coupled to a high‐ SiN cavity.

42 ENGINEERING↗

Tailored Frequency Conversion in Nonlinear Subwavelength Grating Metaslabs

The use of metaslab waveguides consisting of subwavelength gratings (SWGs) is proposed as a highly flexible approach to obtain phase‐matched frequency conversion in nonlinear integrated devices. Control over the phase matching (PM) condition is achieved through modification of the linear modal propagation constants by the SWG. The method is experimentally demonstrated by fabrication of 0.4 mm‐long GaAs metaslabs designed for second harmonic generation of near‐infrared light. A controlled spectral shift of up to 100 nm in the PM condition, relative to a uniform slab is observed. Excellent agreement is found between the experimental results and finite‐difference eigenmode simulations, which predict a tunable frequency conversion bandwidth of ∼300 nm for the design parameters. The different regimes of PM control supported by the structure are discussed. This work highlights the potential of SWGs as building blocks for tailored frequency conversion in nonlinear integrated devices.

Form birefringence↗

Looking at extremal black holes from very far away

Near-extremal black holes are subject to large quantum effects, which modify their low-temperature thermodynamic behavior. Hitherto, these quantum effects were analyzed by separating the geometry into the near-horizon region and its exterior. It is desirable to understand and reproduce such corrections from the full higher-dimensional asymptotically flat or AdS geometry’s perspective. We address this question in this article and fill this gap. Specifically, we find off-shell eigenmodes of the quadratic fluctuation operator of the Euclidean gravitational dynamics, with eigenvalues that vanish linearly with temperature. We illustrate this for BTZ and neutral black holes with hyperbolic horizons in AdS in Einstein-Hilbert theory, and for the charged black holes in Einstein-Maxwell theory. The linear scaling with Matsubara frequency, which is a distinctive feature of the modes, together with the fact that their wavefunctions localize close to the horizon as we approach extremality, identifies them as responsible for the aforementioned quantum effects. We provide a contour prescription to deal with the sign indefiniteness of the Euclidean Einstein-Maxwell action, which we derive to aid our analysis. We also resolve a technical puzzle regarding modes associated with rotational isometries in stationary black hole spacetimes.

AdS-CFT Correspondence↗

Boson star normal modes

Boson stars are gravitationally bound objects that arise in ultralight dark matter models and form in the centers of galactic halos or axion miniclusters. We systematically study the excitations of a boson star, taking into account the mixing between positive and negative frequencies introduced by gravity. We show that the spectrum contains zero-energy modes in the monopole and dipole sectors resulting from spontaneous symmetry breaking by the boson star background. We analyze the general properties of the eigenmodes and derive their orthogonality and completeness conditions which have non-standard form due to the positive-negative frequency mixing. The eigenvalue problem is solved numerically for the first few energy levels in different multipole sectors and the results are compared to the solutions of the Schrödinger equation in fixed boson star gravitational potential. The two solutions differ significantly for the lowest modes, but get close for higher levels. We further confirm the normal mode spectrum in 3D wave simulations where we inject perturbations with different multipoles. As an application of the normal mode solutions, we compute the matrix element entering the evaporation rate of a boson star immersed in a hot axion gas. The computation combines the use of exact wavefunctions for the low-lying bound states and of the Schrödinger approximation for the high-energy excitations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Eigenfrequency constrained topology optimization of finite strain hyperelastic structures

This paper incorporates hyperelastic materials, nonlinear kinematics, and preloads in eigenfrequency constrained density–based topology optimization. The formulation allows for initial finite deformations and subsequent small harmonic oscillations. The optimization problem is solved by the method of moving asymptotes, and the gradients are calculated using the adjoint method. Both simple and degenerate eigenfrequencies are considered in the sensitivity analysis. A well-posed topology optimization problem is formulated by filtering the volume fraction field. Numerical issues associated with excessive distortion and spurious eigenmodes in void regions are reduced by removing low volume fraction elements. The optimization objective is to maximize stiffness subject to a lower bound on the fundamental eigenfrequency. Numerical examples show that the eigenfrequencies drastically change with the load magnitude, and that the optimization is able to produce designs with the desired fundamental eigenfrequency.

42 ENGINEERING↗

Explicit dispersion relations for warm fluid waves in a uniform plasma (invited)

Classical dispersion relations for waves in a fluid plasma are expressed as implicit functions of wave frequency and wave- number. The implicit dispersion relation must in general be solved numerically. This work introduces an astute method for obtaining an explicit dispersion relation for waves in a fluid plasma. The explicit expression has an advantage of providing eigenmodes without numerical computations and giving a dispersion relation more detailed than the implicit expression. As in the case of cold waves, the wavenumber for an arbitrary frequency can be directly obtained from the explicit formula. The explicit formula also enables an accurate evaluation of finite-temperature effects on a dispersion relation even near the cyclotron resonance, where a numerical analysis of the implicit relation is impractical. As a result, the analytic formula can be used to investigate temperature effects on electromagnetic waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗

Nonlinear dynamics of combustor azimuthal modes: Experiments and modeling

This paper provides an experimental and theoretical investigation of azimuthal instabilities in a lean, premixed combustor. It presents data from a high Reynolds number, industrial-scale facility with multiple fuel stages that enable variation of radial and azimuthal symmetries. The combustor exhibits a variety of behaviors – including standing, spinning, mixed waves, and intermittent behaviors, as well as periodic limit cycles and quasi-periodic oscillations – depending upon varying power and azimuthal symmetry. A second key objective of this paper is to present a weakly nonlinear framework that is used to understand and interpret these results, showing the influence of eigenmode degeneracy, background noise, and azimuthal non-uniformities in mean quantities such as temperature and flame-acoustic coupling. In particular, this analysis shows the importance of the linear instability growth/damping rate and azimuthal non-uniformities in controlling the relative prominence of standing, spinning, or mixed waves at the limit cycle, as well as the preference in the spinning direction. It also shows the conditions under which single frequency or quasi-periodic oscillations are dominant, and how these are more fundamentally influenced by azimuthal non-uniformities and background noise.

42 ENGINEERING↗

Re-splitting δ f method for electro-magnetic gyrokinetic particle-in-cell (PIC) simulation of tokamak plasmas

A new method based on a re-splitting δ f technique is proposed to mitigate the cancellation problem in electro-magnetic gyrokinetic PIC simulation. This new scheme, together with the pull-back mitigation (PBM) scheme (Mishchenko et al., 2014), have recently been implemented in the GEM code (Chen et al., 2003). These two algorithms are systematically compared with the original GEM algorithm for the simulation of the ion temperature gradient (ITG) mode, the kinetic ballooning mode (KBM), the toroidal Alfvén eigenmode (TAE) and the energetic particle mode (EPM). The three algorithms agree well for all of these modes. Moreover, both the re-splitting method and PBM can mitigate the cancellation problem and improve the computational efficiency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗