Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle method”

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 37 records · Page 2

Computation of axisymmetric and ionized flows using particle and continuum methods

Comparisons between particle and continuum simulations of hypersonic near-continuum flows are presented. The particle approach employs the direct simulation Monte Carlo method (DSMC), and the continuum approach solves the Euler equations. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code which is efficiently vectorized is developed for this study. In this DSMC code, particular attention is paid to matching the relaxation rates employed in the continuum approach. This investigation represents a continuation of a previous study which considered thermochemical relaxation in one-dimensional shock waves of nitrogen. Comparison of the particle and continuum methods is first made for an axisymmetric blunt-body flow of air at 7 km/s. Very good agreement is obtained for the two solutions. The two techniques also compare well for a one-dimensional shock wave in air at 10 km/s. In both applications, the results are found to be sensitive to various aspects of the chemistry models employed.

Boyd, Iain D.↗

A Review of Discrete Element Method (DEM) Particle Shapes and Size Distributions for Lunar Soil

As part of ongoing efforts to develop models of lunar soil mechanics, this report reviews two topics that are important to discrete element method (DEM) modeling the behavior of soils (such as lunar soils): (1) methods of modeling particle shapes and (2) analytical representations of particle size distribution. The choice of particle shape complexity is driven primarily by opposing tradeoffs with total number of particles, computer memory, and total simulation computer processing time. The choice is also dependent on available DEM software capabilities. For example, PFC2D/PFC3D and EDEM support clustering of spheres; MIMES incorporates superquadric particle shapes; and BLOKS3D provides polyhedra shapes. Most commercial and custom DEM software supports some type of complex particle shape beyond the standard sphere. Convex polyhedra, clusters of spheres and single parametric particle shapes such as the ellipsoid, polyellipsoid, and superquadric, are all motivated by the desire to introduce asymmetry into the particle shape, as well as edges and corners, in order to better simulate actual granular particle shapes and behavior. An empirical particle size distribution (PSD) formula is shown to fit desert sand data from Bagnold. Particle size data of JSC-1a obtained from a fine particle analyzer at the NASA Kennedy Space Center is also fitted to a similar empirical PSD function.

Lane, John E.↗

Computation of axisymmetric and ionized hypersonic flows using particle and continuum methods

Comparisons between particle and continuum simulations of hypersonic near-continuum flows are presented. The particle approach employs the direct simulation Monte Carlo (DSMC) method, and the continuum approach solves the appropriate equations of fluid flow. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code that is efficiently vectorized is developed for this study. In this DSMC code, particular attention is paid to matching the relaxation rates employed in the continuum approach. This investigation represents a continuum of a previous study that considered thermochemical relaxation in one-dimensional shock waves of nitrogen. Comparison of the particle and continuum methods is first made for an axisymmetric blunt-body flow of air at 7 km/s. Very good agreement is obtained for the two solutions. The two techniques also compare well for a one-dimensional shock wave in air at 10 km/s. In both applications, the results are found to be sensitive to various aspects of the chemistry model employed.

Boyd, Iain D.↗

A collision-selection rule for a particle simulation method suited to vector computers

A theory is developed for a selection rule governing collisions in a particle simulation of rarefied gas-dynamic flows. The selection rule leads to an algorithmic form highly compatible with fine grain parallel decomposition, allowing for efficient utilization of supercomputers having vector or massively parallel single instruction multiple data architectures. A comparison of shock-wave profiles obtained using both the selection rule and Bird's direct simulation Monte Carlo (DSMC) method show excellent agreement. The equation on which the selection rule is based is shown to be directly related to the time-counter procedure in the DSMC method. The results of several example simulations of representative rarefied flows are presented, for which the number of particles used ranged from 10 to the 6th to 10 to the 7th demonstrating the greatly improved computational efficiency of the method.

Baganoff, D.↗

An Automated Refinement Process for Particle Trajectory Methods in GlennICE

Computational methods for ice accretion can simulate the impact of water drops and ice crystals on an aircraft surface in a Lagrangian reference frame or in the Eulerian reference frame. In the Eulerian reference frame, particles are considered a continuous fluid while in the Lagrangian frame individual particle trajectories are calculated. Methods that use the Eulerian reference frame are typically easier to develop as established modules used for continuum mechanics can be leveraged. The Eulerian systems can also be faster since the user does not have to simulate millions of particles in order to achieve good results. It is imperative therefore that a Lagrangian method optimize the release points of trajectories such that accurate solutions can be obtained while minimizing as much as possible the number of trajectories computed. This paper will present a methodology for this refinement process and demonstrate its effectiveness on sample three dimensional test cases.

William B Wright↗

An Automated Refinement Process for Particle Trajectory Methods in GlennICE

