Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “toroidal magnetic fields”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Solar Cycle 24 and the Solar Dynamo

We will discuss the polar field precursor method for solar activity prediction, which predicts cycle 24 will be significantly lower than recent activity cycles, and some new ideas rejuvenating Babcock's shallow surface dynamo. The polar field precursor method is based on Babcock and Leighton's dynamo models wherein the polar field at solar minimum plays a major role in generating the next cycle's toroidal field and sunspots. Thus, by examining the polar fields of the Sun near solar minimum, a forecast for the next cycle's activity is obtained. With the current low value for the Sun's polar fields, this method predicts solar cycle 24 will be one of the lowest in recent times, with smoothed F10.7 radio flux values peaking near 135 plus or minus 35 (2 sigma), in the 2012-2013 timeframe (equivalent to smoothed Rz near 80 plus or minus 35 [2 sigma]). One may have to consider solar activity as far back as the early 20th century to find a cycle of comparable magnitude. We discuss unusual behavior in the Sun's polar fields that support this prediction. Normally, the solar precursor method is consistent with the geomagnetic precursor method, wherein geomagnetic variations are thought to be a good measure of the Sun's polar field strength. Because of the unusual polar field, the Earth does not appear to be currently bathed in the Sun's extended polar field (the interplanetary field), hence negating the primal cause behind the geomagnetic precursor technique. We also discuss how percolation may support Babcock's original shallow solar dynamo. In this process ephemeral regions from the solar magnetic carpet, guided by shallow surface fields, may collect to form pores and sunspots.

Pesnell, W. D.↗
DOE OSTI2024

Ripple-induced neoclassical toroidal viscous torque in Augmented-First Plasma operation phase in ITER

Abstract A systematic calculation is performed on the ripple-induced neoclassical toroidal viscous (NTV) torque for new ITER scenarios designed for the Augmented-First Plasma (A-FP) operation phase with the full tungsten wall, where the plasma-wall gap is varied in view of mitigating the impact of tungsten wall-plasma interactions. The torque calculation includes drift kinetic response of the plasma thermal and energetic particles to the n = 18 ( n is the toroidal harmonic number) ripple field. For the plasma scenario with ~45 cm plasma-wall gap at the outboard mid-plane and considering the corrected ripple level of 0.17% by the ferritic steel inserts, the computed net NTV torque acting on the plasma column is in the sub-Nm level. However, with decreasing the plasma-wall gap, the computed net NTV torque can reach a level comparable to that produced by the neutral-beam momentum injection in ITER. Ripple correction by ferritic inserts reduces the net torque by a factor of 3.3 for all the three A-FP scenarios considered. The n ω d = l ω b (with ω d and ω b being the toroidal precession and bounce frequencies of trapped particles, respectively, and l an integer number) type of resonance-enhancement of the NTV torque, due to thermal particles, is found to be weak in ITER despite high- n of 18. The same also holds for the ITER 10 MA steady state scenario from the D-T operation phase, where the aforementioned resonance associated with fusion-born alphas is also included. The ripple-induced NTV torque is well below that produced by the resonant magnetic perturbation applied for controlling the type-I edge-localized mode in ITER.

ITER↗
DOE OSTI2024

Characterizing impurity sourcing and transport in the high temperature boundary of DIII-D wide pedestal QH-mode plasmas

Wide Pedestal QH-mode (WPQH) plasmas in the DIII-D tokamak show a sheath limited SOL, where electron density and temperature remain nearly constant along field lines from the midplane to the divertor. Consequently, parallel gradients are weak. The first Langmuir probe measurements in this regime point to a high sheath temperature with target electron temperature T e up to ~150 eV, leading to high carbon self-sputtering. Midplane carbon densities from SOLPS-ITER modeling of these plasmas in a double-null configuration fall far below experimental measurements unless full drifts are activated. In the drift-dependent SOLPS-ITER modeling, anomalous poloidally uniform, radially varying transport is adjusted to match measured radial electron temperature and density profiles in the pedestal and SOL region. The resulting carbon density (C6+) matches measured carbon densities just inside the separatrix. The C2+ density near the outer strike point is also consistent with spectroscopic imaging. With the ion B×∇B drift towards the X-point, carbon in the lower divertors is redistributed from the private flux region to the high field side (HFS) and pushed upstream in the SOL by poloidal E×B drifts. The B × ∇B drift dominates the radial flow of carbon as it moves upstream. Simulations with reversed toroidal field (ion B × ∇B drift away from the X-point) show a radically different behavior, where carbon accumulates on the low field side in lower divertors, and is pushed towards the HFS in the upper divertors, indicating a strong effect of particle drifts on impurity distribution. In double-null configurations which are usually used in WPQH plasmas, these drift effects tend to counterbalance, so that the carbon density and Z eff are reduced at the outer mid-plane only modestly, by around 7% in the modeling. Further predictive modeling indicates that the carbon density can be significantly reduced by an order of magnitude by implementing a single-null shape with reversed toroidal field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

