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 73 records · Page 4

Gas-Particle Interaction Model Development in Plume Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The GCD PSI Project, planned to be completed over a four year period, contains a Predictive Simulation Capability (PSC) Element focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the PSC Element, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling portion of the GCD PSI project. The objectives of this task are to investigate and construct models of particle drag and gas-particle cloud interactions leading to what is known as particle turbulent kinetic energy (PTKE). While drag on a lone particle suspended in a flow gas of gas is well-known, the drag and drag-induced dynamics of a cloud of particles in a gas flow are not. The PSC element of the PSI project requires models of gas-particle cloud interactions to implement in the Gas Granular Flow Solver (GGFS) in order to arrive at a predictive simulation capability for PSI-induced soil erosion, cratering and ejecta transport. Experiments of dilute mixtures of soil particles in supersonic gas flow are being conducted at Johns Hopkins University in order to provide a dataset describing gas-particle cloud interactions. Augmented by a separate funding source are efforts to perform small-scale particle-resolved direct numerical simulations (PR-DNS) and larger scale Eulerian-Lagrangian simulations. Together, these experiments and simulations are being used to devise models of particle cloud drag and particle cloud drag-induced dynamics. The final product of these efforts will be particle drag models suitable for implementation into the GGFS application and a PTKE model suitable for the same purpose.

Jeff West↗

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

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

To understand the behavior of charged lunar regolith when perturbed by lunar landers, it is important to couple the grain dynamics with mechanical and electrical particle interactions. To accomplish this, improvements have been made to a discrete element method (DEM) software package to include both short- and long-range interactions between spherical particles. Short-range interactions rely on contact between the particles, such as electrical conduction and triboelectric charge transfer driven by work functions. Long-range interactions act at a distance between every pairing of particles, such as electrostatic forces and gravitational forces. Results from simulations between a few particles are compared with theory to verify these added behaviors prior to scaling up to more complex scenarios. The radii, initial charges, electrical conductivities, work functions, and separation of the particles are varied and the resultant charges as well as the time required to reach the final state are determined.

Electrostatics↗

Models for dynamic surface temperatures during rarefied aeropass maneuvers

New models for calculating dynamic surface temperatures in a direct simulation Monte Carlo particle method during atmospheric entry are described. The models take into account the radiative emissivity, thermal conductivity, and heat capacity of the material. They are applied to thin panels in a highly rarefied flows and yield the transient thermal response of individual surface facets. Simulation results for radiation panels with internal conductivity are verified by comparing them with free molecule theory.

Haas, Brian L.↗

Enhanced Flight Dynamics Models with Aerodynamic Interference for Real-Time Simulation of VTOL Concept Vehicles

In an effort to establish common models to support the Advanced Air Mobility (AAM) community, NASA has designed a fleet of VTOL configuration reference models. Advanced Rotorcraft Technologies (ART) in cooperation with the NASA Ames Aeromechanics branch has developed FLIGHTLAB simulation models for several of these AAM VTOL concept vehicles. These simulation models are real-time capable while maintaining accurate flight dynamic characteristics with enhanced interference simulation by extracting important modeling parameters from Viscous Vortex Particle Method (VVPM). The focus of this paper is threefold. First, it introduces the simulation models developed and describes their modeling characteristics important to control and flight dynamics simulation. Second, the method for deriving the interference coefficients between model components such as rotors and wings is outlined. Implementation of the method for the lift+cruise and tiltwing models are discussed in detail. Third, the accuracy of the interference enhancement is assessed, with emphasis placed on the impacts to trim and flight dynamic characteristics. The interference enhancement method is further used to evaluate configuration design decisions by quantifying performance aspects of interference. This research effort culminated in simulation models for the AAM VTOL concept vehicles, with proper interference and real-time capability.

Flight Dynamics↗

The Effect of Stochastic Acceleration on the Dynamics of Solar Energetic Particles in the Heliosphere

