Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

OceanWATERS Lander Robotic Arm Operation

Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) is an open-source simulator for developing onboard autonomy software for robotic exploration of ocean worlds, such as Europa, Enceladus, and Titan, built on the Robot Operating System (ROS) and Gazebo simulation environment. Inevitable ground communication delays increase demand for a high degree of autonomy during excavation, collection and transfer of samples to scientific instruments for in-situ analysis. This paper offers a detailed discussion of the robotic arm design and operation for such autonomous surface exploration, taking as reference the Europa Lander mission. The lander arm, which is designed primarily to acquire icy surface and subsurface samples within the arm’s workspace, is a 6-degree-of-freedom manipulator with two end effectors: a sample excavation tool and a trenching end-effector. The robotic arm’s modes and operations can be summarized as follows: stowed arm, intended as the lander arm default configuration characterized by zero-power consumption; un-stowed arm, target arm configuration after its first deployment; selection and deployment of the end-effector to use next; guarded move, to detect ground level at the desired trenching location; drill ice using the grinder; dig trench at a particular location using the scoop; deliver sample to the sample transfer dock; discard redundant samples. The motion planning tool used for the lander arm is MoveIt, a ROS package. MoveIt uses sampling-based planning and collision checking libraries to determine safe paths. The Rapidly Exploring Random Trees* (RRT*) has been chosen as default planning algorithm as it provides optimal plans with an exponential speed and is guaranteed to find a solution, if feasible solutions exist. Furthermore, this work quantifies and discusses the energy requirements for excavating and collecting samples. In OceanWATERS, force feedback from the terrain, which influences the arm dynamics, is modelled using a discrete element method (DEM) simulation. The DEM and Gazebo software run in parallel and communicate through a co-simulation plugin. This paper presents an analysis and comparison of three DEM open source software (YADE, ESyS-Particle, Project Chrono) for implementation in OceanWATERS and motivates the choice of YADE as most suitable candidate.

Damiana Catanoso↗
NASA NTRS2008

John F. Kennedy Space Center's Technology Development and Application 2006-2007 Report

Topics covered include: Reversible Chemochromic Hydrogen Detectors; Determining Trajectory of Triboelectrically Charged Particles, Using Discrete Element Modeling; Using Indium Tin Oxide To Mitigate Dust on Viewing Ports; High-Performance Polyimide Powder Coatings; Controlled-Release Microcapsules for Smart Coatings for Corrosion Applications; Aerocoat 7 Replacement Coatings; Photocatalytic Coatings for Exploration and Spaceport Design; New Materials for the Repair of Polyimide Electrical Wire Insulation; Commodity-Free Calibration; Novel Ice Mitigation Methods; Crack Offset Measurement With the Projected Laser Target Device; New Materials for Structural Composites and Protective Coatings; Fire Chemistry Testing of Spray-On Foam Insulation (SOFI); Using Aerogel-Based Insulation Material To Prevent Foam Loss on the Liquid-Hydrogen Intertank; Particle Ejection and Levitation Technology (PELT); Electrostatic Characterization of Lunar Dust; Numerical Analysis of Rocket Exhaust Cratering; RESOLVE Projects: Lunar Water Resource Demonstration and Regolith Volatile Characterization; Tribocharging Lunar Soil for Electrostatic Beneficiation; Numerically Modeling the Erosion of Lunar Soil by Rocket Exhaust Plumes; Trajectory Model of Lunar Dust Particles; Using Lunar Module Shadows To Scale the Effects of Rocket Exhaust Plumes; Predicting the Acoustic Environment Induced by the Launch of the Ares I Vehicle; Measuring Ultrasonic Acoustic Velocity in a Thin Sheet of Graphite Epoxy Composite; Hail Size Distribution Mapping; Launch Pad 39 Hail Monitor Array System; Autonomous Flight Safety System - Phase III; The Photogrammetry Cube; Bird Vision System; Automating Range Surveillance Through Radio Interferometry and Field Strength Mapping Techniques; Next-Generation Telemetry Workstation; GPS Metric Tracking Unit; and Space-Based Range.

Source record↗
NASA NTRS1978

Energy deposition in a gaseous mixture

Existing models in planetary aeronomy involve a variety of possible gaseous mixtures subjected to charged particle bombardment. In the present paper, a phenomenological approach to electron degradation in such mixtures is proposed to facilitate analysis of such problems. Existing parameterized yield spectra for a variety of gases, computed by a discrete energy bin method, are used as a basis for constructing a parameterized composite yield spectrum for an arbitrary mixture. The approach proposed is an improvement over the continuous slow-down approximation in that discreteness is taken into account phenomenologically. The advantages are in the ease of application to aeronomical problems involving variable mixtures.

