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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

The CN molecule is an important contributor to radiative heat flux in shock layers around vehicles entering Titan’s atmosphere. Current data for heavy particle (de)excitation rate coefficients of CN leads to uncertainties in the population of CN in its first and second excited states. This in turn leads to uncertainties in the radiative heat flux predicted by Computational Fluid Dynamics (CFD) simulations of Titan atmospheric entry. This work performs ab initio electronic structure calculations of the CNN complex to create Potential Energy Surfaces (PESs) that correlate to the ground and first and second excited states of CN. Specifically, the state combinations of CN(X,A,B) + N(4S𝑜) correlate to six states of CNN (three Quintet A” and three Triplet A”). Initial calculations of these states suggest that heavy particle (de)excitation of CN by N atoms is likely to proceed through collinear geometries on triplet surfaces. Complete PESs will show all of the reaction pathways in detail, and will be used in nonadiabatic dynamics calculations to evaluate improved rate coefficients and reduce uncertainty in the radiative heat flux during Titan entry.

Eric C Geistfeld↗

NASA NTRSDate not supplied

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

Eric Geistfeld↗

NASA NTRS1993

Polar wind ion dynamics in the magnetotail

The circulation of polar wind ions from the high-latitude ionosphere to the plasma sheet is investigated. Single-particle trajectory simulations for the geomagnetic tail show, in addition to the breaking of adiabaticity, the possible creation of new high-altitude mirror points. This trajectory feature results from an intense parallel deceleration imparted by the magnetic field rotation during fast ExB transport. This centrifugal deceleration yields a critical parallel velocity which depends on the magnitude of the convection electric field and below which ions remain trapped inside the neutral sheet.

Delcourt, D. C.↗

NASA NTRS1994

The effect of wave-particle interactions on the polar winds O(+)

