SEARCH · Search NASA

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

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model

The global dynamics of the ionized and neutral gases in the environment of Io plays an important role in the interaction of Jupiter s corotating magnetospheric plasma with Io. Stationary simulations of this problem have already been done using the magnetohydrodynamics (MHD) and the electrodynamics approaches. One of the major results of recent simplified two-fluid model simulations [Saur, J., Neubauer, F.M., Strobel, D.F., Summers, M.E., 2002. J. Geophys. Res. 107 (SMP5), 1-18] was the production of the structure of the double-peak in the magnetic field signature of the Io flyby. These could not be explained before by standard MHD models. In this paper, we present a hybrid simulation for Io with kinetic ions and fluid electrons. This method employs a fluid description for electrons and neutrals, whereas for ions a particle approach is used. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. Our model may provide a much more accurate description for the ion dynamics than previous approaches and allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. The first results of such a simulation of the dynamics of ions in Io s environment are discussed in this paper. Comparison with the Galileo IO flyby results shows that this approach provides an accurate physical basis for the interaction and can therefore naturally reproduce all the observed salient features.

Lipatov, Alexander S.↗

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 NTRS1986

Preparation, analysis and release of simulated interplanetary grains into low Earth orbit

Astronomical observations which reflect the optical and dynamical properties of interstellar and interplanetary grains are the primary means of identifying the shape, size, and the chemistry of extraterrestrial grain materials. Except for recent samplings of extraterrestrial particles in near-Earth orbit and in the stratosphere observations were the only method of deducing the properties of extraterrestrial particles. In order to elucidate the detailed characteristics of observed dust, the observations must be compared with theoretical studies, some of which are discussed in this volume, or compared with terrestrial laboratory experiments. The formation and optical characterization of simulated interstellar and interplanetary dust with particular emphasis on studying the properties on irregularly shaped particles were discussed. Efforts to develop the techniques to allow dust experiments to be carried out in low-Earth orbit were discussed, thus extending the conditions under which dust experiments may be performed.

Stephens, J. R.↗

NASA NTRSDate not supplied

Performance of Coupled Physics Solvers for Multidisciplinary Hypersonic Flow Simulations on Several Classes of Computer Architectures

The application of hypersonic flow simulation tools to realistic flight scenarios will require the coupling of multiple physical effects to the baseline fluid dynamics. Such multiphysics effects can include the aerooelastic response of the airframe or engine components, dynamic transport of atmospheric particles, the deformation of solid-fluid interfaces that can ablate, pyrolyze, or erode, as well as a host of other processes, all of which are governed by unique sets of physical equations and models. Coupling multiple (and potentially disparate) physics solvers to a robust compressible flow solver poses additional challenges related to the stability, performance and scalability of the combined solver. The choices made during the software design process can therefore lead to a variation in simulation efficiency across different computer architectures. In this paper, we will consider two representative multiphysics hypersonic flow scenarios: the interaction of solid particulates with the flow field created by a hypersonic lifting body and the aerooelastic deformation of a model airframe under high-Mach-number flow conditions. For these simulations we explore the behavior of several hypersonic simulation tools, including Kestrel, FUN3D, US3D, and JENRE multiphysics framework, on several high performance computing systems containing various CPU and GPU architectures.

architecture↗

NASA NTRS2013

Comparison of Computed and Measured Vortex Evolution for a UH-60A Rotor in Forward Flight

A Computational Fluid Dynamics (CFD) simulation using the Navier-Stokes equations was performed to determine the evolutionary and dynamical characteristics of the vortex flowfield for a highly flexible aeroelastic UH-60A rotor in forward flight. The experimental wake data were acquired using Particle Image Velocimetry (PIV) during a test of the fullscale UH-60A rotor in the National Full-Scale Aerodynamics Complex 40- by 80-Foot Wind Tunnel. The PIV measurements were made in a stationary cross-flow plane at 90 deg rotor azimuth. The CFD simulation was performed using the OVERFLOW CFD solver loosely coupled with the rotorcraft comprehensive code CAMRAD II. Characteristics of vortices captured in the PIV plane from different blades are compared with CFD calculations. The blade airloads were calculated using two different turbulence models. A limited spatial, temporal, and CFD/comprehensive-code coupling sensitivity analysis was performed in order to verify the unsteady helicopter simulations with a moving rotor grid system.

Evolution↗

NASA NTRS1994

Momentum transport near a magnetic X line in collisionless reconnection

