Search NASA⌕ Search

SEARCH · Search NASA

Results for “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.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

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

Physics of Boundaries and their Interactions in Space Plasmas

This final report describes a brief summary of our accomplishments during the complete contract period. Traditionally, due to computational limitations, it has been impossible to obtain a global view of the magnetosphere on ion time and spatial scales. As a result, kinetic-simulations have concentrated on the local structure of different magnetospheric discontinuities and boundaries. However, due to the emergence of low cost supercomputers, as well as by taking full advantage of latest advances in data mining and visualization technology, we were able to bypass our planned (proposed) regional simulations and proceed to large-scale 3-D and 2-D global hybrid simulations of the magnetosphere. As a result, although we are only finishing the second year of the proposed activity, much of the original scientific objectives have been surpassed and new avenues of investigation have been opened. Such simulations have led us to possible explanations of some long-standing issues in magnetospheric physics. They have also enables us to make a number of important discoveries predictions, which need to be looked for in satellite data. Examples include the finding that the bow shock can become unstable to the Kelvin-Helmholtz (KH), (2) the discovery of a mechanism for intermittent reconnection due to ion physics which may be relevant to the explanation of the recurrence rate of flux transfer events (FTEs), and (3) this finding that the current sheet in the near-Earth magnetotail region can become unstable to KH with detectable, unique ionospheric signatures. Further, we demonstrated a viable mechanism for the onset of reconnection at the magnetopause, examined the detailed structure of the boundary layer incorporating curvature effects, and provided an explanation for the large core fields observed within FTEs as well as flux ropes in the magnetotail.

Omidi, Nojan↗
NASA NTRS1998

Physics of Boundaries and their Interactions in Space Plasmas

This final report describes a brief summary of our accomplishments during the complete contract period. Traditionally, due to computational limitations, it has been impossible to obtain a global view of the magnetosphere on ion time and spatial scales. As a result, kinetic simulations have concentrated on the local structure of different magnetospheric discontinuities and boundaries. However, due to the emergence of low cost desktop superconductors, as well as by taking full advantage of latest advances in data mining and visualization technology, we were able to bypass our planned (proposed) regional simulations and proceed to large-scale 3-D and 2-D global hybrid simulations of the magnetosphere. As a result, although we are only finishing the second year of the proposed activity, much of the original scientific objectives have been surpassed and new avenues of investigation have been opened. Such simulations have led us to possible explanations of some long-standing issues in magnetospheric physics. They have also enabled us to make a number of important discoveries/predictions, which need to be looked for in satellite data. Examples include: (1) the finding that the bow shock can become unstable to the Kelvin-Helmholtz (KH;) (2) the discovery of a mechanism for intermittent reconnection due to ion physics which may be relevant to the explanation of the recurrence rate of flux transfer events (FTEs;) and (3) the finding that the current sheet in the near-Earth magnetotail region can become unstable to KH with detectable, unique ionospheric signatures. Further, we demonstrated a viable mechanism for the onset of reconnection at the magnetopause, examined the detailed structure of the boundary layer incorporating curvature effects, and provided an explanation for the large core fields observed within FTEs as well as flux ropes in the magnetotail.

Omidi, Nojan↗
NASA NTRS2000

Investigating the Role of the Earth's Ionosphere in Space Weather: Modeling and Observations of High-Latitude Ionospheric Outflows

We will present what we believe to be the first reported observations of up- and downflows of topside ionospheric thermal plasmas from multiple near-simultaneous tracks through the high-latitude topside ionosphere. From several Southern polar passes, it has been possible to construct plots of field-aligned flows of 0+ observed by the Thermal Ion Dynamics Experiment(TIDE) on the POLAR spacecraft near 5000 km altitude together with vertical ion flow observations from one or more DMSP spacecraft near 800 km altitude. We also will present simultaneous observations of POLAR auroral LTVI images with field-aligned flows. We will further show new fluid-kinetic simulations of auroral ionospheric flows in the altitude range 120 km to several R-E as subjected to the synergistic auroral effects of soft electron precipitation, transverse wave-driven ion heating, and hot plasma-driven electric potentials.

Horwitz, J. L.↗
NASA NTRS2002

Mechanisms for the Dissipation of Alfven Waves in Near-Earth Space Plasma