On reflection of Alfven waves in the solar wind

We have revisited the problem of propagation of toroidal and linear Alfven waves formulated by Heinemann and Olbert (1980) to compare WKB and non-WKB waves and their effects on the solar wind. They considered two solar wind models and showed that reflection is important for Alfven waves with periods of the order of one day and longer, and that non-WKB Alfven waves are no more effective in accelerating the solar wind than WKB waves. There are several recently published papers which seem to indicate that Alfven waves with periods of the order of several minutes should be treated as non-WKB waves and that these non-WKB waves exert a stronger acceleration force than WKB waves. The purpose of this paper is to study the origin of these discrepancies by performing parametric studies of the behavior of the waves under a variety of different conditions. In addition, we want to investigate two problems that have not been addressed by Heinemann and Olbert, namely, calculate the efficiency of Alfven wave reflection by using the reflection coefficient and identify the region of strongest wave reflection in different wind models. To achieve these goals, we investigated the influence of temperature, electron density distribution, wind velocity and magnetic field strength on the waves. The obtained results clearly demonstrate that Alfven wave reflection is strongly model dependent and that the strongest reflection can be expected in models with the base temperatures higher than 10(exp 6) K and with the base densities lower than 7 x 10(exp 7) cm(exp -3). In these models as well as in the models with lower temperatures and higher densities, Alfven waves with periods as short as several minutes have negligible reflection so that they can be treated as WKB waves; however, for Alfven waves with periods of the order of one hour or longer reflection is significant, requiring a non-WKB treatment. We also show that non-WKB, linear Alfven waves are always less effective in accelerating the plasma than WKB Alfven waves. Finally, it is evident from our results that the region of strongest wave reflection is usually located at the base of the models, and hence that interpretation of wave reflection based solely on the reflection coefficient can be misleading.

Krogulec, M.↗
NASA NTRS1994

Reflection of Alfven waves in the solar wind

We have revisited the problem of propagation of toroidal and linear Alfven waves formulated by Heinemann and Olbert (1980) to compare Wentzel-Kramers-Brillouin (WKB) and non-WKB waves and their effects on the solar wind. They considered two solar wind models and showed that reflection is important for Alfven waves with periods of the order of one day and longer and that non-WKB Alfven waves are no more effective in accelerating the solar wind than in WKB waves. There are several recently published papers that seem to indicate that Alfven waves with periods of the order of several minutes should be treated as non-WKB waves and that these non-WKB waves exert a stronger acceleration force than WKB waves. The purposse of this paper is to study the origin of these discrepancies by performing parametric studies of the behavior of the waves under a variety of different conditions. In addition, we want to investigate two problems that have not been addressed by Heinimann and Olbert, namely, calculate the efficieny of Alfven wave reflection by using the reflection coefficient and identfy the region of strongest wave reflection in different wind models. To achieve these goals, we investigate the influence of temperature, electron desity distribution, wind velocity, and magnetic field strength on te waves. The obtained results clearly demonstrate that Alfven wave reflection is strongly model dependent and that the strongest reflection can be expected in models with the base temperatures higher than 10(exp 6) K and with the base densities lower than 7 x 10(exp 7)/cu cm. In these models as well as in the models with lower temperatures and higher densities Alfven waves with periods as short as several minutes have negligible reflection so that they can be treated as WKB waves; however, for Alfven waves with periods of the order of one hour or longer reflection is significant, requiring a non-WKB treatment. We also show that non-WKB, linear Alfven waves are always less effective in accelerating the plasma than WKB Alfven waves. Finally, it is evident from our results that the region of strongest wave reflection is usually located at the base of the models and hence that interpretation of wave reflection based soley on the reflection coefficient can be misleading.

