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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Particle drift in the Earth's plasma sheet

We generalize the derivation of the average gradient/curvature-drift for a flux tube filled with an isotropic distribution of particles at specified kinetic energy. The present treatment is restricted to a two-dimensional magnetic field with zero electric field, but it includes all chaotic and Speiser orbits, which do not correspond to the simple picture of gradient/curvature drift. We assume that particles are evenly distributed throughout the regions of phase space allowed by their energy and canonical momentum. This assumption is closely related but not exactly equivalent to the assumption of isotropic pitch-angle distribution. Our derivation assumes that the maximum Larmor radius is small compared to the scale length for equatorial variations in the flux tube volume, but it does not involve any restrictions on the curvature of the field line. The resulting expression for the drift rate is valid for situations where the particle drift velocity is comparable to the thermal speed in some regions. The apparent implication of this generalized treatment is that the existence of very complex non-adiabatic particle trajectories in the plasma sheet may not invalidate previous estimates of the average rate of particle drift out the sides of the tail, estimates that were made under the assumption of simple guiding-center drifts.

Wolf, R. A.↗
NASA NTRS1993

Long-term containment of energetic particles in coronal loops

Recent observation from the Compton Gamma-Ray Observatory shows that gamma-ray emission after a solar flare can last for as long as 8 hours. There is also evidence that electrons and protons are accelerated only during the impulsive phase of the flare and are subsequently mirror trapped in coronal magnetic loops. This poses the following dilemma: if the magnetic field lines in the loop are simple plane arches, the protons will drift across the cross section of the loop in seconds to minutes, rather than hours. To solve the dilemma, we use guiding center theory to show that long-term containment of energetic protons in a coronal loop is possible if magnetic field lines have enough twist. We also find that in the trapped region of the loop, the twist angle of field lines between the mirror points of a bounce orbit is approximately 2 pi.

Lau, Yun-Tung↗
NASA NTRS1993

Energy content of stormtime ring current from phase space mapping simulations

We perform a phase space mapping study to estimate the enhancement in energy content that results from stormtime particle transport in the equatorial magnetosphere. Our pre-storm phase space distribution is based on a steady-state transport model. Using results from guiding-center simulations of ion transport during model storms having main phases of 3 hr, 6 hr, and 12 hr, we map phase space distributions of ring current protons from the pre-storm distribution in accordance with Liouville's theorem. We find that transport can account for the entire ten to twenty-fold increase in magnetospheric particle energy content typical of a major storm if a realistic stormtime enhancement of the phase space density f is imposed at the nightside tail plasma sheet (represented by an enhancement of f at the neutral line in our model).

Chen, Margaret W.↗
NASA NTRS1993

Stormtime transport of ring current and radiation belt ions

This is an investigation of stormtime particle transport that leads to formation of the ring current. Our method is to trace the guiding-center motion of representative ions (having selected first adiabatic invariants mu) in response to model substorm-associated impulses in the convection electric field. We compare our simulation results qualitatively with existing analytically tractable idealizations of particle transport (direct convective access and radial diffusion) in order to assess the limits of validity of these approximations. For mu approximately less than 10 MeV/G (E approximately less than 10 keV at L equivalent to 3) the ion drift period on the final (ring-current) drift shell of interest (L equivalent to 3) exceeds the duration of the main phase of our model storm, and we find that the transport of ions to this drift shell is appropriately idealized as direct convective access, typically from open drift paths. Ion transport to a final closed drift path from an open (plasma-sheet) drift trajectory is possible for those portions of that drift path that lie outside the mean stormtime separatrix between closed and open drift trajectories, For mu approximately 10-25 MeV/G (110 keV approximately less than E approximately less than 280 keV at L equivalent to 3) the drift period at L equivalent to 3 is comparable to the postulated 3-hr duration of the storm, and the mode of transport is transitional between direct convective access and transport that resembles radial diffusion. (This particle population is transitional between the ring current and radiation belt). For mu approximately greater than 25 MeV/G (radiation-belt ions having E approximately greater than 280 keV at L equivalent to 3) the ion drift period is considerably shorter than the main phase of a typical storm, and ions gain access to the ring-current region essentially via radial diffusion. By computing the mean and mean-square cumulative changes in 1/L among (in this case) 12 representative ions equally spaced in drift time around the steady-state drift shell of interest (L equivalent to 3), we have estimated (from both our forward and our time-reversed simulations) the time-integrated radial-diffusion coefficients D(sup sim)(sub LL) for particles having selected values of mu approximately greater than 15 MeV/G. The results agree surprisingly well with the predictions (D(sup ql)(sub LL)) of quasilinear radial diffusion theory, despite the rather brief duration (approximately 3 hrs) of our model storm and despite the extreme variability (with frequency) of the spectral-density function that characterizes the applied electric field during our model storm. As expected, the values of D(sup sim)(sub LL) deduced (respectively) from our forward and time-reversed simulations agree even better with each other and with D(sup sim)(sub LL) when the impulse amplitudes which characterize the individual substorms of our model storm are systematically reduced.

