Noise and error analysis and optimization in particle-based kinetic plasma simulations
Explore the source record for details and available documents.
E.G. Evstatiev↗
SEARCH · Engineering Papers
Results for “Plasma kinetic simulation”
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.
Kinetic simulation of plasma flows in the inner magnetosphere
A one-dimensional hybrid particle code is used to study the interactions between upflowing thermal ions from conjugate ionospheres. The simulation model allows for multiple species, convection of plasmaspheric flux tubes, and Coulomb self-collisions which conserve momentum and energy locally. The model incorporates a variable-flux boundary condition where the flux, at the boundaries, approaches zero as the plasmasphere fills and equilibrium conditions are reached. The effects of two important processes on plasmaspheric refilling have been considered. The first includes convection of the plasmaspheric flux tube. The second is the interaction of ionospheric thermal plasma and particle injection from an external source. Particle injection seems to play an important role in the evolution of the total particle distribution on the early timescales (t less than 1 hour); however, for late timescales (t larger than 8 days) the thermal plasma from the ionosphere dominates the particle distribution.
Miller, Ronald H.↗
Plasma Wake Simulations and Object Charging in a Shadowed Lunar Crater During a Solar Storm
Within a permanently shadowed lunar crater the horizontal flow of solar wind is obstructed by upstream topography, forming a plasma wake that electrostatically diverts ions toward the crater floor and generates a surface potential that can reach kilovolts. In the present work kinetic plasma simulations are employed to investigate the morphology of a lunar crater wake during passage of a solar storm. Results are cast in terms of leading dimensionless ratios including the ion Mach number, ratio of crater depth to plasma Debye length, peak secondary electron yield, and electron temperature vs. electron impact energy at peak secondary yield. This small set of ratios allows generalization to a much wider range of scenarios. The kinetic simulation results are fed forward into an equivalent-circuit model of a roving astronaut. In very low-plasma-current environments triboelectric charging of the astronaut suit becomes effectively perpetual, representing a critical engineering concern for roving within shadowed lunar regions. Finally, simulated ion fluxes are used to explore sputtering and implantation processes within an idealized crater. It is suggested that the physics of plasma mini-wakes formed in the vicinity of permanently shadowed topography may play a critical role in modulating the enigmatic spatial distribution of volatiles at the lunar poles.
Zimmerman, Michael I.↗
Comment on 'Comparison of observed and calculated implanted ion distributions outside Comet Halley's bow shock' by T. I. Gombosi, M. Neugebauer, A. D. Johnstone, A. J. Coates, and D. E. Huddleston
It is argued that the formulation used by Gombosi (1991) does not distinguish between two fundamental types of resonant interactions (those between cometary pickup ions and LF electromagnetic waves). Two problems arise from using their overly simplified expressions. The wave polarization of the solar wind Alfven wave is not purely right-handed, and values from alpha sub com and alpha sub sw are not well known, nor are the 'polarizations' of the 'backward'-going waves in either case. Gombosi et al. discuss the antisunward vs the sunward propagating solar wind Alfven wave, and assume a reasonable value of 0.2 for alpha sub sw. However, alpha sub com is totally unknown. It is concluded that a more correct analysis of wave polarization can lead to significant changes in the formulation and results from Gombosi et al. In his reply Gombosi argues that his simplified wave-particle interaction model is consistent with Giotto observations and plasma kinetic simulations and it leads to a reasonable approximation of the implanted ion transport coefficients. The more complicated approach suggested by Tsurutani & Thorne (1993) requires the introduction of several ad hoc parameters and it does not lead to any change of the energy diffusion coefficient.
Tsurutani, Bruce T.↗
Kinetic Simulation of Inhomogeneous Plasma with a Variable Sized Grid System
Space plasmas are usually inhomogeneous and irregular. Particle in cell (PIC) plasma simulations, however, often consider the plasma to be homogeneous and periodic. Here non-periodic PIC simulations with inhomogeneities in the density, temperature and background magnetic field are considered. Boundary conditions and a variable sized grid system are discussed, along with an application to a magnetospheric plasma physics problem in the auroral zone.
Schriver, David↗
Kinetic simulation technique for plasma flow in strong external magnetic field
Explore the source record for details and available documents.
Frans H. Ebersohn↗
UCLA IGPP Space Plasma Simulation Group
During the past 10 years the UCLA IGPP Space Plasma Simulation Group has pursued its theoretical effort to develop a Mission Oriented Theory (MOT) for the International Solar Terrestrial Physics (ISTP) program. This effort has been based on a combination of approaches: analytical theory, large scale kinetic (LSK) calculations, global magnetohydrodynamic (MHD) simulations and self-consistent plasma kinetic (SCK) simulations. These models have been used to formulate a global interpretation of local measurements made by the ISTP spacecraft. The regions of applications of the MOT cover most of the magnetosphere: the solar wind, the low- and high-latitude magnetospheric boundary, the near-Earth and distant magnetotail, and the auroral region. Most recent investigations include: plasma processes in the electron foreshock, response of the magnetospheric cusp, particle entry in the magnetosphere, sources of observed distribution functions in the magnetotail, transport of oxygen ions, self-consistent evolution of the magnetotail, substorm studies, effects of explosive reconnection, and auroral acceleration simulations.
Source record↗
Theoretical Technology Research for the International Solar Terrestrial Physics (ISTP) Program
During the last four years the UCLA (University of California, Los Angeles) IGPP (Institute of Geophysics and Planetary Physics) Space Plasma Simulation Group has continued its theoretical effort to develop a Mission Oriented Theory (MOT) for the International Solar Terrestrial Physics (ISTP) program. This effort has been based on a combination of approaches: analytical theory, large-scale kinetic (LSK) calculations, global magnetohydrodynamic (MHD) simulations and self-consistent plasma kinetic (SCK) simulations. These models have been used to formulate a global interpretation of local measurements made by the ISTP spacecraft. The regions of applications of the MOT cover most of the magnetosphere: solar wind, low- and high- latitude magnetospheric boundary, near-Earth and distant magnetotail, and auroral region. Most recent investigations include: plasma processes in the electron foreshock, response of the magnetospheric cusp, particle entry in the magnetosphere, sources of observed distribution functions in the magnetotail, transport of oxygen ions, self-consistent evolution of the magnetotail, substorm studies, effects of explosive reconnection, and auroral acceleration simulations. A complete list of the activities completed under the grant follow.
Ashour-Abdalla, Maha↗
Quantifying the Agyrotropy of Proton and Electron Heating in Turbulent Plasmas
An important aspect of energy dissipation in weakly collisional plasmas is that of energy partitioning between different species (e.g., protons and electrons) and between different energy channels. Here we analyse pressure–strain interaction to quantify the fractions of isotropic compressive, gyrotropic, and nongyrotropic heating for each species. An analysis of kinetic turbulence simulations is compared and contrasted with corresponding observational results from Magnetospheric Multiscale Mission data in the magnetosheath. In assessing how protons and electrons respond to different ingredients of the pressure–strain interaction, we find that compressive heating is stronger than incompressive heating in the magnetosheath for both electrons and protons, while incompressive heating is stronger in kinetic plasma turbulence simulations. Concerning incompressive heating, the gyrotropic contribution for electrons is dominant over the nongyrotropic contribution, while for protons nongyrotropic heating is enhanced in both simulations and observations. Variations with plasma β are also discussed, and protons tend to gain more heating with increasing β.
Yan Yang↗
Simulation studies of plasma waves in the electron foreshock - The transition from reactive to kinetic instability
Particle simulation experiments were used to analyze the electron beam-plasma instability. It is shown that there is a transition from the reactive state of the electron beam-plasma instability to the kinetic instability of Langmuir waves. Quantitative tests, which include an evaluation of the dispersion relation for the evolving non-Maxwellian beam distribution, show that a quasi-linear theory describes the onset of this transition and applies again fully to the kinetic stage. This stage is practically identical to the late stage seen in simulations of plasma waves in the electron foreshock described by Dum (1990).
Dum, C. T.↗
Scaling of Magnetic Reconnection in Relativistic Collisionless Pair Plasmas
Using fully kinetic simulations, we study the scaling of the inflow speed of collisionless magnetic reconnection in electron-positron plasmas from the non-relativistic to ultra-relativistic limit. In the anti-parallel configuration, the inflow speed increases with the upstream magnetization parameter sigma and approaches the speed of light when sigma is greater than O(100), leading to an enhanced reconnection rate. In all regimes, the divergence of the pressure tensor is the dominant term responsible for breaking the frozen-in condition at the x-line. The observed scaling agrees well with a simple model that accounts for the Lorentz contraction of the plasma passing through the diffusion region. The results demonstrate that the aspect ratio of the diffusion region, modified by the compression factor of proper density, remains approximately 0.1 in both the non-relativistic and relativistic limits.
plasma↗
Hybrid simulations of the effects of interstellar pickup hydrogen on the solar wind termination shock
Hybrid (kinetic ions/fluid electrons) plasma simulations are used to study the effects of a population of energetic interstellar pickup hydrogen ions on the solar wind termination shock. The pickup hydrogen is treated as a second ion species in the simulations, and thus the effects of the pick-ups on the shock, as well as the effects of the shock on the pickups, are treated in a fully self-consistent manner. For quasi-perpendicular shocks with 10-20 percent pickup hydrogen the pickup ions manifest themselves in a small foot ahead of the shock ramp caused by pickup ion reflection. For oblique shocks with smaller angles between the field and the shock normal, a large fraction of the pickup ions are reflected and move back upstream where they excite large amplitude magnetosonic waves which steepen into shocklets. These backstreaming pickup ions may provide advance warning of a spacecraft encounter with the termination shock.
Liewer, P. C.↗
Jovian Plasma Torus Interaction with Europa: 3D Hybrid Kinetic Simulation. First results
The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa-moon-magnetosphere system with respect to variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo orbiter mission, and for planning flyby and orbital measurements, (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa (Cassidy etal.,2007;Shematovichetal.,2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyro radius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions).Non-thermal distributions of upstream plasma will be addressed in future work. Photoionization,electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider two models for background plasma:(a) with O(++) ions; (b) with O(++) and S(++) ions. The majority of O2 atmosphere is thermal with an extended cold population (Cassidyetal.,2007). A few first simulations already include an induced magnetic dipole; however, several important effects of induced magnetic fields arising from oceanic shell conductivity will be addressed in later work.
Lipatov, A. S.↗
Effects of equatorially trapped ions on refilling of the plasmasphere
The kinetic aspects of shock formation in response to the equatorial heating of ions are reported by means of small-scale kinetic simulations of countersteaming plasma flows along a spatially varying magnetic field having the features of the dipolar magnetic field lines. Attention is given to the trapping of the field-aligned flowing ions in response to the simulated equatorial heating of ions and the associated self-consistent structures in the electric potentials and the f-plasma flow pattern. Trapping is found to lead to the formation of electrostatic shocks, which propagate away from the 'equator' as seen in the hydrodynamic model of Singh and Torr (1990). The shocks form near the mirror points of the 'equatorially' heated ions. The second stage of the refilling with equatorially trapped ions is reported. A variety of structures in density and temperature distributions as the refilling proceeds slowly is found.
Singh, Nagendra↗
The formation of the wall region - Consequences in the near earth magnetotail
This paper discusses important new findings obtained from global kinetic simulations of magnetotail plasma. A region of strongly nonadiabatic ion acceleration (known as the wall region) exists in the near earth tail and demarcates two very different regimes of ion motion: adiabatic and quasi-adiabatic. A strong enhancement of the cross-tail current occurs on the tailward side of the wall. A comparison of numerical and adiabatic pressure profiles indicates that nonadiabatic processes operating in this region may contribute significantly to a pressure balance relief in the course of quasi-steady magnetospheric convection.
Ashour-Abdalla, Maha↗
Effects of near-Earth stochastic acceleration and reflections of magnetotail ions on the formation of auroral arcs
Findings obtained from global kinetic simulations of magnetotail plasma are discussed. A region of strongly nonadiabatic ion acceleration (known as the 'wall' region) exists in the near earth tail and demarcates two very different regimes of ion motion: adiabatic and quasi-adiabatic. After convection through the wall, ion distributions rapidly become isotropized and thermalized. A strong enhancement of the cross tail current occurs on the tailward side of the wall. Comparison of numerical and adiabatic pressure profiles indicates that nonadiabatic processes operating in this region may contribute significantly to a pressure balance relief in the course of quasi-steady magnetospheric convection.
Ashour-Abdalla, Maha↗
Doubly-Shadowed Regions in Lunar Polar Craters: Hydrogen Accumulation in the Presence of Recursive Plasma Wakes
Permanently shadowed regions (PSRs) of the Moon have been identified as unique environments of extreme cold and comprise a natural cold trap for sequestering volatiles [Paige et al. 2010]. The diverse chemical composition of the LCROSS impact plume provided evidence for a volatile-rich and chemically-complex PSR environment [Cola prete et al. 2010, Schultz et al. 2010]. Additionally, the polar electrostatic environment is highly complex, with the possibility of strong, localized electric fields that divert solar wind ions directly into polar cold traps [Farrell et al. 2010, Zimmerman et al. 2011]. Thus, regional plasma physics processes couple directly with volatile sequestration. In the present work, kinetic simulations show that recursive plasma wake structure arises in the presence of step-like topographic features (Le. doubly-shadowed craters). Combining the plasma code with a numerical sputtering model demonstrates that solar wind protons can be either a hydrogen source via implantation or a volatile loss mechanism via sputtering, depending on properties of the regolith and solar wind. The present model provides a novel theoretical pathway toward understanding the lunar surface/solar wind physical and chemical interactions for complex topography near the poles.
Zimmerman, Michael I.↗
Heating and scattering of ring-beam distributions by turbulence
Pickup ions in the solar wind are initially are born in ring-beam distributions, i.e. f(v) varies as delta(v(sub perpendicular) - V(sub sw)sin(Theta)) delta(v(sub parallel) - V(sub sw)cos(Theta)), where Theta is the angle between the solar wind velocity and the IMF(Interplanetary Magnetic Field), and V(sub sw) is the solar wind speed. Often the distribution has been presumed to relax to a distribution that is isotropic in Theta and essentially mono-energetic, a shell or a 'bi-spherical distribution.' However solar wind turbulence is capable of heating the ring distribution on the timescale of a few tens of gyroperiods, a timescale not greatly distinct from that required for pitch angle scattering to a shell. To describe this effect, we have performed test-particle studies of the heating/scattering of the ring beam distribution by MHD turbulence, adopting various models for the MHD fluctuations, including slab and fully dynamic 2D and 3D incompressible turbulence. Furthermore, a system composed of a cold ion ring and a background plasma is unstable to several kinetic plasma instabilities. We carried out kinetic simulations of the ring beam distribution, showing that plasma instabilities also rapidly energize and scatter particles. Results will be presented comparing relaxation and heating rates of the ring-beam distribution by the various mechanisms.
Gray, P. C.↗