Alfven waves are a major mechanism for the transport of electromagnetic energy from the distant part of the magnetosphere to the near-Earth space. This is especially true for the auroral and polar regions of the Earth. However, the mechanisms for their dissipation have remained illusive. One of the mechanisms is the formation of double layers when the current associated with Alfven waves in the inertial regime interact with density cavities, which either are generated nonlinearly by the waves themselves or are a part of the ambient plasma turbulence. Depending on the strength of the cavities, weak and strong double layers could form. Such double layers are transient; their lifetimes depend on that of the cavities. Thus they impulsively accelerate ions and electrons. Another mechanism is the resonant absorption of broadband Alfven- wave noise by the ions at the ion cyclotron frequencies. But this resonant absorption may not be possible for the very low frequency waves, and it may be more suited for electromagnetic ion cyclotron waves. A third mechanism is the excitation of secondary waves by the drifts of electrons and ions in the Alfven wave fields. It is found that under suitable conditions, the relative drifts between different ion species and/or between electrons and ions are large enough to drive lower hybrid waves, which could cause transverse accelerations of ions and parallel accelerations of electrons. This mechanism is being further studied by means of kinetic simulations using 2.5- and 3-D particle-in-cell codes. The ongoing modeling efforts on space weather require quantitative estimates of energy inputs of various kinds, including the electromagnetic energy. Our studies described here contribute to the methods of determining the estimates of the input from ubiquitous Alfven waves.

Singh, Nagendra↗
NASA NTRS2008

Magnetotail Current Sheet Thinning and Magnetic Reconnection Dynamics in Global Modeling of Substorms

Magnetotail current sheet thinning and magnetic reconnection are key elements of magnetospheric substorms. We utilized the global MHD model BATS-R-US with Adaptive Mesh Refinement developed at the University of Michigan to investigate the formation and dynamic evolution of the magnetotail thin current sheet. The BATSRUS adaptive grid structure allows resolving magnetotail regions with increased current density up to ion kinetic scales. We investigated dynamics of magnetotail current sheet thinning in response to southwards IMF turning. Gradual slow current sheet thinning during the early growth phase become exponentially fast during the last few minutes prior to nightside reconnection onset. The later stage of current sheet thinning is accompanied by earthward flows and rapid suppression of normal magnetic field component $B-z$. Current sheet thinning set the stage for near-earth magnetic reconnection. In collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. One of the major challenges in global MHD modeling of the magnetotail magnetic reconnection is to reproduce fast reconnection rates typically observed in smallscale kinetic simulations. Bursts of fast reconnection cause fast magnetic field reconfiguration typical for magnetospheric substorms. To incorporate nongyritropic effects in diffusion regions we developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated self-consistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites and spatial scales of the diffusion region are updated during the simulations. Such an approach allows quantifying the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions localized near reconnection sites. To clarify the role of smallscale non-MHD effects in diffusion region on the global magnetospheric dynamic and to test different models of dissipation we perform simulations with steady southward IMF driving of the magnetosphere.

Kuznetsova, M. M.↗
NASA NTRS2008

How Ionospheric Ions Populate the Magnetosphere during a Magnetic Storm

Ionospheric oxygen ions have been observed throughout the magnetosphere, from the plasma sheet to the ring current region. I t has been found that the O+ /H+ density ratio in the magnetosphere increases with geomagnetic activity and varies with storm phases. During the magnetic storm in late September to earIy October 2002, Cluster was orbiting in the plasma sheet and ring current regions. At prestorm time, Cluster observed high H+ density and low O+ density in the plasma sheet and lobes. During the storm main phase, 0+ density has increased by 10 times over the pre-storm level. Strong field-aligned beams of O+ were observed in the lobes. O+ fluxes were significantly reduced in the central plasma sheet during the storm recovery. However, 0+ was still evident on the boundaries of the plasma sheet and in the lobes. In order to interpret the Cluster observations and to understand how O+ ions populate the magnetosphere during a magnetic storm, we model the storm in early October 2002 using our global ion kinetic simulation (GIK). We use the LFN global simulation model to produce electric and magnetic fields in the outer magnetosphere, the Strangeway outflow scaling with Delcourt ion trajectories to include ionospheric outflows, and the Fok inner magnetospheric model for the plasmaspheric and ring current response to all particle populations. We find that the observed composition features are qualitatively reproduced by the simulations, with some quantitative differences that point to future improvements in the models.