Solar Energetic Particles (SEPs) are crucial in space weather phenomena, impacting space-based technologies and human activities in space. Understanding the dynamics of SEPs, including their acceleration and transport mechanisms, is essential for improving predictive models and mitigating adverse effects. This study investigates the impact of stochastic acceleration (SA) on the dynamics of SEPs within the heliosphere. By integrating stochastic acceleration into a SEP transport model, we assess its impact relative to diffusive shock acceleration (DSA). The investigation utilizes the Adaptive Mesh Particle Simulator (AMPS) within the Space Weather Modeling Framework (SWMF) to simulate SEP dynamics, incorporating solar wind parameters and turbulence from the Alfven Wave Solar Model (AWSoM). To do this, we solve the focused transport equation to model the propagation of SEPs along magnetic field lines, extending these lines to 5 AU to incorporate the potential effects of pitch angle scattering beyond 1 AU. This aspect could notably impact the decay phase of a SEP event. Modeling the transport of SEPs is coupled with the simulation of solar wind dynamics, the interplanetary magnetic field, and the characteristics of Alfven wave turbulence. This presentation covers the modeling approach used in this research, the characterization of the impact of the stochastic acceleration, and the role of pitch angle scattering at various heliocentric distances on the dynamics of the SEP events' decay phase.

SEP↗

A Global, Multi-Waveband Model for the Zodiacal Cloud

This recently completed three-year project was undertaken by the PI at the University of Florida, NASA Goddard and JPL, and by the Co-I and Collaborator at the University of Florida. The funding was used to support a continuation of research conducted at the University of Florida over the last decade which focuses on the dynamics of dust particles in the interplanetary environment. The main objectives of this proposal were: To produce improved dynamical models of the zodiacal cloud by performing numerical simulations of the orbital evolution of asteroidal and cometary dust particles. To provide visualizations of the results using our visualization software package, SIMUL, simulating the viewing geometries of IRAS and COBE and comparing the model results with archived data. To use the results to provide a more accurate model of the brightness distribution of the zodiacal cloud than existing empirical models. In addition, our dynamical approach can provide insight into fundamental properties of the cloud, including but not limited to the total mass and surface area of dust, the size-frequency distribution of dust, and the relative contributions of asteroidal and cometary material. The model can also be used to provide constraints on trace signals from other sources, such as dust associated with the "Plutinos" , objects captured in the 2:3 resonance with Neptune.

Grogan, Keith↗

Internal Flow Simulation of Enhanced Performance Solid Rocket Booster for the Space Transportation System

An enhanced performance solid rocket booster concept for the space shuttle system has been proposed. The concept booster will have strong commonality with the existing, proven, reliable four-segment Space Shuttle Reusable Solid Rocket Motors (RSRM) with individual component design (nozzle, insulator, etc.) optimized for a five-segment configuration. Increased performance is desirable to further enhance safety/reliability and/or increase payload capability. Performance increase will be achieved by adding a fifth propellant segment to the current four-segment booster and opening the throat to accommodate the increased mass flow while maintaining current pressure levels. One development concept under consideration is the static test of a "standard" RSRM with a fifth propellant segment inserted and appropriate minimum motor modifications. Feasibility studies are being conducted to assess the potential for any significant departure in component performance/loading from the well-characterized RSRM. An area of concern is the aft motor (submerged nozzle inlet, aft dome, etc.) where the altered internal flow resulting from the performance enhancing features (25% increase in mass flow rate, higher Mach numbers, modified subsonic nozzle contour) may result in increased component erosion and char. To assess this issue and to define the minimum design changes required to successfully static test a fifth segment RSRM engineering test motor, internal flow studies have been initiated. Internal aero-thermal environments were quantified in terms of conventional convective heating and discrete phase alumina particle impact/concentration and accretion calculations via Computational Fluid Dynamics (CFD) simulation. Two sets of comparative CFD simulations of the RSRM and the five-segment (IBM) concept motor were conducted with CFD commercial code FLUENT. The first simulation involved a two-dimensional axi-symmetric model of the full motor, initial grain RSRM. The second set of analyses included three-dimensional models of the RSRM and FSM aft motors with four-degree vectored nozzles.

Ahmad, Rashid A.↗

Building a better leapfrog