Chen, Margaret W.↗
NASA NTRS1992

Inferences concerning the magnetospheric source region for auroral breakup

Inferences concerning the magnetospheric source region for auroral arcs obtained from particle measurements on polar orbiting satellites are presented and contrasted with other ideas. An argument that the magnetospheric source region for auroral arc breakup and substorm initiation is along Boundary Plasma Sheet (BPS) magnetic field lines is given. This source region lies beyond a distinct central plasma sheet region and sufficiently far from the Earth that energetic ion motion violates the guiding center approximation (is chaotic). The source region is not constrained to any particular range of distances from the Earth, and substorm initiation may be possible over a wide range of distances from near synchronous orbit to the distant tail. It is also argued that the layer of low energy electrons and velocity dispersed ion beams observed at low altitudes on Aureol 3 is not a different region from the region of auroral arcs. Both comprise the BPS. The two regions occasionally appear distinct at low altitudes because of the effects of arc field aligned potential drops on precipitating particles.

Lyons, L. R.↗
NASA NTRS1994

Simulations of phase space distributions of storm time proton ring current

We use results of guiding-center simulations of ion transport to map phase space densities of the stormtime proton ring current. We model a storm as a sequence of substorm-associated enhancements in the convection electric field. Our pre-storm phase space distribution is an analytical solution to a steady-state transport model in which quiet-time radial diffusion balances charge exchange. This pre-storm phase space spectra at L approximately 2 to 4 reproduce many of the features found in observed quiet-time spectra. Using results from simulations of ion transport during model storms having main phases of 3, 6, and 12 hr, we map phase space distributions from the pre-storm distribution in accordance with Liouville's theorem. We find stormtime enhancements in the phase space densities at energies E approximately 30-160 keV for L approximately 2.5 to 4. These enhancements agree well with the observed stormtime ring current. For storms with shorter main phases (approximately 3 hr), the enhancements are caused mainly by the trapping of ions injected from open night side trajectories, and diffusive transport of higher-energy (greater than or approximately 160 keV) ions contributes little to the stormtime ring current. However, the stormtime ring current is augmented also by the diffusive transport of higher-energy ions (E greater than or approximately 160 keV) durinng stroms having longer main phases (greater than or approximately 6 hr). In order to account for the increase in Dst associated with the formation of the stormtime ring current, we estimate the enhancement in particle-energy content that results from stormtime ion transport in the equatorial magnetosphere. We find that transport alone cannot account for the entire increase in absolute value of Dst typical of a major storm. However, we can account for the entire increase in absolute value of Dst by realistically increasing the stormtime outer boundary value of the phase space density relative to the quiet-time value. We compute the magnetic field produced by the ring current itself and find that radial profiles of the magnetic field depression resemble those obtained from observational data.

Chen, Margaret W.↗
NASA NTRS1993

Simulation of the prompt energization and transport of radiation belt particles during the March 24, 1991 SSC

We model the rapid (about 1 min) formation of a new electron radiation belt at L about or = 2.5 that resulted from the Storm Sudden Commencement (SSC) of March 24, 1991 as observed by the Combined Release and Radiation Effects Satellite (CRRES) satellite. Guided by the observed electric and magnetic fields, we represent the time-dependent magnetospheric electric field during the SSC by an asymmetric bipolar pulse that is associated with the compression and relaxation of the Earth's magnetic field. We follow the electrons using a relativistic guiding center code. The test-particle simulations show that electrons with energies of a few MeV at L greater than 6 were energized up to 40 MeV and transported to L about or = 2.5 during a fraction of their drift period. The energization process conserves the first adiabatic invariant and is enhanced due to resonance of the electron drift motion with the time-varying electric field. Our simulation results, with an initial W(exp -8) energy flux spectra, reproduce the observed electron drift echoes and show that the interplanetary shock impacted the magnetosphere between 1500 and 1800 MLT.

