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 775 records · Page 43

Use of the thin sheath approximation for obtaining ion temperatures from the ISEE 1 limited aperture RPA

A procedure for analyzing low-energy (less than approximately 100 eV) ion data from the plasma composition experiment on ISEE 1 is set forth. The method is based on a derived analytic expression for particle flux to a limited aperture retarding potential analyzer (RPA) in the thin sheath approximation, which makes allowance for some effects of a charged spacecraft on plasma particle trajectories. Calculations using simulated data are employed in testing the efficacy and accuracy of the technique. On the basis of an analysis of these calculation results and the mathematical model, the method is seen as being able to provide accurate ion temperatures from all good plasmaspheric RPA data. It is noted that corresponding densities and spacecraft potentials should be accurate when spacecraft potentials are negative but that they are subject to error for positive spacecraft potentials, particularly when ion Mach numbers are much less than 1. An analysis of data from a representative ISEE 1 pass produces a plasmasphere temperature profile that is consistent in overall structure with previous observations.

Comfort, R. H.↗

Non-local transport in turbulent MHD convection

The nonlocal non-diffusive transport of passive scalars in turbulent magnetohydrodynamic (MHD) convection is investigated using transilient matrices. These matrices describe the probability that a tracer particle beginning at one position in a flow will be advected to another position after some time. A method for the calculation of these matrices from simulation data which involves following the trajectories of passive tracer particles and calculating their transport statistics, is presented. The method is applied to study the transport in several simulations of turbulent, rotating, three dimensional compressible, penetrative MDH convection. Transport coefficients and other diagnostics are used to quantify the transport, which is found to resemble advection more closely than diffusion. Some of the results are found to have direct relevance to other physical problems, such as the light element depletion in sun-type stars. The large kurtosis found for downward moving particles at the base of the convection zone implies several extreme events.

Miesch, Mark↗

Soil Moisture Data Assimilation to Estimate Irrigation Water Use

Knowledge of irrigation is essential to support food security, manage depleting water resources, and comprehensively understand the global water and energy cycles. Despite the importance of understanding irrigation, little consistent information exists on the amount of water that is applied for irrigation. In this study, we develop and evaluate a new method to predict daily to seasonal irrigation magnitude using a particle batch smoother data assimilation approach, where land surface model soil moisture is applied in different configurations to understand how characteristics of remotely sensed soil moisture may impact the performance of the method. The study employs a suite of synthetic data assimilation experiments, allowing for systematic diagnosis of known error sources. Assimilation of daily synthetic soil moisture observations with zero noise produces irrigation estimates with a seasonal bias of 0.66% and a correlation of 0.95 relative to a known truth irrigation. When synthetic observations were subjected to an irregular overpass interval and random noise similar to the Soil Moisture Active Passive satellite (0.04 cm(exp 3) cm(exp -3)), irrigation estimates produced a median seasonal bias of <1% and a correlation of 0.69. When systematic biases commensurate with those between NLDAS‐2 land surface models and Soil Moisture Active Passive are imposed, irrigation estimates show larger biases. In this application, the particle batch smoother outperformed the particle filter. The presented framework has the potential to provide new information into irrigation magnitude over spatially continuous domains, yet its broad applicability is contingent upon identifying new method(s) of determining irrigation schedule and correcting biases between observed and simulated soil moisture, as these errors markedly degraded performance.

R Abolafia-Rosenzweig↗

A generalized Zel'dovich approximation to gravitational instability

The orbits of particles undergoing gravitational instability are parameterized by generalizing the Zel'dovich approximation to a series expansion of arbitrary accuracy in the nonlinear regime. The coefficients of this series are determined from an action principle, or, more generally, from moments of the equation of motion. It is shown that the series is more rapidly convergent than previous nonlinear approximations. The method therefore provides a practical means of determining particle orbits, even for highly nonlinear perturbations. As an alternative, we also outline how the nonlinear dynamics may be computed as a field theory in which the evolution of the density and the velocity is determined in fixed comoving Eulerian coordinates.

