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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Snowfall Rate Retrieval Using Passive Microwave Measurements and Its Applications in Weather Forecast and Hydrology

(AMSU), Microwave Humidity Sounder (MHS) and Advance Technology Microwave Sounder (ATMS). ATMS is the follow‐on sensor to AMSU and MHS. Currently, an AMSU and MHS based land snowfall rate (SFR) product is running operationally at NOAA/NESDIS. Based on the AMSU/MHS SFR, an ATMS SFR algorithm has also been developed. The algorithm performs retrieval in three steps: snowfall detection, retrieval of cloud properties, and estimation of snow particle terminal velocity and snowfall rate. The snowfall detection component utilizes principal component analysis and a logistic regression model. It employs a combination of temperature and water vapor sounding channels to detect the scattering signal from falling snow and derives the probability of snowfall. Cloud properties are retrieved using an inversion method with an iteration algorithm and a two‐stream radiative transfer model. A method adopted to calculate snow particle terminal velocity. Finally, snowfall rate is computed by numerically solving a complex integral. The SFR products are being used mainly in two communities: hydrology and weather forecast. Global blended precipitation products traditionally do not include snowfall derived from satellites because such products were not available operationally in the past. The ATMS and AMSU/MHS SFR now provide the winter precipitation information for these blended precipitation products. Weather forecasters mainly rely on radar and station observations for snowfall forecast. The SFR products can fill in gaps where no conventional snowfall data are available to forecasters. The products can also be used to confirm radar and gauge snowfall data and increase forecasters' confidence in their prediction.

Meng, Huan↗
NASA NTRS1985

Geometry of the diffusive propagation region in the August 14, 1982 solar electron event

On August 14, 1982, relativistic electrons arrived promptly after an impulsive gamma ray flare, indicating that very little scattering was taking place in interplanetary space. By ignoring anisotropy data the time profile of the event is well described by interplanetary diffusion except for the derived particle injection time. This discrepancy provides independent evidence that the particles are diffusing in a volume close to the Sun rather than in interplanetary space. The flux at maximum method of determining the number of particles produced is still a good approximation when appropriately applied.

Evenson, P. A.↗
NASA NTRSDate not supplied

Transient Resonances in EMRIs: A Scalar Toy Model

We present calculations of the scalar self-force for a scalar point-charge following different $r\theta$-resonant geodesics in Kerr spacetime. Using this scalar model, we present two methods for constructing the self-force for resonant sources using two different phase parameterizations of the point-particle motion. We relate these two methods via a simple shifting relation that holds for both the scalar and gravitational self-forces. We discuss how this shifting relation allows us to compare self-force results both at and away from resonances. Finally, we provide new evidence that conservative scalar perturbations can drive the leading-order evolution of the orbital dynamics during resonances (i.e., the orbit-averaged rate of change of the Carter constant depends on contributions from the conservative self-force). This evidence conflicts with a commonly assumed conjecture that these integrated conservative contributions should vanish during both resonant and non-resonant motion.

Zachary Nasipak↗
NASA NTRS1992

Computational studies of metal-metal and metal-ligand interactions

The geometric structure of Cr(CO)6 is optimized at the modified coupled-pair functional (MCPF), single and double excitation coupled-cluster (CCSD) and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), and the force constants for the totally symmetric representation are determined. The geometry of Cr(CO)5 is partially optimized at the MCPF, CCSD and CCSD(T) levels of theory. Comparison with experimental data shows that the CCSD(T) method gives the best results for the structures and force constants, and that remaining errors are probably due to deficiencies in the one-particle basis sets used for CO. A detailed comparison of the properties of free CO is therefore given, at both the MCPF and CCSD/CCSD(T) levels of treatment, using a variety of basis sets. With very large one-particle basis sets, the SSCD(T) method gives excellent results for the bond distance, dipole moment and harmonic frequency of free CO. The total binding energies of Cr(CO)6 and Cr(CO)5 are also determined at the MCPF, CCSD and CCSD(T) levels of theory. The CCSD(T) method gives a much larger total binding energy than either the MCPF or CCSD methods. An analysis of the basis set superposition error (BSSE) at the MCPF level of treatment points out limitations in the one-particle basis used here and in a previous study. Calculations using larger basis sets reduced the BSSE, but the total binding energy of Cr(CO)6 is still significantly smaller than the experimental value, although the first CO bond dissociation energy of Cr(CO)6 is well described. An investigation of 3s3p correlation reveals only a small effect. The remaining discrepancy between the experimental and theoretical total binding energy of Cr(CO)6 is probably due to limitations in the one-particle basis, rather than limitations in the correlation treatment. In particular an additional d function and an f function on each C and O are needed to obtain quantitative results. This is underscored by the fact that even using a very large primitive se (1042 primitive functions contracted to 300 basis functions), the superposition error for the total binding energy of Cr(CO)6 is 22 kcal/mol at the MCPF level of treatment.