Peterson, L. R.↗
NASA NTRS2019

Comprehensive Thematic T-Matrix Reference Database: a 2017-2019 Update

Since its inception in the mid-1960s, the T-matrix method has been one of the most versatile and efficient numerically exact computer solvers of the time-harmonic macroscopic Maxwell equations. It has been widely used for the computation of electromagnetic scattering by single and composite particles, discrete random media, periodic structures (including metamaterials), and particles in the vicinity of plane or rough interfaces separating media with different refractive indices. This compilation is the ninth update to the comprehensive thematic database of peer-reviewed T-matrix publications initiated in 2004 and lists relevant publications that have appeared since 2017. It also includes a few earlier publications that have so far been overlooked.

Michael I Mishchenko↗
NASA NTRS1995

A deterministic particle method for one-dimensional reaction-diffusion equations

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Mascagni, Michael↗
NASA NTRS2002

Particle Swarm Optimization

The purpose of this paper is to show how the search algorithm known as particle swarm optimization performs. Here, particle swarm optimization is applied to structural design problems, but the method has a much wider range of possible applications. The paper's new contributions are improvements to the particle swarm optimization algorithm and conclusions and recommendations as to the utility of the algorithm, Results of numerical experiments for both continuous and discrete applications are presented in the paper. The results indicate that the particle swarm optimization algorithm does locate the constrained minimum design in continuous applications with very good precision, albeit at a much higher computational cost than that of a typical gradient based optimizer. However, the true potential of particle swarm optimization is primarily in applications with discrete and/or discontinuous functions and variables. Additionally, particle swarm optimization has the potential of efficient computation with very large numbers of concurrently operating processors.

Venter, Gerhard↗
NASA NTRS1988

Light scattering of large rough particles application to cometary grains

While the electromagnetic field scattered by a spherical particle is classically obtained by the Helmholtz equation, the general case of an arbitrary particle may be investigated in the general framework of the interaction of a wave with a scattering potential. A wave function then satisfies the Schroedinger equation. The general solution of the Schroedinger equation is given. The main disadvantage of this approach are its restriction to large particles and its scalar nature preventing the calculation of the polarization. However, Perrin and Lamy have shown how to avoid the second limitation and retrieve a vectorial description. They proved that in the case of large spheres when the ad hoc assumptions are satisfied, the expression of the scattering amplitude may be approximated by an expansion series in partial waves, i.e., on a discrete basis. The analogy may be generalized, and the ratio of the two components for a rough particle obtained by taking the ratio of the reflectivities for the two directions of polarization. These reflectivities involve the simple and double reflections calculated following the method developed by Wolff for rough surfaces. The theory is further detailed.

Lamy, Philippe L.↗
NASA NTRS2018

An Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. II: Approach Based on the Dyson and Bethe-Salpeter Equations

In this paper, the vector radiative transfer equation is derived by means of the vector integral Foldy equations describing the electromagnetic scattering by a group of particles. By assuming that in a discrete random medium the positions of the particles are statistically independent and by applying the Twersky approximation to the order-of-scattering expansion of the total field, we derive the Dyson equation for the coherent field and the ladder approximated Bethe–Salpeter equation for the dyadic correlation function. Then, under the far-field assumption for sparsely distributed particles, the Dyson equation is reduced to the Foldy integral equation for the coherent field, while the iterated solution of the Bethe–Salpeter equation ultimately yields the vector radiative transfer equation.

Electromagnetic scattering↗
NASA NTRS1973

Calculations of neutron flux spectra induced in the earth's atmosphere by galactic cosmic rays.

Calculations have been carried out to determine the neutron flux induced in the earth's atmosphere by galactic protons and alpha particles at solar minimum for a geomagnetic latitude of 42 N. Neutron flux spectra in the energy range from about 10 to the minus 8th to about 100,000 MeV at various depths in the atmosphere were calculated by using Monte Carlo and discrete ordinates methods, and various comparisons with experimental data are presented. The magnitude and shape of the calculated neutron leakage spectrum at the particular latitude considered support the theory that the cosmic ray albedo neutron decay (Crand) mechanism is the source of the protons trapped in the inner radiation belt.

Armstrong, T. W.↗
NASA NTRS1985

Real-time measurement of sodium chloride in individual aerosol particles by mass spectrometry