In stellar dynamical computer simulations, as well as other types of simulations using particles, time step size is often held constant in order to guarantee a high degree of energy conservation. In many applications, allowing the time step size to change in time can offer a great saving in computational cost, but variable-size time steps usually imply a substantial degradation in energy conservation. We present a meta-algorithm' for choosing time steps in such a way as to guarantee time symmetry in any integration scheme, thus allowing vastly improved energy conservation for orbital calculations with variable time steps. We apply the algorithm to the familiar leapfrog scheme, and generalize to higher order integration schemes, showing how the stability properties of the fixed-step leapfrog scheme can be extended to higher order, variable-step integrators such as the Hermite method. We illustrate the remarkable properties of these time-symmetric integrators for the case of a highly eccentric elliptical Kepler orbit and discuss applications to more complex problems.

Hut, Piet↗

High-Density Plasma Reactors: Simulations for Design

The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.

Hash, David B.↗

Equilibrium properties of a one-dimensional kinetic system.

One-dimensional systems of N = 500 and 250 particles in equilibrium are numerically simulated utilizing the method of molecular dynamics. Periodic boundary conditions are imposed. The classical two-body interaction potential is short range, repulsive and has a corresponding finite force. The equations of state are determined for densities both less and greater than one. Corresponding theoretical isochores are determined from models based on nearest-neighbor interactions and on a truncated virial expansion, and a comparison is made with the experimental isochores. Time independent radial distributions are constructed numerically and discussed. A change of state from a solidlike state to a fluid-gas state based on the penetrability of the particles is predicted. The transition temperatures are estimated from the radial distribution functions and the nearest-neighbor model. Self-diffusion is observed and the corresponding constants are determined from the velocity autocorrelation functions.

Williams, J. H.↗

Coordinated Microanalyses of Seven Particles of Probable Interstellar Origin from the Stardust Mission

Stardust, a NASA Discovery-class mission, was the first sample-return mission to return solid samples from beyond the Moon. Stardust was effectively two missions in one spacecraft: it returned the first materials from a known primitive solar system body, the Jupiter-family comet Wild 2; Stardust also returned a collector that was exposed to the contemporary interstellar dust stream for 200 days during the interplanetary cruise. Both collections present severe technical challenges in sample preparation and in analysis. By far the largest collection is the cometary one: approximately 300 micro g of material was returned from Wild 2, mostly consisting of approx. 1 ng particles embedded in aerogel or captured as residues in craters on aluminum foils. Because of their relatively large size, identification of the impacts of cometary particles in the collection media is straightforward. Reliable techniques have been developed for the extraction of these particles from aerogel. Coordinated analyses are also relatively straightforward, often beginning with synchrotron-based x-ray fluorescence (S-XRF), X-ray Absorption Near-Edge Spectoscopy (XANES) and x-ray diffraction (S-XRD) analyses of particles while still embedded in small extracted wedges of aerogel called ``keystones'', followed by ultramicrotomy and TEM, Scanning Transmission X-ray Microscopy (STXM) and ion microprobe analyses (e.g., Ogliore et al., 2010). Impacts in foils can be readily analyzed by SEM-EDX, and TEM analysis after FIB liftout sample preparation. In contrast, the interstellar dust collection is vastly more challenging. The sample size is approximately six orders of magnitude smaller in total mass. The largest particles are only a few pg in mass, of which there may be only approx.10 in the entire collection. The technical challenges, however, are matched by the scientific importance of the collection. We formed a consortium carry out the Stardust Interstellar Preliminary Examination (ISPE) to carry out an assessment of this collection, partly in order to characterize the collection in sufficient detail so that future investigators could make well-informed sample requests. The ISPE is the sixth PE on extraterrestrial collections carried out with NASA support. Some of the basic questions that we asked were: how many impacts are there in the collector, and what fraction of them have characteristics consistent with extraterrestrial materials? What is the elemental composition of the rock-forming elements? Is there crystalline material? Are there organics? Here we present coordinated microanalyses of particles captured in aerogel, using S-FTIR, S-XRF, STXM, S-XRD; and coordinated microanalyses of residues in aluminum foil, using SEMEDX, Auger spectroscopy, STEM, and ion microprobe. We discuss a novel approach that we employed for identification of tracks in aerogel, and new sample preparation techniques developed during the ISPE. We have identified seven particles - three in aerogel and four in foils - that are most consistent with an interstellar origin. The seven particles exhibit a large diversity in elemental composition. Dynamical evidence, supported supported by laboratory simulations of interstellar dust impacts in aerogel and foils, and numerical modeling of interstellar dust propagation in the heliosphere, suggests that at least some of the particles have high optical cross-section, perhaps due to an aggregate structure. However, the observations are most consistent with a variety of morphologies