Krogulec, M.↗
DOE OSTI2021

Microtearing modes as the source of magnetic fluctuations in the JET pedestal

Abstract We report on a detailed study of magnetic fluctuations in the JET pedestal, employing basic theoretical considerations, gyrokinetic simulations, and experimental fluctuation data to establish the physical basis for their origin, role, and distinctive characteristics. We demonstrate quantitative agreement between gyrokinetic simulations of microtearing modes (MTMs) and two magnetic frequency bands with corresponding toroidal mode numbers n = 4 and 8. Such disparate fluctuation scales, with substantial gaps between toroidal mode numbers, are commonly observed in pedestal fluctuations. Here we provide a clear explanation, namely the alignment of the relevant rational surfaces (and not others) with the peak in the ω * profile, which is localized in the steep gradient region of the pedestal. We demonstrate that a global treatment is required to capture this effect. Nonlinear simulations suggest that the MTM fluctuations produce experimentally-relevant transport levels and saturate by relaxing the background electron temperature gradient, slightly downshifting the fluctuation frequencies from the linear predictions. Scans in collisionality are compared with a simple MTM dispersion relation. At the experimental points considered, MTM growth rates can either increase or decrease with collision frequency depending on the parameters thus defying any simple characterization of collisionality dependence.

Physics↗
DOE OSTI2022

Multi-harmonic electron cyclotron heating and current drive scenarios for non-inductive start-up and ramp-up in high field ST-40 spherical tokamak

We report non-inductive start-up and ramp-up is an important topic for spherical tokamak reactor design as the central solenoid implementation is highly restrictive particularly for the low-aspect-ratio tokamak configuration. In the high field spherical tokamak (ST), ST-40 with B T0 ≤ 3 T, a preparation is underway for high power ECH and ECCD current start-up/ramp-up experiments utilizing two MW-class 140/105 GHz gyrotrons. Here, we explored various ECH/ECCD scenarios for a low-field-side (LFS) launch-angle steerable waveguide launcher placed near the mid-plane region. Due to the large toroidal field variation of ST configuration, multiple cyclotron harmonic resonance layers could exist within the plasma. In this start-up and ramp-up regime, both fundamental and second harmonic ECH resonances must be considered. We find that even with the presence of X-II resonance layer in the plasma, an efficient X-I ECH and ECCD regime can be accessed for the low electron temperature T e0 as low as 200 eV which is a typical starting temperature of ECH heated plasmas in an open-field-line configuration. The presence of X-II resonance could become significant at higher T e0 as X-II absorption increases with T e0 which could reduce the current ramp-up efficiency as the power reaching X-I is reduced. Finally for the pure X-I regime where the 2Ω e resonance is moved outside the plasma with B T0 ~ 3.4 T, we find that it is possible to reach the full current of I p ~ 1 MA fully non-inductively with the ECH power of ~1 MW at n e0 ~ 1.0 × 10 19 m -3 using 105 GHz frequency gyrotron. By reducing the outer limiter position R L ~ 78 cm to 70 cm, the pure X-I regime is recovered at the rated ST-40 magnetic field of B T0 ~ 3.0 T. This X-I regime is accessible with a relatively broad range of launched n II or the launching angles. A survey of X-mode X-II ECH and ECCD at higher density regimes is also shown for completeness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulations of plasma current induced by toroidal field ramping down on tokamaks

In tokamak experiments, the BT ramping down can drive off-axis and parallel inductive current. This approach leads to a decrease in q95 and an increase in normalized beta, βN. The off-axis inductive current also broadens the current profile. Typically, the 1.5D transport code is used to simulate the time evolution of the plasma current profile, which is based on the flux-surface-averaged Faraday’s law. The ONETWO code is one of those transport codes. However, this code cannot simulate the situation of the evolving toroidal field BT. This article proposes a new time-dependent model to take into account the BT ramping down situation. A modified formula of flux-surface-averaged Faraday’s law was derived to consider the effect of the BT ramping rate on the current evolution, and it was implemented in the ONETWO code. Then, the modified ONETWO code was used to simulate the BT ramping down experiment on DIII-D. The simulation result of the plasma current evolution with BT ramping down shows a broader current profile with smaller ohmic current induced by the poloidal field, compared to that without BT ramping down.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Asymptotic solutions and convergence studies of the resistive inner region equations

