Search NASA⌕ Search

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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Motion Imagery and Robotics Application (MIRA): Standards-Based Robotics

This technology development originated from the need to assess the debris threat resulting from soil material erosion induced by landing spacecraft rocket plume impingement on extraterrestrial planetary surfaces. The impact of soil debris was observed to be highly detrimental during NASA s Apollo lunar missions and will pose a threat for any future landings on the Moon, Mars, and other exploration targets. The innovation developed under this program provides a simulation tool that combines modeling of the diverse disciplines of rocket plume impingement gas dynamics, granular soil material liberation, and soil debris particle kinetics into one unified simulation system. The Unified Flow Solver (UFS) developed by CFDRC enabled the efficient, seamless simulation of mixed continuum and rarefied rocket plume flow utilizing a novel direct numerical simulation technique of the Boltzmann gas dynamics equation. The characteristics of the soil granular material response and modeling of the erosion and liberation processes were enabled through novel first principle-based granular mechanics models developed by the University of Florida specifically for the highly irregularly shaped and cohesive lunar regolith material. These tools were integrated into a unique simulation system that accounts for all relevant physics aspects: (1) Modeling of spacecraft rocket plume impingement flow under lunar vacuum environment resulting in a mixed continuum and rarefied flow; (2) Modeling of lunar soil characteristics to capture soil-specific effects of particle size and shape composition, soil layer cohesion and granular flow physics; and (3) Accurate tracking of soil-borne debris particles beginning with aerodynamically driven motion inside the plume to purely ballistic motion in lunar far field conditions.

Martinez, Lindolfo↗
NASA NTRSDate not supplied

Experimental and Computational Analysis of a Rotating Detonation Combustor

Rotating detonation combustor (RDC) research has progressed along parallel experimental and computational paths with limited opportunities for validation. Validation is a crucial step to ensure RDC simulations provide an accurate representation of the physical phenomena, and thus, can be used with confidence for design optimization in different applications. In this study, data from a quasi-two-dimensional computational fluid dynamic simulation in an annular RDC are compared with experimental measurements including particle image velocimetry at 100 kHz. Results show good agreement between measured and simulated dynamic pressures, Chapman-Jouguet speeds, and static pressure distributions. Velocity measurements at the annulus exit show good temporal and quantitative agreement with simulations including spikes from the passing shock.

Detonation↗
NASA NTRS1993

Energy storage and dissipation in the magnetotail during substorms. I - Particle simulations. II - MHD simulations

2D electromagnetic particle simulations are used to investigate the dynamics of the tail during development of substorms under the influence of the pressure in the magnetospheric boundary layer and the dawn-to-dusk electric field. It is shown that pressure pulses result in thinning of the tail current sheet as the magnetic field becomes pinched near the region where the pressure pulse is applied. The pinching leads to the tailward flow of the current sheet plasma and the eventual formation and injection of a plasmoid. Surges in the dawn-to-dusk electric field cause plasma on the flanks to convect into the center of the current sheet, thereby thinning the current sheet. The pressure in the magnetospheric boundary laser is coupled to the dawn-to-dusk electric field through the conductivity of the tail. Changes in the predicted evolution of the magnetosphere during substorms due to changes in the resistivity are investigated under the assumption that MHD theory provides a suitable representation of the global or large-scale evolution of the magnetotail to changes in the solar wind and to reconnection at the dayside magnetopause. It is shown that the overall evolution of the magnetosphere is about the same for three different resistivity distributions with plasmoid formation and ejection in each case.

Winglee, R. M.↗
NASA NTRS1994

Powder agglomeration in a microgravity environment

This is the final report for NASA Grant NAG3-755 entitled 'Powder Agglomeration in a Microgravity Environment.' The research program included both two types of numerical models and two types of experiments. The numerical modeling included the use of Monte Carlo type simulations of agglomerate growth including hydrodynamic screening and molecular dynamics type simulations of the rearrangement of particles within an agglomerate under a gravitational field. Experiments included direct observation of the agglomeration of submicron alumina and indirect observation, using small angle light scattering, of the agglomeration of colloidal silica and aluminum monohydroxide. In the former class of experiments, the powders were constrained to move on a two-dimensional surface oriented to minimize the effect of gravity. In the latter, some experiments involved mixture of suspensions containing particles of opposite charge which resulted in agglomeration on a very short time scale relative to settling under gravity.

Cawley, James D.↗
NASA NTRS1993

Magnetosheath-ionspheric plasma interactions in the cusp/cleft. 2: Mesoscale particle simulations