Giavalisco, M.↗

Global Aspects of Charged Particle Motion in Axially Symmetric Multipole Magnetic Fields

The motion of a single charged particle in the space outside of a compact region of steady currents is investigated. The charged particle is assumed to produce negligible electromagnetic radiation, so that its energy is conserved. The source of the magnetic field is represented as a point multipole. After a general description, attention is focused on magnetic fields with axial symmetry. Lagrangian dynamical theory is utilized to identify constants of the motion as well as the equations of motion themselves. The qualitative method of Stonner is used to examine charged particle motion in axisymmetric multipole fields of all orders. Although the equations of motion generally have no analytical solutions and must be integrated numerically to produce a specific orbit, a topological examination of dynamics is possible, and can be used, d la Stonner, to completely describe the global aspects of the motion of a single charged particle in a space with an axisymmetric multipole magnetic field.

Shebalin, John V.↗

Evaluation of thermochemical models for particle and continuum simulations of hypersonic flow

Computations are presented for one-dimensional, strong shock waves that are typical of those that form in front of a reentering spacecraft. The fluid mechanics and thermochemistry are modeled using two different approaches. The first employs traditional continuum techniques in solving the Navier-Stokes equations. The second approach employs a particle simulation technique (the direct simulation Monte Carlo method, DSMC). The thermochemical models employed in these two techniques are quite different. The present investigation presents an evaluation of thermochemical models for nitrogen under hypersonic flow conditions. Four separate cases are considered that are dominated in turn by vibrational relaxation, weak dissociation, strong dissociation and weak ionization. In near-continuum, hypersonic flow, the nonequilibrium thermochemical models employed in continuum and particle simulations produce nearly identical solutions. Further, the two approaches are evaluated successfully against available experimental data for weakly and strongly dissociating flows.

Boyd, Iain D.↗

Radiation Shielding Materials Containing Hydrogen, Boron, and Nitrogen: Systematic Computational and Experimental Study

The key objectives of this study are to investigate, both computationally and experimentally, which forms, compositions, and layerings of hydrogen, boron, and nitrogen containing materials will offer the greatest shielding in the most structurally robust combination against galactic cosmic radiation (GCR), secondary neutrons, and solar energetic particles (SEP). The objectives and expected significance of this research are to develop a space radiation shielding materials system that has high efficacy for shielding radiation and that also has high strength for load bearing primary structures. Such a materials system does not yet exist. The boron nitride nanotube (BNNT) can theoretically be processed into structural BNNT and used for load bearing structures. Furthermore, the BNNT can be incorporated into high hydrogen polymers and the combination used as matrix reinforcement for structural composites. BNNT's molecular structure is attractive for hydrogen storage and hydrogenation. There are two methods or techniques for introducing hydrogen into BNNT: (1) hydrogen storage in BNNT, and (2) hydrogenation of BNNT (hydrogenated BNNT). In the hydrogen storage method, nanotubes are favored to store hydrogen over particles and sheets because they have much larger surface areas and higher hydrogen binding energy. The carbon nanotube (CNT) and BNNT have been studied as potentially outstanding hydrogen storage materials since 1997. Our study of hydrogen storage in BNNT - as a function of temperature, pressure, and hydrogen gas concentration - will be performed with a hydrogen storage chamber equipped with a hydrogen generator. The second method of introducing hydrogen into BNNT is hydrogenation of BNNT, where hydrogen is covalently bonded onto boron, nitrogen, or both. Hydrogenation of BN and BNNT has been studied theoretically. Hyper-hydrogenated BNNT has been theoretically predicted with hydrogen coverage up to 100% of the individual atoms. This is a higher hydrogen content than possible with hydrogen storage; however, a systematic experimental hydrogenation study has not been reported. A combination of the two approaches may be explored to provide yet higher hydrogen content. The hydrogen containing BNNT produced in our study will be characterized for hydrogen content and thermal stability in simulated space service environments. These new materials systems will be tested for their radiation shielding effectiveness against high energy protons and high energy heavy ions at the HIMAC facility in Japan, or a comparable facility. These high energy particles simulate exposure to SEP and GCR environments. They will also be tested in the LaRC Neutron Exposure Laboratory for their neutron shielding effectiveness, an attribute that determines their capability to shield against the secondary neutrons found inside structures and on lunar and planetary surfaces. The potential significance is to produce a radiation protection enabling technology for future exploration missions. Crew on deep space human exploration missions greater than approximately 90 days cannot remain below current crew Permissible Exposure Limits without shielding and/or biological countermeasures. The intent of this research is to bring the Agency closer to extending space missions beyond the 90-day limit, with 1 year as a long-term goal. We are advocating a systems solution with a structural materials component. Our intent is to develop the best materials system for that materials component. In this Phase I study, we have shown, computationally, that hydrogen containing BNNT is effective for shielding against GCR, SEP, and neutrons over a wide range of energies. This is why we are focusing on hydrogen containing BNNT as an innovative advanced concept. In our future work, we plan to demonstrate, experimentally, that hydrogen, boron, and nitrogen based materials can provide mechanically strong, thermally stable, structural materials with effective radiation shielding against GCR, SEP, and neutrons.