A recent publication presents a detailed prescription for constructing global eigenfunctions and complex growth rates for resistive instabilities in axisymmetric toroidal systems, using the method of matched asymptotic expansions. At each singular surface in the plasma, the small-x asymptotic solutions of the ideal MHD outer-region equations are matched to the large-x asymptotic solutions of the resistive inner-region equations, with x being the distance from the singular surface. The success of this method depends upon accurate asymptotic solutions of the resistive inner region equations. Extensive studies have shown that there are not only regimes of excellent agreement with straight-through methods, without asymptotic matching, but also regimes of strong disagreement. As such, in an effort to find the cause of this discrepancy, we present here a different method of constructing asymptotic solutions of the inner region equations and we study the convergence properties of the new and old solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Edge-localized-mode simulation in CFETR steady-state scenario

We report the EPED1 model and self-consistent core-pedestal coupling in integrated modeling are used to design the pedestal structure of the China Fusion Engineering Testing Reactor (CFETR) steady-state scenario. The key parameters, such as β p and q 95 , are based on the grassy edge-localized-mode (ELM) experimental database. In this work, we use the BOUT++ six-field two-fluid code to simulate the onset of the ELM in the CFETR steady-state scenario. The ELM size is around 0.2% in nonlinear simulations, which is in the experimental range of the grassy ELM discharges, 0.1%–1% observed in multiple tokamak devices. Linear and nonlinear simulations show that the dominant high-n ballooning modes peak around n = 40. Compared to type-I ELM crashing dynamics, grassy ELM crashing has a smaller initial crash and is then followed by three phases of turbulence spreading, which are dominated by multi-modes, a high-n mode of n = 45 and low-n mode of n = 5, respectively. In contras to type-I ELM, the perturbation of the high-n mode has a narrow width around ψ = 0.95, and magnetic island formation and reconnection occur only beyond ψ = 0.95, leading to a small initial crash. Mode–mode interaction in the multi-mode coexistence stage stops the growth of individual modes and reduces the transport of particles and heat, and these are the two reasons why the ELM size is small. In–out asymmetry of transient heat flux with a ratio of E out /E in = 3.5 is found during grassy ELM crash. The rise and delay times of the heat flux match the calculation from the free-streaming model. To evaluate the erosion of the divertor target, the energy fluence at the outer divertor target is calculated, which is 0.029 MJ m -2 , 5.5 times smaller than the tungsten melting limit 0.16 MJ m -2 . The calculated energy fluency still follows the experimental scaling law from type-I ELM experiments. The fluctuation eddies in the toroidal direction show a filament structure at the outer mid-plane. Parallel heat flux patterns with a toroidal mode number n = 10 are found at the outer divertor with an amplitude of 680 MW m -2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modelling the Alfvén eigenmode induced fast-ion flow measured by an imaging neutral particle analyzer

An imaging neutral particle analyzer (INPA) provides energy and radially resolved measurements of the confined fast-ion population ranging from the high-field side to the edge on the midplane of the DIII-D tokamak. In recent experiments, it was used to diagnose fast-ion flow in the INPA-interrogated phase-space driven by multiple, marginally unstable Alfvén eigenmodes (AEs). Here, the key features of this measured fast-ion flow are: (I) a fast-ion flow from q min and the injection energy (81 keV) towards lower energies and plasma periphery.(II) A flow from the same location towards higher energies and the plasma core, (III) a phase-space 'hole' at the injected energy and plasma core and (IV) a pile-up at the plasma core at lower energies (~60 keV). Ad hoc energetic particle diffusivity modelling of TRANSP significantly deviates from the observation. Comparably, a reduced modelling, i.e. a combination of NOVA-K and ASCOT5 code with the measured mode structure and amplitude, generally reproduce some key features of the observed phase-space flow, but largely failed to interpret fast ion depletion near the plasma axis. At last, self-consistent, first-principle multi-phase hybrid simulations that include realistic neutral beam injection and collisions are able to reproduce most features of the time-resolved phase-space flow. During consecutive hybrid phases, an RSAE consistent with the experiment grows and saturates, redistributing the injected fast ions. The resulting synthetic INPA images are in good agreement with the measurement near the injection energy. The simulations track the fast-ion redistribution within the INPA range, confirming that the measured fast-ion flow follows streamlines defined by the intersection of phase-space surfaces of constant magnetic moment μ and constant E' = nE + ωP φ , where n and ω are the instability toroidal mode number and frequency, and E and P φ the ion energy and toroidal canonical momentum. Nonperturbative effects are required to reproduce the depletion of fast ions near the magnetic axis at the injection energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

