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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Stormtime ring current and radiation belt ion transport: Simulations and interpretations

We use a dynamical guiding-center model to investigate the stormtime transport of ring current and radiation-belt ions. We trace the motion of representative ions' guiding centers in response to model substorm-associated impulses in the convection electric field for a range of ion energies. Our simple magnetospheric model allows us to compare our numerical results quantitatively with analytical descriptions of particle transport, (e.g., with the quasilinear theory of radial diffusion). We find that 10-145-keV ions gain access to L approximately 3, where they can form the stormtime ring current, mainly from outside the (trapping) region in which particles execute closed drift paths. Conversely, the transport of higher-energy ions (approximately greater than 145 keV at L approximately 3) turns out to resemble radial diffusion. The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with particle energy, since our impulses occur at specific (although randomly determined) times. Despite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of the transport process for approximately greater than 145-keV ions, even for the case of an individual storm. For 4 different realizations of our model storm, the geometric mean discrepancies between diffusion coefficients D(sup sim, sub LL) obtained from the simulations and the quasilinear diffusion coefficient D(sup ql, sub LL) amount to factors of 2.3, 2.3, 1.5, and 3.0, respectively. We have found that these discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp minima and maxima in D(sup ql, sub LL). The mean of the remaining discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) for the 4 different storms then amount to factors of 1.9, 2.1, 1.5, and 2.7, respectively. We find even better agreement when we reduce the impulse amplitudes systematically in a given model storm (e.g., reduction of all the impulse amplitudes by half reduces the discrepancy factor by at least its square root) and also when we average our results over an ensemble of 20 model storms (agreement is within a factor of 1.2 without impulse-amplitude reduction). We use our simulation results also to map phase-space densities f in accordance with Liouville's theorem. We find that the stormtime transport of approximately greater than 145-keV ions produces little change in f-bar the drift-averaged phase-space density on any drift shell of interest. However, the stormtime transport produces a major enhancement from the pre-storm phase-space density at energies approximately 30-145 keV, which are representative of the stormtime ring current.

Lyons, Larry R.↗
DOE OSTI2023

Measurement of ion temperature and toroidal flow during magnetic reconnection with a large guide field

Here, ion temperature and toroidal flow along the guide field direction are measured using a new ion tomographic diagnostic on the Magnetic Reconnection eXperiment (MRX) during magnetic reconnection with a guide field strength of about 1.4 and 2.1 times the strength of the reconnecting component. Strong toroidal flows, beyond what has been measured in anti-parallel and lower guide field conditions on MRX, are observed. Sustained ion heating with no discernible structure within the measurement region is also observed. Probe measurements including Langmuir and Mach probe measurements are made to support the tomographic inversion of line-integrated measurements, as well as to provide local measurements of plasma parameters. Measurements of toroidal velocity and ion temperature are supported with time series data. Energy flow into and out of the X-line region is estimated using a guiding center framework and presented in the Appendix of this manuscript, suggesting an outsized role played by parallel electric field in energizing ions. The guiding center approximation is not well satisfied in the region of interest; however, the estimates provide a springboard for future, further experimentation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electron energization in reconnection: Eulerian vs Lagrangian perspectives

Particle energization due to magnetic reconnection is an important unsolved problem for myriad space and astrophysical plasmas. Electron energization in magnetic reconnection has traditionally been examined from a particle, or Lagrangian, perspective using particle-in-cell (PIC) simulations. Guiding-center analyses of ensembles of PIC particles have suggested that Fermi (curvature drift) acceleration and direct acceleration via the reconnection electric field are the primary electron energization mechanisms. However, both PIC guiding-center ensemble analyses and spacecraft observations are performed in an Eulerian perspective. For this work, we employ the continuum Vlasov–Maxwell solver within the Gkeyll simulation framework to reexamine electron energization from a kinetic continuum, Eulerian, perspective. We separately examine the contribution of each drift energization component to determine the dominant electron energization mechanisms in a moderate guide-field Gkeyll reconnection simulation. In the Eulerian perspective, we find that the diamagnetic and agyrotropic drifts are the primary electron energization mechanisms away from the reconnection x-point, where direct acceleration dominates. We compare the Eulerian (Vlasov Gkeyll) results with the wisdom gained from Lagrangian (PIC) analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Stormtime transport of ring current and radiation belt ions

