Search NASA⌕ Search

SEARCH · Search NASA

Results for “Guiding center motion”

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 19 records

Hamiltonian formulation of guiding center motion

The nonrelativistic guiding center motion of a charged particle in a static magnetic field is derived using the Hamiltonian formalism. By repeated application of first-order canonical perturbation theory, the first two adiabatic invariants and their averaged Hamiltonians are obtained, including the first-order correction terms. Other features of guiding center theory are also given, including lowest order drifts and the flux invariant.

Stern, D. P.↗

Extensions of guiding center motion to higher order

In a static magnetic field, some well-known guiding-center equations maintain their form when extended to next order in gyroradius. In these cases, it is only necessary to include the next order term in the magnetic moment series. The differential equation for guiding-center motion which describes both the parallel and perpendicular velocities correctly through first order in gyroradius is given. The question of how to define the guiding center position through second order arises and is discussed, and second order drifts are derived for one usual definition. The toroidal canonical angular momentum, P-phi, of the guiding center in an axisymmetric field is shown to be conserved using the guiding center velocity correct through first order. When second-order motion is included, P-phi is no longer a constant. The above extensions of guiding-center theory help to resolve the different tokamak orbits obtained either by using the guiding-center equations of motion or by using conservation of P-phi.

Northrop, T. G.↗

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.↗

Motion of trapped electrons and protons in Saturn's inner magnetosphere

A summary is given of basic formulas for the guiding center motion of energetic charged particles trapped in a dipolar magnetic field. These formulas for longitudinal drift rates, latitudinal bounce periods, equatorial gyroradii, and equatorial gyroperiods are then stated in convenient numerical form for electrons and protons as functions of kinetic energy E, magnetic shell parameter L, and equatorial pitch angle alpha 0 for a slightly simplified model of the observed magnetic field of Saturn. To aid in the study of the interaction of charged particles with the rings and inner satellites of Saturn, additional formulas are given for the time interval between successive encounters of charged particles with a satellite in a circular prograde orbit and for the energy of electrons whose longitudinal angular velocity is resonant, or synchronous, with the Keplerian angular velocity of such a satellite.

Thomsen, M. F.↗

Kinetic and thermodynamic properties of a convecting plasma in a two-dimensional dipole field

Charged particle guiding center motion is considered in the magnetic field of a two-dimensional ('line') dipole on which is superimposed a small, static, perpendicular electric field. The parallel equation of motion is that of a simple harmonic oscillator for cos theta, the cosine of magnetic colatitude theta. Equations for the perpendicular electric and magnetic drifts are derived as well as their bounce-averaged forms. The latter are solved to yield a bounce-averaged guiding center trajectory, which is the same as that obtained from conversation of magnetic moment mu, longitudinal invariant J, and total (kinetic plus electrostatic) energy K. The algebraic simplicity of the trajectory equations is also manifest in the forms of the invariants. An interesting result is that guiding centers drift in such a way that they preserve the values of their equatorial pitch angles and (equivalently) mirror latitudes. The most general Maxwellian form of the equilibrium one-particle distribution function f is constructed from the invariants, and spatially varying density and pressure moments, parallel and perpendicular to the magnetic field, are identified. Much of the paper deals with the more restricted problem in which f is specified as a bi-Maxwellian over a straight line of finite length in the equatorial plane of the dipole and perpendicular to field lines. This might be thought of as specifying a cross-tail ion injection source; our formalism then describes the subsequent spatial development. The distribution away from the source is a scaled bi-Maxwellian but one that is cut off at large and small kinetic energies, which depend on position. Density and pressure components are reduced from the values they would have if the total content of individual flux tubes convected intact. The equatorial and meridional variations of density and pressure components are examined and compared systematically for the isotropic and highly anisotropic situations. There appears to be little qualitative difference due to anisotropy. An anisotropy measure is defined, and its spatial variation determined as a signature of possible MHD instability. Extreme values are found, larger than at the source, but the plasma beta in such regions is probably so low as to render the effect inconsequential energetically. Finally, the possible consequence of 'nonadia- batic' pressure profiles on electrostatic interchanges is considered, and a boundary delineating stabilizing and destabilizing regions determined.

Huang, T. S.↗

A search for the interplanetary quantity controlling the development of geomagnetic storms