Computational methods for ice accretion can simulate the impact of water drops and ice crystals on an aircraft surface in a Lagrangian reference frame or in the Eulerian reference frame. In the Eulerian reference frame, particles are considered a continuous fluid while in the Lagrangian frame individual particle trajectories are calculated. Methods that use the Eulerian reference frame are typically easier to develop as established modules used for continuum mechanics can be leveraged. The Eulerian systems can also be faster since the user does not have to simulate millions of particles in order to achieve good results. It is imperative therefore that a Lagrangian method optimize the release points of trajectories such that accurate solutions can be obtained while minimizing as much as possible the number of trajectories computed. This paper will present a methodology for this refinement process and demonstrate its effectiveness on sample three dimensional test cases.

William Wright↗

Particle Seeding Method for Small-Scale, High-Pressure Nozzles

NASA tests new launch and reentry vehicle configurations in wind tunnels, where flow visualizations and quantitative flowfield measurements are often desired. Some of these vehicles have rocket motors for propulsion, retro-propulsion, or reaction control. High-pressure air is used to supply these rocket motor plumes. However, it is difficult to make off-body measurements in these regions, for several reasons. First, the plumes themselves lack seeding particles for flow diagnostics, and at the low pressures after expansion, Rayleigh scattering or other molecular techniques yield insufficient signal for flow velocity measurements. Second, the plumes displace particle-seeded tunnel air, preventing measurements in the vicinity of the plume. Third, the plumes force shock waves ahead of the vehicle, which melts the ice crystal fog commonly used for visualization and measurement techniques in certain facilities. In the current work, a novel method for seeding the flow in these small-scale, high-pressure nozzles has been devised and initially demonstrated, potentially enabling quantitative and qualitative measurements with particle-based instruments such as Doppler global velocimetry or particle image velocimetry. The method involves a Venturi contraction to draw the seed liquid out of a reservoir and into the nozzle channel, wherein shearing forces atomize the seed into particles. The concept was tested with a laser sheet visualization, which demonstrated that the flow rate of liquid spray was controllable; a valve could be adjusted to drop the flow rate by up to 65%. This relatively inexpensive and simple technique may prove useful in wind tunnel experiments involving particle-based laser diagnostics and small, high-pressure nozzles.

Aditya S. Acharya↗

Particle Seeding Method for Small-Scale, High-Pressure Nozzles

NASA tests new launch and reentry vehicle configurations in wind tunnels, where flow visualizations and quantitative flowfield measurements are often desired. Some of these vehicles have rocket motors for propulsion, retro-propulsion, or reaction control. High-pressure air is used to supply these rocket motor plumes. However, it is difficult to make off-body measurements in these regions, for several reasons. First, the plumes themselves lack seeding particles for flow diagnostics, and at the low pressures after expansion, Rayleigh scattering or other molecular techniques yield insufficient signal for flow velocity measurements. Second, the plumes displace particle-seeded tunnel air, preventing measurements in the vicinity of the plume. Third, the plumes force shock waves ahead of the vehicle, which melts the ice crystal fog commonly used for visualization and measurement techniques in certain facilities. In the current work, a novel method for seeding the flow in these small-scale, high-pressure nozzles has been devised and initially demonstrated, potentially enabling quantitative and qualitative measurements with particle-based instruments such as Doppler global velocimetry or particle image velocimetry. The method involves a Venturi contraction to draw the seed liquid out of a reservoir and into the nozzle channel, wherein shearing forces atomize the seed into particles. The concept was tested with a laser sheet visualization, which demonstrated that the flow rate of liquid spray was controllable; a valve could be adjusted to drop the flow rate by up to 65%. This relatively inexpensive and simple technique may prove useful in wind tunnel experiments involving particle-based laser diagnostics and small, high-pressure nozzles.

lasers↗

A comparative study of turbulence decay using Navier-Stokes and a discrete particle simulation

A comparative study of the two dimensional temporal decay of an initial turbulent state of flow is presented using a direct Navier-Stokes simulation and a particle method, ranging from the near continuum to more rarefied regimes. Various topics related to matching the initial conditions between the two simulations are considered. The determination of the initial velocity distribution function in the particle method was found to play an important role in the comparison. This distribution was first developed by matching the initial Navier-Stokes state of stress, but was found to be inadequate beyond the near continuum regime. An alternative approach of using the Lees two-sided Maxwellian to match the initial strain-rate is discussed. Results of the comparison of the temporal decay of mean kinetic energy are presented for a range of Knudsen numbers. As expected, good agreement was observed for the near continuum regime, but the differences found for the more rarefied conditions were unexpectedly small.

Goswami, A.↗