Barnes, Leslie A.↗
NASA NTRS1990

Determination of the optical thickness and effective particle radius of clouds from reflected solar radiation measurements. I - Theory

A method is presented for determining the optical thickness and effective particle radius of stratiform cloud layers from reflected solar radiation measurements. A detailed study is presented which shows that the cloud optical thickness (tau c) and effective particle radius (r/e/) of water clouds can be determined solely from reflection function measurements at 0.75 micron and 2.16 microns, provided tau c is not less than 4 and r(e) is not less than 6 microns. For optically thin clouds, the retrieval becomes ambiguous, resulting in two possible solutions for the effective radius and optical thickness. Adding a third channel near 1.65 micron does not improve the situation noticeably, whereas the addition of a channel near 3.70 microns reduces the ambiguity in deriving the effective radius. The effective radius determined by the above procedure corresponds to the droplet radius at some optical depth within the cloud layer.

Nakajima, Teruyuki↗
NASA NTRS1992

Adaptation of multidimensional group particle tracking and particle wall-boundary condition model to the FDNS code

A particulate two-phase flow CFD model was developed based on the FDNS code which is a pressure based predictor plus multi-corrector Navier-Stokes flow solver. Turbulence models with compressibility correction and the wall function models were employed as submodels. A finite-rate chemistry model was used for reacting flow simulation. For particulate two-phase flow simulations, a Eulerian-Lagrangian solution method using an efficient implicit particle trajectory integration scheme was developed in this study. Effects of particle-gas reaction and particle size change to agglomeration or fragmentation were not considered in this investigation. At the onset of the present study, a two-dimensional version of FDNS which had been modified to treat Lagrangian tracking of particles (FDNS-2DEL) had already been written and was operational. The FDNS-2DEL code was too slow for practical use, mainly because it had not been written in a form amenable to vectorization on the Cray, nor was the full three-dimensional form of FDNS utilized. The specific objective of this study was to reorder to calculations into long single arrays for automatic vectorization on the Cray and to implement the full three-dimensional version of FDNS to produce the FDNS-3DEL code. Since the FDNS-2DEL code was slow, a very limited number of test cases had been run with it. This study was also intended to increase the number of cases simulated to verify and improve, as necessary, the particle tracking methodology coded in FDNS.

Chen, Y. S.↗
NASA NTRS2015

Methodology for Producing a Uniform Distribution of UO2 in a Tungsten Matrix

Current work at NASA's Marshall Space Flight Center (MSFC) is focused on the development CERMET fuel materials for Nuclear Thermal Propulsion (NTP). The CERMETs consist of uranium dioxide (UO2) fuel particles embedded in a tungsten (W) metal matrix. Initial testing of W-UO2 samples fabricated from fine angular powders performed reasonably well, but suffered from significant fuel loss during repeated thermal cycling due to agglomeration of the UO2 (1). The blended powder mixtures resulted in a non-uniform dispersion of the UO2 particles in the tungsten matrix, which allows rapid vaporization of the interconnected UO2 from the sample edges into the bulk material. Also, the angular powders create areas of stress concentrations due to thermal expansion mismatch, which eventually cracks the tungsten matrix. Evenly coating spherical UO2 particles with chemical vapor deposited (CVD) tungsten prior to consolidation was previously demonstrated to provide improved performance. However, the CVD processing technology is expensive and not currently available. In order to reduce cost and enhance performance, a powder coating process has been developed at MSFC to produce a uniform distribution of the spherical UO2 particles in a tungsten matrix. The method involves utilization of a polyethylene binder during mixing which leads to fine tungsten powders clinging to the larger UO2 spherical particles. This process was developed using HfO2 as a surrogate for UO2. Enough powder was mixed to make 8 discs (2cm diameter x 8mm thickness) using spark plasma sintering. A uniaxial pressure of 50 MPa was used at four different temperatures (2 samples at each temperature). The first two samples were heated to 1400C and 1500C respectively for 5 minutes. Densities for these samples were less than 85% of theoretical, so the time at temperature was increased to 20 minutes for the remaining samples. The highest densities were achieved for the two samples sintered at 1700C (approx. 92% of theoretical). Scanning electron microscopy (SEM) of the mixed powders and the sintered samples along with energy dispersive x-ray analysis was obtained. The SEM of the powders clearly show the fine W powder adhered to the larger HfO2 particles and a uniform distribution of HfO2 particles in a tungsten matrix upon densification. Vicker's Microhardness testing was also performed on all samples using 0.5, 1.0 and 2.0 kg loads. Five indents were made at each load level. All indents were placed in the tungsten matrix to assist as a proxy in measuring densification. The highest hardness value was obtained for the 1700C specimens. The hardness average for these samples was 312.14 MPa. This powder processing method has been applied to W/UO2 powders with the SEM of the powders appearing similar to the W/HfO2 powder images.