A dynamical guiding-center simulation model is used to study the stormtime ion transport which leads to the formation of the ring current and diffusion in the radiation belts. Representative ions guiding-center motion in response to model storm-associated impulses in the convection electric field is traced for a range of ion mu values. The present numerical results are compared with previously formulated limiting idealization of particle transport in order to assess the limits of validity of these approximations. For ions having drift periods that exceed the duration of the main phase of the storm, their inward transport to form the stormtime ring current is appropriately described as direct convective access. For ions having drift periods comparable to the duration of the main phase of the storm, there is a transition between direct convective access and transport that resembles radial diffusion. Lower-energy ring-current ions at L of about 3 are freshly injected there from open adiabatic trajectories, whereas the higher-energy ring-current population consists of a mixture of freshly injected and previously trapped ions.

Chen, Margaret W.↗
NASA NTRS2021

East–West Proton Flux Anisotropy Observed with the PAMELA Mission

We present a study of the east–west anisotropy of trapped-proton fluxes in low-Earth orbit based on the measurements of the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics(PAMELA)experiment. The differential intensities of eastward- and westward-traveling protons detected in the South Atlantic Anomaly region were estimated as a function of equatorial pitch angle and drift shell, for six energy bins between80 MeV and 2 GeV. We found that, as a consequence of the strong atmospheric gradient coupled with the large gyro radius in this energy range, the intensities of eastward fluxes exceed those of westward fluxes by a factor of∼10–20. However, the reported directional asymmetry also depends on the sign of the local flux gradient, resulting in more intense westward fluxes beyond the radial distances where the inner belt peaks. PAMELA observations can be used to improve the description of the near-Earth radiation environment at lowest altitudes and highest trapping energies, where current theoretical and empirical models are affected by the largest uncertainties. Unified Astronomy Thesaurus concepts: Cosmic rays(329);Van Allen radiation belts(1758)1. Introduction Low-altitude inner-belt protons are strongly influenced by the density distribution of Earth’s atmosphere, mostly through interactions with its neutral constituents, which induce significant flux anisotropies. In prim is, the atmospheric loss cone results in a steep pitch-angle distribution, which becomes narrower for lower drift shells. A further, azimuthal anisotropy originates from finite gyro radius effects at proton energies in excess of a few tens of MeV(Haerendel1962; Lenchek & Singer1962). In fact, for a given spacecraft position, protons with the same pitch angle but different gyro phase—the azimuth angle associated with the gyration motion—have their guiding centers on different drift shells. In particular, protons from the west and from the east gyrate around magnetic field lines located at higher and lower altitudes, respectively. The guiding-center separationΔhincreaseswith increasing energy, so that protons moving eastward will encounter progressively lower drift-averaged densities, thus experiencing less atmospheric absorption; the opposite situation will occur for protons traveling westward, resulting in an east–west asymmetry of flux intensities whenΔhbecomes comparable to or larger than the flux scale height(Garmire1963; Heckman &Nakano1963).The trapped-flux anisotropy is a relevant aspect of the modeling of the low Earth orbit(LEO)radiation environment, given the significant engineering implications, especially for The Astrophysical Journal,919:114(6pp), 2021 October 1https://doi.org/10.3847/1538-4357/ac1677© 2021. The American Astronomical Society.

A. Bruno↗
NASA NTRS1982

Charged particle motions in the distended magnetospheres of Jupiter and Saturn

Charged particle motion in the guiding center approximation is analyzed for models of the Jovian and Saturnian magnetospheric magnetic fields based on Voyager magnetometer observations. Field lines are traced and exhibit the distention which arises from azimuthally circulating magnetospheric currents. The spatial dependencies of the guiding center bounce period and azimuthal drift rate are investigated for the model fields. Non-dipolar effects in the gradient-curvature drift rate are most important at the equator and affect particles with all mirror latitudes. The effect is a factor of 10-15 for Jupiter with its strong magnetodisc current and 1-2 for Saturn with its more moderate ring current. Limits of adiabaticity, where particle gyroradii become comparable with magnetic scale lengths, are discussed and are shown to occur at quite modest kinetic energies for protons and heavier ions.

Birmingham, T. J.↗
NASA NTRS1981

Diffusion of cold magnetospheric ions

The paper investigates conditions under which a diffusion equation is a valid description of the transport across magnetic field lines of low energy magnetospheric plasma. The case in which the guiding center drift and its fluctuating component are comparable in size to the particle thermal speed is considered. It is concluded that a two-dimensional diffusion equation is valid in the alpha-beta space of the Euler potentials, provided that the fluctuations have a short enough correlation time so that the guiding center drifts remain negligible, and ions in Jupiter's Io plasma torus can thus be described by a diffusion equation if the correlation time is no more than a few rotation periods of the planet.

Birmingham, T. J.↗
NASA NTRS1974

Energy and pitch angle distributions for auroral ions using the current sheet acceleration model