Radiation↗

A Melting Layer Model for Passive/Active Microwave Remote Sensing Applications: Model Formulation and Comparison with Observations - Part 1

In this study, a 1-D steady-state microphysical model which describes the vertical distribution of melting precipitation particles is developed. The model is driven by the ice-phase precipitation distributions just above the freezing level at applicable gridpoints of "parent" 3-D cloud-resolving model (CRM) simulations. It extends these simulations by providing the number density and meltwater fraction of each particle in finely separated size categories through the melting layer. The depth of the modeled melting layer is primarily determined by the initial material density of the ice-phase precipitation. The radiative properties of melting precipitation at microwave frequencies are calculated based upon different methods for describing the dielectric properties of mixed phase particles. Particle absorption and scattering efficiencies at the Tropical Rainfall Measuring Mission Microwave Imager frequencies (10.65 to 85.5 GHz) are enhanced greatly for relatively small (approx. 0.1) meltwater fractions. The relatively large number of partially-melted particles just below the freezing level in stratiform regions leads to significant microwave absorption, well-exceeding the absorption by rain at the base of the melting layer. Calculated precipitation backscatter efficiencies at the Precipitation Radar frequency (13.8 GHz) increase in proportion to the particle meltwater fraction, leading to a "bright-band" of enhanced radar reflectivities in agreement with previous studies. The radiative properties of the melting layer are determined by the choice of dielectric models and the initial water contents and material densities of the "seeding" ice-phase precipitation particles. Simulated melting layer profiles based upon snow described by the Fabry-Szyrmer core-shell dielectric model and graupel described by the Maxwell-Garnett water matrix dielectric model lead to reasonable agreement with radar-derived melting layer optical depth distributions. Moreover, control profiles that do not contain mixed-phase precipitation particles yield optical depths that are systematically lower than those observed. Therefore, the use of the melting layer model to extend 3-D CRM simulations appears justified, at least until more realistic spectral methods for describing melting precipitation in high-resolution, 3-D CRM's are implemented.

Olson, William S.↗

Development and application of computational aerothermodynamics flowfield computer codes

Computations are presented for one-dimensional, strong shock waves that are typical of those that form in front of a reentering spacecraft. The fluid mechanics and thermochemistry are modeled using two different approaches. The first employs traditional continuum techniques in solving the Navier-Stokes equations. The second-approach employs a particle simulation technique (the direct simulation Monte Carlo method, DSMC). The thermochemical models employed in these two techniques are quite different. The present investigation presents an evaluation of thermochemical models for nitrogen under hypersonic flow conditions. Four separate cases are considered. The cases are governed, respectively, by the following: vibrational relaxation; weak dissociation; strong dissociation; and weak ionization. In near-continuum, hypersonic flow, the nonequilibrium thermochemical models employed in continuum and particle simulations produce nearly identical solutions. Further, the two approaches are evaluated successfully against available experimental data for weakly and strongly dissociating flows.

