SEARCH · Engineering Papers

Results for “Particle representation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1987

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

NASA NTRS1990

Interactive numerical flow visualization using stream surfaces

Particle traces and ribbons are often used to depict the structure of three-dimensional flowfields, but images produced using these models can be ambiguous. Stream surfaces offer a more visually intuitive method for the depiction of flowfields, but interactive response is needed to allow the user to place surfaces which reveal the essential features of a given flowfield. FLORA, a software package which supports the interactive calculation and display of stream surfaces on silicon graphics workstations, is described. Alternative methods for the integration of particle traces are examined, and calculation through computational space is found to provide rapid results with accuracy adequate for most purposes. Rapid calculation of traces is teamed with progressive refinement of appoximated surfaces. An initial approximation provides immediate user feedback, and subsequent improvement of the surface ensures that the final image is an accurate representation of the flowfield.

Hultquist, J. P. M.↗

NASA NTRS2015

Rao-Blackwellization for Adaptive Gaussian Sum Nonlinear Model Propagation

When dealing with imperfect data and general models of dynamic systems, the best estimate is always sought in the presence of uncertainty or unknown parameters. In many cases, as the first attempt, the Extended Kalman filter (EKF) provides sufficient solutions to handling issues arising from nonlinear and non-Gaussian estimation problems. But these issues may lead unacceptable performance and even divergence. In order to accurately capture the nonlinearities of most real-world dynamic systems, advanced filtering methods have been created to reduce filter divergence while enhancing performance. Approaches, such as Gaussian sum filtering, grid based Bayesian methods and particle filters are well-known examples of advanced methods used to represent and recursively reproduce an approximation to the state probability density function (pdf). Some of these filtering methods were conceptually developed years before their widespread uses were realized. Advanced nonlinear filtering methods currently benefit from the computing advancements in computational speeds, memory, and parallel processing. Grid based methods, multiple-model approaches and Gaussian sum filtering are numerical solutions that take advantage of different state coordinates or multiple-model methods that reduced the amount of approximations used. Choosing an efficient grid is very difficult for multi-dimensional state spaces, and oftentimes expensive computations must be done at each point. For the original Gaussian sum filter, a weighted sum of Gaussian density functions approximates the pdf but suffers at the update step for the individual component weight selections. In order to improve upon the original Gaussian sum filter, Ref. [2] introduces a weight update approach at the filter propagation stage instead of the measurement update stage. This weight update is performed by minimizing the integral square difference between the true forecast pdf and its Gaussian sum approximation. By adaptively updating each component weight during the nonlinear propagation stage an approximation of the true pdf can be successfully reconstructed. Particle filtering (PF) methods have gained popularity recently for solving nonlinear estimation problems due to their straightforward approach and the processing capabilities mentioned above. The basic concept behind PF is to represent any pdf as a set of random samples. As the number of samples increases, they will theoretically converge to the exact, equivalent representation of the desired pdf. When the estimated qth moment is needed, the samples are used for its construction allowing further analysis of the pdf characteristics. However, filter performance deteriorates as the dimension of the state vector increases. To overcome this problem Ref. [5] applies a marginalization technique for PF methods, decreasing complexity of the system to one linear and another nonlinear state estimation problem. The marginalization theory was originally developed by Rao and Blackwell independently. According to Ref. [6] it improves any given estimator under every convex loss function. The improvement comes from calculating a conditional expected value, often involving integrating out a supportive statistic. In other words, Rao-Blackwellization allows for smaller but separate computations to be carried out while reaching the main objective of the estimator. In the case of improving an estimator's variance, any supporting statistic can be removed and its variance determined. Next, any other information that dependents on the supporting statistic is found along with its respective variance. A new approach is developed here by utilizing the strengths of the adaptive Gaussian sum propagation in Ref. [2] and a marginalization approach used for PF methods found in Ref. [7]. In the following sections a modified filtering approach is presented based on a special state-space model within nonlinear systems to reduce the dimensionality of the optimization problem in Ref. [2]. First, the adaptive Gaussian sum propagation is explained and then the new marginalized adaptive Gaussian sum propagation is derived. Finally, an example simulation is presented.

state estimation↗

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 NTRS1992

