Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “guiding center”

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

An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields

Here we present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit motion in the small time-step limit, but also recovers all the first-order guiding center drifts as well as the correct gyroradius when stepping over the gyration time-scale. In contrast to previous efforts in this direction, the new scheme also features exact energy conservation. In the derivation of the scheme, we find that a new numerical time-scale is introduced. This scale is analyzed and the resulting restrictions on time-step are derived. Based on this analysis, we develop an adaptive time-stepping strategy the respects these constraints while stepping over the gyration scale when physically justified. It is shown through numerical tests on single-particle motion that the scheme's energy conservation property results in tremendous improvements in accuracy, and that the scheme is able to transition smoothly between magnetized and unmagnetized regimes as a result of the adaptive time-stepping.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Modeling the Effects of Drift Orbit Bifurcation on Radiation Belt Electrons

Abstract Drift orbit bifurcation (DOB) has been suggested to play a major role in the loss and transport of radiation belt electrons since it violates the second adiabatic invariants of particles and makes the third invariant undefined. Results from our guiding‐center test particle simulations using the Tsyganenko‐1989c magnetic field model show that the DOB could affect a broad region of the outer radiation belt, which can penetrate inside the geosynchronous orbit at Kp ≥ 3, and its effects are more significant further away from the Earth, at a higher Kp and for higher electron energies. Specifically, the short‐term simulation results after one electron drift show both traditional and nontraditional DOB transport of electrons, with the nontraditional DOB, caused by a third minimum of the magnetic field strength near the equator, reported here for the first time. Moreover, our results show large ballistic jumps in the second invariant and radial distance for electrons at high equatorial pitch angles after one drift. In addition, long‐term DOB transport coefficients of electrons over many drifts are calculated based on our simulation results. We find that the pitch angle and radial diffusion coefficients of electrons due to DOB could be comparable to or even larger than those caused by electron interactions with chorus and Ultra‐Low‐Frequency waves, respectively. In sum, our results demonstrate that DOB could cause effective loss and transport of radiation belt electrons even in the absence of waves.

Astronomy & Astrophysics↗
DOE OSTI2020

Runaway electron transport in stochastic toroidal magnetic fields

Here, we study the transport and confinement properties of runaway electrons (RE) in the presence of magnetic fields with perturbations producing different levels of stochasticity. We use Kinetic Orbit Runaway Electron Code (KORC) [Carbajal et al. , Phys. Plasmas 24, 042512 (2017) and del-Castillo-Negrete et al. , Phys. Plasmas 25, 056104 (2018)] for simulating the full-orbit (FO) and guiding-center (GC) dynamics of RE in perturbed magnetic fields that exhibit magnetic islands. We extend previous works on this problem [Wingen et al. , Nucl. Fusion 46, 941 (2006); Izzo et al. , Nucl. Fusion 51, 063032 (2011); Papp et al. , Nucl. Fusion 51, 043004 (2011); V. Izzo and P. Parks, Phys. Plasmas 24, 060705 (2017); and Sommariva et al. , Nucl. Fusion 58, 016043 (2018)] by studying in detail full-orbit effects on the RE dynamics. We quantify FO effects on RE transport by performing one-to-one comparisons between FO and GC simulations. It is found that, for the magnetic field configurations considered, GC simulations predict twice the RE losses of FO simulations for 1 MeV and four times the RE losses of FO simulations for 25 MeV. Similarly, we show how different GC and FO dynamics of RE moving around magnetic islands can be, especially in the scenario where the RE Larmor radius is on the order of the size of the magnetic island. We also study the role of rotation of the magnetic islands on RE confinement, and we find that low-frequency toroidal rotation has no observable effect on RE transport in the cases considered. These results shed some light into the potential of avoidance or mitigation mechanisms based on magnetic perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Some mathematics for quasi-symmetry

Quasi-symmetry of a steady magnetic field means integrability of first-order guiding-center motion by a spatial symmetry. Here, we derive many restrictions on the possibilities for a quasi-symmetry. We also derive an analog of the Grad–Shafranov equation for the flux function in a quasi-symmetric magnetohydrostatic field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Spatially dependent modeling and simulation of runaway electron mitigation in DIII-D