First Demonstration of Resonant Pitch-Angle Scattering of Relativistic Electrons by Externally Launched Helicon Waves

Helicon waves (a.k.a whistler waves) satisfying the normal wave-particle cyclotron resonance are observed to limit the growth and maximum energy of relativistic electrons (REs) in low-density Ohmic DIII-D tokamak plasmas. Following the application of helicon waves, pitch-angle scattering of high-energy REs causes an increase in both synchrotron and electron-cyclotron emissions. The hard x-ray emission, a proxy for the RE population, ceases to grow. Energy-resolved hard x-ray measurements also show a striking decrease in the number of high-energy REs (above the resonance at approximately 8MeV) to below the noise floor and an increase in low-energy (∼ 4 MeV) REs. This occurs despite the toroidal electric field remaining high enough to drive exponential RE growth in the absence of helicon waves. Furthermore, these results open new directions for limiting the maximum energy of RE populations in laboratory and fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Alfvén mode spectrum in the advanced steady-state tokamak scenario with negative triangularity

In the earlier studies [Zheng et al. Fundam. Plasma Phys. 10, 100051 (2024)], the negative triangularity tokamak was shown to have a unique and outstanding feature for the steady-state confinement. In addition to relieving the divertor heat load concern, it can reach a much higher beta limit in the nearly steady-state confinement, about twice the Troyon limit in the positive triangularity case. In this work, we extend the investigation to the Alfvén spectrum and toroidal Alfvén eigenmodes (TAEs). This helps to evaluate the property of energetic particle confinement in this advanced scenario of a negative triangularity tokamak. We found that, because of the large safety factor and reversed magnetic shear in the core plasma region near the magnetic axis, the TAE gap tends to be closed or have only a short radial extension between the continua. in the on-axis core region and open primarily in the off-axis region, where TAEs are shown to be present. This indicates that energetic particle confinement in the core region near the magnetic axis, where nuclear fusion reactions prevail, is not seriously affected by the excitation of TAEs in this scenario. The configuration with off-axis TAEs helps the confinement of energetic particles as compared to the configuration with near-axis TAEs. The results show that the advanced steady-state tokamak scenario with negative triangularity, because of the rapid variation of the safety factor, also has a favorable feature for energetic particle confinement.

Fourier analysis↗
DOE OSTI2021

Further modeling of q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

An improved resonant plasma response model that more accurately captures the physics of the interaction between a tokamak plasma and a resonant magnetic perturbation (RMP) is developed. The model interpolates between the linear and the nonlinear response regimes and takes into account the fact that the slip-frequency is non-zero in the nonlinear regime. The improved model is incorporated into the extended perturbed equilibrium code (EPEC) toroidal asymptotic matching code. The modified EPEC code is used to investigate RMP-induced edge-localized-mode (ELM) suppression in DIII-D H-mode discharge #145380. Somewhat surprisingly, allowing for a finite slip-frequency (i.e., relaxing the so-called no-slip constraint) is found to only slightly facilitate the locking of driven magnetic island chains to the RMP, and, hence, to only slightly facilitate RMP-induced ELM suppression. This is true despite the fact that the nature of non-locked island solutions is radically different when the no-slip constraint is imposed compared to when it is relaxed (in the first case, the widths of the island chains driven at the rational surfaces pulsate, and in the second case, they remain steady). The previously obtained conclusion that the response of a typical H-mode tokamak plasma to an RMP cannot be accurately modeled by linear theory is confirmed. Here, the previously obtained conclusion that the best agreement between theory and observations is achieved by assuming that the natural frequencies of tearing modes, in the absence of the RMP, are determined by the local equilibrium E × B velocity is also confirmed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of Magnetic Field Configuration on Heat Transport in Stellarators and Heliotrons