Tucker, Dennis S.↗
NASA NTRS2009

Modeling Techniques for Evaluation the Effectiveness of Particle Damping in Turbomachinery

High power turbopumps are frequently used to supply propellants to the combustion chambers of rocket engines. Due to the high pressures and flow-rates required, turbopump components are subjected to harsh environments which include dynamic excitation due to random, sine, and acoustic vibration. Additionally, fluid-induced forces can couple with the dynamics of the structure resulting in flow induced instabilities (flutter). Structural response to these forms of excitation results in reduced fatigue life and increases the likelihood of an operational failure. Particle damping has been used successfully on vibration problems in the past by increasing the damping and therefore reducing the response to acceptable levels. Empirical methods have typically been employed to evaluate the performance of the particles in reducing the structural response. This report explores the use of finite element methods to estimate the effectiveness of particle damping in a typical non-rotating turbopump component. Axisymmetric harmonic models are used to estimate the increase in modal damping produced by the addition of particles in the cavity of an axisymmetric seal. Target modes of vibration are evaluated to quantify how the effective particle damping is altered by geometry changes in the seal design. A new method to predict the performance of particle dampers is developed and shown to provide more reasonable estimates of damping.

Ehrgott, R.↗
NASA NTRS2014

Combined Uncertainty and A-Posteriori Error Bound Estimates for General CFD Calculations: Theory and Software Implementation

This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.

CFD↗
NASA NTRS1984

Measurement of mass distribution of chemical species in aerosol particles

Aerosols may be generated through the nebulizing of solutions and the evaporation of their solvent, leaving the dry solute particles. Attention is presently given to a method for the direct determination of the masses of chemical species in individual aerosol particles on a continuous, real-time basis, using mass spectrometry. After the aerosol particles are introduced into the ion source of a quadrupole mass spectrometer, the particles impinge on a hot rhenium filament in the mass spectrometer's ion source. The resulting vapor plume is ionized by electron bombardment, and a pulse of ions is generated by each particle. The intensities of different masses in the ion pulses can then be measured by the mass spectrometer.

Sinha, M. P.↗
NASA NTRS1986

Preparation, analysis and release of simulated interplanetary grains into low Earth orbit

Astronomical observations which reflect the optical and dynamical properties of interstellar and interplanetary grains are the primary means of identifying the shape, size, and the chemistry of extraterrestrial grain materials. Except for recent samplings of extraterrestrial particles in near-Earth orbit and in the stratosphere observations were the only method of deducing the properties of extraterrestrial particles. In order to elucidate the detailed characteristics of observed dust, the observations must be compared with theoretical studies, some of which are discussed in this volume, or compared with terrestrial laboratory experiments. The formation and optical characterization of simulated interstellar and interplanetary dust with particular emphasis on studying the properties on irregularly shaped particles were discussed. Efforts to develop the techniques to allow dust experiments to be carried out in low-Earth orbit were discussed, thus extending the conditions under which dust experiments may be performed.

Stephens, J. R.↗
NASA NTRS1990

Potential, velocity, and density fields from sparse and noisy redshift-distance samples - Method

A method for recovering the three-dimensional potential, velocity, and density fields from large-scale redshift-distance samples is described. Galaxies are taken as tracers of the velocity field, not of the mass. The density field and the initial conditions are calculated using an iterative procedure that applies the no-vorticity assumption at an initial time and uses the Zel'dovich approximation to relate initial and final positions of particles on a grid. The method is tested using a cosmological N-body simulation 'observed' at the positions of real galaxies in a redshift-distance sample, taking into account their distance measurement errors. Malmquist bias and other systematic and statistical errors are extensively explored using both analytical techniques and Monte Carlo simulations.