New simulations with the Kinetic Orbit Runaway electron (RE) Code (KORC) show RE deconfinement losses to the wall during plasma scrape off are the primary current dissipation mechanism in DIII-D experiments with high-Z impurity injection, and not collisional slowing down. The majority of simulations also exhibit an increase in the RE beam energy due to acceleration by the induced toroidal electric field, even while the RE beam current is decreasing. In this study, KORC integrates RE orbits using the relativistic guiding center equations of motion and incorporates time-sequenced, experimental reconstructions of the magnetic and electric fields and line integrated electron density to construct spatiotemporal models of electron and partially ionized impurity transport in the companion plasma. Comparisons of experimental current evolution and KORC results demonstrate the importance of including Coulomb collisions with partially ionized impurity physics, initial RE energy, pitch angle, and spatial distributions, and spatiotemporal electron and partially ionized impurity transport. The research presented here provides an initial quantification of the efficacy of RE mitigation via injected impurities and identification of the critical role played by loss of confinement due to plasma scrape off on the inner wall as compared to the relatively slow collisional damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Alpha particle channeling in ITER

Alpha particle channeling, the attempt to extract energy from high-energy fusion alpha particles using a combination of high- and low-frequency waves introduced in the plasma with antennae, is studied using a guiding center code. We examine the effect of the combination of cooling caused by a cyclotron mode, particle diffusion due to Alfvén modes introduced for ash removal, the inclusion of microturbulent diffusion, and the slowing down of the alpha particles on electrons. Induced cooling is achieved with a cyclotron resonance, producing additional losses of particles above 3.5 MeV in an amount about equal to the prompt loss. The removal of the low-energy cooled particles by additional modes or by microturbulent diffusion is the hardest challenge for alpha channeling. The cooled alpha particles remain deep in the plasma center, so for this removal, the method used must be capable of providing transport from the plasma center to the edge without also causing losses of high-energy particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Gyrokinetic and extended-MHD simulations of a flow shear stabilized Z-pinch experiment

Axisymmetric (m = 0) gyrokinetic and extended-MHD simulations of a sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The extended-MHD model is electromagnetic and includes the effects of the gyro-viscous pressure tensor, diamagnetic electron and ion heat fluxes, and generalized Ohm's law. Additionally, a prominent feature of this work is that the radial profiles for the plasma density and temperature are taken from the fusion Z-pinch experiment (FuZE), and the magnetic field profile is obtained as a solution of the MHD force balance equation. Such an approach allows to address realistic plasma parameters and provide insights into the current and planned experiments. In particular, it is demonstrated that the radial profiles play an important role in stabilization, as the embedded guiding center (E×B ) drift has a strong radial shear, which can contribute to the Z-pinch stabilization even in the absence of the fluid flow shear. The results of simulations for the FuZE plasma parameters show a decrease in the linear growth rate with an increase in the flow shear; however, full stabilization in the linear regime is not observed even for large (comparable to the Alfvén velocity) radial variations of the axial flow. Nonlinear stability properties of the FuZE plasmas are also studied, and it is found that profile broadening can have a pronounced stabilizing effect in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Normal stability of slow manifolds in nearly periodic Hamiltonian systems

Kruskal [J. Math. Phys. 3, 806 (1962)] showed that each nearly periodic dynamical system admits a formal U(1) symmetry, generated by the so-called roto-rate. We prove that such systems also admit nearly invariant manifolds of each order, near which rapid oscillations are suppressed. We study the nonlinear normal stability of these slow manifolds for nearly periodic Hamiltonian systems on barely symplectic manifolds—manifolds equipped with closed, non-degenerate 2-forms that may be degenerate to leading order. In particular, we establish a sufficient condition for long-term normal stability based on second derivatives of the well-known adiabatic invariant. We use these results to investigate the problem of embedding guiding center dynamics of a magnetized charged particle as a slow manifold in a nearly periodic system. Here, we prove that one previous embedding and two new embeddings enjoy long-term normal stability and thereby strengthen the theoretical justification for these models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Polarized imaging of visible synchrotron emission from runaway electron plateaus in DIII-D

