Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle simulations”

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 217 records · Page 12

Energetic Ions Downtail of the Reconnection Site

Using combined MHD/test particle simulations, we explore characteristics of ion (proton) acceleration tailward of a near-Earth reconnection site. We present spatial distributions and explore acceleration mechanisms and sources of accelerated ions. Acceleration is due primarily due simple crossings of the enhanced electric field near the x-line or in the departing plasmoid. The energetic particle distributions show the expected energy dispersed tailward streaming at the plasma sheet boundary, while equatorial distributions are more complicated, resulting from different acceleration sites within the moving plasmoid. Sources are mostly inside the central plasma sheet dawnward of the plasmoid.

Joachim Birn↗

Electron acceleration at nearly perpendicular collisionless shocks. II - Reflection at curved shocks

Test particle simulations by Krauss-Varban et al. (1989), carried out for plane shocks, have confirmed previous results of Wu (1984) and Leroy and Mangeney (1984) that electrons can be effectively accelerated at nearly perpendicular shocks. This paper investigates the reflection and acceleration of electrons at a nearly perpendicular shock, using two-dimensional test-particle calculations which account for the effect of shock curvature. The computations show that reflected electron fluxes are of the order of observed fluxes. For several reasons, the combined effects of shock curvature are far less severe than anticipated.

Krauss-Varban, D.↗

Lunar Simulant Deposition Technique for Dust Tolerance Studies

A renewed interest in lunar exploration has spawned an array of development efforts for lunar surface assets. These systems depend on the reliable operation of mechanisms and components that may be susceptible to performance degradations or failure due to dust. The Uniform Dust Deposition System was developed at the NASA Glenn Research Center to provide repeatable, uniform, and automated deposition of simulants on surfaces of interest for dust mitigation testing. The system is capable of depositing simulants on test articles up to 60 cm in diameter and 15 cm high in a dry air environment with less than 1 percent relative humidity while keeping users safe from aerosolized dust. The automation of the system allows for high testing throughput while not sacrificing test quality and allows the user to reduce data in parallel. The additional development of a simulant preparation technique complements the repeatability of the deposition physics during testing. The system includes an imaging subsystem that leverages the power of machine learning to count simulant particles and measure their size, thereby allowing for accurate predictions of surface deposition densities (coefficient of determination R^(2) = 0.93) from images alone. The coverage of dust on a surface was shown to be uniform (coefficient of variation CV < 0.11), allowing developers to accurately evaluate the performance of their technology with a prescribed amount of lunar simulant, information that can be used to develop and refine models. The accuracy of the system is currently less than desired for a single deposition run, with a standard deviation (SD) ranging from 18 to 24 mg, or 0.839 to 1.184 mg/sq. cm , for a 5-cm-diameter area. However, the accuracy can be improved by performing multiple deposition runs to build dust to a desired level. Testing has shown that a SD of 0.2 to 0.6 mg, or 0.076 to 0.227 mg/sq. cm, can be achieved for a 5-cm-diameter area using this technique.

Stephen Gerdts↗

Simulation of electrostatic turbulence in the plasma sheet boundary layer with electron currents and bean-shaped ion beams

Plasma data from ISEE-1 show the presence of electron currents as well as energetic ion beams in the plasma sheet boundary layer. Broadband electrostatic noise and low-frequency electromagnetic bursts are detected in the plasma sheet boundary layer, especially in the presence of strong ion flows, currents, and steep spacial gradients in the fluxes of few-keV electrons and ions. Particle simulations have been performed to investigate electrostatic turbulence driven by a cold electron beam and/or ion beams with a bean-shaped velocity distribution. The simulation results show that the counterstreaming ion beams as well as the counterstreaming of the cold electron beam and the ion beam excite ion acoustic waves with a given Doppler-shifted real frequency. However, the effect of the bean-shaped ion velocity distributions reduces the growth rates of ion acoustic instability. The simulation results also show that the slowing down of the ion bean is larger at the larger perpendicular velocity. The wave spectra of the electric fields at some points of the simulations show turbulence generated by growing waves.