Venkatapathy, Ethiraj↗

A New Optical Aerosol Spectrometer

An optical particle spectrometer capable of measuring aerosol particle size distributions from 0.02 to 100 micrometers has been developed. This instrument combines several optical methods in one, in-situ configuration; it can provide continuous data collection to encompass the wide dynamic size ranges and concentrations found in studies of modeled planetary atmospheres as well as terrestrial air quality research. Currently, the system is incorporated into an eight liter capacity spherical pressure vessel that is appropriate both for flowthrough and for in-situ particle generation. The optical sizing methods include polarization ratio, The scattering, and forward scattering detectors, with illumination from a fiber-coupled, Argon-ion laser. As particle sizes increase above 0.1 micrometer, a customized electronics and software system automatically shifts from polarization to diffraction-based measurements as the angular scattering detectors attain acceptable signal-to-noise ratios. The number concentration detection limits are estimated to be in the part-per-trillion (ppT by volume) range, or roughly 1000 submicron particles per cubic centimeter. Results from static experiments using HFC134A (approved light scattering gas standard), flow-through experiments using sodium chloride (NaCl) and carbon particles, and dynamic 'Tholin' (photochemical produced particles from ultraviolet (UV)-irradiated acetylene and nitrogen) experiments have been obtained. The optical spectrometer data obtained with particles have compared well with particle sizes determined by electron microscopy. The 'Tholin' tests provided real-time size and concentration data as the particles grew from about 30 nanometers to about 0.8 micrometers, with concentrations ranging from ppT to ppB, by volume. Tests are still underway, to better define sizing accuracy and concentration limits, these results will be reported.

Fonda, Mark↗

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↗

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↗

Problems in calibrating a polar nephelometer

A calibration procedure for nephelometers which is based on the method developed by Pritchard and Elliott (1960) is outlined. Basically, the method consists of comparing the intensity of the light scattered from the particles in the common volume with the intensity reflected from a calibration standard. A proof of the calibration method is given.

Holland, A.↗

New Method Developed To Purify Single Wall Carbon Nanotubes for Aerospace Applications

Single wall carbon nanotubes have attracted considerable attention because of their remarkable mechanical properties and electrical and thermal conductivities. Use of these materials as primary or secondary reinforcements in polymers or ceramics could lead to new materials with significantly enhanced mechanical strength and electrical and thermal conductivity. Use of carbon-nanotube-reinforced materials in aerospace components will enable substantial reductions in component weight and improvements in durability and safety. Potential applications for single wall carbon nanotubes include lightweight components for vehicle structures and propulsion systems, fuel cell components (bipolar plates and electrodes) and battery electrodes, and ultra-lightweight materials for use in solar sails. A major barrier to the successful use of carbon nanotubes in these components is the need for methods to economically produce pure carbon nanotubes in large enough quantities to not only evaluate their suitability for certain applications but also produce actual components. Most carbon nanotube synthesis methods, including the HiPCO (high pressure carbon monoxide) method developed by Smalley and others, employ metal catalysts that remain trapped in the final product. These catalyst impurities can affect nanotube properties and accelerate their decomposition. The development of techniques to remove most, if not all, of these impurities is essential to their successful use in practical applications. A new method has been developed at the NASA Glenn Research Center to purify gram-scale quantities of single wall carbon nanotubes. This method, a modification of a gas phase purification technique previously reported by Smalley and others, uses a combination of high-temperature oxidations and repeated extractions with nitric and hydrochloric acid. This improved procedure significantly reduces the amount of impurities (catalyst and nonnanotube forms of carbon) within the nanotubes, increasing their stability significantly. The onset of decomposition of the purified nanotubes (determined by thermal gravimetric analysis in air) is more than 300 C higher than that of the crude nanotubes. Transmission electron microscopy analysis of nanotubes purified by this method reveals near complete removal of iron catalyst particles. Analysis of the nanotubes using inductively coupled plasma spectroscopy revealed that the iron content of the nanotubes was reduced from 22.7 wt% in the crude nanotubes to less than 0.02 wt%. X-ray photoelectron spectroscopy revealed a decrease in iron content after purification as well as an increase in oxygen content due to the formation of carboxylic acid groups on the surface of the nanotubes. Nanotubes purified by this improved method can be readily dispersed in common organic solvents, in particular N,N-dimethylformamide, using prolonged ultrasonic treatment. These dispersions can then be used to incorporate single wall carbon nanotubes into polymer films.

