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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

The dynamics of charged dust in the tail of Comet Giacobini-Zinner

A numerical simulation has been performed to calculate the distributions of dust particles of various sizes down the tail of comet P/Giacobini-Zinner. When the electrostatic charging of the grains in the plasma and radiative environment of the comet is taken into account, it is found that the distribution of the grains (particularly at the lowest end of the mass spectrum) in a plane normal to the orbital plane is entirely different from what is expected had the grains been uncharged. Although the NASA-ICE spacecraft, which will fly through the tail of this comet almost normal to its orbital plane on September 11, 1985, has no dedicated dust experiments, it is expected that the plasma wave instrument will serve as an indirect detector, even of the smallest grains, via the plasma clouds created by the high-velocity dust impacts. Knowledge of the spatial variation of the grain sizes encountered along the flight path of the spacecraft will provide us with the information necessary to calculate the electrostatic potential of the grains, which in turn will lead to an estimation of the role of the electromagnetic forces on the dynamics of such grains.

Horanyi, M.↗
NASA NTRS1999

Chain Dynamics in Magnetorheological Suspensions

Magnetorheological (MR) suspensions are composed of colloidal particles which acquire dipole moments when subjected to an external magnetic field. At sufficient field strengths and concentrations, the dipolar particles rapidly aggregate to form long chains. Subsequent lateral cross-linking of the dipolar chains is responsible for a rapid liquid-to-solid-like rheological transition. The unique, magnetically-activated rheological properties of MR suspensions make them ideal for interfacing mechanical systems to electronic controls. Additionally, the ability to experimentally probe colloidal suspensions interacting through tunable anisotropic potentials is of fundamental interest. Our current experimental work has focused on understanding the fluctuations of dipolar chains. It has been proposed by Halsey and Toor (HT) that the strong Landau-Peierls thermal fluctuations of dipolar chains could be responsible for long-range attractions between chains. Such interactions will govern the long-time relaxation of MR suspensions. We have synthesized monodisperse neutrally buoyant MR suspensions by density matching stabilized ferrofluid emulsion droplets with D2O. This allows us to probe the dynamics of the dipolar chains using light scattering without gravitational, interfacial, and polydispersity effects to resolve the short-wavelength dynamics of the dipolar chains. We used diffusing wave spectroscopy to measure these dynamics. The particle displacements at short times that show an independence to the field strength, but at long times exhibit a constrained, sub-diffusive motion that slows as the dipole strength is increased. The experiments are in good qualitative agreement with Brownian dynamics simulations of dipolar chains. Although there have been several important and detailed studies of the structure and interactions in MR suspensions, there has not been conclusive evidence that supports or contradicts the HT model prediction that long-range interactions exist between fluctuating chains of dipolar particles. Resolving this issue would contribute greatly to the understanding of these interesting and important materials. We have begun to test the predictions of the HT model by both examining the dynamics of individual chains and by measuring the forces between dipolar chains directly to accurately and quantitatively assess the interactions that they experience. To do so, we employ optical trapping techniques and video-microscopy to manipulate and observe our samples on the microscopic level. With these techniques, it is possible to observe chains that are fluctuating freely in three-dimensions, independent of interfacial effects. More importantly, we are able to controllably observe the interactions of two chains at various separations to measure the force-distance profile. The techniques also allow us to study the mechanical properties of individual chains and chain clusters. Our work to this point has focused on reversibly-formed dipolar chains due to field induced dipoles where the combination of this chaining, the dipolar forces, and the hydrodynamic interactions that dictate the rheology of the suspensions. One can envision, however, many situations where optical, electronic, or rheological behavior may be optimized with magneto-responsive anisotropic particles. Chains of polarizable particles may have the best properties as they can coil and flex in the absence of a field and stiffen and orient when a field is applied. We have recently demonstrated a synthesis of stable, permanent paramagnetic chains by both covalently and physically linking paramagnetic colloidal particles. The method employed allows us to create monodisperse chains of controlled length. We observed the stability, field-alignment, and rigidity of this new class of materials. The chains may exhibit unique rheological properties in an applied magnetic field over isotropic suspensions of paramagnetic particles. They are also useful rheological models as bead-spring systems. These chains form the basis for our current experiments with optical traps.

Gast, A. P.↗
NASA NTRS1978

Similar explosive eruptions of lunar and terrestrial volcanoes