Ionospheric plasma flowing out from the cusp can be an important source of plasma to the magnetosphere. One source of free energy that can drive this outflow is the injection of magnetosheath plasma into the cusp. Two-dimensional (three velocity) mesoscale particle simulations are used to investigate the particle dynamics in the cusp during southward interplanetary magnetic field. This mesoscale model self-consistently incorporates (1) global influences such as the convection of plasma across the cusp, the action of the mirror force, and the injection of the magnetosheath plasma, and (2) wave-particle interactions which produce the actual coupling between the magnetosheath and ionospheric plasmas. It is shown that, because the thermal speed of the electrons is higher than the bulk motion of the magnetosheath plasma, an upward current is formed on the equatorward edge of the injection region with return currents on either side. However, the poleward return currents are the stronger due to the convection and mirroring of many of the magnetosheath electrons. The electron distribution in this latter region evolves from upward directed streams to single-sided loss cones or possibly electron conics. The ion distribution also shows a variety of distinct features that are produced by spatial and/or temporal effects associated with varying convection patterns and wave-particle interactions. On the equatorward edge the distribution has a downflowing magnetosheath component and an upflowing cold ionospheric component due to continuous convection of ionospheric plasma into the region. In the center of the magnetosheath region, heating from the development of an ion-ion streaming instability causes the suppression of the cold ionospheric component and the formation of downward ionospheric streams. Further poleward there is velocity filtering of ions with low pitch angles, so that the magnetosheath ions develop a ring-beam distribution and the ensuing wave instabilities generate downward ionospheric conics. These downward ionospheric components are eventually turned by the mirror force, leading to the production of upward conics at elevated energies throughout the region.

Winglee, R. M.↗
NASA NTRS2014

Onset of Reconnection in the near Magnetotail: PIC Simulations

Using 2.5-dimensional particle-in-cell (PIC) simulations of magnetotail dynamics, we investigate the onset of reconnection in two-dimensional tail configurations with finite Bz. Reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. We found a clear distinction between stable and unstable cases, dependent on deformation amplitude and ion/electron mass ratio. The threshold appears consistent with electron tearing. The evolution prior to onset, as well as the evolution of stable cases, are largely independent of the mass ratio, governed by integral flux tube entropy conservation as imposed in MHD (magnetohydrodynamics). This suggests that ballooning instability in the tail should not be expected prior to the onset of tearing and reconnection. The onset time and other onset properties depend on the mass ratio, consistent with expectations for electron tearing. At onset,we found electron anisotropies T⊥⁄ T∥ (bottom tail divided by parallel tail) equals 1.1-1.3, raising growth rates and wavenumbers. Our simulations have provided a quantitative onset criterion that is easily evaluated in MHD simulations, provided the spatial resolution is sufficient. The evolution prior to onset and after the formation of a neutral line does not depend on the electron physics, which should permit an approximation by MHD simulations with appropriate dissipation terms.

PIC↗
NASA NTRS2019

Technology Development and Infusion by NASA's Entry Systems Modeling Project

This paper describes recent development of modeling and simulation technologies for entry systems and their infusion into NASA's exploration missions. Technology development is organized and prioritized using a system-level perspective, resulting in four broad technical areas of investment: (1) Thermal protection material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. The paper will highlight key contributions from each of these areas, their impacts from a spacecraft and mission design perspective, and discuss planned future investment. Aspects of each technical area are only briefly summarized here. Thermal protection material modeling is geared toward high-fidelity, predictive models capable of optimizing design performance, post-flight reconstruction, and quantifying thermal protection system reliability. New computational tools and experimental techniques have been applied to Orion, MSL/Mars 2020, Mars InSight, and Mars Sample Return missions. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. In both cases, substantial improvements in model uncertainty have directly impacted the development of mission margin policies, flight instrumentation design and analysis (Orion and Mars 2020), and have even revealed the importance of neglected phenomena like mid-wave infrared radiation of CO2. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics affecting Orion, Commercial Crew, and Mars programs; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate augmentation of aerothermal environments at Mars during dust storms; and studies of roughness-induced heating augmentation relevant to tiled (Orion, Mars 2020) and woven (Mars Sample Return) thermal protection systems. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent, with precision landing requirements driven by Mars human exploration goals.

Barnhardt, Michael D.↗
NASA NTRS1991

Dynamics of charged particles in a Paul radio-frequency quadrupole trap

A molecular-dynamics simulation of hundreds of ions confined in a Paul trap has been performed. The simulation includes the trapped particles' micromotion and interparticle Coulomb interactions. A random walk in velocity was implemented to bring the secular motion to a given temperature which was numerically measured. When the coupling Gamma is large the ions from concentric shells which undergo a quadrupole oscillation at the RF frequency, while the ions within a shell form a 2D hexagonal lattice. Ion clouds at 5 mK show no RF heating for q(z) less than about 0.6, whereas rapid heating is seen for qz = 0.8.

Prestage, J. D.↗
NASA NTRS1979