Neutral-particle wake method for measuring the atmospheric temperature from a satellite.

Description of a method that would permit a satellite-borne neutral mass spectrometer to measure the atmospheric temperature. The spectrometer examines the partial pressure variations that occur as the wake of a small rectangular baffle is swept across the entrance orifice of the spectrometer. For a given baffle size and for a mounting distance from the orifice, the depth of the resulting pressure minimum depends only on the thermal velocity or temperature of the observed species. The validity of the method can be checked by measuring the wake characteristics of more than one species and/or by employing each of several baffle sizes. The theory includes the effect of a finite orifice size, finite baffle length, and the backscattering of particles from the baffle into the orifice. It is found that a suitable baffle arrangement can be achieved that will permit the temperature to be measured over at least the range normally encountered in the thermosphere (200 to 2000 K) and, depending on the sensitivity and background pressure of the spectrometer, over an altitude range of about 140 to 600 km.

Brace, L. H.↗

Particle-Surface Interaction Model and Method of Determining Particle-Surface Interactions

A method and model of predicting particle-surface interactions with a surface, such as the surface of a spacecraft. The method includes the steps of: determining a trajectory path of a plurality of moving particles; predicting whether any of the moving particles will intersect a surface; predicting whether any of the particles will be captured by the surface and/or; predicting a reflected trajectory and velocity of particles reflected from the surface.

Hughes, David W.↗

A neutral particle wake method for measuring the atmospheric temperature from a satellite

A method is described that permits a satellite-borne neutral mass spectrometer to perform measurements of the atmospheric temperature. The method employs the spectrometer to examine the partial pressure variations that occur as the wake of a small rectangular baffle is swept across the entrance orifice of the spectrometer. For a given baffle size and mounting distance from the orifice the depth of the resulting pressure minimum depends only upon the thermal velocity, or temperature, of the observed species. A check upon the validity of the method is obtained by measuring the wake characteristics of more than one species and/or by employing each of several baffle sizes.

Brace, L. H.↗

Laser Doppler spectrometer method of particle sizing

A spectrometer for the detection of airborne particulate pollution in the submicron size range is described. In this device, airborne particles are accelerated through a supersonic nozzle, with different sizes achieving different velocities in the gas flow. Information about the velocities of the accelerated particles is obtained with a laser-heterodyne optical system through the Doppler shift of light scattered from the particles. Detection is accomplished by means of a photomultiplier. Nozzle design and signal processing techniques are also discussed.

Weber, F. N.↗

Calculating Charge Transport in Semiconductors

Method, called flux method, improvement over conventional methods, easy to use, and applies in operating regime not accessible to conventional macroscopic methods. Key concept of flux method follows behavior of particle flux rather than concentration of particles. Method not subject to limitations that follow from local or infinitesimal analysis.

Byvik, C. E.↗

Models of thermal relaxation mechanics for particle simulation methods

A comparative study is conducted of the simulation of molecular vibrational relaxation in rarefied flows by means of existing models and a novel phenomenological microscopic model. This improved technique iterates between translation-rotation and rotation-vibration exchanges; equilibrium is in this way promoted if the model for the latter process is compatible with quantized oscillators. The iteration-equipartition model retains computational simplicity and promotes thermal equilibrium, even in the case of nondegenerate anharmonic quantized oscillators and multispecies gas mixtures.

Haas, Brian L.↗

Determining Wheel-Soil Interaction Loads Using a Meshfree Finite Element Approach Assisting Future Missions with Rover Wheel Design

A wheel experiencing sinkage and slippage events poses a high risk to rover missions as evidenced by recent mobility challenges on the Mars Exploration Rover (MER) project. Because several factors contribute to wheel sinkage and slippage conditions such as soil composition, large deformation soil behavior, wheel geometry, nonlinear contact forces, terrain irregularity, etc., there are significant benefits to modeling these events to a sufficient degree of complexity. For the purposes of modeling wheel sinkage and slippage at an engineering scale, meshfree finite element approaches enable simulations that capture sufficient detail of wheel-soil interaction while remaining computationally feasible. This study demonstrates some of the large deformation modeling capability of meshfree methods and the realistic solutions obtained by accounting for the soil material properties. A benchmark wheel-soil interaction problem is developed and analyzed using a specific class of meshfree methods called Reproducing Kernel Particle Method (RKPM). The benchmark problem is also analyzed using a commercially available finite element approach with Lagrangian meshing for comparison. RKPM results are comparable to classical pressure-sinkage terramechanics relationships proposed by Bekker-Wong. Pending experimental calibration by future work, the meshfree modeling technique will be a viable simulation tool for trade studies assisting rover wheel design.

Contreras, Michael T.↗

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