Westphal, Andrew J.↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗

Dynamics of planetary rings

The modeling of the dynamics of particle collisions within planetary rings is discussed. Particles in the rings collide with one another because they have small random motions in addition to their orbital velocity. The orbital speed is roughly 10 km/s, while the random motions have an average speed of about a tenth of a millimeter per second. As a result, the particle collisions are very gentle. Numerical analysis and simulation of the ring dynamics, performed with the aid of a supercomputer, is outlined.

Araki, Suguru↗

Higher order moments of the matter distribution in scale-free cosmological simulations with large dynamic range

We calculate reduced moments (xi bar)(sub q) of the matter density fluctuations, up to order q = 5, from counts in cells produced by particle-mesh numerical simulations with scale-free Gaussian initial conditions. We use power-law spectra P(k) proportional to k(exp n) with indices n = -3, -2, -1, 0, 1. Due to the supposed absence of characteristic times or scales in our models, all quantities are expected to depend on a single scaling variable. For each model, the moments at all times can be expressed in terms of the variance (xi bar)(sub 2), alone. We look for agreement with the hierarchical scaling ansatz, according to which ((xi bar)(sub q)) proportional to ((xi bar)(sub 2))(exp (q - 1)). For n less than or equal to -2 models, we find strong deviations from the hierarchy, which are mostly due to the presence of boundary problems in the simulations. A small, residual signal of deviation from the hierarchical scaling is however also found in n greater than or equal to -1 models. The wide range of spectra considered and the large dynamic range, with careful checks of scaling and shot-noise effects, allows us to reliably detect evolution away from the perturbation theory result.

Lucchin, Francesco↗

High Performance Simulations of Accretion Disk Dynamics and Jet Formations Around Kerr Black Holes

We investigate jet formation in black-hole systems using 3-D General Relativistic Particle-In-Cell (GRPIC) and 3-D GRMHD simulations. GRPIC simulations, which allow charge separations in a collisionless plasma, do not need to invoke the frozen condition as in GRMHD simulations. 3-D GRPIC simulations show that jets are launched from Kerr black holes as in 3-D GRMHD simulations, but jet formation in the two cases may not be identical. Comparative study of black hole systems with GRPIC and GRMHD simulations with the inclusion of radiate transfer will further clarify the mechanisms that drive the evolution of disk-jet systems.

Nishikawa, Ken-Ichi↗

Numerical Model of the Mars Electrostatic Precipitator

NASA's future human exploration missions will require chemical processing plants to convert local resources into consumables to support astronaut activities. The thin and mostly carbon dioxide atmosphere of Mars is estimated to have 1 - 10 particles/cu.cm with diameters of 1 - 10 m and up to 1000 particles/cu.cm during storms. The dust in the Martian atmosphere can foul chemical reactors and pose a risk to life support systems. Electrostatic precipitation (ESP) removes dust particles from the Martian atmosphere. The Electrostatics and Surface Physics Laboratory at NASA's Kennedy Space Center has developed a COMSOL Multiphysics(Registered Trademark) model of an ESP for dust filtration on Mars. The fundamental principles of an ESP can be simulated by four physics modules: plasma, AC/DC electromagnetics, computational fluid dynamics (CFD), and particle tracing. In the ESP model presented here, the plasma module was solved to estimate particle charge. The AC/DC and CFD module were solved for the electrostatic force and fluid force. The particle-tracing module was solved for particle collection efficiency.

Wang, Jerry J.↗