A mathematical model of one-dimensional, steady duct flow of a mixture of gas and small solid particles (rock), has been applied to the earth and the moon under geometrically and dynamically similar conditions to simulate explosive eruption conditions. Numerical results for equilibrium two-phase flows for lunar and terrestrial explosive eruptions under similar conditions indicate that: (1) a lunar explosive vent is much larger than the corresponding terrestrial vent, (2) the exit velocity from the lunar explosive flow may be higher than the lunar escape velocity, but the exit velocity of terrestrial explosive flow is much less than that of the lunar case (this result supports the hypothesis that Australian tektites came from the moon as a stream of a mixture of rock and gas of extremely high speed), and (3) the thermal effects on the lunar explosive flows are much larger than those of the terrestrial case.

Pai, S. I.↗
NASA NTRS1995

A comparison of two algorithms for simulating collisionless systems

Two completely different simulation algorithms are compared by applying them to the same stellar dynamical problems: one is a self-consistent field (SCF) method for solving Poisson's equation and the other is a phase-space method for integrating the collisionless Boltzmann equation. We consider simulations of spherical stellar systems which are initially far from equilibrium and relax to their final states by gravitational collapse. The initial conditions consist of either uniform-density spheres or nonequilibrium models having Plummer density profiles, in which velocity dispersions are assigned according to given virial ratios. If a few tens of radial expansion terms with hundreds of thousands of particles are used in the SCF code, excellent agreement is found between the results it generates and those obtained with the phase-space solver, provided that a sufficiently large number of grid cells are employed with the latter. These findings imply that for simulating collisionless systmes over many dynamical times, the SCF approach based on sampling phase space is competitive with the approach treating phase space as a continuous fluid. The results of our tests make it possible to estimate the number of particles and basis functions required in situations like those modeled. Limitations of the SCF method and the choice of an optimal set of basis functions are also discussed.

Hozumi, Shunsuke↗
NASA NTRS2004

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

The global dynamics of the ionized and neutral components in the environment of Io plays an important role in the interaction of Jupiter's corotating magnetospheric plasma with Io. The stationary simulation of this problem was done in the MHD and the electrodynamics approaches. One of the main significant results from the simplified two-fluid model simulations was a production of the structure of the double-peak in the magnetic field signature of the I0 flyby that could not be explained by standard MHD models. In this paper, we develop a method of kinetic ion simulation. This method employs the fluid description for electrons and neutrals whereas for ions multilevel, drift-kinetic and particle, approaches are used. We also take into account charge-exchange and photoionization processes. Our model provides much more accurate description for ion dynamics and allows us to take into account the realistic anisotropic ion distribution that cannot be done in fluid simulations. The first results of such simulation of the dynamics of ions in the Io's environment are discussed in this paper.

Lipatov, Alexander S.↗
NASA NTRSDate not supplied

Semi-Empirical Interatomic Potential for Large-Scale Molecular Dynamics Simulation of Metal-Oxide Systems

Classical molecular dynamics (MD) simulation can be applied to systems containing billions of atoms during times up to microseconds. However, utilization of a large-scale MD simulation requires reliable but computationally cheap interatomic potentials. In the case of metallic systems, embedded atom method (EAM) and Finnis-Sinclair (FS) potentials are the reasonable choices, but their development for multicomponent alloys is a challenge. Addition of oxygen atoms makes the problem of developing semi-empirical potentials even more difficult because the long-range Coulomb interaction cannot be neglected, and different atoms have different charges which vary during MD relaxation such that a charge equilibration procedure should be applied every MD step. As results researchers tend to avoid simulating metal-oxide systems. In this talk, a new Finnis-Sinclair (FS) Ni-Co-Cr potential will be presented. Special attention will be paid to reproducing of the element melting temperatures and clustering. Next, a charge transfer ionic potential (CTIP) developed to simulate the NiO properties will be presented. This potential utilizes the Ni FS potential and accounts for the Coulomb interaction in the vicinity of oxygen atoms. Fitting of the CTIP parameters and FS potential functions describing the interaction with oxygen atoms to the available experimental and ab initio data will be discussed. Results of MD simulation of interaction of NiO particles with dislocations in Ni will be shown.

molecular dynamics simulation↗
NASA NTRS1994

A 'skewed' lognormal approximation to the probablility distribution function of the large-scale density field

I propose a method to fit the probability distribution function (PDF) of the large-scale density field rho, motivated by a Lagrangian version of the continuity equation. It consists in applying the Edgeworth expansion to the quantity Phi identical with log rho - mean value of log rho. The method is tested on the matter particle distribution in two cold dark matter N-body simulations of different physical sizes to cover a large dynamic range. It is seen to be very efficient, even in the nonlinear regime, and may thus be used as an analytical tool to study the effect on the PDF of the transition between the weakly nonlinear regime and the highly nonlinear regime.