We assess the magnetic field configuration in modern fusion devices by comparing experiments with the same heating power, between a stellarator and a heliotron. The key role of turbulence is evident in the optimized stellarator, while neoclassical processes largely determine the transport in the heliotron device. Gyrokinetic simulations elucidate the underlying mechanisms promoting stronger ion scale turbulence in the stellarator. Similar plasma performances in these experiments suggests that neoclassical and turbulent transport should both be optimized in next step reactor designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Initial Results from High-Field-Side Transient CHI Start-Up on QUEST

Transient coaxial helicity injection (t-CHI) current start-up using a new design simple electrode configuration has been implemented on the QUEST. Discharges injected from the low field side (LFS) and from the high field side (HFS) were examined. Compared to the LFS injection case, the HFS injection has the advantages of providing access to a higher toroidal field and better controlling the location of the injector flux footprint location. Although the present PF coils on QUEST are not well positioned to form the injector flux on the HFS injector region and there has been a frequent occurrence of the spurious arcs, known as absorber arcs, HFS injection has shown flux evolution in a shape that is suitable for the formation of closed flux surfaces. The discharges were improved by installing an in-vessel-coil and adding a new cylindrical electrode to the existing CHI electrode. Overall, the results show that the new cylindrical electrode allowed the flux to evolve stably while allowing both the inner and the outer injector flux footprint to remain in the vicinity of the cylindrical electrode. This configuration which inherently generates a narrow injector flux footprint width resulted in discharges that strongly suggested the persistence of the CHI generated plasma after the injector current was reduced to zero. These studies have informed us of the need to improve the CHI gas injection system so that the absorber arcs could be better controlled in the HFS injection configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Interplay between beam-driven chirping modes and plasma confinement transitions in spherical tokamak ST40

Experiments on the high field spherical tokamak ST40 have led to the recent observation of interplay between beam-driven modes of sweeping frequency (chirping modes) and transitions to the enhanced global confinement regime (H-mode) and back to the low confinement regime (L-mode). The H-modes of plasma confinement are identified from decreased intensity of D α signal and from clear distinctions in the edge gradients of the visible plasma boundary (observed as a sharp plasma edge in camera images). The beam-driven chirping modes, identified as ideal magnetohydrodynamics beta-induced Alfvén acoustic eigenmodes modes, are observed in Mirnov coil signals, interferometry, and soft x-ray diagnostics. A moderate amplitude ‘primer’ chirping mode usually precedes an H–L transition. This is followed by a ‘dominant’ chirping mode with higher amplitude during the L-mode. The L–H transition back to the improved confinement occurs on a longer time scale of tens of ms, consistent with the slowing down time scale of fast beam ions. A dramatic decrease in toroidal plasma rotation is systematically observed associated with chirping modes sweeping down to zero frequency. Resonance maps built for the beam-driven chirping modes with the ASCOT (accelerated simulation of charged particle orbits in torodoial devices) code show that the resonant beam ions have orbits near the trapped-passing boundary. The ASCOT modelling assesses how losses of the resonant fast ions caused by the chirping modes with high enough amplitude modify the torque, potentially affecting the plasma rotation.

magnetic confinement fusion↗
DOE OSTI2024

Sparsified time-dependent Fourier neural operators for fusion simulations

This paper presents a sparsified Fourier neural operator for coupled time-dependent partial differential equations (ST-FNO) as an efficient machine learning surrogate for fluid and particle-based fusion codes such as NIMROD (Non-Ideal Magnetohydrodynamics with Rotation - Open Discussion) and GTC (Gyrokinetic Toroidal Code). ST-FNO leverages the structures in the governing equations and utilizes neural operators to represent Green's function-like numerical operators in the corresponding numerical solvers. Once trained, ST-FNO can rapidly and accurately predict dynamics in fusion devices compared with first-principle numerical algorithms. In general, ST-FNO represents an efficient and accurate machine learning surrogate for numerical simulators for multi-variable nonlinear time-dependent partial differential equations, with the proposed architectures and loss functions. The efficacy of ST-FNO has been demonstrated using quiescent H-mode simulation data from NIMROD and kink-mode simulation data from GTC. The ST-FNO H-mode results show orders of magnitude reduction in memory and central processing unit usage in comparison with the numerical solvers in NIMROD when computing fields over a selected poloidal plane. The ST-FNO kink-mode results achieve a factor of 2 reduction in the number of parameters compared to baseline FNO models without accuracy loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