Nishikawa, K.-I.↗

Anomalous ion mixing within a Kelvin-Helmholtz vortex in a collisionless plasma

Anomalously fast ion mixing is observed in a hybrid code simulation (particle ions and a massless electron fluid) of the Kelvin-Helmholtz instability in a collisionless plasma. While the traditional view predicts that the ion mixing occurs on a time scale longer than the roll-up time T(v) of the Kelvin-Helmholtz vortex by a factor of (velocity shear scale length)/(ion Larmor radius), the simulation results show that the mixing occurs within a time scale comparable to T(v). It is shown that this fast mixing is due to the scattering of ions by unsteady electromagnetic fields.

Terasawa, T.↗

Simulations and Observations of Heating of Auroral Ion Beams

In the auroral zone, quasi-static parallel electric fields produce beams of ionospheric ions (e.g., H(+), He(+) and 0(+)), which flow outward into the magnetosphere, providing a significant source of ions for the ring current and plasma sheet. Because the velocities to which these beams are accelerated is dependent on the mass of the ions, differential flows between the various ion species can develop which are unstable to an ion-ion streaming instability. Particle simulations and observations from DE 1 are used to investigate the heating of the ion beams produced by this instability. It is shown that there is net transfer of energy from the light ions to the heavy ions, with the heavy ions reaching maximum velocities near the beam velocity of the light ions. Bulk heating of the heavy ions occurs when their relative density is low while high-energy tails are produced when their relative density is high. The heating is primarily parallel to the magnetic field if the difference in the heavy and light ion beam velocities is subsonic while both perpendicular and parallel heating can occur if it is supersonic. In the latter case, very strong heating of an intermediate ion's species such as He(+) can also occur. Comparison with observations shows features consistent with heating via the ion-ion instability including perpendicular heating in the supersonic regime and parallel heating in the subsonic regime and a change in the heating between these regimes as the ratio of the H(+) beam speed to the local sound speed is observed to decrease. This heating is, however, not always observed in association with enhanced wave emissions. This lack of waves is attributed to reabsorption of the waves as the ions become heated.

Winglee, R. M.↗

Experience in using SIMD and MIMD parallelism for computational fluid dynamics

One of the key objectives of the Applied Research Branch in the Numerical Aerodynamic Simulation (NAS) Systems Division at NASA Ames Research Center is the accelerated introduction of highly parallel machines into a fully operational environment. In this report we summarize some of the experiences with the parallel testbed machines at the NAS Applied Research Branch. We discuss the performance results obtained from the implementation of two computational fluid dynamics (CFD) applications, an unstructured grid solver and a particle simulation, on the Connection Machine CM-2 and the Intel iPSC/860.

Simon, Horst D.↗

Magnetic Reconnection in a Quasi‐Parallel Shock: Two‐Dimensional Local Particle‐in‐Cell Simulation

Magnetic reconnection in a quasi‐parallel bow shock is investigated with two‐dimensional local particle‐in‐cell simulations. In the shock transition and downstream regions, large amplitude magnetic fluctuations exist, and abundant current sheets form. In some current sheets, reconnection occurs, and ion‐scale and electron‐scale magnetic islands are generated. In electron‐scale island regions, only electron outflow jets are observed, producing a quadrupolar out‐of‐plane magnetic field pattern, while in ion‐scale islands, both ions and electrons are involved and energized in reconnection. Normalized reconnection rates are obtained to be between around 0.1 to 0.2, and particle acceleration signatures are seen in distribution functions.

Bessho, N.↗

Protocol Dependent Frictional Granular Jamming Simulations: Cyclical, Compression and Expansion

Granular matter takes many paths to pack in natural and industrial processes. The path influences the packing microstructure, particularly for frictional grains. We perform discrete element modeling simulations of different paths to construct packings of frictional spheres. Specifically, we explore four stress-controlled protocols implementing packing expansions and compressions in various combinations thereof. We characterize the eventual packed states through their dependence of the packing fraction and coordination number on packing pressure, identifying non-monotonicities with pressure that correlate with the fraction of frictional contacts. These stress-controlled, bulk-like particle simulations access very low-pressure packings, namely, the marginally stable limit, and demonstrate the strong protocol dependence of frictional granular matter.

