Test particle method in kinetic theory of a plasma
Many body correlation and conditional probability functions of plasma shield clouds surrounding test particles
SEARCH · Search NASA
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.
Many body correlation and conditional probability functions of plasma shield clouds surrounding test particles
The present evaluation of the basic techniques and illustrative results of charged particle-modeling numerical codes suitable for particle acceleration at oblique, fast-mode collisionless shocks emphasizes the treatment of ions as test particles, calculating particle dynamics through numerical integration along exact phase-space orbits. Attention is given to the acceleration of particles at planar, infinitessimally thin shocks, as well as to plasma simulations in which low-energy ions are injected and accelerated at quasi-perpendicular shocks with internal structure.
A study is made of the mathematical solution of the differential equation of motion of a test particle in the equatorial plane of the Kerr gravitational field, using S (Schwarzschild-like) coordinates. A qualitative solution of this equation leads to the conclusion that there can only be 25 different types of orbits. For each value of a, the results are presented in a master diagram for which h and e are the parameters. A master diagram divides the h, e parameter space into regions such that at each point within one of these regions the types of admissible orbits are qualitatively the same. A pictorial representation of the physical orbits in the r, phi plane is also given.
An approximation method for statistical mechanics is presented and applied to a class of problems which contains a test particle propagation problem. All of the available basic equations used in statistical mechanics are cast in the form of a single equation which is integrodifferential in time and which is then used as the starting point for the construction of the local approximation method. Simplification of the integrodifferential equation is achieved through approximation to the Laplace transform of its kernel. The approximation is valid near the origin in the Laplace space and is based on the assumption of small Laplace variable. No other small parameter is necessary for the construction of this approximation method. The n'th level of approximation is constructed formally, and the first five levels of approximation are calculated explicitly. It is shown that each level of approximation is governed by an inhomogeneous partial differential equation in time with time independent operator coefficients. The order in time of these partial differential equations is found to increase as n does. At n = 0 the most local first order partial differential equation which governs the Markovian limit is regained.
In magnetically confined fusion plasmas, the role played by zonal 𝐸 x 𝐵 flow shear layers in the suppression of turbulent transport is relatively well understood. However, less is understood about the role played by the weak shear regions that arise in the non-monotonic radial electric field profiles often associated with these shear layers. In electrostatic simulations from the global total-𝑓 gyrokinetic particle-in-cell code XGC, we demonstrate how shearless regions with non-zero flow curvature form zonal ‘jets’ that, in conjunction with neighbouring regions of shear, can act as robust barriers to particle transport and turbulence spreading. By isolating quasi-coherent fluctuations radially localised to the zonal jets, we construct a map model for the Lagrangian dynamics of gyrokinetic test particles in the presence of drift waves. We identify the presence of shearless invariant tori in this model and verify that these tori act as partial phase-space transport barriers in the simulations. We also demonstrate how avalanches impinging on these shearless tori cause eddy detachment events that form ‘cold/warm core ring’ structures analogous to those found in oceanic jets, facilitating transport across the barriers without destroying them completely. We discuss how shearless tori may generically arise from tertiary instabilities or other types of discrete eigenmodes, suggesting their potential relevance to broader classes of turbulent fluctuations.
Nonlinear interactions of water group ions with large-amplitude whistler wave packets detected at the leading edge of steepened magnetosonic waves observed near Comet Giacobini-Zinner (GZ) are studied using test particle simulations of water-ion interactions with a model wave based on GZ data. Some of the water ions are found to be decelerated in the steepened portion of the magnetosonic wave to the resonance velocity with the whistler wave packets. Through resonance and related nonlinear interaction with the large-amplitude whistler waves, the water ions become trapped by the packet. An energy balance calculation demonstrates that the trapped ions lose their kinetic energy during the trapped motion in the packet. Thus, the nonlinear trapping motion in the wave structure leads to effective energy transfer from the water group ions to the whistler wave packets in the leading edge of the steepened MHD waves.
Ionization of cometary neutral molecules produces ions which are picked-up by the solar wind. The cometary ion pick-up process for Comet Giacobini-Zinner is studied using a test particle method in which trajectories are numerically calculated for several thousand ions whose initial locations were chosen randomly with probability proportional to the neutral density. The cometary ion distribution function was calculated at several locations upstream of the bow shock. Several types of magnetic fluctuations were included. Circularly and linearly polarized monochromatic waves result in less pitch-angle scattering of pick-up ions than turbulent fluctuations which have about the same amplitude. Some anisotropy remains in the distribution function just upstream of the shock, even for the turbulent fluctuations.
Quasi-linear theory supposes that the energy in resonant fluctuations is small compared to that in the mean magnetic field. This is evident in the fact that the zeroth-order particle trajectories are helices about a mean field B(sub o) that is spatially uniform over many correlation lengths. However, in the solar wind it is often the case that the fluctuating part of the field is comparable in magnitude to the mean part. It is generally expected that quasi-linear theory remains viable for particles that are in resonance with a region of the fluctuation spectrum having only small energy density, but even so, care must be taken when comparing simulations to theoretical predictions. We have performed a series of test-particle simulations to explore the evolution of ion distributions in turbulent situations with varying levels of magnetic fluctuations. As delta-B/B(sub o) is increased the distinctions among absolute pitch angle (defined relative to B(sub o)), local pitch angle (defined relative to B(x)), and magnetic moment become important, some of them exhibiting periodic sloshing unrelated to the nonadiabatic processes of interest. Comparing and contrasting the various runs illustrates the phenomena that must be considered when the premise underlying quasi-linear theory are relaxed.
The particle impact ignition mechanism is considered to be the most effective means for igniting metals in gaseous oxygen environments and is a concern for every oxygen system. For this reason, the National Aeronautics and Space Administration (NASA) commits a large amount of resources attempting to understand and defend against this phenomenon. The most practical means of accomplishing this is through extensive of testing. To this end, over the years NASA has developed and built a facility and various apparatuses to test material or component vulnerability to a particle impact ignition. In testing materials, NASA uses sub-sonic and supersonic nozzles to replicate the conditions that a material is expected to experience in service. However, challenges include obtaining an understanding of the gas dynamics and of the particle behavior through and up to the point of impact at the material under test. One means of addressing these challenges is the use of computational fluid dynamics (CFD). If a CFD model can be built and validated, the model can then be used to confidently interrogate changes to inlet conditions and particle geometries, as well as changes to nozzle geometries to better simulate a material’s in-situ conditions and greatly reduce or eliminate trial-and-error approaches in testing. This paper will show NASA’s approach to and results of using CFD to better understand their particle impact nozzles.
While we have been utilizing standard fluence levels such as those listed in the JESD57 document, we have begun revisiting what an appropriate test fluence is when it comes to qualifying a device for single events. Instead of a fixed fluence level or until a specific number of events occurs, a different thought process is required.
While we have been utilizing standard fluence levels such as those listed in the JESD57 document, we have begun revisiting what an appropriate test fluence is when it comes to qualifying a device for single events. Instead of a fixed fluence level or until a specific number of events occurs, a different thought process is required.
Explore the source record for details and available documents.
Rostoker particle theorem for multispecies plasma in framework of Klimontovich formalism
Atmospheric model for high energy particle acceleration in plasma turbulence
The equation which governs the quasi-linear approximation to the ensemble and gyro-phase averaged one-body probability distribution function is constructed from first principles. This derived equation is subjected to a thorough investigation in order to calculate the possible limitations of the quasi-linear approximation. It is shown that the reduction of this equation to a standard diffusion equation in the Markovian limit can be accomplished through the application of the adiabatic approximation. A numerical solution of the standard diffusion equation in the Markovian limit is obtained for the narrow parallel beam injection. Comparison of the diabatic and adiabatic results explicitly demonstrates the failure of the Markovian description of the probability distribution function. Through the use of a linear time-scale extension the failure of the adiabatic approximation, which leads to the Markovian limit, is shown to be due to mixing of the relaxation and interaction time scales in the presence of the strong mean field.
The asymptotic behavior, for large time, of the quasi-linear diabatic solutions and their local approximations is considered. A time averaging procedure is introduced which yields the averages of these solutions over time intervals which contain only large time values. A discussion of the quasi-linear diabatic solutions which is limited to those solutions that are bounded from below as functions of time is given. It is shown that as the upper limit of the time averaging interval is allowed to approach infinity the time averaged quasi-linear diabatic solutions must approach isotropy (mu-independence). The first derivative with respect to mu of these solutions is also considered. This discussion is limited to first derivatives which are bounded functions of time. It is shown that as the upper limit of the time averaging interval is allowed to approach infinity, the time averaged first derivative must approach zero everywhere in mu except at mu = 0 where it must approach a large value which is calculated. The impact of this large derivative on the quasi-linear expansion scheme is discussed. An H-theorem for the first local approximation to the quasi-linear diabatic solutions is constructed. Without time averaging, the H-theorem is used to determine sufficient conditions for the first local approximate solutions to asymptote, with increasing time, to exactly the same final state which the time averaged quasi-linear diabatic solutions must approach as discussed above.
The motion of a satellite of negligible mass about a massive spherically symmetric object in a space of Schwarzschild geometry is calculated by using ordinary Newtonian methods for an inverse-square gravitational field acted upon by an inverse-cube disturbing function. The disturbing function is expressed in terms of the Keplerian elements of the orbit and substituted in the Lagrange planetary equations. The equations can be integrated to calculate the displacement in position of the satellite due to the relativistic potential. For example, the Beacon Explorer C satellite is found to be displaced by about 17.4 cm after one revolution. For Mercury the advance of the perihelion is about 85 km after one orbit, and the maximum periodic displacement is about 13 km.
Explore the source record for details and available documents.