Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle methods”

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 199 records · Page 11

Particle kinetic simulation of high altitude hypervelocity flight

A new model for simulating the transfer of energy between the translational and rotational modes was derived for a homogeneous gas of diatomic molecules. The model was developed specifically for use in discrete particle simulation methods where molecular motion and inter-molecular collisions are treated at the molecular level. In such methods it is normal to assume a constant rotational collision number for the entire flowfield. The new model differs in that a temperature dependence is introduced which has been predicted by theory and observed in experiment. The new model is applied to the relaxation of rotational temperature, and is found to produce significant differences in comparison with the model normally employed at both high and low temperatures. Calculations have also been performed for a Mach 7 normal shock wave. Large differences in the solutions are again observed, with the new model offering an improved correspondence to the available experimental data.

Boyd, Iain D.↗

Rotational-translational energy transfer in rarefied nonequilibrium flows

A new model for simulating the transfer of energy between the translational and rotational modes is derived for a homogeneous gas of diatomic molecules. The model has been developed specifically for use in discrete particle simulation methods where molecular motion and intermolecular collisions are treated at the molecular level. A temperature dependence is introduced which has been predicted by theory and observed in experiment. The new model is applied to the relaxation of rotational temperature, and is found to produce significant differences in comparison with the model normally employed at both high and low temperatures. Calculations have also been performed for a Mach 7 normal shock wave.

Boyd, Iain D.↗

Fundamentals of chemistry modeling applicable to a vectorized particle simulation

This paper describes the fundamentals of extending the vectorized particle simulation method derived by Baganoff and McDonald (1990), McDonald and Baganoff (1988), and McDonald (1989) for modeling chemically reacting flows. Details of reaction mechanics per reaction are presented, with particular attention given to the quantum nature of the vibrational mode. The models of reactive flows developed here were verified through a simulation of a superheated diatomic gas relaxing thermochemically to equilibrium in a reservoir.

Haas, Brian L.↗

Three dimensional Discrete Particle Simulation about the AFE geometry

The Discrete Particle Simulation method, due to Baganoff, has recently been extended to allow representation of gases composed of multiple species, to general power-law molecular interactions and to permit flows in thermal non-equilibrium. Particular attention has been paid to the implementation of this physics while retaining the efficiency of the original algorithm. Here, the enhanced algorithm is applied to the simulation of the flow field about the Aeroassisted Flight Experiment (AFE) vehicle with the same flight parameters as in a previous paper. The enhancements to the algorithm are introduced and comparisons are made to the previous calculation.

Feiereisen, William J.↗

A simulation of the collapse and fragmentation of cooling molecular clouds

The application of the Smoothed Particle Hydrodynamics method to the fragmentation of rotating cloud and disk systems is described, allowing for molecular cooling due to H2 and CO. A novel approach to solving Poisson's equation for disklike structures which exploits the multigrid algorithm is also described. Numerical studies are presented which investigate the evolution of both rotating clouds and Maclaurin disks, in each case with both an isothermal equation of state and with molecular cooling. The results establish the influence of molecular cooling on the fragmentation of molecular clouds. The isothermal sequences, if they fragment at all, do so into far fewer lumps than the cooling sequences. This is not due to a cooling instability as such, but rather to the reduced thermal support. One of the sequences shows a remarkable similarity to the W49A star-forming region.

Monaghan, Joe J.↗

Particle simulation in a multiprocessor environment

A parallel implementation of a particle simulation method that is portable between a wide class of multiprocessor computers is presented. A fine grain spatial decomposition is utilized where several subdomains having a regular structure are computed at each processing node. This leads directly to an efficient and straightforward load balancing scheme if the number of subdomains at each processor is permitted to vary in an appropriate manner. Three dimensional simulations incorporating full thermochemical nonequilibrium are possible using the resulting code. Vectorizable algorithms are retained from earlier work allowing efficient use of deeply pipelined node processors where available. Performance results are presented from three different machine architectures demonstrating the portability of the code. On a 128-node Intel iPSC/860, performance is twice that of a single Cray-Y/MP CPU running a highly vectorized simulation code. Speedup is linear over the full range of number of processors on all target machines, indicating scalability of the method to higher degrees of parallelism.

Mcdonald, Jeffrey D.↗

The origin of the moon and the single impact hypothesis. IV