Using a dipole-plus-tail magnetic field model, H(+), He(+2), and O(+6) ions are followed numerically, backward in time, from an output plane perpendicular to the axis of the geomagnetic tail, to their point of entrance to the magnetosphere as solar wind particles in the magnetosheath. An adiabatic or guiding center approximation is used in regions where the particles do not interact directly with the current sheet. A Maxwellian distribution with bulk flow is assumed for solar wind particles in the magnetosheath. Bulk velocity, density, and temperature along the magnetopause are taken from the fluid calculations of Spreiter. Using Liouville's theorem, and varying initial conditions at the output plane, the distribution function is found as a function of energy and pitch angle at the output plane. These results are then mapped to the auroral ionosphere using guiding center theory.

Jaeger, E. F.↗
NASA NTRS1977

Adiabatic particle motion in a nearly drift-free magnetic field: Application to the geomagnetic tail

The guiding center motion of particles in a nearly drift free magnetic field is analyzed in order to investigate the dependence of mean drift velocity on equatorial pitch angle, the variation of local drift velocity along the trajectory, and other properties. The mean drift for adiabatic particles is expressed by means of elliptic integrals. Approximations to the twice-averaged Hamiltonian W near z = O are derived, permitting simple representation of drift paths if an electric potential also exists. In addition, the use of W or of expressions for the longitudinal invariant allows the derivation of the twice averaged Liouville equation and of the corresponding Vlasov equation. Bounce times are calculated (using the drift-free approximation), as are instantaneous guiding center drift velocities, which are then used to provide a numerical check on the formulas for the mean drift.

Stern, D. P.↗
NASA NTRS1992

Ion radial diffusion in an electrostatic impulse model for stormtime ring current formation

Guiding-center simulations of stormtime transport of ring-current and radiation-belt ions having first adiabatic invariants mu is approximately greater than 15 MeV/G (E is approximately greater than 165 keV at L is approximately 3) are surprisingly well described (typically within a factor of approximately less than 4) by the quasilinear theory of radial diffusion. This holds even for the case of an individual model storm characterized by substorm-associated impulses in the convection electric field, provided that the actual spectrum of the electric field is incorporated in the quasilinear theory. Correction of the quasilinear diffusion coefficient D(sub LL)(sup ql) for drift-resonance broadening (so as to define D(sub LL)(sup ql)) reduced the typical discrepancy with the diffusion coefficients D(sub LL)(sup sim) deduced from guiding-center simulations of representative-particle trajectories to a factor of approximately 3. The typical discrepancy was reduced to a factor of approximately 1.4 by averaging D(sub LL)(sup sim), D(sub LL)(sup ql), and D(sub LL)(sup rb) over an ensemble of model storms characterized by different (but statistically equivalent) sets of substorm-onset times.

Chen, Margaret W.↗
NASA NTRS1995

Bounce-Averaged Hamiltonian for Charged Particles in an Axisymmetric but Nondipolar Model Magnetosphere

In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M, J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. The L value of a field line is (by definition) inversely proportional to the flux enclosed by the corresponding magnetic shell of equatorial radius r(sub 0), and the L value at the neutral line (r(sub 0) = b) is denoted L*. The azimuthal coordinate phi measures magnetic local time. The best functional representation found for the normalized difference (L(exp 3)a(exp 3)/mu(sub E))(B(sub m) - B(sub 0)) between mirror-point field B(sub m) and equatorial field B(sub 0) along any field line is a 5-term expansion in powers (2/3 through 6/3) of the quantity X equivalent to (La/mu(sub E))(exp 1/2)K, where K equivalent to (J(exp 2)/8m(sub 0)M)(exp 1/2) is an adiabatically conserved quantity independent of particle energy, m(sub 0) is the rest mass of the particle, and a is the radius of the Earth. This functional form is motivated by results for limiting cases in which particles mirror very near and very far from the magnetic equator. Expansion coefficients corresponding to various powers of X are obtained from least squares fits to numerically computed results for X as a function of L and B(sub m). These are accurately expressible as fourth-order polynomials in (r(sub 0)/b)(exp 3), hence indirectly as functions of L/L* = 3La/2b. This representation, which leads (except for a manageably small region of parameter space) to better than 1% accuracy in the specification of B(sub m) as a function of K and L, allows bounce-averaged guiding center simulations to be performed without actually tracing the bounce motions of individual particles. Bounce-averaged drifts L' (meridional) and phi' (azimuthal) are proportional to derivatives of the Hamiltonian H (sum of kinetic and potential energies) with respect to phi and L, respectively. Our formulation thus provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Schulz, Michael↗
NASA NTRS1993

Evolution of molecular clouds