Plasma dynamics and momentum transport near an X line during time-dependent magnetic reconnection in a collisionless plasma are investigated based on two-dimensional particle simulations. We find that a weakly skewed velocity distribution is formed near the magnetic X line, leading to the presence of off-diagonal elements of the plasma pressure tensor. Let the reconnection electric field be in the y direction. The gradients of the off-diagonal elements of the pressure tensor can provide a transport of the y momentum. During the normal magnetic reconnection, the momentum transport associated with the off-diagonal terms of the pressure tensor mediates a transfer of the y momentum from the region near the X line to regions outside the X line. A period of 'reverse magnetic reconnection,' during which the plasma kinetic energy is converted into magnetic energy, is also observed in the simulation. When reverse reconnection occurs, the gradients of the off-diagonal pressure tensor elements can mediate a transfer of y momentum into the X line. It is found that the inertial term also plays a significant role in the force balance near the magnetic X line. An explanation for the origin of the off-diagonal pressure terms is also given in this paper.

Cai, H. J.↗

NASA NTRS1998

Virtual and Experimental Visualization of Flows in Packed Beds of Spheres Simulating Porous Media Flows

A videotape presentation of flow in a packed bed of spheres is provided. The flow experiment consisted of three principal elements: (1) an oil tunnel 76.2 mm by 76.2 mm in cross section, (2) a packed bed of spheres in regular and irregular arrays, and (3) a flow characterization methodology, either (a) full flow field tracking (FFFT) or (b) computational fluid dynamic (CFD) simulation. The refraction indices of the oil and the test array of spheres were closely matched, and the flow was seeded with aluminum oxide particles. Planar laser light provided a two-dimensional projection of the flow field, and a traverse simulated a three-dimensional image of the entire flow field. Light focusing and reflection rendered the spheres black, permitting visualization of the planar circular interfaces in both the axial and transverse directions. Flows were observed near the wall-sphere interface and within the set of spheres. The CFD model required that a representative section of a packed bed be formed and gridded, enclosing and cutting six spheres so that symmetry conditions could be imposed at all cross-boundaries. Simulations had to be made with the flow direction at right angles to that used in the experiments, however, to take advantage of flow symmetry. Careful attention to detail was required for proper gridding. The flow field was three-dimensional and complex to describe, yet the most prominent finding was flow threads, as computed in the representative 'cube' of spheres with face symmetry and conclusively demonstrated experimentally herein. Random packing and bed voids tended to disrupt the laminar flow, creating vortices.

Hendricks, R. C.↗

NASA NTRSDate not supplied

A Hybrid CFD/Engineering Model Tool for Lunar Lander Surface Erosion Prediction

Plume-Surface Interaction (PSI) between lander engine plumes and landing area regolith poses risks to the landers and space exploration missions through view obscuration and high-energy ejecta impacts from eroded particle clouds and vehicle tilting from plume-induced craters. To address these risks, NASA MSFC Fluid Dynamics Branch (ER42) has developed a cascade of predictive simulation capabilities for PSI, including the recent development of a hybrid CFD/engineering model tool that allows for rapid PSI simulations with a descending/ascending vehicle and an eroding surface. The development of this fully-coupled, moving body/eroding surface simulation capability is detailed along with a demonstration of the tool’s capabilities on the Apollo 12 Lunar Module landing. Predictions of eroded mass flow rates as a function of time are verified against Apollo 12 flight data, from which the developed, viscous erosion model has been calibrated.

Plume Surface Interaction↗

NASA NTRSDate not supplied

A Hybrid CFD/Engineering Model Tool for Lunar Lander Surface Erosion Prediction

Plume Surface Interaction↗

NASA NTRS1987

Collisional dynamics of perturbed particle disks in the solar system

Investigations of the collisional evolution of particulate disks subject to the gravitational perturbation of a more massive particle orbiting within the disk are underway. Both numerical N-body simulations using a novel collision algorithm and analytical kinetic theory are being employed to extend our understanding of perturbed disks in planetary rings and during the formation of the solar system. Particular problems proposed for investigation are: (1) The development and testing of general criteria for a small moonlet to clear a gap and produce observable morphological features in planetary rings; (2) The development of detailed models of collisional damping of the wavy edges observed on the Encke division of Saturn's A ring; and (3) The determination of the extent of runaway growth of the few largest planetesimals during the early stages of planetary accretion.

Roberts, W. W.↗

NASA NTRS1994

Charged dust in planetary magnetospheres: Hamiltonian dynamics and numerical simulations for highly charged grains