The smoothed-particle hydrodynamics method developed to simulate the conditions under which a major planetary collision may have led to the formation of the moon is presently used in 41 simulations exploring the relevant parameter space. Two possible scenarios emerge: one in which the moon is formed directly in the collision or only after dissipation and evolution of the disk. Attention is given to the Giant Blowoff hypothesis, in which not only the primitive atmosphere but large amounts of mantle material flow over the surface of the orbiting disk and blow off in the forward direction. Large amounts of angular momentum can also be lost.

Cameron, A. G. W.↗

Particle simulation of rarefied aeropass maneuvers of the Magellan spacecraft

Rarefied and free molecular flowfields corresponding to proposed aeropass maneuvers through the atmosphere of Venus by the Magellan spacecraft are computed with a vectorized particle simulation method. The significance of rarefaction and reaction effects are assessed and surface heat flux, drag coefficients, and flowfield properties are computed. A simple surface heat transfer model is coupled directly to the simulation and assumes that each surface element is in radiative equilibrium with deep space. This allows direct computation of surface temperature distributions rather than requiring prescribed isothermal boundary conditions. Uncoupled dual-node heat transfer models, which account for heat capacity and thermal conductivity within each structural component, permit more accurate determination of surface temperatures. Such temperatures restrict the allowable altitudes and entry speeds of the aeropass maneuvers, and must therefore be estimated accurately. Simulations with an entry velocity of 8600 m/s at altitudes between 125 and 140 km reveal that allowable surface temperatures occur only at the highest altitudes where the flowfield is appropriately modeled as free-molecular. Excessive surface heat flux occurs at lower altitudes where molecular collisions are significant.

Haas, Brian L.↗

High voltage spheres in an unmagnetized plasma - Fluid and PIC simulations

The basic physics involving the interaction of a high voltage sphere with a LEO-Plasma was investigated via computer experiments. Two approaches were used: the fluid description and the more general particle pushing method. Both techniques displayed qualitatively similar features. Thus, the initial time-dependent response of the plasma predicted by the fluid model, including the initial current surge and sheath formation, the formation of ion and electron density shells about the sphere with a double layer potential structure, and the subsequent propagation of the density shells away from the sphere. The PIC results also displayed an oscillatory character and slightly different time constants for the growth and decay mechanisms involved.

Thiemann, H.↗

Computational fluid dynamics combustion analysis evaluation

This study involves the development of numerical modelling in spray combustion. These modelling efforts are mainly motivated to improve the computational efficiency in the stochastic particle tracking method as well as to incorporate the physical submodels of turbulence, combustion, vaporization, and dense spray effects. The present mathematical formulation and numerical methodologies can be casted in any time-marching pressure correction methodologies (PCM) such as FDNS code and MAST code. A sequence of validation cases involving steady burning sprays and transient evaporating sprays will be included.

Kim, Y. M.↗

Nonaxisymmetric evolution in protostellar disks

We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.

Laughlin, Gregory↗

Beta Absorption Mass Monitoring of Particulates - A Review

The theory and application of beta-radiation absorption for the measurement and monitoring of airborne particulates are discussed. The use of this technique, both for source testing and for ambient air quality monitoring is reviewed. Various particle collection methods used in conjunction with beta absorption sensing configurations are considered. State of the art and current developments of instrumentation approaches for the automated measurement of mass concentration and size distribution of aerosols by beta absorption are discussed. Methods for electronic signal processing and recording are presented. The Beta absorption technique appears as a powerful tool for the unattended measurement of the mass of particulate pollution, compatible with telemetry and central data processing methods.

Lilienfeld, Pedro↗

Interference of Backscatter from Two Droplets in a Focused Continuous Wave CO2 Doppler Lidar Beam