On solar cosmic rays

The report summarizes the major results obtained over the last four years in the following areas: (1) transport in interplanetary space which includes solar particle observations combined with theoretical efforts to develop models which simulate the dynamics of the event; (2) transport in the solar corona and release of particles in the interplanetary medium; and (3) characteristics of source spectra and acceleration models. Modelling the propagation of solar cosmic rays through interplanetary space and establishing the magnitude and the radial and energy dependence of the diffusion coefficient K are considered. An important result due to statistical observation is that the coronal propagation must be independent of or weakly dependent on both rigidity and energy. Increasing evidence is found for prolonged injection at the sun.

Van Hollebeke, M. A. I.↗
NASA NTRS2013

Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars' Force-Free Magnetospheres: Space Charge Limited Flow

We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.

Timokhin, A. N.↗
NASA NTRS1993

Advanced Multi-Phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

It is known that the simulations of solid rocket motor internal flow field with AL-based propellants require complex multi-phase turbulent flow model. The objective of this study is to develop an advanced particulate multi-phase flow model which includes the effects of particle dynamics, chemical reaction and hot gas flow turbulence. The inclusion of particle agglomeration, particle/gas reaction and mass transfer, particle collision, coalescence and breakup mechanisms in modeling the particle dynamics will allow the proposed model to realistically simulate the flowfield inside a solid rocket motor. The Finite Difference Navier-Stokes numerical code FDNS is used to simulate the steady-state multi-phase particulate flow field for a 3-zone 2-D axisymmetric ASRM model and a 6-zone 3-D ASRM model at launch conditions. The 2-D model includes aft-end cavity and submerged nozzle. The 3-D model represents the whole ASRM geometry, including additional grain port area in the gas cavity and two inhibitors. FDNS is a pressure based finite difference Navier-Stokes flow solver with time-accurate adaptive second-order upwind schemes, standard and extended k-epsilon models with compressibility corrections, multi zone body-fitted formulations, and turbulence particle interaction model. Eulerian/Lagrangian multi-phase solution method is applied for multi-zone mesh. To simulate the chemical reaction, penalty function corrected efficient finite-rate chemistry integration method is used in FDNS. For the AL particle combustion rate, the Hermsen correlation is employed. To simulate the turbulent dispersion of particles, the Gaussian probability distribution with standard deviation equal to (2k/3)(exp 1/2) is used for the random turbulent velocity components. The computational results reveal that the flow field near the juncture of aft-end cavity and the submerged nozzle is very complex. The effects of the turbulent particles affect the flow field significantly and provide better prediction of the ASRM performance. The multi-phase flow analysis using the FDNS code in the present research can be used as a design tool for solid rocket motor applications.

Liaw, Paul↗
NASA NTRS1991

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

Tracking Debris Shed by a Space-Shuttle Launch Vehicle

The DEBRIS software predicts the trajectories of debris particles shed by a space-shuttle launch vehicle during ascent, to aid in assessing potential harm to the space-shuttle orbiter and crew. The user specifies the location of release and other initial conditions for a debris particle. DEBRIS tracks the particle within an overset grid system by means of a computational fluid dynamics (CFD) simulation of the local flow field and a ballistic simulation that takes account of the mass of the particle and its aerodynamic properties in the flow field. The computed particle trajectory is stored in a file to be post-processed by other software for viewing and analyzing the trajectory. DEBRIS supplants a prior debris tracking code that took .15 minutes to calculate a single particle trajectory: DEBRIS can calculate 1,000 trajectories in .20 seconds on a desktop computer. Other improvements over the prior code include adaptive time-stepping to ensure accuracy, forcing at least one step per grid cell to ensure resolution of all CFD-resolved flow features, ability to simulate rebound of debris from surfaces, extensive error checking, a builtin suite of test cases, and dynamic allocation of memory.

Stuart, Phillip C.↗
NASA NTRS1994

Cosmic bullets as particle accelerators and radio sources