Fok, Mei-Ching↗
NASA NTRS2011

Effects of Solar Wind Conditions on the Plasma Wake Within a Polar Crater: Preliminary Results

As the solar wind sweeps horizontally past a shadowed lunar crater it simultaneously diffuses toward the surface through an ambipolar process, forming a plasma wake (e.g., Figure 1). Importantly, the resulting electric field structure diverts solar wind protons toward the cold crater floor where they may represent a source of surficial hydrogen. We present a handful of two-dimensional kinetic simulations exploring the range of wake structures and surface particle fluxes possible under various background plasma conditions.

Zimmerman, M. I.↗
NASA NTRS2011

Solar Wind Access to Lunar Polar Craters: Feedback Between Surface Charging and Plasma Expansion

Determining the plasma environment within permanently shadowed lunar craters is critical to understanding local processes such as surface charging, electrostatic dust transport, volatile sequestration, and space weathering. In order to investigate the nature of this plasma environment, the first two-dimensional kinetic simulations of solar wind expansion into a lunar crater with a self-consistent plasma-surface interaction have been undertaken. The present results reveal how the plasma expansion into a crater couples with the electrically-charged lunar surface to produce a quasi-steady wake structure. In particular, there is a negative feedback between surface charging and ambipolar wake potential that allows an equilibrium to be achieved, with secondary electron emission strongly moderating the process. A range of secondary electron yields is explored, and two distinct limits are highlighted in which either surface charging or ambipoiar expansion is responsible for determining the overall wake structure.

Zimmerman, M. I.↗
NASA NTRS2012

Lunar Crater Mini-Wakes: Structure, Variability, and Volatiles

Within a permanently shadowed lunar crater the horizontal flow of solar wind is obstructed by upstream topography, forming a regional plasma mini-wake. In the present work kinetic simulations are utilized to investigate how the most prominent structural aspects of a crater mini-wake are modulated during passage of a solar storm. In addition, the simulated particle fluxes are coupled into an equivalent-circuit model of a roving astronaut,. including triboelectric charging due to frictional contact with the lunar regolith, to characterize charging of the astronaut suit during the various stages of the storm. In some cases, triboelectric charging of the astronaut suit becomes effectively perpetual, representing a critical engineering concern for roving within shadowed lunar regions. Finally, the present results suggest that wake structure plays a critical role in modulating the spatial distribution of volatiles at the lunar poles.

Zimmerman, Michael I.↗
NASA NTRS2012

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

Ringing After a High-Energy Collision: Ambipolar Oscillations During Impact Plasma Expansion

High-velocity impacts on the Moon and other airless bodies deliver energy and material to the lunar surface and exosphere. The target and i mpactor material may become vaporized and ionized to form a collision al plasma that expands outward and eventually becomes collisionless. In the present work, kinetic simulations of the later collision less stage of impact plasma expansion are performed. Attention is paid to characterizing "ambipolar oscillations" in which thermodynamic distur bances propagate outward to generate "ringing" within the expanding e lectron cloud, which could radiate an electromagnetic signature of lo cal plasma conditions. The process is not unlike a beam-plasma intera ction, with the perturbing electron population in the present case ac ting as a highly thermal "beam" that resonates along the expanding de nsity gradient. Understanding the electromagnetic aspects of impact p lasma expansion could provide insight into the lasting effects of nat ural, impact-generated currents on airless surfaces and charging haza rds to human exploration infrastructure and instrumentation.

Zimmerman, M. I.↗
NASA NTRS2015

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↗
NASA NTRS2014

The Link Between Shocks, Turbulence, and Magnetic Reconnection in Collisionless Plasmas