Li, Xinlin↗
NASA NTRS1995

Simulation of proton radiation belt formation during the March 24, 1991 SSC

The rapid formation of a new proton radiation belt at L approximately = 2.5 following the March 24, 1991 Storm Sudden Commencement (SSC) observed at the Combined Release and Radiation Effects Satellite (CRRES) satellite is modeled using a relativistic guiding center test particle code. The SSC is modeled by a bipolar electric field and associated compression and relaxation in the magnetic field, superimposed on a dipole magnetic field. The source population consists of both solar and trapped inner zone protons. The simulations show that while both populations contribute to drift echoes in the 20-80 MeV range, primary conditions is from the solar protons. Proton acceleration by the SSC differs from relativistic electron acceleration in that different source populations contribute and nonrelativistic conservation of the first adiabatic invariation leads to greater energization of protons for a given decrease in L. Model drift echoes and flux distribution in L at the time of injection compare well with CRRES observations.

Hudson, M. K.↗
NASA NTRS1994

Consequences of magnetotail ion dynamics

The trajectories of a large ensemble of particles are calculated in a modified Tsyganenko magnetic field model with a uniform cross-tail electric field. The model magnetotail can be divided into several distinct dynamical regimes of ion motion. Near Earth, where the field lines are dipolar the adiabatic formalism is adequate. In the mid-tail and distant tail, guiding-center theory breaks down and must be replaced by a quasi-adiabatic formalism. There is an important transition region between the adiabatic and quasi-adiabatic regions where ion trajectories become more complicated and no simple analytical description holds. This wall region is characterized by rapid ion acceleration and a major loss of particles to the dusk flank. The moments of the ion distribution function are constructued from the ion trajectories, including density, temperature, and pressure in the x-z and x-y planes. In the noon-midnight meridian plane, parameters are relatively constant except near the Earth, while the x-y plots show strong gradients across the magnetotail. Magnetotail plasma convects earthward, drifts toward dusk, and is squeezed out of the tail in the near-Earth region. A thin current sheet forms in the quasi-adiabatic region, and the pressure tensor has significant off-diagonal terms at its edges. These terms are the result of quasi-adiabatic ion trajectories which lead to azimuthally asymmetric distribution functions capable of maintaining approximate stress balance across the current sheet. Simplified analytical descriptions provide further physical insight into ion dynamics that are observed.

Ashour-Abdalla, Maha↗
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 vector-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. Our formulation 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 NTRS1994

Charged dust in planetary magnetospheres: Hamiltonian dynamics and numerical simulations for highly charged grains

We use a combination of analytical and numerical methods to investigate the dynamics of charged dust grains in planetary magnetospheres. Our emphasis is on obtaining results valid for particles that are not necessarily dominated either by gravitational or electromagnetic forces. A Hamiltonian formulation of the problem yields exact results, for all values of charge-to-mass ratio, when we introduce two constraints: particles remain in the equatorial plane and the magnetic field is taken as axially symmetric. In particular, we obtain locations of equilibrium points, the frequencies of stable periodic orbits, the topology of separatrices in phase space, and the rate of longitudinal drift. These results are significant for specific applications: motion in the nearly aligned dipolar field of Saturn, and the trajectories of arbitrarily charged particles in complex magnetic fields for limited periods of time after ejection from parent bodies. Since the model is restrictive, we also use numerical integrations of the full three-dimensional equations of motion and illustrate under what conditions the constrained problem yields reasonable results. We show that a large fraction of the intermediately charged and highly charged (gyrating) particles will always be lost to a planet's atmosphere within a few hundred hours, for motion through tilted-dipole magnetic fields. We find that grains must have a very high charge-to-mass ratio in order to be mirrored back to the ring plane. Thus, except perhaps at Saturn where the dipole tilt is very small, the likely inhabitants of the dusty ring systems are those particles that are either nearly Keplerian (weakly charged) grains or grains whose charges place them in the lower end of the intermediate charge zone. Fianlly, we demonstrate the effect of plasma drag on the orbits of gyrating particles to be a rapid decrease in gyroradius followed by a slow radial evolution of the guiding center.

Schaffer, L.↗
NASA NTRS1997

Force Balance and Substorm Effects in the Magnetotail