We have simulated in two dimensions the dynamical evolution of dense gas clouds(`cosmic bullets') moving supersonically through a uniform low-density medium. The diffusive shock acceleration of relativistic protons (cosmic rays) and their dynamical feedback on the background flow are included by the two-fluid model for this process. The acceleration of relativistic electrons is approximated by a test-particle model, and a passive magnetic field is followed by a simple advection scheme. Strong bow shocks, with Mach numbers similar to that of a bullet's motion, are the most important particle accelerators in the flow, while tail shocks and shocks inside the bullets do not play generally significant roles in this regard. For our simulation parameters, approximately greater than 10% of the initial bullet kinetic energy is converted to a combination of internal energy of gas and cosmic-ray protons by the time the bullets begin to be disrupted. Characteristically, the cosmic rays gain several percent of the available kinetic energy. Bullet destruction on timescales only a little larger than the ram pressure bullet crushing time begins in response to Kelvin-Helmholtz and especially to Rayleigh-Taylor instabilities along the forward bullet surface. For dense bullets this happens before the bullet is stopped by ram pressure. According to our simple model for synchrotron emission from relativistic electrons accelerated and transported within the flows, that emission increases rapidly as the bullet begins to fragment, when it is strongly dominated by field enhancement in sheared flows. Synchrotron emission from the acceleration region within the bow shock is, by contrast, much weaker.

Jones, T. W.↗
NASA NTRS1987

Macroscale particle simulation of kinetic Alfven waves

Two types of simulations of the kinetic Alfven wave are presented using a macroscale particle simulation code (Tanaka and Sato, 1986) which enables individual particle dynamics to be followed in the MHD scales. In this code, low frequency electromagnetic fields are solved by eliminating high frequency oscillations such as the light modes, and the scalar potential electric field is solved by eliminating Lagrangian oscillations. The dependences of the frequency and the Landau damping on the perpendicular wavenumber were studied, and good agreement was found between simulation and theoretical predictions. Some fundamental nonlinear interactions of the kinetic Alfven wave with the particles (parallel acceleration of the electrons) were also noted.

Tanaka, Motohiko↗
NASA NTRS2011

NASA Tech Briefs, September 2011

Topics covered include: Fused Reality for Enhanced Flight Test Capabilities; Thermography to Inspect Insulation of Large Cryogenic Tanks; Crush Test Abuse Stand; Test Generator for MATLAB Simulations; Dynamic Monitoring of Cleanroom Fallout Using an Air Particle Counter; Enhancement to Non-Contacting Stress Measurement of Blade Vibration Frequency; Positively Verifying Mating of Previously Unverifiable Flight Connectors; Radiation-Tolerant Intelligent Memory Stack - RTIMS; Ultra-Low-Dropout Linear Regulator; Excitation of a Parallel Plate Waveguide by an Array of Rectangular Waveguides; FPGA for Power Control of MSL Avionics; UAVSAR Active Electronically Scanned Array; Lockout/Tagout (LOTO) Simulator; Silicon Carbide Mounts for Fabry-Perot Interferometers; Measuring the In-Process Figure, Final Prescription, and System Alignment of Large; Optics and Segmented Mirrors Using Lidar Metrology; Fiber-Reinforced Reactive Nano-Epoxy Composites; Polymerization Initiated at the Sidewalls of Carbon Nanotubes; Metal-Matrix/Hollow-Ceramic-Sphere Composites; Piezoelectrically Enhanced Photocathodes; Iridium-Doped Ruthenium Oxide Catalyst for Oxygen Evolution; Improved Mo-Re VPS Alloys for High-Temperature Uses; Data Service Provider Cost Estimation Tool; Hybrid Power Management-Based Vehicle Architecture; Force Limit System; Levitated Duct Fan (LDF) Aircraft Auxiliary Generator; Compact, Two-Sided Structural Cold Plate Configuration; AN Fitting Reconditioning Tool; Active Response Gravity Offload System; Method and Apparatus for Forming Nanodroplets; Rapid Detection of the Varicella Zoster Virus in Saliva; Improved Devices for Collecting Sweat for Chemical Analysis; Phase-Controlled Magnetic Mirror for Wavefront Correction; and Frame-Transfer Gating Raman Spectroscopy for Time-Resolved Multiscalar Combustion Diagnostics.

Source record↗
NASA NTRS1984

Particle acceleration by turbulent magnetohydro-dynamic reconnection

Test particles in a two dimensional, turbulent MHD simulation are found to undergo significant acceleration. The magnetic field configuration is a periodic sheet pinch which undergoes reconnection. The test particles are trapped in the reconnection region for times of order an Alfven transit time in the large electric fields that characterize the turbulent reconnection process at the relatively large magnetic Reynolds number used in the simulation. The maximum speed attained by these particles is consistent with an analytic estimate which depends on the reconnection electric field, the Alfven speed, and the ratio of Larmor period to the Alfven transit time.

Matthaeus, W. H.↗
NASA NTRS1990

A phase-space fluid simulation of a two-component narrow planetary ring - Particle size segregation, edge formation, and spreading rates

The Krook kinetic equation for identical planetary ring particles is presently generalized for the case of two-component systems, and the equations are numerically solved on the basis of Brophy and Esposito's (1989) phase-space CFD method. Attention is given to the simulation results obtained for a two-component narrow ring, in which the large particles are eight times as massive as the small particles. This ring's unconstrained edge dynamics are resolved by the simulation, and are found to exhibit a sharpening that would not have been expected in single-component rings.

Brophy, Thomas G.↗