The method of particle analysis by mass spectrometry has been applied to the quantitative measurement of sodium chloride in individual particles on a real-time basis. Particles of known masses are individually introduced, in the form of a beam, into a miniature Knudsen cell oven (1600 K). The oven is fabricated from rhenium metal sheet (0.018 mm thick) and is situated in the ion source of a quadrupole mass spectrometer. A particle once inside the oven is trapped and completely volatilized; this overcomes the problem of partial volatilization due to particles bouncing from the filament surface. Individual particles are thermally volatilized and ionized inside the rhenium oven, and produce discrete sodium ion pulses whose intensities are measured with the quadrupole mass spectrometer. An ion pulse width of several milliseconds (4-12 ms) is found for particles in the mass range 1.3 x 10 to the -13th to 5.4 x 10 to the -11th g. The sodium ion intensity is found to be proportional to the particle mass to the 0.86-power. The intensity distribution for monodisperse aerosol particles possesses a geometric standard deviation of 1.09, showing that the method can be used for the determination of the mass distribution function with good resolution in a polydisperse aerosol.

Sinha, M. P.↗
NASA NTRS2021

Evaluation of Higher-order Quadrature Schemes in Improving Computational Efficiency for Orientation-averaged Single-Scattering Properties of Nonspherical Ice Particles

We evaluate several high-order quadrature schemes for accuracy and efficacy in obtaining orientation-averaged single-scattering properties (SSPs). We use the highly efficient MIDAS to perform electromagnetic scattering calculations to evaluate the gain in efficiency from these schemes. MIDAS is shown to be superior to DDSCAT, a popular discrete dipole approximation (DDA) method. This study is motivated by the fact that quality physical precipitation retrievals rely on using accurate orientation-averaged SSPs derived from realistic hydrometeors as input to radiative transfer Models (RTMs). The DDA has been a popular choice for single-scattering calculations, due to its versatility with respect to target geometry. However, being iterative-solver-based (ISB), the most used DDA codes, e.g. DDSCAT and ADDA, must solve the scattering problem for each orientation of the target separately. As the size parameter and geometric anisotropy of the hydrometeor increase, the number of orientations needed to obtain accurate orientation-averages can increase drastically and so does the computation cost incurred by the ISB-DDA methods. MIDAS is a Direct-Solver-Based (DSB) code, its decomposition of the original large matrix with a high rank into multiple more manageable smaller matrices of lower ranks makes it much more computationally efficient while maintaining excellent accuracy. In addition, direct solvers consider all requested orientations at once, giving MIDAS further advantage over popular ISB-DDA methods. MIDAS, when combined with high-order quadrature for orientation averaging, can be greater than three orders of magnitude more efficient in obtaining RTM-ready SSPs of complex-shaped hydrometeors than existing ISB-DDA methods, with the native quadrature schemes they offer.

Ines Fenni↗
NASA NTRS1998

Measurements of Heat Transfer, Flow, and Pressures in a Simulated Turbine Blade Internal Cooling Passage

An experimental study was made to obtain quantitative information on heat transfer, flow, and pressure distribution in a branched duct test section that had several significant features of an internal cooling passage of a turbine blade. The objective of this study was to generate a set of experimental data that could be used for validation of computer codes that would be used to model internal cooling. Surface heat transfer coefficients and entrance flow conditions were measured at nominal entrance Reynolds numbers of 45,000, 335,000, and 726,000. Heat transfer data were obtained by using a steady-state technique in which an Inconel heater sheet is attached to the surface and coated with liquid crystals. Visual and quantitative flow-field data from particle image velocimetry measurements for a plane at midchannel height for a Reynolds number of 45,000 were also obtained. The flow was seeded with polystyrene particles and illuminated by a laser light sheet. Pressure distribution measurements were made both on the surface with discrete holes and in the flow field with a total pressure probe. The flow-field measurements yielded flow-field velocities at selected locations. A relatively new method, pressure sensitive paint, was also used to measure surface pressure distribution. The pressure paint data obtained at Reynolds numbers of 335,000 and 726,000 compared well with the more standard method of measuring pressures by using discrete holes.

Russell, Louis M.↗
NASA NTRS2022

Multi-Scale Modelling of the Bound Metal Deposition Manufacturing of Ti6Al4V

Nonlinear shrinkage of the metal part during manufacturing by bound metal deposition, both on the ground and under microgravity, is considered. A multi-scale physics-based approach is developed to address the problem. It spans timescales from atomistic dynamics on the order of nanoseconds to full-part shrinkage on the order of hours. This approach enables estimation of the key parameters of the problem, including the widths of grain boundaries, the coefficient of surface diffusion, the initial redistribution of particles during the debinding stage, the evolution of the microstructure from round particles to densely-packed grains, the corresponding changes in the total and chemical free energies, and the sintering stress. The method has been used to predict shrinkage at the levels of two particles, of the filament cross-section, of the sub-model, and of the whole green, brown, and metal parts.

Nonlinear shrinkage↗
NASA NTRS1993