The structure and energetics of Cr(CO)6 and Cr(CO)5

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. 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. Calculations using larger basis sets reduce 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. In the largest basis set, the total CO binding energy of Cr(CO)6 is estimated to be 140 kcal/mol at the CCSD(T) level of theory, or about 86 percent of the experimental value. The remaining discrepancy between the experimental and theoretical value 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 set (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 NTRS1967

Rotating emission ring in binary systems.

Representation of rotating gaseous emission rings in eclipsing binary systems by periodic orbits around more massive component in restricted three-body problem

PERIODIC ORBIT↗

NASA NTRS1979

Corotating energetic particle and fast plasma streams in the inner and outer solar system: Radial dependence and energy spectra

Interplanetary acceleration processes are shown as the most plausible explanation for the observed corotating energetic particle events. The relation between the energetic particle events and the properties of the high speed solar wind streams observed at 1 AU were investigated along with the form of the energy spectrum of the corotating energetic particle streams and its variation with respect to CIR boundaries and with radial distance. It is shown that: (1) at 1 AU a correlation exists between the j particle intensity and the solar wind velocity measured during the rising part of the event, of the form I is proportional to exp (V sub sw/V sub o); and (2) the energy spectra from .5 to 20 MeV are well represented by an exponential in momentum of the form dJ/dP = C exp (-P/P sub o). This representation is found to apply from .45 AU to beyond 5 AU. The variation of P sub o with respect to the CIR boundaries was studied using a method of superposed epoch analysis. It is shown that at 1 AU the spectrum remains constant during the first two days and then progressively flattens; between 3-4AU.

Vanhollebeke, M. A. I.↗

NASA NTRS1978

Characteristics of the unsteady motion on transversely sheared mean flows

An explicit representation for the unsteady motion on a transversely sheared mean flow is obtained which corresponds to the gustline motion on a uniform mean flow. The important features of this motion are discussed. It is shown that its velocity, pressure and vorticity are all induced by a certain disturbance field that is a linear combination of the vorticity and particle-displacement fields and is everywhere frozen in the mean flow. The general ideas are illustrated by considering the scattering of a gust by a half-plane embedded in a shear flow.

Goldstein, M. E.↗

NASA NTRS2004

Plasma Dispersion Function for the Kappa Distribution

The plasma dispersion function is computed for a homogeneous isotropic plasma in which the particle velocities are distributed according to a Kappa distribution. An ordinary differential equation is derived for the plasma dispersion function and it is shown that the solution can be written in terms of Gauss' hypergeometric function. Using the extensive theory of the hypergeometric function, various mathematical properties of the plasma dispersion function are derived including symmetry relations, series expansions, integral representations, and closed form expressions for integer and half-integer values of K.

Podesta, John J.↗

NASA NTRS2016

Understanding Cirrus Ice Crystal Number Variability for Different Heterogeneous Ice Nucleation Spectra

Along with minimizing parameter uncertainty, understanding the cause of temporal and spatial variability of the nucleated ice crystal number, Ni, is key to improving the representation of cirrus clouds in climate models. To this end, sensitivities of Ni to input variables like aerosol number and diameter provide valuable information about nucleation regime and efficiency for a given model formulation. Here we use the adjoint model of the adjoint of a cirrus formation parameterization (Barahona and Nenes, 2009b) to understand Ni variability for various ice-nucleating particle (INP) spectra. Inputs are generated with the Community Atmosphere Model version 5, and simulations are done with a theoretically derived spectrum, an empirical lab-based spectrum and two field-based empirical spectra that differ in the nucleation threshold for black carbon particles and in the active site density for dust. The magnitude and sign of Ni sensitivity to insoluble aerosol number can be directly linked to nucleation regime and efficiency of various INP. The lab-based spectrum calculates much higher INP efficiencies than field-based ones, which reveals a disparity in aerosol surface properties. Ni sensitivity to temperature tends to be low, due to the compensating effects of temperature on INP spectrum parameters; this low temperature sensitivity regime has been experimentally reported before but never deconstructed as done here.

Cirrus↗

NASA NTRS2016

Scattering and Diffraction of Electromagnetic Radiation: An Effective Probe to Material Structure

Scattered electromagnetic waves from material bodies of different forms contain, in an intricate way, precise information on the intrinsic, geometrical and physical properties of the objects. Scattering theories, ever deepening, aim to provide dependable interpretation and prediction to the complicated interaction of electromagnetic radiation with matter. There are well-established multiple-scattering formulations based on classical electromagnetic theories. An example is the Generalized Multi-particle Mie-solution (GMM), which has recently been extended to a special version ̶ the GMM-PA approach, applicable to finite periodic arrays consisting of a huge number (e.g., >>106) of identical scattering centers [1]. The framework of the GMM-PA is nearly complete. When the size of the constituent unit scatterers becomes considerably small in comparison with incident wavelength, an appropriate array of such small element volumes may well be a satisfactory representation of a material entity having an arbitrary structure. X-ray diffraction is a powerful characterization tool used in a variety of scientific and technical fields, including material science. A diffraction pattern is nothing more than the spatial distribution of scattered intensity, determined by the distribution of scattering matter by way of its Fourier transform [1]. Since all linear dimensions entered into Maxwell's equations are normalized by wavelength, an analogy exists between optical and X-ray diffraction patterns. A large set of optical diffraction patterns experimentally obtained can be found in the literature [e.g., 2,3]. Theoretical results from the GMM-PA have been scrutinized using a large collection of publically accessible, experimentally obtained Fraunhofer diffraction patterns. As far as characteristic structures of the patterns are concerned, theoretical and experimental results are in uniform agreement; no exception has been found so far. Closely connected with the spatial distribution of scattered intensities are cross sections, such as for extinction, scattering, absorption, and radiation pressure, as a critical type of key quantity addressed in most theoretical and experimental studies of radiative scattering. Cross sections predicted from different scattering theories are supposed to be in general agreement. For objects of irregular shape, the GMM-PA solutions can be compared with the highly flexible Discrete Dipole Approximation (DDA) [4,5] when dividing a target to no more than ~106 unit cells. Also, there are different ways to calculate the cross sections in the GMM-PA, providing an additional means to examine the accuracy of the numerical solutions and to unveil potential issues concerning the theoretical formulations and numerical aspects. To solve multiple scattering by an assembly of material volumes through classical theories such as the GMM-PA, the radiative properties of the component scatterers, the complex refractive index in particular, must be provided as input parameters. When using a PA to characterize a material body, this involves the use of an adequate theoretical tool, an effective medium theory, to connect Maxwell's phenomenogical theory with the atomistic theory of matter. In the atomic theory, one regards matter as composed of interacting particles (atoms and molecules) embedded in the vacuum [6]. However, the radiative properties of atomic-scaled particles are known to be substantially different from bulk materials. Intensive research efforts in the fields of cluster science and nanoscience attempt to bridge the gap between bulk and atom and to understand the transition from classical to quantum physics. The GMM-PA calculations, which place virtually no restriction on the component-particle size, might help to gain certain insight into the transition.

Xu, Yu-Lin↗

NASA NTRS2013

A Simple Model of Global Aerosol Indirect Effects

Most estimates of the global mean indirect effect of anthropogenic aerosol on the Earth's energy balance are from simulations by global models of the aerosol lifecycle coupled with global models of clouds and the hydrologic cycle. Extremely simple models have been developed for integrated assessment models, but lack the flexibility to distinguish between primary and secondary sources of aerosol. Here a simple but more physically based model expresses the aerosol indirect effect (AIE) using analytic representations of cloud and aerosol distributions and processes. Although the simple model is able to produce estimates of AIEs that are comparable to those from some global aerosol models using the same global mean aerosol properties, the estimates by the simple model are sensitive to preindustrial cloud condensation nuclei concentration, preindustrial accumulation mode radius, width of the accumulation mode, size of primary particles, cloud thickness, primary and secondary anthropogenic emissions, the fraction of the secondary anthropogenic emissions that accumulates on the coarse mode, the fraction of the secondary mass that forms new particles, and the sensitivity of liquid water path to droplet number concentration. Estimates of present-day AIEs as low as 5 W/sq m and as high as 0.3 W/sq m are obtained for plausible sets of parameter values. Estimates are surprisingly linear in emissions. The estimates depend on parameter values in ways that are consistent with results from detailed global aerosol-climate simulation models, which adds to understanding of the dependence on AIE uncertainty on uncertainty in parameter values.

Clouds↗

NASA NTRS1976

New fine structure cooling rate

One of the dominant electron cooling processes in the ionosphere is caused by electron impact induced fine structure transitions among the ground state levels of atomic oxygen. This fine structure cooling rate is based on theoretical cross sections. Recent advances in the numerical cross section determinations to include polarization effects and more accurate representations of the atomic target result in new lower values. These cross sections are employed in this paper to derive a new fine structure cooling rate which is between 40% and 60% of the currently used rate. A new generalized formula is presented for the cooling rate (from which the fine structure cooling rate is derived), valid for arbitrary mass and temperature difference of the colliding particles and arbitrary inelastic energy difference.

Hoegy, W. R.↗

NASA NTRS1976

New fine structure cooling rate

Hoegy, W. R.↗

NASA NTRS1983

The nonlocal theory of periodic density drift instabilities

The nonlocal theory of electrostatic density drift instabilities is developed for an arbitrary, periodic variation in plasma density. Linear Vlasov theory is applied to Fourier series representations of the plasma density, and the resulting equations are solved numerically. The growth rates of the universal drift instability are compared by using local and nonlocal theories. It is demonstrated that the local theory is inadequate if the nonlocal eigenfunctions of the electrostatic modes have significant amplitude over spatial regions wider than the regions of large plasma density gradients. Nonlocal effects can enhance wave-particle diffusion of a triangular-shaped variation and reduce wave-particle dissipation of a square-shaped irregularity.

Bernhardt, P. A.↗

NASA NTRS1992

r-process nucleosynthesis in the high-entropy supernova bubble

We show that the high-temperature, high-entropy evacuated region outside the recent neutron star in a core-collapse supernova may be an ideal r-process site. In this high-entropy environment it is possible that most nucleons are in the form of free neutrons or bound into alpha particles. Thus, there can be many neutrons per seed nucleus even though the material is not particularly neutron rich. The predicted amount of r-process material ejected per event from this environment agrees well with that required by simple galactic evolution arguments. When averaged over regions of different neutron excess in the supernova ejecta, the calculated r-process abundance curve can give a good representation of the solar-system r-process abundances as long as the entropy per baryon is sufficiently high. Neutrino irradiation may aid in smoothing the final abundance distribution.

Meyer, B. S.↗

NASA NTRS2006

NASA Tech Briefs, November 2006

Topics include: Simulator for Testing Spacecraft Separation Devices; Apparatus for Hot Impact Testing of Material Specimens; Instrument for Aircraft-Icing and Cloud-Physics Measurements; Advances in Measurement of Skin Friction in Airflow; Improved Apparatus for Testing Monoball Bearings; High-Speed Laser Scanner Maps a Surface in Three Dimensions; Electro-Optical Imaging Fourier-Transform Spectrometer; Infrared Instrument for Detecting Hydrogen Fires; Modified Coaxial Probe Feeds for Layered Antennas; Detecting Negative Obstacles by Use of Radar; Cryogenic Pound Circuits for Cryogenic Sapphire Oscillators; PixelLearn; New Software for Predicting Charging of Spacecraft; Conversion Between Osculating and Mean Orbital Elements; Generating a 2D Representation of a Complex Data Structure; Making Activated Carbon by Wet Pressurized Pyrolysis; Composite Solid Electrolyte Containing Li+- Conducting Fibers; Electrically Conductive Anodized Aluminum Surfaces; Rapid-Chill Cryogenic Coaxial Direct-Acting Solenoid Valve; Variable-Tension-Cord Suspension/Vibration- Isolation System; Techniques for Connecting Superconducting Thin Films; Versatile Friction Stir Welding/Friction Plug Welding System; Thermal Spore Exposure Vessels; Enumerating Spore-Forming Bacteria Airborne with Particles; Miniature Oxidizer Ionizer for a Fuel Cell; Miniature Ion-Array Spectrometer; Promoted-Combustion Chamber with Induction Heating Coil; Miniature Ion-Mobility Spectrometer; Mixed-Salt/Ester Electrolytes for Low-Temperature Li+ Cells; Miniature Free-Space Electrostatic Ion Thrusters; Miniature Bipolar Electrostatic Ion Thruster; Holographic Plossl Retroreflectors; Miniature Electrostatic Ion Thruster With Magnet; Using Apex To Construct CPM-GOMS Models; Sequence Detection for PPM Optical Communication With ISI; Algorithm for Rapid Searching Among Star-Catalog Entries; Expectation-Based Control of Noise and Chaos; Radio Heating of Lunar Soil to Release Gases; Using Electrostriction to Manipulate Ullage in Microgravity; Equations for Scoring Rules When Data Are Missing; Insulating Material for Next-Generation Spacecraft; and Pseudorandom Switching for Adding Radar to the AFF Sensor.

Source record↗

NASA NTRS2005

Nonlinear Attitude Filtering Methods

The extended Kalman filter (EKF) is the workhorse of real-time spacecraft attitude estimation. Since the group SO3 of rotation matrices has dimension three, most attitude determination EKFs use lower- dimensional attitude parameterizations than the nine-parameter attitude matrix itself. The fact that all three- parameter representations of SO3 are singular or discontinuous for certain attitudes has led to extended discussions of constraints and attitude representations in EKFs. The most successful EKF uses a nonsingular parameterization for the global attitude, which necessarily has more than three parameters, while employing a three-component representation for the attitude errors. This filter has become known as the Multiplicative Extended Kalman Filter. These issues are now well understood, however, and the EKF has performed admirably in the vast majority of attitude determination applications. Nevertheless, poor performance or even divergence arising from the linearization implicit in the EKF has led to the development of nonlinear filters, most recently sigma point or unscented filters and particle filters.

Markley, F. Landis↗