A new approach to estimate the pitch angle of highly energetic postdisruption runaway electrons is presented by imaging the polarized runaway electron synchrotron emission. The vertical to horizontal polarization ratio is found to range between 3 and 14 in image-average studies looking at repeated shots as well as the polarization in a single discharge. Numerical studies with guiding center and full orbit simulations performed with the SOFT and KORC codes are given and compared to the experiment. The pitch angle evolution analyzed using a 1D impurity diffusion model in combination with kinetic simulations is also presented. Furthermore, these studies find agreement in the magnitude of the pitch angle, but not in the temporal rate of change of pitch angle. The measured pitch angle evolution timescale of order 50 ms falls between a modeled timescale of 20 ms in an equilibrium model and >>100 ms in a nonequilibrium model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Modification of Lie's transform perturbation theory for charged particle motion in a magnetic field

It is pointed out that the conventional Lie transform perturbation theory for the guiding center motion of charged particles in a magnetic field needs to be modified for ordering inconsistency. There are two reasons. First, the ordering difference between the temporal variation of gyrophase and that of the other phase space coordinates needs to be taken into account. Second, it is also important to note that the parametric limit of the derivative of a function is not equivalent to the derivative of the limit function. When these facts are taken into account, the near identity transformation rule for one form related to the Lagrangian is modified. With the modified near identity transformation rule, the drift motion of charged particles can be described in the first order, instead of the second order and beyond through a tedious expansion process as in the conventional formulation. Here, this resolves the discrepancy between the direct and Lie transform treatments in the Lagrangian perturbation theory for charged particle motion in a magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Alfvén mode induced particle loss in LHD

Herein losses of beam ions due to Alfvén modes are simulated for discharges in the Large Helical Device at Toki, Japan, using a guiding center code. High frequency beam particle resonances are found for the equilibria for passing particles, providing locations for the destablization of Alfvén modes observed in the device. Losses due to even small modes due to these resonances are significant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau

We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.

Algorithms and data structure↗
DOE OSTI2026

Origins of anisotropic linear magnetoresistance with isotropic mobility in Cd 3 As 2 films on GaAs$[110]$

Thin film synthesis allows for the potential to orient crystals in different orientations, permitting measurement of orientation-dependent material aspects such as band structures and transport anisotropy. Here, Dirac semimetal Cd 3 As 2 films are epitaxially grown on GaAs$[110]$ substrates, which has a $[001]$ orientation in-plane. Films contain domains of two different c-axis orientations resulting from an aligned a-axis in-plane. Magnetoresistance measurements performed along both $[1\bar{1}0]$ and $[001]$ substrate directions reveal similar mobility and carrier concentration, but much larger magnetoresistance along the $[001]$ direction, which can be explained by the guiding center diffusion model as arising from anisotropic disorder and different atomic spacings.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Reduced fast-ion transport calculations of m = n = 1 fishbone-like instabilities in MAST-U

Fast-ion transport associated with an m = n = 1 fishbone-like burst in MAST-U discharge 47128 is investigated using a reduced guiding-center-based transport model (ORBIT-Kick) constrained by multi-diagnostic measurements. The two-dimensional beam-emission spectroscopy system provides measurements of the core poloidal mode structure and fluctuation amplitude, while EFIT++ reconstructions constrained by the motional Stark effect diagnostic indicate a flat q-profile with q 0 > 1⁠, indicating the absence of a resonant q = 1 surface and supporting a pressure-driven infernal-mode interpretation. Analytic m = n = 1 displacement profiles consistent with the measured core mode structure and equilibrium constraints are used as the mode structure inputs to ORBIT-Kick. The calculations show that the dominant resonances occur between the mode and co-passing fast ions, producing redistribution localized near the magnetic axis. Synthetic neutron camera signals from TRANSP-Kick recover up to 90% of the experimentally observed neutron deficit at the time of peak mode amplitude, indicating that the measured m = n = 1 mode is a dominant contributor to core fast-ion transport. However, the synthetic neutron signals recover rapidly, whereas the measured neutron emission continues to decrease after the peak amplitude. In conclusion, the remaining discrepancy may arise from contributions not included in the present single-harmonic model, including higher-m and higher-n harmonics, multi-harmonic interactions, and additional transport mechanisms, motivating future diagnostic development and modeling efforts to resolve and incorporate these additional contributions.