Lebron, Marisabel↗

Satellite determination of nature and microstructure of atmospheric aerosols

A method is presented for the determination of aerosol physical parameters on the basis of scattered radiance measurements. The reconstruction of aerosol particle size distribution in the case of a spectral forward scattering method is considered, taking into account an inverse diffraction integral expression, a model cloud, the effect of lower wavenumber cut-off, the effect of higher wavenumber cut-off, the effect of spectral resolution, and the effects of multiple scattering and background noise. The determination of the aerosol complex refractive index with the aid of a minimization search method is also discussed.

Fymat, A. L.↗

Extinction-sedimentation inversion technique for measuring size distribution of artificial fogs

In measuring the size distribution of artificial fog particles, it is important that the natural state of the particles not be disturbed by the measuring device, such as occurs when samples are drawn through tubes. This paper describes a method for carrying out such a measurement by allowing the fog particles to settle in quiet air inside an enclosure through which traverses a parallel beam of light for measuring the optical depth as a function of time. An analytic function fit to the optical depth time decay curve can be directly inverted to yield the size distribution. Results of one such experiment performed on artificial fogs are shown as an example. The forwardscattering corrections to the measured extinction coefficient are also discussed with the aim of optimizing the experimental design so that the error due to forwardscattering is minimized.

Deepak, A.↗

Petrography and Mineralogy of Basalt Fragments in Apollo 17 Double Drive Tube 73002 and 73001

The ANGSA program was formed, in part, to examine previously unstudied Apollo 17 double drive tube samples. As a part of this project, the Moon United team received eight basalt fragments from 73001 and 73002. Here we present the petrography and mineralogy of these basalt fragments. Methods: We analyzed the petrography and mineralogy of the 73001/73002 particles using Electron Probe Microanalysis (EMPA) at both Carnegie Institution for Science and the Smithsonian Institution. We made quantitative compositional maps for each sample following the methods of [1,2].

S. N. Valencia↗

The 3-D world modeling with updating capability based on combinatorial geometry

A 3-D world modeling technique using range data is discribed. Range data quantify the distances from the sensor focal plane to the object surface, i.e., the 3-D coordinates of discrete points on the object surface are known. The approach proposed herein for 3-D world modeling is based on the Combinatorial Geometry (CG) method which is widely used in Monte Carlo particle transport calculations. First, each measured point on the object surface is surrounded by a small sphere with a radius determined by the range to that point. Then, the 3-D shapes of the visible surfaces are obtained by taking the (Boolean) union of all the spheres. The result is an unambiguous representation of the object's boundary surfaces. The pre-learned partial knowledge of the environment can be also represented using the CG Method with a relatively small amount of data. Using the CG type of representation, distances in desired directions to boundary surfaces of various objects are efficiently calculated. This feature is particularly useful for continuously verifying the world model against the data provided by a range finder, and for integrating range data from successive locations of the robot during motion. The efficiency of the proposed approach is illustrated by simulations of a spherical robot in a 3-D room in the presence of moving obstacles and inadequate prelearned partial knowledge of the environment.

Goldstein, M.↗