A model of the quiet time middle magnetotail is developed using a consistent orbit tracing technique. The momentum equation is used to calculate geocentric solar magnetospheric components of the particle and electromagnetic forces throughout the current sheet. Ions generate the dominant x and z force components. Electron and ion forces almost cancel in the y direction because the two species drift earthward at comparable speeds. The force viewpoint is applied to a study of some substorm processes. Generation of the rapid flows seen during substorm injection and bursty bulk flow events implies substantial force imbalances. The formation of a substorm diversion loop is one cause of changes in the magnetic field and therefore in the electromagnetic force. It is found that larger forces are produced when the cross-tail current is diverted to the ionosphere than would be produced if the entire tail current system simply decreased. Plasma is accelerated while the forces are unbalanced resulting in field lines within a diversion loop becoming more dipolar. Field lines become more stretched and the plasma sheet becomes thinner outside a diversion loop. Mechanisms that require thin current sheets to produce current disruption then can create additional diversion loops in the newly thinned regions. This process may be important during multiple expansion substorms and in differentiating pseudoexpansions from full substorms. It is found that the tail field model used here can be generated by a variety of particle distribution functions. However, for a given energy distribution the mixture of particle mirror or reflection points is constrained by the consistency requirement. The study of uniqueness also leads to the development of a technique to select guiding center electrons that will produce charge neutrality all along a flux tube containing nonguiding center ions without the imposition of a parallel electric field.

Kaufmann, Richard L.↗
NASA NTRS1996

Electron Flow to a Satellite at High Positive Potential

The Tethered Satellite System (TSS) is designed to deploy a 1.6 m diameter spherical satellite a distance of 20 km above the space shuttle orbiter on an insulated conducting tether. Because of the passage of the conducting tether through the earth's magnetic field, an emf is generated producing a positive satellite potential of about 5000 V. Electron flow under the influence of this high positive potential is the focus of the present analysis. The ionospheric parameters at TSS orbit altitude are; thermal velocity of electrons, 1.9 x 10(exp 5) M/S, thermal velocity of the ions, 1.1 x 10(exp 3) m/s, velocity of the satellite 8 x 10(exp 3) m/s. The electrons, with a Debye length, lambda(D) = 0.49 cm, spiral about the earth's magnetic field lines (0.4 Gauss) with a radius of about 3 cm and the ions spiral with a radius of 5 m. Under these conditions, the electron thermal energy, kT is 0.17 eV. The TSS satellite radius, r(p) is 163 Debye lengths. There is an extensive literature on the interaction of satellites with the near-earth ionospheric plasma. The space charge limitation to the electron current collected by a sphere at positive electrical potential was calculated by Langmuir and Blodgett (1924). Parker and Murphy (1967) recognized the importance of the influence of the earth's magnetic field and used the guiding center approximation to calculate the electron current collected by a positive charged satellite. More recently Ma and Schunk (1989) have calculated the time dependent flow of electrons to a spherical satellite at positive potential utilizing numerical methods and Sheldon (1994) used similar methods to solve this problem for the steady state. In order to analyze some of the phenomena that occurred in the ionosphere during the TSS flights, it would be useful to have analytic expressions for these electron flows. The governing equations are very complex and an exact analytical solution is not likely. An approximate analytical solution is feasible however, and the results of one attempt are presented herein.

Sheldon, John W.↗
NASA NTRS2002

Wave-Particle Interactions Associated with Nongyrotropic Distribution Functions: A Hybrid Simulation Study

Nongyrotropic plasma distribution functions can be formed in regions of space where guiding center motion breaks down as a result of strongly curved and weak ambient magnetic fields. Such are the conditions near the current sheet in the Earth's middle and distant magnetotail, where observations of nongyrotropic ion distributions have been made. Here a systematic parameter study of nongyrotropic proton distributions using electromagnetic hybrid simulations is made. We model the observed nongyrotropic distributions by removing a number of arc length segments from a cold ring distribution and find significant differences with the results of simulations that initially have a gyrotropic ring distribution. Model nongyrotropic distributions with initially small perpendicular thermalization produce growing fluctuations that diffuse the ions into a stable Maxwellian-like distribution within a few proton gyro periods. The growing waves produced by nongyrotropic distributions are similar to the electromagnetic proton cyclotron waves produced by a gyrotropic proton ring distribution in that they propagate parallel to the background magnetic field and occur at frequencies on the order of the proton gyrofrequency, The maximum energy of the fluctuating magnetic field increases as the initial proton distribution is made more nongyrotropic, that is, more highly bunched in perpendicular velocity space. This increase can be as much as twice the energy produced in the gyrotropic case.