DEM↗

High-performance multiprocessor architecture for a 3-D lattice gas model

The lattice gas method has recently emerged as a promising discrete particle simulation method in areas such as fluid dynamics. We present a very high-performance scalable multiprocessor architecture, called ALGE, proposed for the simulation of a realistic 3-D lattice gas model, Henon's 24-bit FCHC isometric model. Each of these VLSI processors is as powerful as a CRAY-2 for this application. ALGE is scalable in the sense that it achieves linear speedup for both fixed and increasing problem sizes with more processors. The core computation of a lattice gas model consists of many repetitions of two alternating phases: particle collision and propagation. Functional decomposition by symmetry group and virtual move are the respective keys to efficient implementation of collision and propagation.

Lee, F.↗

Cold plasma heating in the plasma sheet boundary layer - Theory and simulations

Satellite observations in recent years have confirmed that the plasma sheet boundary layer is a permanent feature of the earth's magnetotail located between the lobe and central plasma sheet during both quiet and active magnetic periods. Distinct features of the boundary layer include field aligned ion beams and intense electrostatic emissions known as broadband electrostatic noise. Since the plasma sheet boundary layer is a spatial feature of the magnetotail, within it will occur thermal mixing of the resident warm boundary layer plasma with inflowing (convecting) cold ionospheric plasma. A theoretical study involving linear theory and nonlinear numerical particle simulations is presented which examines ion beam instabilities in the presence of a thermally mixed hot and cold background plasma. It is found that the free energy in the ion beams can heat the cool ionospheric plasma to ambient plasma sheet boundary layer temperatures via broadband electrostatic noise. These results, along with recent observational reports that ionospheric outflow can account for measured plasma sheet densities, suggest that the ionospheric role in plasma sheet dynamics and content may be as large as the solar wind.

Schriver, David↗

SPH modelling of energy partitioning during impacts on Venus

Impact cratering of the Venusian planetary surface by meteorites was investigated numerically using the Smoothed Particle Hydrodynamics (SPH) method. Venus presently has a dense atmosphere. Vigorous transfer of energy between impacting meteorites, the planetary surface, and the atmosphere is expected during impact events. The investigation concentrated on the effects of the atmosphere on energy partitioning and the flow of ejecta and gas. The SPH method is particularly suitable for studying complex motion, especially because of its ability to be extended to three dimensions. In our simulations, particles representing impactors and targets are initially set to a uniform density, and those of atmosphere are set to be in hydrostatic equilibrium. Target, impactor, and atmosphere are represented by 9800, 80, and 4200 particles, respectively. A Tillotson equation of state for granite is assumed for the target and impactor, and an ideal gas with constant specific heat ratio is used for the atmosphere. Two dimensional axisymmetric geometry was assumed and normal impacts of 10km diameter projectiles with velocities of 5, 10, 20, and 40 km/s, both with and without an atmosphere present were modeled.

Takata, T.↗

Energetic ion acceleration at collisionless shocks

An example is presented from a test particle simulation designed to study ion acceleration at oblique turbulent shocks. For conditions appropriate at interplanetary shocks near 1 AU, it is found that a shock with theta sub B n = 60 deg is capable of producing an energy spectrum extending from 10 keV to approx. 1 MeV in approx 1 hour. In this case total energy gains result primarily from several separate episodes of shock drift acceleration, each of which occurs when particles are scattered back to the shock by magnetic fluctuations in the shock vicinity.

Decker, R. B.↗

Ion Velocity Distributions in Dipolarization Events: Beams in the Vicinity of the Plasma Sheet Boundary