Colombi, S.↗
NASA NTRS2011

Background and Pickup Ion Velocity Distribution Dynamics in Titan's Plasma Environment: 3D Hybrid Simulation and Comparison with CAPS T9 Observations

In this report we discuss the ion velocity distribution dynamics from the 3D hybrid simulation. In our model the background, pickup, and ionospheric ions are considered as a particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. The current simulation shows that mass loading by pickup ions H(+); H2(+), CH4(+) and N2(+) is stronger than in the previous simulations when O+ ions are introduced into the background plasma. In our hybrid simulations we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M = 28 amu ions that were generated inside the ionosphere. The moon is considered as a weakly conducting body. Special attention will be paid to comparing the simulated pickup ion velocity distribution with CAPS T9 observations. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of the Alfve n wing-like structures. The simulation also shows that the ring-like velocity distribution for pickup ions relaxes to a Maxwellian core and a shell-like halo.

Lipatov, A. S.↗
NASA NTRSDate not supplied

Impacts of Improved Process Representation of Particle Dry Deposition on PM Pollution in a Global Chemistry-Climate Model: Differences Across Regions, Seasons, and PM Sizes

Dry deposition (DD) is a primary removal pathway of particulate matter (PM). The aerosol DD schemes in most global models do not reflect current mechanistic understanding gleaned from observations. The NASA GISS global chemistry-climate model has a new and more dynamic DD scheme that largely captures observed changes in deposition velocities with particle size. We quantify the response of simulated PM to changes in the DD scheme for the Northeast US, Central Europe, North China Plain, and Punjab. Relative to the widely used old scheme, the new scheme shows higher annual PM2.5 for all regions (up to +14%) except C. Europe where there are very small decreases. For PM1, annual increases occur over all regions (up to +20%). For PM10, there are decreases in C. Europe (-8%) and very small decreases in the NE US yet increases (up to +12%) in Punjab and N. China Plain. While there are always increases across seasons for Punjab and N. China Plain, there are both seasonal increases and decreases for the NE US and C. Europe. Given incomplete understanding of observed variations in deposition velocities for a given particle size, we perform sensitivity simulations that perturb the magnitude of the deposition velocities simulated by the new scheme. The annual PM response to increasing DD is similar in magnitude to decreasing DD, implying linearity in the PM sensitivity to DD. The relative annual response to perturbing the DD magnitude is weaker over Punjab and sometimes N. China Plain than the NE US and C. Europe. Higher PM over Punjab and N. China Plain implies a stronger sensitivity to DD when aerosol abundances are low. More mechanistic representation of aerosol DD can sometimes improve or worsen existing model PM biases, which suggests that PM biases due to other processes can be confounded or compounded by biases in DD. Further improvements to DD parameterizations require not only more observational constraints on aerosol deposition velocities but also an advanced understanding of the processes controlling observed variability.

Environmental pollution↗
NASA NTRSDate not supplied

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

The development of a predictive simulation capability for extra-terrestrial Plume Surface Interaction (PSI) environments is undertaken by the Fluid Dynamics Branch at NASA/MSFC under the Game Changing Development (GCD) Program funded by NASA’s Space Technology Mission Directorate (STMD). Predictive simulation capability for propulsive Martian and Lunar landing PSI requires accurate modeling of the complex high-speed plume impingement and resulting gas/particle-cloud and ejecta stream formation. Physics modeling gap analysis during formulation of the PSI project scope identified two particular phenomena of first order importance in gas-particle interactions that lacked existing modeling formulations. The first arises from the lack of models capturing the unsteady drag forces imparted on particles by the rapidly expanding supersonic compressible flow. The second addresses the fluctuating forces and dispersions imparted on both the particle and the fluid resulting from the interference and wake turbulence generated by close proximity particle clouds, dubbed the Particle Turbulent Kinetic Energy(PTKE). Their first order significance has been identified in experiments, but simulation models currently do not exist for either effect. The development of models and the eventual implementation into the Eulerian Gas-Granular Flow Solver (GGFS) simulation tools was constructed as one of four major tasks of the PSI project. In this process, data on particle kinetics and gas-particle interactions are collected from carefully designed experiments of particles embedded in jets. The effects observed in the experiments are then replicated in high-fidelity particle-resolved CFD simulations to inform the formulation of improvements to particle phase drag models for implementations in the more efficient Eulerian-Lagrangian CFD simulations. The resulting models are ultimately ported to the Eulerian-Eulerian models applied for most efficient simulations in PSI production application tools. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling task of the PSI project.