We use a combination of analytical and numerical methods to investigate the dynamics of charged dust grains in planetary magnetospheres. Our emphasis is on obtaining results valid for particles that are not necessarily dominated either by gravitational or electromagnetic forces. A Hamiltonian formulation of the problem yields exact results, for all values of charge-to-mass ratio, when we introduce two constraints: particles remain in the equatorial plane and the magnetic field is taken as axially symmetric. In particular, we obtain locations of equilibrium points, the frequencies of stable periodic orbits, the topology of separatrices in phase space, and the rate of longitudinal drift. These results are significant for specific applications: motion in the nearly aligned dipolar field of Saturn, and the trajectories of arbitrarily charged particles in complex magnetic fields for limited periods of time after ejection from parent bodies. Since the model is restrictive, we also use numerical integrations of the full three-dimensional equations of motion and illustrate under what conditions the constrained problem yields reasonable results. We show that a large fraction of the intermediately charged and highly charged (gyrating) particles will always be lost to a planet's atmosphere within a few hundred hours, for motion through tilted-dipole magnetic fields. We find that grains must have a very high charge-to-mass ratio in order to be mirrored back to the ring plane. Thus, except perhaps at Saturn where the dipole tilt is very small, the likely inhabitants of the dusty ring systems are those particles that are either nearly Keplerian (weakly charged) grains or grains whose charges place them in the lower end of the intermediate charge zone. Fianlly, we demonstrate the effect of plasma drag on the orbits of gyrating particles to be a rapid decrease in gyroradius followed by a slow radial evolution of the guiding center.

Schaffer, L.↗

NASA NTRS1995

Fast Multipole Methods for Three-Dimensional N-body Problems

We are developing computational tools for the simulations of three-dimensional flows past bodies undergoing arbitrary motions. High resolution viscous vortex methods have been developed that allow for extended simulations of two-dimensional configurations such as vortex generators. Our objective is to extend this methodology to three dimensions and develop a robust computational scheme for the simulation of such flows. A fundamental issue in the use of vortex methods is the ability of employing efficiently large numbers of computational elements to resolve the large range of scales that exist in complex flows. The traditional cost of the method scales as Omicron (N(sup 2)) as the N computational elements/particles induce velocities at each other, making the method unacceptable for simulations involving more than a few tens of thousands of particles. In the last decade fast methods have been developed that have operation counts of Omicron (N log N) or Omicron (N) (referred to as BH and GR respectively) depending on the details of the algorithm. These methods are based on the observation that the effect of a cluster of particles at a certain distance may be approximated by a finite series expansion. In order to exploit this observation we need to decompose the element population spatially into clusters of particles and build a hierarchy of clusters (a tree data structure) - smaller neighboring clusters combine to form a cluster of the next size up in the hierarchy and so on. This hierarchy of clusters allows one to determine efficiently when the approximation is valid. This algorithm is an N-body solver that appears in many fields of engineering and science. Some examples of its diverse use are in astrophysics, molecular dynamics, micro-magnetics, boundary element simulations of electromagnetic problems, and computer animation. More recently these N-body solvers have been implemented and applied in simulations involving vortex methods. Koumoutsakos and Leonard (1995) implemented the GR scheme in two dimensions for vector computer architectures allowing for simulations of bluff body flows using millions of particles. Winckelmans presented three-dimensional, viscous simulations of interacting vortex rings, using vortons and an implementation of a BH scheme for parallel computer architectures. Bhatt presented a vortex filament method to perform inviscid vortex ring interactions, with an alternative implementation of a BH scheme for a Connection Machine parallel computer architecture.

Koumoutsakos, P.↗

NASA NTRS2019

A new environment to simulate the dynamics in the close proximity of rubble-pile asteroids

This paper presents a new environment to simulate close-proximity dynamics around rubble-pile asteroids. The code provides methods for modeling the asteroid’s gravity field and surface through granular dynamics. It implements stateof-the-art techniques to model both gravity and contact interaction between particles: 1) mutual gravity as either direct N2 or Barnes-Hut GPU-parallel octree and 2) contact dynamics with a soft-body (force-based, smooth dynamics), hard-body (constraint-based, non-smooth dynamics), or hybrid (constraint-based with compliance and damping) approach. A very relevant feature of the code is its ability to handle complex-shaped rigid bodies and their full 6D motion. Examples of spacecraft close-proximity scenarios and their numerical simulations are shown.

Ferrari, Fabio↗

NASA NTRS1980

Development of test methods for scale model simulation of aerial applications in the NASA Langley Vortex Facility

Methods have been developed at the Langley Vortex Research Facility to simulate and measure the deposition patterns of aerially applied sprays and granular materials by means of tests with small-scale models of agricultural and dynamically scaled test particles. Interactions between the aircraft wake and the dispersed particles are studied with the aim of modifying wake characteristics and dispersal techniques in order to increase swath width, improve deposition pattern uniformity, and minimize drift. This paper examines the particle sizing analysis, test methods for particle dispersal from the model aircraft, and measurement and computer analysis of test deposition patterns. Results that indicate improved control of chemical drift by winglets are presented to demonstrate test methods.

Jordan, F. L., Jr.↗