An historical account is presented concerning the evolution of our present concept of geomagnetic storms. The present concept was formulated by Chapman (1927) in his magnetic data statistical studies of 'the initial rise' (now termed the initial phase) 'and subsequent larger decrease' (now termed the main phase) in H, followed by 'slow recovery'. The concept introduced by Alfven in 1940 of guiding center motions of a charged particle in a nonuniform magnetic field (ring currents) is also discussed. By 1963 it became quite certain that the ring current, namely a storm-time Van Allen belt, is formed in the magnetosphere during the storm's main phase. The search then began for the solar wind quantity controlling the development of the main phase. The author then gives a personal account of how our concept of geomagnetic storms has advanced and how new findings based on satellite and ground-based observations have made it possible to arrive at a first-approximation expression for the interplanetary quantity controlling the development of geomagnetic storms. Since a geomagnetic storm is a magnetic manifestation of a magnetospheric storm, which is a nonlinear superposition of intense magnetospheric substorms, the main emphasis is shifted toward the understanding of magnetospheric substorms in order to arrive at the parameters controlling the development of geomagnetic storms.

Akasofu, S.-I.↗

Evidence for specularly reflected ions upstream from the quasi-parallel bow shock

Ion velocity distributions in the form of bunches of gyrating particles traveling along helical paths have been observed moving sunward immediately upstream from quasi-parallel parts of the earth's bow shock using Los Alamos/Garching instruments on ISEE-1 and -2. These distributions have characteristics which indicate that they are produced by the nearly specular reflection at the shock of a portion of the incident solar wind ions. In particular, the guiding center motion and the gyrospeeds of the gyrating ions are quantitatively consistent with simple geometrical considerations for specular reflection. These considerations reveal that specularly reflected ions can escape upstream when the angle between the upstream magnetic field and the local shock normal is less than 45 deg but not when the angle is greater than 45 deg. These upstream gyrating ions are an important signature of one of the processes by which solar wind streaming energy is dissipated into other forms of energy at the shock.

Gosling, J. T.↗

Ions upstream of the earth's bow shock - A theoretical comparison of alternative source populations

The trajectories of ions reflected or leaked upstream from the earth's bow shock and subject solely to the Lorentz force in a steady interplanetary magnetic field B and the V x B electric field are studied theoretically. Expressions are obtained for the guiding center motion and gyromotion in a frame (the Hoffman-Teller frame) moving parallel to the shock surface with sufficient speed to transform the incident solar wind velocity into motion entirely along the interplanetary magnetic field. Equations are derived which transform these motions back to the observer's frame. The predicted upstream motions for four different source models for upstream ions are compared using these expressions: magnetic moment-conserving reflection of solar wind ions, specular reflection of solar wind ions, magnetic moment-conserving leakage of magnetosheath ions, and leakage of magnetosheath ions parallel to the shock normal.

Schwartz, S. J.↗

Thermal equilibrium of a cryogenic magnetized pure electron plasma

The thermal equilibrium correlation properties of a magnetically confined pure electron plasma (McPEP) are related to those of a one-component plasma (OCP). The N-particle spatial distribution rho sub s and the Helmholtz free energy F are evaluated for the McPEP to O(lambda sub d-squared/a-squared), where lambda sub d is the thermal de Broglie wavelength and is an interparticle spacing. The electron gyromotion is allowed to be fully quantized while the guiding center motion is quasi-classical. The distribution rho sub s is shown to be identical to that of a classical OCP with a slightly modified potential. To O(lambda sub d-squared/a-squared) this modification does not affect that part of F that is caused by correlations, as long as certain requirements concerning the size of the plasma are met. This theory is motivated by a current series of experiments that involve the cooling of a magnetically confined pure electron plasma to the cryogenic temperature range.

Dubin, D. H. E.↗

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.↗

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.↗

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.↗

Adiabatic charged particle motion in rapidly rotating magnetospheres

Some aspects of adiabatic drift theory are investigated in the regime where the E x B drift velocity is comparable with or larger than the gyro velocity. Particles undergo five drifts in addition to the E x B gradient, and line curvature drifts have three more terms in the parallel equation of motion. The case of the rapidly rotating rigid magnetic field configuration is found to be an exception due to the corotating particle with a nearly periodic motion. The guiding center drift velocity and parallel differential equation of motion are obtained, and the drift velocity in the rotating frame is found to consist of the expected field gradient and line curvature drift, plus centrifugal and coriolis force induced drifts. It is shown that the second invariant is conserved in the rotating frame by the four drifts, so that a particle slowly drifts around on its drift shell and returns to its original field line. Thus, there is no long term energy change, and any energy change is periodic on the bounce and drift time scale.

Northrop, T. G.↗

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.↗

Guiding center model to interpret neutral particle analyzer results

An attempt is made to model the essential features of ion motion in a plasma. A guiding center model is hypothesized for a plasma heating device whose radial electric fields and axial magnetic fields are crossed. The guiding center model is shown to give results which agree with experimental data from four different plasma heating devices and are helpful for the interpretation of neutral particle analyzer results.

Englert, G. W.↗

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.↗