Jeff West↗
NASA NTRS1972

A simulation of saltation dynamics with application to terrestrial and Martian model surfaces.

The feasibility of predicting the saltation threshold from a dynamic model of the saltation process is demonstrated. The model is based on mechanical equations that describe the trajectories of particles (sand grains) near the ground-atmosphere interface. The equations of particle motion are formulated in terms of forces and moments. The forces employed are based on theoretical and experimental results obtained for continuum flow. Surface structuring analytics and simulated fluid turbulence further modify the particle behavior. The model is used to predict the threshold velocities on Mars.

Andres, R. M.↗
NASA NTRS1992

Computer simulations of comet- and asteroidlike bodies passing through the Venusian atmosphere: Preliminary results on atmospheric and ground shock effects

We have completed computer simulations that model shock effects in the venusian atmosphere caused during the passage of two cometlike bodies 100 m and 1000 m in diameter and an asteroidlike body 10 km in diameter. Our objective is to examine hypervelocity-generated shock effects in the venusian atmosphere for bodies of different types and sizes in order to understand the following: (1) their deceleration and depth of penetration through the atmosphere; and (2) the onset of possible ground-surface shock effects such as splotches, craters, and ejecta formations. The three bodies were chosen to include both a range of general conditions applicable to Venus as well as three specific cases of current interest. These calculations use a new multiphase computer code (DICE-MAZ) designed by California Research & Technology for shock-dynamics simulations in complex environments. The code was tested and calibrated in large-scale explosion, cratering, and ejecta research. It treats a wide range of different multiphase conditions, including material types (vapor, melt, solid), particle-size distributions, and shock-induced dynamic changes in velocities, pressures, temperatures (internal energies), densities, and other related parameters, all of which were recorded in our calculations.

Roddy, D.↗
NASA NTRS1978

NASA Agricultural Aircraft Research Program in the Langley Vortex Research Facility and the Langley Full Scale Wind Tunnel

The current status of aerial applications technology research at the Langley's Vortex Research Facility and Full-Scale Wind Tunnel is reviewed. Efforts have been directed mainly toward developing and validating the required experimental and theoretical research tools. A capability to simulate aerial dispersal of materials from agricultural airplanes with small-scale airplane models, numerical methods, and dynamically scaled test particles was demonstrated. Tests on wake modification concepts have proved the feasibility of tailoring wake properties aerodynamically to produce favorable changes in deposition and to provide drift control. An aerodynamic evaluation of the Thrush Commander 800 agricultural airplane with various dispersal systems installed is described. A number of modifications intended to provide system improvement to both airplane and dispersal system are examined, and a technique for documenting near-field spray characteristics is evaluated.

Jordan, F. L., Jr.↗
NASA NTRS1984

A stochastic model for photon noise induced by charged particles in multiplier phototubes of the space telescope fine guidance sensors

The Space Telescope (ST) is subjected to charged particle strikes in its space environment. ST's onboard fine guidance sensors utilize multiplier phototubes (PMT) for attitude determination. These tubes, when subjected to charged particle strikes, generate spurious photons in the form of Cerenkov radiation and fluorescence which give rise to unwanted disturbances in the pointing of the telescope. A stochastic model for the number of these spurious photons which strike the photocathode of the multiplier phototube which in turn produce the unwanted photon noise are presented. The model is applicable to both galactic cosmic rays and charged particles trapped in the Earth's radiation belts. The model which was programmed allows for easy adaption to a wide range of particles and different parameters for the phototube of the multiplier. The probability density functions for photons noise caused by protons, alpha particles, and carbon nuclei were using thousands of simulated strikes. These distributions are used as part of an overall ST dynamics simulation. The sensitivity of the density function to changes in the window parameters was also investigated.

Howell, L. W.↗
NASA NTRS1986

Stochastic model for photon noise induced by charged particles in multiplier phototubes of the Hubble Space Telescope fine guidance sensors

The Space Telescope (ST) is subjected to charged particle strikes in its space environment. ST's onboard fine guidance sensors utilize multiplier phototubes (PMT) for attitude determination. These tubes, when subjected to charged particle strikes, generate spurious photons in the form of Cerenkov radiation and fluorescence which give rise to unwanted disturbances in the pointing of the telescope. A stochastic model for the number of these spurious photons which strike the photocathodes of the multiplier phototube which in turn produce the unwanted photon noise are presented. The model is applicable to both galactic cosmic rays and charged particles trapped in the earth's radiation belts. The model which was programmed allows for easy adaption to a wide range of particles and different parameters for the phototube of the multiplier. The probability density functions for photons noise caused by protons, alpha particles, and carbon nuclei were using thousands of simulated strikes. These distributions are used as part of an overall ST dynamics simulation. The sensitivity of the density function to changes in the window parameters was also investigated.