Interference due to the superposition of backscatter, beta, from two micron-sized droplets using a NASA/MSFC cw (continuous wave) CO2 Doppler lidar at 9.1 micrometer wavelength was detected for the first time. The resultant single beta signal from both droplets contained an interference structure with a well-defined periodicity which was accurately measured and compared with cw lidar theory. The agreement between measurements and theory is excellent, indicating that the interference arises because the droplets are moving at different speeds and, therefore, the relative droplet separation is not constant. This gives the superimposed beta signal from both droplets in the lidar beam moving in and out of constructive and destructive interference with a well-defined periodic structure. A measurement of a time-resolved signal pulse with an oscilloscope of the combined from two approx. 14.12 micrometers in diameter silicone oil droplets responding to the Gaussian lidar beam intensity at the lidar beam focus is shown. Full details of this laboratory experimental setup, particle generation method, measurement technique, and the cw lidar can be found elsewhere. The stream of silicone oil droplets resided at a Doppler-shift center frequency of f (sub D) approx. (3.4 plus or minus 0.2) MHz, giving droplet speed v approx.(21.9 plus or minus 1.3) ms (exp. -1). Also shown on a separate channel is the corresponding signature using an amplitude demodulator circuit designed to detect the amplitude envelope of f(sub D) within the pulse profile. beta from simultaneous droplet events show a complete cyclic interference structure of maximum and minimum. The average period T of the complete cycle of interference is 13.02 plus or minus.39 microseconds. Toward the right edge of the profile, the interference disappears because one of the droplets is leaving the lidar beam while the other one remains in the beam, thus, giving beta for a single droplet.

Jarzembski, Maurice A.↗

Modeling the Compression of Merged Compact Toroids by Multiple Plasma Jets

A fusion propulsion scheme has been proposed that makes use of the merging of a spherical distribution of plasma jets to dynamically form a gaseous liner. The gaseous liner is used to implode a magnetized target to produce the fusion reaction in a standoff manner. In this paper, the merging of the plasma jets to form the gaseous liner is investigated numerically. The Los Alamos SPHINX code, based on the smoothed particle hydrodynamics method is used to model the interaction of the jets. 2-D and 3-D simulations have been performed to study the characteristics of the resulting flow when these jets collide. The results show that the jets merge to form a plasma liner that converge radially which may be used to compress the central plasma to fusion conditions. Details of the computational model and the SPH numerical methods will be presented together with the numerical results.

Thio, Y. C. Francis↗

The New MODIS-Terra, and the Proposed COBRA Mission: First Global Aerosol Distribution and Properties Over Land and Ocean, and Plans to Measure Global Black Carbon Absorption Over the Ocean Glint

The MODIS instrument was launched on the NASA Terra satellite in Dec. 1999. Since last Oct, the sensor and the aerosol algorithm reached maturity and provide global daily retrievals of aerosol optical thickness and properties. MODIS has 36 spectral channels in the visible to IR with resolution down to 250 m. This allows accurate cloud screening and multi-spectral aerosol retrievals. We derive the aerosol optical thickness over the ocean and most of the land areas, distinguishing between fine (mainly man-made aerosol) and coarse (mainly natural) aerosol particles. New methods to derive the aerosol absorption of sunlight are also being developed. These measurements are use to track different aerosol sources, transport and the radiative forcing at the top and bottom of the atmosphere. However MODIS or any present satellite sensor cannot measure absorption by Black Carbon over the oceans, a critical component in studying climate change and human health. For this purpose we propose the COBRA mission that observes the ocean at glint and off glint simultaneously measuring the spectral polarized light and deriving precisely the aerosol absorption.

Kaufman, Yoram J.↗

Excess chemical potential of small solutes across water--membrane and water--hexane interfaces

The excess chemical potentials of five small, structurally related solutes, CH4, CH3F, CH2F2, CHF3, and CF4, across the water-glycerol 1-monooleate bilayer and water-hexane interfaces were calculated at 300, 310, and 340 K using the particle insertion method. The excess chemical potentials of nonpolar molecules (CH4 and CF4) decrease monotonically or nearly monotonically from water to a nonpolar phase. In contrast, for molecules that possess permanent dipole moments (CH3F, CH2F, and CHF3), the excess chemical potentials exhibit an interfacial minimum that arises from superposition of two monotonically and oppositely changing contributions: electrostatic and nonelectrostatic. The nonelectrostatic term, dominated by the reversible work of creating a cavity that accommodates the solute, decreases, whereas the electrostatic term increases across the interface from water to the membrane interior. In water, the dependence of this term on the dipole moment is accurately described by second order perturbation theory. To achieve the same accuracy at the interface, third order terms must also be included. In the interfacial region, the molecular structure of the solvent influences both the excess chemical potential and solute orientations. The excess chemical potential across the interface increases with temperature, but this effect is rather small. Our analysis indicates that a broad range of small, moderately polar molecules should be surface active at the water-membrane and water-oil interfaces. The biological and medical significance of this result, especially in relation to the mechanism of anesthetic action, is discussed.

NASA Center ARC↗