NASA NTRS1988

Computer modeling of test particle acceleration at oblique shocks

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.

Decker, Robert B.↗

NASA NTRS1991

SEPAC data analysis in support of the environmental interaction program

Data analyses of the Space Experiments with Particle Accelerators (SEPAC) data and computer modeling were conducted to investigate spacecraft environmental effects associated with injection of electron beams, plasma clouds, and neutral gas clouds from the Shuttle orbiter. The data analysis indicates that Extremely Low Frequency oscillations from 150 to 200 Hz were seen in the Langmuir probe current when the beam was fired in a continuous mode. The strongest oscillations occurred when the ambient pressure was augmented by neutral gas releases from the SEPAC plasma accelerator magnetoplasma-dynamic (MPD) arcjet. To understand the dependence of spacecraft charging potential on beam density and other plasma parameters, a two-dimensional electrostatic particle code was used to simulate the injection of electron beams from an infinite conductor into a plasma. The simulations show that the conductor charging potential depends critically on the reflection coefficient of the conductor surface, which is defined as the percentage of incident particles reflected by the conductor. The ionization effects on spacecraft charging were examined by including interactions of electrons with neutral gas. The simulations show that the conductor charging potential decreases with increasing neutral background density due to the production of secondary electrons near the conductor surface. The simulations also indicate that the beam radius is generally proportional to the beam electron gyroradius when the conductor is charged to a large potential. It appears that the charge buildup at the beam stagnation point causes the beam radial expansion. A survey of the simulation results suggests that the ratio of the beam radius to the beam electron gyroradius increases with the square root of beam density and decreases inversely with beam injection velocity. These results are useful for explaining the spacecraft charging phenomena observed during SEPAC experiments from Spacelab 1.

Lin, Chin S.↗

NASA NTRS2010

On the Role of Global Magnetic Field Configuration in Affecting Ring Current Dynamics

Plasma and field interaction is one important aspect of inner magnetospheric physics. The magnetic field controls particle motion through gradient, curvature drifts and E cross B drift. In this presentation, we show how the global magnetic field affects dynamics of the ring current through simulations of two moderate geomagnetic storms (20 November 2007 and 8-9 March 2008). Preliminary results of coupling the Comprehensive Ring Current Model (CRCM) with a three-dimensional plasma force balance code (to achieve self-consistency in both E and B fields) indicate that inclusion of self-consistency in B tends to mitigate the intensification of the ring current as other similar coupling efforts have shown. In our approach, self-consistency in the electric field is already an existing capability of the CRCM. The magnetic self-consistency is achieved by computing the three-dimensional magnetic field in force balance with anisotropic ring current ion distributions. We discuss the coupling methodology and its further improvement. In addition, comparative studies by using various magnetic field models will be shown. Simulation results will be put into a global context by analyzing the morphology of the ring current, its anisotropy and characteristics ofthe interconnected region 2 field-aligned currents.

Zheng, Y.↗

NASA NTRS2004

A Massively Parallel Hybrid Dusty-Gasdynamics and Kinetic Direct Simulation Monte Carlo Model for Planetary Applications

In order to understand the global structure, dynamics, and physical and chemical processes occurring in the upper atmospheres, exospheres, and ionospheres of the Earth, the other planets, comets and planetary satellites and their interactions with their outer particles and fields environs, it is often necessary to address the fundamentally non-equilibrium aspects of the physical environment. These are regions where complex chemistry, energetics, and electromagnetic field influences are important. Traditional approaches are based largely on hydrodynamic or magnetohydrodynamic (MHD) formulations and are very important and highly useful. However, these methods often have limitations in rarefied physical regimes where the molecular collision rates and ion gyrofrequencies are small and where interactions with ionospheres and upper neutral atmospheres are important. At the University of Michigan we have an established base of experience and expertise in numerical simulations based on particle codes which address these physical regimes. The Principal Investigator, Dr. Michael Combi, has over 20 years of experience in the development of particle-kinetic and hybrid kinetichydrodynamics models and their direct use in data analysis. He has also worked in ground-based and space-based remote observational work and on spacecraft instrument teams. His research has involved studies of cometary atmospheres and ionospheres and their interaction with the solar wind, the neutral gas clouds escaping from Jupiter s moon Io, the interaction of the atmospheres/ionospheres of Io and Europa with Jupiter s corotating magnetosphere, as well as Earth s ionosphere. This report describes our progress during the year. The contained in section 2 of this report will serve as the basis of a paper describing the method and its application to the cometary coma that will be continued under a research and analysis grant that supports various applications of theoretical comet models to understanding the inner comae of comets (grant NAGS- 13239 from the Planetary Atmospheres program).

Combi, Michael R.↗