Howell, L. W.↗
NASA NTRS1994

Particle acceleration in the dynamic magnetotail: Orbits in self-consistent three-dimensional MHD fields

The acceleration of protons in a dynamically evolving magnetotail is investigated by tracing particles in the fields obtained from a three-dimensional resistive magnetohydrodynamic (MHD) simulation. The MHD simulation, representing plasmoid formation and ejection through a near-Earth reconnection process, leads to cross-tail electric fields of up to approximately 4 mV/m with integrated voltages across the tail of up to approximately 200 kV. Energization of particles takes place over a wide range along the tail, due to the large spatial extent of the increased electric field together with the finite cross-tail extent of the electric field region. Such accelerated particles appear earthward of the neutral line over a significant portion of the closed field line region inside of the separatrix, not just in the vicinity of the separatrix. Two different acceleration processes are identified: a 'quasi-potential' acceleration, due to particle motion in the direction of the cross-tail electric field, and a 'quasi-betatron' effect, which consists of multiple energy gains from repeated crossings of the acceleration region, mostly on Speiser-type orbits, in the spatially varying induced electric field. The major source region for accelerated particles in the hundreds of keV range is the central plasma sheet at the dawn flank outside the reconnection site. Since this source plasma is already hot and dense, its moderate energization by a factor of approximately 2 may be sufficient to explain the observed increases in the energetic particle fluxes. Particles from the tail are the source of beams at the plasma sheet/lobe boundary. The temporal increase in the energetic particle fluxes, estimated from the increase in energy gain, occurs on a fast timescale of a few minutes, coincident with a strong increase in B(sub z), despite the fact that the inner boundary ('injection boundary') of the distribution of energized particles is fairly smooth.

Birn, Joachim↗
NASA NTRS2004

Dynamics of Aerosol Particles in Stationary, Isotropic Turbulence

A detailed study of the dynamics of sub-Kolmogorov-size aerosol particles in stationary isotropic turbulence has been performed. The study combined direct numerical simulations (DNS; directed by Prof. Collins) and high-resolution experimental measurements (directed by Prof. Meng) under conditions of nearly perfect geometric and parametric overlap. The goal was to measure the accumulation of particles in low-vorticity regions of the flow that arises from the effect commonly referred to as preferential concentration. The grant technically was initiated on June 13, 2000; however, funding was not available until July 11, 2000. The grant was originally awarded to Penn State University (numerical simulations) and SUNY-Buffalo (experiments); however, Prof. Collins effort was moved to Cornell University on January 2002 when he joined that university. He completed the study there. A list of the specific tasks that were completed under this study is presented.

Collins, Lance R.↗
NASA NTRSDate not supplied

Advancing xEMU Lunar Dust Mitigation Devices

Fine, electrically charged, glass like dust particles caused significant damage to the Apollo EMU9 during lunar EVAs, identified as one of the greatest challenges to future exploration. Passive Lunar Dust Mitigation Devices (LDMD) were developed, within a SBIR Phase II, to prohibit this dust from interrupting venting space suit component operation. A Computational Fluid-Dynamics and Discrete Element Method (CFD-DEM) Simulation Tool was developed at the University of Colorado, Boulder to predict venting gas flow ability to self-clean adhered dust particles from LDMD surfaces. Lunar dust properties (i.e., adhesion and cohesion strengths) required to complete Simulation Tool analysis are relatively unknown due to considerable differences between the Earth and the Moon (i.e., gravity, humidity) and due to an absence of dust particles in their native state. Analytically determining gas flow velocity, density and direction within the fluid Boundary Layer, microns from LDMD surfaces presented a second challenge. Dusty Plasma Chamber testing is being performed at Auburn University, Auburn to observe electrostatically charged dust behavior as it adheres to LDMD prototypes and specific geometry and is then blown away by metered gas flow. Test articles were developed to offer insight into the impact of different flow geometries, surface roughness and dust removal within the gas Boundary-Layer. Observed dust behavior is currently being developed to support the CFD-DEM analysis. Many Simulation Tool analytical cases have been processed to support the intention of completing sensitivity studies to assess how different dust adherence values and Boundary Layer fluid properties impact LDMD self-cleaning effectivity.

Thomas J Stapleton↗