Preliminary design of an intermittent smoke flow visualization system

A prototype intermittent flow visualization system that was designed to study vortex flow field dynamics has been constructed and tested through its ground test phase. It produces discrete pulses of dense white smoke consisting of particles of terephthalic acid by the pulsing action of a fast-acting three-way valve. The trajectories of the smoke pulses can be tracked by a video imaging system without intruding in the flow around in flight. Two methods of pulsing the smoke were examined. The simplest and safest approach is to simply divert the smoke between the two outlet ports on the valve; this approach should be particularly effective if it were desired to inject smoke at two locations during the same test event. The second approach involves closing off one of the outlet ports to momentarily block the flow. The second approach requires careful control of valve dwell times to avoid excessive pressure buildup within the cartridge container. This method also increases the velocity of the smoke injected into the flow. The flow of the smoke has been blocked for periods ranging from 30 to 80 milliseconds, depending on the system volume and the length of time the valve is allowed to remain open between valve closings.

Ward, Donald T.↗
NASA NTRS1984

An atomistic model of slip formation

The results of an atomistic model for the simulation of the early stages of crack initiation in a two-dimensional triangular lattice are presented. In the current model, each particle in the system is treated discretely and assumed to be interacting with the surrounding particles via Lennard-Jones potentials. A uniaxial load (in incremental elongations) is applied to the rectangular two-dimensional slab in either the x or the y direction. After each incremental elongation the system is equilibrated using a static method. Initially, elastic behavior in the x and y directions is observed. Continued elongation results in plastic deformation. In lattices with point defects, the defects first move to the surface, creating vacancies which trigger plastic deformation.

Halicioglu, T.↗
NASA NTRS1983

Helicopter rotor loads using discretized matched asymptotic expansions

The numerical practicality of a matched asymptotic expansion approach for the computation of unsteady three dimensional airloads on a helicopter rotor was improved. This effort utilizes a discretized repesentation of the doublet strength distribution and helical streamlines to decrease the computational requirements of the original analysis. The continuous variation of the doublet strength was approximated by piecewise constant or piecewise quadratic distributions, and the helical trajectory of a fluid particle was approximated by connected straight line segments. As a direct result of these simplified representations the computational time required for the execution of a typical flight condition was reduced by an order of magnitude with respect to the requirements of the original analysis. Airloads which were computed using the discretized method for a two bladed model rotor and a full scale four bladed rotor are in close agreement with measured results and airloads from the original asymptotic analysis. For conditions characterized by significant rotor/wake interaction the piecewise constant representation requires a reduced azimuth spacing to maintain acceptable accuracy.

Pierce, G. A.↗
NASA NTRS1992

Numerical modeling for dilute and dense sprays

We have successfully implemented a numerical model for spray-combustion calculations. In this model, the governing gas-phase equations in Eulerian coordinate are solved by a time-marching multiple pressure correction procedure based on the operator-splitting technique. The droplet-phase equations in Lagrangian coordinate are solved by a stochastic discrete particle technique. In order to simplify the calculation procedure for the circulating droplets, the effective conductivity model is utilized. The k-epsilon models are utilized to characterize the time and length scales of the gas phase in conjunction with turbulent modulation by droplets and droplet dispersion by turbulence. This method entails random sampling of instantaneous gas flow properties and the stochastic process requires a large number of computational parcels to produce the satisfactory dispersion distributions even for rather dilute sprays. Two major improvements in spray combustion modelings were made. Firstly, we have developed a probability density function approach in multidimensional space to represent a specific computational particle. Secondly, we incorporate the Taylor Analogy Breakup (TAB) model for handling the dense spray effects. This breakup model is based on the reasonable assumption that atomization and drop breakup are indistinguishable processes within a dense spray near the nozzle exit. Accordingly, atomization is prescribed by injecting drops which have a characteristic size equal to the nozzle exit diameter. Example problems include the nearly homogeneous and inhomogeneous turbulent particle dispersion, and the non-evaporating, evaporating, and burning dense sprays. Comparison with experimental data will be discussed in detail.

Chen, C. P.↗
NASA NTRS1991

Chaotic jumps in the generalized first adiabatic invariant in current sheets

The present study examines how the changes in the generalized first adiabatic invariant J derived from the separatrix crossing theory can be incorporated into the drift variable approach to generating distribution functions. A method is proposed for determining distribution functions for an ensemble of particles following interaction with the tail current sheet by treating the interaction as a scattering problem characterized by changes in the invariant. Generalized drift velocities are obtained for a 1D tail configuration by using the generalized first invariant. The invariant remained constant except for the discrete changes caused by chaotic scattering as the particles cross the separatrix.

Brittnacher, M. J.↗