The evolution of interstellar molecular hydrogen was studied, with a special interest for the formation and evolution of molecular clouds and star formation within them, by a two-dimensional hydrodynamical simulation performed on a rectangular grid of physical sizes on the order of 100 pc. It is filled with an initial density of approx. 1 cm(exp -3), except for one cell (approx. 1 pc(exp 2)) at the center of the grid where an accretion core of 1-10(exp 3) solar masses is placed. The grid is co-moving with the gridcenter that is on a circular orbit around the Galactic center and that also is the guiding center of epicyclic approximation of orbits of the matter surrounding it. The initial radial velocity is zero; to account for differential rotation the initial tangential velocity (i.e. the movement around the galactic center) is proportional to the radial distance to the grid center. The rate is comparable to the rotation rate at the Local Standard of Rest. The influence of galactic rotation is noticed by spiral or elliptical forms, but on much longer time scales than self gravitation and cooling processes. Density and temperature are kept constant at the boundaries and no inflow is allowed along the tangential boundaries.

Sevenster, M.↗
DOE OSTI2023

Response to “Comment on ‘Modification of Lie's transform perturbation theory for charged particle motion in a magnetic field’” [Phys. Plasmas 30 , 104701 (2023)]

Dr. Brizard's comment on my work is based on a conceived procedure that does not come from my work. The defense of his claim that the modification of the so-called standard Lie's transform theory is unnecessary is also unsupported. This response reveals in detail the inconsistency issues in the so-called standard Lie's transform theory by analyzing both its results and root causes. The problem in the so-called standard Lie's transform theory is beyond the issue to take into account the ordering difference between the guiding center motion and gyromotion. Here, the inconsistent commutation of derivative and limit causes another issue. In addition, the so-called standard Lie's transform formulation leads to an unnecessarily lengthy and tedious derivation process for a one or two page task under the singular (or renormalized) formalism described in my paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Minimizing separatrix crossings through isoprominence

Abstract A simple property of magnetic fields that minimizes bouncing to passing type transitions of guiding center orbits is defined and discussed. This property, called isoprominence, is explored through the framework of a near-axis expansion. It is shown that isoprominent magnetic fields for a toroidal configuration exist to all orders in a formal expansion about a magnetic axis. Some key geometric features of these fields are described.

confinement↗
DOE OSTI2025

High order interpolation of magnetic fields with vector potential reconstruction for particle simulations

We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary differential equation (ODE) integrators, such as the Dormand-Prince method. The method provides C(m) continuity and achieves high-order accuracy, making it particularly suited for particle trajectory integration and Poincaré section analysis under optimal integration order and timestep adjustments. Through numerical experiments, we demonstrate that the Hermite interpolation method preserves volume and continuity, which are critical for conserving toroidal canonical momentum and magnetic moment in guiding center simulations, especially over long-term trajectory integration. Furthermore, we analyze the impact of insufficient derivative continuity on Runge-Kutta schemes and show how it degrades accuracy at low error tolerances, introducing discontinuity-induced truncation errors. Lastly, we demonstrate performant Poincaré section analysis in two relevant settings of field data collocated from finite element meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

Electric field corrections in the guidance-center plasma

Electric field autocorrelations for the two-dimensional electrostatic guiding-center plasma are calculated numerically. It is concluded that the autocorrelation, averaged over a thermal equilibrium ensemble, is damped in an approximately exponential fashion, as predicted by Taylor and McNamara. Oscillatory behavior of the type predicted by Taylor and Thompson is not observed.

Joyce, G.↗
NASA NTRS1986

Vibroacoustic payload environment prediction system (VAPEPS): Data base management center remote access guide

A Vibroacoustic Data Base Management Center has been established at the Jet Propulsion Laboratory (JPL). The center utilizes the Vibroacoustic Payload Environment Prediction System (VAPEPS) software package to manage a data base of shuttle and expendable launch vehicle flight and ground test data. Remote terminal access over telephone lines to a dedicated VAPEPS computer system has been established to provide the payload community a convenient means of querying the global VAPEPS data base. This guide describes the functions of the JPL Data Base Management Center and contains instructions for utilizing the resources of the center.

Thomas, V. C.↗
NASA NTRS1988

Optical power splitter and polarization splitter

A 3.times.3 optical guiding arrangement is disclosed for providing both power splitting between a pair of output guides in the 3.times.3 guide structure and polarization splitting (into the TE and TM modes) between the central guide and the pair of outer guides. In accordance with the present invention, the mutual coupling length L and separation distance d between adjacent guides are chosen to provide the desired polarization splitting and total intensity (power) in the outer guides. For example, an input signal of intensity I comprising both the TE and TM modes may be launched in the center guide and with correct choices for L and d, will result in output signals of I.sub.TM /2 each of the outer pair of guides and I.sub.TE from the central guide.

Bogert, Gail A.↗