Dekel, Avishai↗
NASA NTRS1992

Mean free paths of energetic particles at very large heliodistances (Pioneer 11 at 20 AU)

The parallel mean free path and the diffusion coefficient parallel to the magnetic field line are derived from magnetic field data at 20 AU to characterize heliospheric modulation and energetic-particle/magnetic-field interaction. The computational method of Moussas et al. (1975, 1982) is employed, and the values of the parallel mean free path are shown to be significantly larger than the values estimated in studies of up to 6 AU. The distance dependence of the parallel diffusion mean free path is found to follow a power law, and the diffusion coefficient dependence upon energy is determined by a constant mean free path and the velocity of the particle. The contribution of the diffusion coefficient perpendicular to the magnetic field is expected to dominate the radial diffusion coefficient of cosmic rays, although the contribution of the diffusion parallel to the field is important with respect to the small-scale structure of intensity gradients.

Moussas, X.↗
NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven↗
NASA NTRS1999

Investigation of Next-Generation Earth Radiation Budget Radiometry

The current effort addresses two issues important to the research conducted by the Thermal Radiation Group at Virginia Tech. The first research topic involves the development of a method which can properly model the diffraction of radiation as it enters an instrument aperture. The second topic involves the study of a potential next-generation space-borne radiometric instrument concept. Presented are multiple modeling efforts to describe the diffraction of monochromatic radiant energy passing through an aperture for use in the Monte-Carlo ray-trace environment. Described in detail is a deterministic model based upon Heisenberg's uncertainty principle and the particle theory of light. This method is applicable to either Fraunhofer or Fresnel diffraction situations, but is incapable of predicting the secondary fringes in a diffraction pattern. Also presented is a second diffraction model, based on the Huygens-Fresnel principle with a correcting obliquity factor. This model is useful for predicting Fraunhofer diffraction, and can predict the secondary fringes because it keeps track of phase. NASA is planning for the next-generation of instruments to follow CERES (Clouds and the Earth's Radiant Energy System), an instrument which measures components of the Earth's radiant energy budget in three spectral bands. A potential next-generation concept involves modification of the current CERES instrument to measure in a larger number of wavelength bands. This increased spectral partitioning would be achieved by the addition of filters and detectors to the current CERES geometry. The capacity of the CERES telescope to serve for this purpose is addressed in this thesis.

Coffey, Katherine L.↗
NASA NTRS1985

Pitch Angle Distributions of Solar Energetic Particles and the Local Scattering Properties of the Interplanetary Magnetic Field

An approximate solution of the Fokker-Planck equation containing pitch angle scattering and adiabatic focusing is discussed. For modest focusing effects the omnidirectional density obeys an ordinary diffusion equation with a modified diffusion coefficient. The anisotropic part of the distribution function is properly normalized and split into an even and an odd part. The even part is determined by the ratio between the scattering mean free path and the focusing length and by the degree of polarization of the magnetic field fluctuations. The odd part is determined by the deviation of the pitch angle scattering from isotropic scattering. The method supplies a powerful tool to obtain the local characteristics of pitch angle scattering. It is insensitive to long lasting solar injections and to moderate radial variations of the mean free path. The method is applied to solar particle events observed on Helio-1 and -2.

Beeck, J.↗
NASA NTRS1983

Retrieval of stratospheric aerosol size distribution from atmospheric extinction of solar radiation at two wavelengths

A parametric study is presented of the effects of the size distribution and composition of stratospheric aerosols on the value R, the ratio of the extinction of solar radiation at wavelength 0.45 micron to the extinction at wavelength 1.0 micron. The aerosol size distributions under study comprise nine analytical expressions, including most of the stratospheric aerosol models used by investigators. The aerosol compositions considered are the supercooled sulfuric-acid droplets, with different weight percentages of H2SO4 in the aerosol. It is found that whereas R is not very sensitive to either the composition or the radii limits of the stratospheric aerosols under consideration, it is quite sensitive to the value of the variable parameter governing the mode radius of the size distribution of aerosol particles. The results of the parametric study underlie the method proposed here of retrieving the size distribution of aerosol particles in the stratosphere from the experimental results of the extinction of solar radiation at two wavelengths.

Yue, G. K.↗