Wong, Henry H. [University of California, Los Ange↗
DOE OSTI2022

Precise stellarator quasi-symmetry can be achieved with electromagnetic coils

Magnetic fields with quasi-symmetry are known to provide good confinement of charged particles and plasmas, but the extent to which quasi-symmetry can be achieved in practice has remained an open question. Recent work [M. Landreman and E. Paul, Phys. Rev. Lett. 128, 035001, 2022] reports the discovery of toroidal magnetic fields that are quasi-symmetric to orders-of-magnitude higher precision than previously known fields. We show that these fields can be accurately produced using electromagnetic coils of only moderate engineering complexity, that is, coils that have low curvature and that are sufficiently separated from each other. Our results demonstrate that these new quasi-symmetric fields are relevant for applications requiring the confinement of energetic charged particles for long time scales, such as nuclear fusion. The coils’ length plays an important role for how well the quasi-symmetric fields can be approximated. For the longest coil set considered and a mean field strength of 1 T, the departure from quasi-symmetry is of the order of Earth’s magnetic field. Additionally, we find that magnetic surfaces extend far outside the plasma boundary used by Landreman and Paul, providing confinement far from the core. Simulations confirm that the magnetic fields generated by the new coils confine particles with high kinetic energy substantially longer than previously known coil configurations. In particular, when scaled to a reactor, the best found configuration loses only 0.04% of energetic particles born at midradius when following guiding center trajectories for 200 ms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Modeling fast ion losses due to tearing and internal kink perturbations in MAST-U

Fast ion (FI) loss properties in the presence of tearing mode and internal kink perturbations are numerically investigated for discharges in the MAST-U spherical tokamak, utilizing the MARS-F magnetohydrodynamic stability code and the REORBIT test particle guiding-center orbit-following module. Here, tracing about 100 000 particle markers sampled from the equilibrium distribution of the neutral-beam injection induced FIs, it is found that about 10% out of the total strike the limiting surface (including the divertor surface) in MAST-U discharge 46943, assuming a maximum perturbation of 100 G inside the plasma (corresponding to ~6 G at the Mirnov probe location at the outboard mid-plane). Detailed particle tracing, assuming a uniform initial distribution in the 2D phase space (at given radial locations), reveals that initially counter-current FIs launched near the plasma edge are subject to significant prompt losses, while almost all initially co-current ions remain well confined at the assumed perturbation level. Most lost FIs strike the lower-half of the limiting surface. Finite gyro-radius effects prevent lost ions from striking the top-outer corner of the super-X divertor chamber. A scan of the perturbation level (based on discharge 45163) reveals, not surprisingly, an approximately linear scaling of the particle loss fraction (for counter-current FIs) with respect to the perturbation amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Construction of quasisymmetric stellarators using a direct coordinate approach

In this work, optimized stellarator configurations and their analytical properties are obtained using a near-axis expansion approach. Such configurations are associated with good confinement as the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak. This makes them appealing as fusion reactor candidates. Using a direct coordinate approach, where the magnetic field and flux surface functions are found explicitly in terms of the position vector at successive orders in the distance to the axis, the set of ordinary differential equations for first and second order quasisymmetry is derived. Examples of quasi-axisymmetric shapes are constructed using a pseudospectral numerical method. Finally, the direct coordinate approach is used to independently verify two hypotheses commonly associated with quasisymmetric magnetic fields, namely that the number of equations exceeds the number of parameters at third order in the expansion and that the near-axis expansion does not prohibit exact quasisymmetry from being achieved on a single flux surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