Using combined MHD/test particle simulations, we further explore characteristic ion velocity distributions in relation to magnetotail reconnection and dipolarization events, focusing on distributions at and near the plasma sheet boundary layer (PSBL). Simulated distributions right at the boundary are characterized by a single earthward beam, as discussed earlier. However, farther inside, the distributions consist of multiple beams parallel and antiparallel to the magnetic field, remarkably similar to recent Magnetospheric Multiscale observations. The simulations provide insight into the mechanisms: the lowest earthward beam results from direct acceleration at an earthward propagating dipolarization front (DF), with a return beam at somewhat higher energy. A higher-energy earthward beam results from dual acceleration, first near the reconnection site and then at the DF, again with a corresponding return beam resulting from mirroring closer to Earth. Multiple acceleration at the X line or the propagating DF with intermediate bounces may produce even higher-energy beams. Particles contributing to the lower energy beams are found to originate from the PSBL with thermal source energies, increasing with increasing beam energy. In contrast, the highest-energy beams consist mostly of particles that have entered the acceleration region via cross-tail drift with source energies in the suprathermal range.

Birn, J.↗

Electron distributions in quasi-perpendicular shocks: Simulations

One of the outstanding problems of collisionless shock physics is the microphysical description of electron thermalization and associated wave processes. Observations at the Earth's bow shock have shown various features in the electron velocity distribution such as a loss-cone, temperature anisotropies, and beams before a thermalized, flat-topped distribution is reached downstream. It is generally believed that in addition to the macroscopic fields, which determine the overall mapping of the distribution function, electromagnetic whistler waves and lower hybrid waves are responsible for the pitch-angle scattering and thermalization of the electrons. However, a genuine physical explanation requires analysis of the evolution of the electron distribution in the self-consistent wave fields. Particle simulations are ideally suited to do this, and include processes in which both electrons and ions are involved, such as cross-field streaming instabilities. We show results from l-D and 2-D simulations using an implicit code, which allows time steps large compared to the inverse electron plasma frequency. It is shown that in order to capture the correct propagation characteristics of the waves, a large mass ratio between the protons and electrons is generally necessary. Detailed diagnostics of the waves as well as the evolution of the electron velocity distribution will be presented and will be discussed in the frame work of satellite observations.

Pantellini, F.↗

The production of ion conics by oblique double layers

Magnetized test ions are subjected to acceleration through a numerically simulated oblique double layer in order to determine whether they emerge with velocity vectors aligned with or oblique to the ambient magnetic field. A criterion for oblique alignment, depending on the double-layer parameters and on the external magnetization, is obtained. When it is applied to observed and theoretical auroral double layers, this criterion predicts that accelerated heavy ions will be substantially less magnetic field aligned than will accelerated hydrogen ions, thus suggesting auroral double layers as a source of high-energy ion conics. Test particle simulations are also used to investigate the perpendicular heating of ions at low altitudes by the electric fields associated with moving auroral arcs. The rapid motion of small-scale structures in the arcs is suggested as a source of low-energy conical ion distributions, and the slow drifts of the entire arc forms are inferred to heat ionospheric ions.

Borovsky, J. E.↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

Studies of dynamic processes related to active experiments in space plasmas

This is the final report for grant NAGw-2055, 'Studies of Dynamic Processes Related to Active Experiments in Space Plasmas', covering research performed at the University of Michigan. The grant was awarded to study: (1) theoretical and data analysis of data from the CHARGE-2 rocket experiment (1keV; 1-46 mA electron beam ejections) and the Spacelab-2 shuttle experiment (1keV; 100 mA); (2) studies of the interaction of an electron beam, emitted from an ionospheric platform, with the ambient neutral atmosphere and plasma by means of a newly developed computer simulation model, relating model predictions with CHARGE-2 observations of return currents observed during electron beam emissions; and (3) development of a self-consistent model for the charge distribution on a moving conducting tether in a magnetized plasma and for the potential structure in the plasma surrounding the tether. Our main results include: (1) the computer code developed for the interaction of electrons beams with the neutral atmosphere and plasma is able to model observed return fluxes to the CHARGE-2 sounding rocket payload; and (2) a 3-D electromagnetic and relativistic particle simulation code was developed.

Banks, Peter M.↗