Global hybrid (electron fluid, kinetic ions) and fully kinetic simulations of the magnetosphere have been used to show surprising interconnection between shocks, turbulence and magnetic reconnection. In particular collisionless shocks with their reflected ions that can get upstream before retransmission can generate previously unforeseen phenomena in the post shocked flows: (i) formation of reconnecting current sheets and magnetic islands with sizes up to tens of ion inertial length. (ii) Generation of large scale low frequency electromagnetic waves that are compressed and amplified as they cross the shock. These 'wavefronts' maintain their integrity for tens of ion cyclotron times but eventually disrupt and dissipate their energy. (iii) Rippling of the shock front, which can in turn lead to formation of fast collimated jets extending to hundreds of ion inertial lengths downstream of the shock. The jets, which have high dynamical pressure, 'stir' the downstream region, creating large scale disturbances such as vortices, sunward flows, and can trigger flux ropes along the magnetopause. This phenomenology closes the loop between shocks, turbulence and magnetic reconnection in ways previously unrealized. These interconnections appear generic for the collisionless plasmas typical of space, and are expected even at planar shocks, although they will also occur at curved shocks as occur at planets or around ejecta.

MAGNETOSPHERE↗
NASA NTRS2015

Orientation of X Lines in Asymmetric Magnetic Reconnection-Mass Ratio Dependency

Using fully kinetic simulations, we study the X line orientation of magnetic reconnection in an asymmetric configuration. A spatially localized perturbation is employed to induce a single X line, which has sufficient freedom to choose its orientation in three-dimensional systems. The effect of ion to electron mass ratio is investigated, and the X line appears to bisect the magnetic shear angle across the current sheet in the large mass ratio limit. The orientation can generally be deduced by scanning through the corresponding 2-D simulations to find the reconnection plane that maximizes the peak reconnection electric field. The deviation from the bisection angle in the lower mass ratio limit is consistent with the orientation shift of the most unstable linear tearing mode in an electron-scale current sheet.

Reconnection↗
NASA NTRS2016

Suppression of Collisionless Magnetic Reconnection in Asymmetric Current Sheets

Using fully kinetic simulations, we study the suppression of asymmetric reconnection in the limit where the diamagnetic drift speed >> Alfven speed and the magnetic shear angle is moderate. We demonstrate that the slippage between electrons and the magnetic flux mitigates the suppression and can even result in fast reconnection that lacks one of the outflow jets. Through comparing a case where the diamagnetic drift is supported by the temperature gradient with a companion case that has a density gradient instead, we identify a robust suppression mechanism. The drift of the x-line is slowed down locally by the asymmetric nature of the x-line, and then the x-line is run over and swallowed by the faster-moving following flux.

Magnetic↗
NASA NTRS2016

MMS Observations of Large Guide Field Symmetric Reconnection Between Colliding Reconnection Jets at the Center of a Magnetic Flux Rope at the Magnetopause

We report evidence for reconnection between colliding reconnection jets in a compressed current sheet at the center of a magnetic flux rope at Earth's magnetopause. The reconnection involved nearly symmetric Inflow boundary conditions with a strong guide field of two. The thin (2.5 ion-skin depth (d(sub i) width) current sheet (at approximately 12 d(sub i) downstream of the X line) was well resolved by MMS, which revealed large asymmetries in plasma and field structures in the exhaust. Ion perpendicular heating, electron parallel heating, and density compression occurred on one side of the exhaust, while ion parallel heating and density depression were shifted to the other side. The normal electric field and double out-of-plane (bifurcated) currents spanned almost the entire exhaust. These observations are in good agreement with a kinetic simulation for similar boundary conditions, demonstrating in new detail that the structure of large guide field symmetric reconnection is distinctly different from antiparallel reconnection.

Oieroset, M.↗
NASA NTRS2018

Effect of the Reconnection Electric Field on Electron Distribution Functions in the Diffusion Region of Magnetotail Reconnection

Electron distribution functions in the electron diffusion region during symmetric magnetic reconnection are investigated by means of theory and fully kinetic simulations. Crescent‐like striations are formed in distribution functions in the velocity plane perpendicular to the magnetic field. Using an analytical current sheet, we solve the equation of motion for electrons, and derive the shape of a crescent distribution, as a function of the distance from the neutral line, field gradients, and the reconnection electric field. Each crescent is tilted in the velocity plane because of the acceleration by the reconnection electric field, and multiple stripes appear due to multiple meandering bounces. Applying the theory to distribution functions observed in Earth's magnetotail, we deduce the amplitude of the reconnection electric field.

Bessho, N.↗