The escape of the polar wind plasma is an important element in the ionosphere-magnetosphere coupling. Both theory and observations indicate that the wave-particle interactions (WPI) play a significant role in the dynamics of ion outflow along open geomagnetic field lines. A Monte Carlo simulation was developed in order to include the effect of the WPI in addition to the factors that are traditionally included in the 'classical' polar wind (i.e. gravity, electrostatic field, and divergence of geomagnetic field lines). The ion distribution function (f(sub j)), as well as the profiles of its moments (density, drift velocity, temperature, etc.) were found for different levels of WPI, that is, for different values of the normalized diffusion rate in the velocity space (D(tilde) (sub j perpendicular). Although the model included O(=), H(+) and eletrons, we presented only the results related to the O(+) ion. We found that (1) both the density and drift velocity of O(+) increased with the WPI strength, and consequently, the O(+) escape flux was enhanced by a factor of up to 10(exp 5), (2) The O(+) ions could be energized up to a few electron volts; (3) for moderate and high levels of WPI D(tilde) (sub perpendicular) (O(+) greater than (tilde) 1, the distribution function f(O(+)) displayed very pronounced conic features at altitudes around 3 R(sub e). Finally, the interplay between the downward body force, the upward mirror force, and the perpendicular heating resulted in the formation of the 'pressure cooker' effect. This phenomena explained some interesting features of our solution, such as, the peak in the O(=) temperature, and the formation of 'ears' and conics for f(O(+)) around 2.5 R(sub e).

Barakat, A. R.↗

NASA NTRS2008

Comparison of ALE and SPH Simulations of Vertical Drop Tests of a Composite Fuselage Section into Water

Simulation of multi-terrain impact has been identified as an important research area for improved prediction of rotorcraft crashworthiness within the NASA Subsonic Rotary Wing Aeronautics Program on Rotorcraft Crashworthiness. As part of this effort, two vertical drop tests were conducted of a 5-ft-diameter composite fuselage section into water. For the first test, the fuselage section was impacted in a baseline configuration without energy absorbers. For the second test, the fuselage section was retrofitted with a composite honeycomb energy absorber. Both tests were conducted at a nominal velocity of 25-ft/s. A detailed finite element model was developed to represent each test article and water impact was simulated using both Arbitrary Lagrangian Eulerian (ALE) and Smooth Particle Hydrodynamics (SPH) approaches in LS-DYNA, a nonlinear, explicit transient dynamic finite element code. Analytical predictions were correlated with experimental data for both test configurations. In addition, studies were performed to evaluate the influence of mesh density on test-analysis correlation.

Jackson, Karen E.↗

NASA NTRS1983

Cloud-particle galactic gas dynamics and star formation

Galactic gas dynamics, spiral structure, and star formation are discussed in relation to N-body computational studies based on a cloud-particle model of the interstellar medium. On the small scale, the interstellar medium is seen as cloud-dominated and supernova-perturbed. It is noted that the cloud-particle model simulates cloud-cloud collisions, the formation of stellar associations, and supernova explosions as dominant local processes. On the large scale, in response to a spiral galactic gravitational field, global density waves and galactic shocks develop having large-scale characteristics similar to those found in continuum gas dynamical studies. Both the system of gas clouds and the system of young stellar associations forming from the clouds figure in the global spiral structure. However, with the attributes of neither assuming a continuum of gas (as in continuum gas dynamical studies) or requiring a prescribed equation of state (such as the isothermal condition), the cloud-particle picture retains much of the detail lost in earlier work. By detail is meant the small-scale features and structures so important in understanding the local, turbulent state of the interstellar medium as well as the degree of raggedness often seen to be superposed on the global spiral structure.

Roberts, W. W., Jr.↗

NASA NTRS1995

Advanced in Visualization of 3D Time-Dependent CFD Solutions

Numerical simulations of complex 3D time-dependent (unsteady) flows are becoming increasingly feasible because of the progress in computing systems. Unfortunately, many existing flow visualization systems were developed for time-independent (steady) solutions and do not adequately depict solutions from unsteady flow simulations. Furthermore, most systems only handle one time step of the solutions individually and do not consider the time-dependent nature of the solutions. For example, instantaneous streamlines are computed by tracking the particles using one time step of the solution. However, for streaklines and timelines, particles need to be tracked through all time steps. Streaklines can reveal quite different information about the flow than those revealed by instantaneous streamlines. Comparisons of instantaneous streamlines with dynamic streaklines are shown. For a complex 3D flow simulation, it is common to generate a grid system with several millions of grid points and to have tens of thousands of time steps. The disk requirement for storing the flow data can easily be tens of gigabytes. Visualizing solutions of this magnitude is a challenging problem with today's computer hardware technology. Even interactive visualization of one time step of the flow data can be a problem for some existing flow visualization systems because of the size of the grid. Current approaches for visualizing complex 3D time-dependent CFD solutions are described. The flow visualization system developed at NASA Ames Research Center to compute time-dependent particle traces from unsteady CFD solutions is described. The system computes particle traces (streaklines) by integrating through the time steps. This system has been used by several NASA scientists to visualize their CFD time-dependent solutions. The flow visualization capabilities of this system are described, and visualization results are shown.

Lane, David A.↗

NASA NTRSDate not supplied

Framework for Analyzing the Complex Interactions Between Spacecraft Motion and Slosh Dynamics in Low-G Environments

The fuel-to-dry-mass ratio of spacecraft continues to grow as new human spaceflight missions target destinations farther from Earth. Large amounts of liquid propellant can lead to significant coupling between the rigid-body dynamics of the spacecraft and the motion of the fuel within its storage tank. The present work gives an overview of the dynamic features and a flowchart for a method of simulating the motion of a spacecraft with fuel slosh inside a cylindrical, domed tank in a low-g environment. The method involves modeling the liquid propellant as a particle that transfers momentum to the spacecraft through perfectly inelastic collisions with the tank wall. The foundation of the modeling methodology is the approach taken during the Apollo program used to predict the effect of fuel slosh on the complex motion exhibited by the Service Module following separation from the Command Module. This paper discusses the motivation, methodology, and conclusions from the Apollo-era method, then presents corrections to the derivation of the dynamics and fills in the gaps left from the unavailability of the detailed contractor report and simulation code. The results presented in this paper provide an example that demonstrates the effect that fuel slosh can have on the trajectory of a spacecraft in a low-g environment.

William J Elke III↗

NASA NTRS1985

Simulation study of Type 2 counterstreaming electrons along auroral field lines

The production of counterstreaming electrons associated with parallel fields along auroral field lines is examined through the use of computer simulation. A 2 1/2-dimensional (two spatial and three velocity dimensions) electrostatic particle algorithm and auroral boundary conditions are used to set up a self-consistent V potential structure. The simulation produces signatures of counterstreaming electrons resembling those observed by the Dynamics Explorer 1 satellite. The main signatures are as follows: (1) the phase space contours of the electron distribution function are elongated along the V-parallel axis, and (2) the energy of electrons streaming in the upward direction is comparable to the energy of the accelerated electron beam. The simulation indicates that a portion of the accelerated electron beam is trapped by large amplitude electrostatic waves produced through the two-stream instability. Strong wave-particle interactions then thermalize the trapped electrons to produce suprathermal electrons streaming in the direction opposite to that of the accelerated electron beam. These results suggest a possible mechanism of producing counterstreaming electron fluxes through nonlinear processes of the two-stream instability.

Wagner, J. S.↗

NASA NTRS1994

Dynamics of binary-disk interaction. 1: Resonances and disk gap sizes

We investigate the gravitational interaction of a generally eccentric binary star system with circumbinary and circumstellar gaseous disks. The disks are assumed to be coplanar with the binary, geometrically thin, and primarily governed by gas pressure and (turbulent) viscosity but not self-gravity. Both ordinary and eccentric Lindblad resonances are primarily responsible for truncating the disks in binaries with arbitrary eccentricity and nonextreme mass ratio. Starting from a smooth disk configuration, after the gravitational field of the binary truncates the disk on the dynamical timescale, a quasi-equilibrium is achieved, in which the resonant and viscous torques balance each other and any changes in the structure of the disk (e.g., due to global viscous evolution) occur slowly, preserving the average size of the gap. We analytically compute the approximate sizes of disks (or disk gaps) as a function of binary mass ratio and eccentricity in this quasi-equilibrium. Comparing the gap sizes with results of direct simulations using the smoothed particle hydrodynamics (SPH), we obtain a good agreement. As a by-product of the computations, we verify that standard SPH codes can adequately represent the dynamics of disks with moderate viscosity, Reynolds number R approximately 10(exp 3). For typical viscous disk parameters, and with a denoting the binary semimajor axis, the inner edge location of a circumbinary disk varies from 1.8a to 2.6a with binary eccentricity increasing from 0 to 0.25. For eccentricities 0 less than e less than 0.75, the minimum separation between a component star and the circumbinary disk inner edge is greater than a. Our calculations are relevant, among others, to protobinary stars and the recently discovered T Tau pre-main-sequence binaries. We briefly examine the case of a pre-main-sequence spectroscopic binary GW Ori and conclude that circumbinary disk truncation to the size required by one proposed spectroscopic model cannot be due to Linblad resonances, even if the disk is nonviscous.

Artymowicz, Pawel↗

NASA NTRSDate not supplied

Framework for Analyzing the Complex Interactions Between Spacecraft Motion and Slosh Dynamics in Low-G Environments

The fuel-to-dry-mass ratio of spacecraft continues to grow as new human spaceflight missions target destinations from Earth. Large amounts of liquid propellant can lead to significant coupling between the rigid-body dynamics of the spacecraft and the motion of the fuel within its storage tank. The present work gives an overview of the dynamic features and a flowchart for a method of simulating the motion of a spacecraft with fuel slosh inside a cylindrical, domed tank in a low-g environment. The method involves modeling the liquid propellant as a particle that transfers momentum to the spacecraft through perfectly inelastic collisions with the tank wall. The foundation of the modeling methodology is the approach taken during the Apollo program used to predict the effect of fuel slosh on the complex motion exhibited by the Service Module following separation from the Crew Module. This paper discusses the motivation, methodology, and conclusions from the Apollo-era method, and subsequently builds upon it by incorporating corrections to the derivation of the dynamics and filling in the gaps due to unavailability of the detailed contractor report and simulation code. The results presented in this paper provide an example that demonstrates the effect that fuel slosh can have on the trajectory of a spacecraft in a low-g environment.

William J Elke III↗

NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗

NASA NTRS2009

Shock Structure and Magnetic Fields Generation Associated with Relativistic Jets Unmagnetized Pair Plasma

Using 3D and 2D particle-in-cell simulations we investigate a shock structure, magnetic field generation, and particle acceleration associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized pair plasma. The simulations use long computational grids which allow to study the formation and dynamics of the system in a spatial and temporal way. We find for the first time a relativistic shock system comparable to a predicted magnetohydrodynamic shock structure consisting of leading and trailing shocks separated by a contact discontinuity. Strong electromagnetic fields resulting from the Weibel two-stream instability are generated in the trailing shock where jet matter is thermalized and decelerated. We analyze the formation and nonlinear development through saturation and dissipation of those fields and associated particle acceleration. In the AGN context the trailing shock corresponds to the jet shock at the head of a relativistic astrophysical jet. In the GRB context this trailing shock can be identified with the bow shock driven by relativistic ejecta. The strong electromagnetic field region in the trailing shock provides the emission site for the hot spot at the leading edge of AGN jets and for afterglow emission from GRBs.

Niemiec, J.↗

NASA NTRSDate not supplied

Computational Fluid Dynamics Models of the NASA White Sands Test Facility Particle Impact Test Nozzles

The particle impact ignition mechanism is considered to be the most effective means for igniting metals in gaseous oxygen environments and is a concern for every oxygen system. For this reason, the National Aeronautics and Space Administration (NASA) commits a large amount of resources attempting to understand and defend against this phenomenon. The most practical means of accomplishing this is through extensive of testing. To this end, over the years NASA has developed and built a facility and various apparatuses to test material or component vulnerability to a particle impact ignition. In testing materials, NASA uses sub-sonic and supersonic nozzles to replicate the conditions that a material is expected to experience in service. However, challenges include obtaining an understanding of the gas dynamics and of the particle behavior through and up to the point of impact at the material under test. One means of addressing these challenges is the use of computational fluid dynamics (CFD). If a CFD model can be built and validated, the model can then be used to confidently interrogate changes to inlet conditions and particle geometries, as well as changes to nozzle geometries to better simulate a material’s in-situ conditions and greatly reduce or eliminate trial-and-error approaches in testing. This paper will show NASA’s approach to and results of using CFD to better understand their particle impact nozzles.

Computational Fluid Dynamics↗

NASA NTRS2000

Interactive Soil Dust Aerosol Model in the GISS GCM: Sensitivity of the Soil Dust Cycle to Radiative Properties of Soil Dust Aerosols - Part 1

The sensitivity of the soil dust aerosol cycle to the radiative forcing by soil dust aerosols is studied. Four experiments with the NASA/GISS atmospheric general circulation model, which includes a soil dust aerosol model, are compared, all using a prescribed climatological sea surface temperature as lower boundary condition. In one experiment, dust is included as dynamic tracer only (without interacting with radiation), whereas dust interacts with radiation in the other simulations. Although the single scattering albedo of dust particles is prescribed to be globally uniform in the experiments with radiatively active dust, a different single scattering albedo is used in those experiments to estimate whether regional variations in dust optical properties, corresponding to variations in mineralogical composition among different source regions, are important for the soil dust cycle and the climate state. On a global scale, the radiative forcing by dust generally causes a reduction in the atmospheric dust load corresponding to a decreased dust source flux. That is, there is a negative feedback in the climate system due to the radiative effect of dust. The dust source flux and its changes were analyzed in more detail for the main dust source regions. This analysis shows that the reduction varies both with the season and with the single scattering albedo of the dust particles. By examining the correlation with the surface wind, it was found that the dust emission from the Saharan/Sahelian source region and from the Arabian peninsula, along with the sensitivity of the emission to the single scattering albedo of dust particles, are related to large scale circulation patterns, in particular to the trade winds during Northern Hemisphere winter and to the Indian monsoon circulation during summer. In the other regions, such relations to the large scale circulation were not found. There, the dependence of dust deflation to radiative forcing by dust particles is probably dominated by physical processes with short time scales. The experiments show that dust radiative forcing can lead to significant changes both in the soil dust cycle and in the climate state. To estimate dust concentration and radiative forcing by dust more accurately, dust size distributions and dust single scattering albedo in the model should be a function of the source region, because dust concentration and climate response to dust radiative forcing are sensitive to dust radiative parameters.

Perlwitz, Jan↗

NASA NTRS2013

Energetic Particles Dynamics in Mercury's Magnetosphere

We investigate the drift paths of energetic particles in Mercury's magnetosphere by tracing their motion through a model magnetic field. Test particle simulations solving the full Lorentz force show a quasi-trapped energetic particle population that gradient and curvature drift around the planet via "Shabansky" orbits, passing though high latitudes in the compressed dayside by equatorial latitudes on the nightside. Due to their large gyroradii, energetic H+ and Na+ ions will typically collide with the planet or the magnetopause and will not be able to complete a full drift orbit. These simulations provide direct comparison for recent spacecraft measurements from MESSENGER. Mercury's offset dipole results in an asymmetric loss cone and therefore an asymmetry in particle precipitation with more particles precipitating in the southern hemisphere. Since the planet lacks an atmosphere, precipitating particles will collide directly with the surface of the planet. The incident charged particles can kick up neutrals from the surface and have implications for the formation of the exosphere and weathering of the surface

Energetic particles↗

NASA NTRS1990

Simulations of gas clouds in interacting galaxies

A companion can induce a variety of morphological changes in a galaxy. The author uses N-body simulations to study the effects of different kinds of perturbations on the dynamics of a disk galaxy. The model is two-dimensional, with a disk consisting of about 60,000 particles. Most of the particles (80%) represent the old stellar population with a high velocity dispersion, while the rest (20%) represent gas clouds with a low velocity dispersion. Initially, the velocity dispersion corresponds to Q = 1 for the star particles, and Q = O for the gas particles, where Q is Toomre's (1964) stability parameter. The gas clouds can collide inelastically. The disk is stabilized by a rigid halo potential, and by the random motions of the old star particles. To simulate the effect of an encounter on the disk, a companion galaxy, modelled as a point mass, can move in a co-planar orbit around the disk. A complete description of the N-body code is found in Thomasson (1989). The spiral structures caused by a companion in first a direct and then a retrograde (with respect to the rotation of the disk) parabolic orbit are presented. The associated velocity fields suggest a way to observationally distinguish between leading and trailing spiral arms. The stability of the gas component in a disk in which tidally triggered infall of gas to the center occurs is studied. Finally, the author shows how a ring of gas can form in a disk as a result of a co-planar encounter with another galaxy.

Thomasson, Magnus↗

NASA NTRS2013

A Particle-in-Cell Simulation for the Traveling Wave Direct Energy Converter (TWDEC) for Fusion Propulsion

A Particle-in-cell simulation model has been developed to study the physics of the Traveling Wave Direct Energy Converter (TWDEC) applied to the conversion of charged fusion products into electricity. In this model the availability of a beam of collimated fusion products is assumed; the simulation is focused on the conversion of the beam kinetic energy into alternating current (AC) electric power. The model is electrostatic, as the electro-dynamics of the relatively slow ions can be treated in the quasistatic approximation. A two-dimensional, axisymmetric (radial-axial coordinates) geometry is considered. Ion beam particles are injected on one end and travel along the axis through ring-shaped electrodes with externally applied time-varying voltages, thus modulating the beam by forming a sinusoidal pattern in the beam density. Further downstream, the modulated beam passes through another set of ring electrodes, now electrically oating. The modulated beam induces a time alternating potential di erence between adjacent electrodes. Power can be drawn from the electrodes by connecting a resistive load. As energy is dissipated in the load, a corresponding drop in beam energy is measured. The simulation encapsulates the TWDEC process by reproducing the time-dependent transfer of energy and the particle deceleration due to the electric eld phase time variations.

Chap, Andrew↗