Convery, P. D.↗
NASA NTRS2018

Modeling the Proton Radiation Belt With Van Allen Probes Relativistic Electron-Proton Telescope Data

An empirical model of the proton radiation belt is constructed from data taken during 2013-2017 by the Relativistic Electron-Proton Telescopes on the Van Allen Probes satellites. The model intensity is a function of time, kinetic energy in the range 18-600 megaelectronvolts, equatorial pitch angle, and L shell of proton guiding centers. Data are selected, on the basis of energy deposits in each of the nine silicon detectors, to reduce background caused by hard proton energy spectra at low L. Instrument response functions are computed by Monte Carlo integration, using simulated proton paths through a simplified structural model, to account for energy loss in shielding material for protons outside the nominal field of view. Overlap of energy channels, their wide angular response, and changing satellite orientation require the model dependencies on all three independent variables be determined simultaneously. This is done by least squares minimization with a customized steepest descent algorithm. Model uncertainty accounts for statistical data error and systematic error in the simulated instrument response. A proton energy spectrum is also computed from data taken during the 8 January 2014 solar event, to illustrate methods for the simpler case of an isotropic and homogeneous model distribution. Radiation belt and solar proton results are compared to intensities computed with a simplified, on-axis response that can provide a good approximation under limited circumstances.

Selesnick, R. S.↗
NASA NTRS2023

Parameter Dependencies of Early-Stage Tangential Discontinuity-Driven Foreshock Bubbles in Local Hybrid Simulations

Foreshock bubbles (FBs) are significant foreshock transients that can accelerate particles and disturb the magnetosphere-ionosphere system. In the kinetic formation model, foreshock ions interact with the discontinuity by performing partial gyrations to generate currents that change the magnetic field topology around the discontinuity. However, how different foreshock ion properties affect the growth of the field variations is not well understood. Therefore, we use 2-D local hybrid simulations to study the effects of different foreshock ion distributions and properties on the growth of tangential discontinuity (TD)-driven FBs. We discover that for a gyrophase-bunched distribution with an initial phase where the guiding center is on the other side of the TD, the foreshock ions gyrate together across the TD, causing more foreshock ions to cross the TD and leading to a faster expansion of the structure than for a Maxwellian distribution. A ring distribution also yields higher expansion speeds because of the higher projected velocity into the new perpendicular direction. For Maxwellian distributions, there are positive and linear correlations of the FB expansion speeds with the initial foreshock ion densities, thermal speeds, parallel speeds, and sine of the TD magnetic shear angles. These parameter dependencies grow in strength as the structures evolve with time. The foreshock ion distributions and properties that lead to stronger currents produce more significant magnetic field variations and higher expansion speeds. Our study helps quantify the formation and expansion of FBs to forecast their space weather effects and contribution to shock acceleration.

Andrew Vu↗
NASA NTRS2022

Preferential Acceleration of Heavy Ions in a Spontaneously Fragmenting Flare Current Sheet

We study the ion acceleration in a mesoscale, spontaneously fragmenting flare current sheet (SFCS) characterized by the presence of a plasmoid cascade. The main subject of our investigation is to determine whether and how plasmoid cascades at intermediate scales in a fragmented current sheet of a solar flare can impact the (preferential) acceleration of specific ions. The time evolution of the SFCS is obtained from high-resolution 2.5D MHD simulations. The ion trajectories (in the background fields resulting from the MHD model), energies, and pitch angles are calculated using a relativistic test-particle code based on the half-acceleration–rotation–half-acceleration method. For light ions, the main acceleration effects of electromagnetic fields within the SFCS are analyzed using the guiding center approximation. We identify regions with the most-efficient ion acceleration within the SFCS, the accelerator efficiency, and spectra of the accelerated ions. The influence of the charge-to-mass ratio on ion behavior is also studied and resulting ion abundances are compared with observational data. The main ion acceleration takes place in the regions with a strong polarization term, which is part of the first-order Fermi acceleration. Because the term is mass dependent, heavier ions undergo preferential acceleration. The ion energy spectra, abundance-enhancement factors, and differential fluxes, obtained from the model, exhibit power-law profiles, in agreement with observed solar energetic particle events. Nonetheless, the obtained slopes for the abundance-enhancement factor do not exactly match the observed data. The computed slopes and profiles are not sensitive to changes in the initial plasma temperature.

David Kramoliš↗