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 1,207 records · Page 67

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↗

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

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↗

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

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

Signatures and detection prospects for sub-GeV dark matter with superfluid helium

We explore the possibility of using superfluid helium for direct detection of sub-GeV dark matter (DM). We discuss the relevant phenomenology resulting from the scattering of an incident dark matter particle on a Helium nucleus. Rather than directly exciting quasi-particles, DM in this mass range will interact with a single He atom, triggering an atomic cascade which eventually also includes emission and thermalization of quasi-particles. We present in detail the analytical framework needed for modeling these processes and determining the resulting flux of quasi-particles. We propose a novel method for detecting this flux with modern force-sensitive devices, such as nanoelectro-mechanical system (NEMS) oscillators, and derive the sensitivity projections for a generic sub-GeV DM detection experiment using such sensors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Three-dimensional phase field sintering simulations accounting for the rigid-body motion of individual grains

Sintering is a widely used powder processing technique in industrial applications. During sintering, atoms migrate to decrease the energy of the system via two main mechanisms: coarsening and densification, both of which lead to significant morphological variation of the sintered microstructure. When simulating sintering dynamics, the phase-field method has been broadly utilized because of its convenience in tracking morphology evolution. When a large number of grains is involved, it is common to use the same order parameter to describe multiple grains that are not in direct contact with one another (in order to reduce the computational memory demands). However, with this treatment it is difficult to handle the rigid-body motion of individual grains during densification. In this work, an implementation scheme is introduced to overcome the challenge of calculating individual particle motion based on existing equations. It uses a grouping algorithm and sets a cutoff radius on each grain for calculating the particle velocity during densification. This method allows for the incorporation of the densification mechanism, which has been commonly ignored in previous work, into phase-field sintering models in three-dimensional simulations with a large number of particles/grains. Moreover, through combination with the smoothed boundary method, material properties of sintered microstructures, such as the effective diffusivity and Young’s modulus, can be calculated during the sintering processes.

36 MATERIALS SCIENCE↗

Determination of Seven Organic Impurities in FD&C Yellow No. 6 by Ultra-High-Performance Liquid Chromatography

The U.S. Food and Drug Administration (FDA) batch certifies FD&C Yellow No. 6 (Y6) to ensure that the color additive meets requirements published in the Code of Federal Regulations (CFR), including specifications for seven organic manufacturing impurities consisting of two intermediates, a reaction by-product, an impurity originating from an intermediate, and three subsidiary colors. An ultra-high-performance liquid chromatography (UHPLC) method was developed and validated for determining seven organic impurities in Y6. For comparison, the currently used high-performance liquid chromatography (HPLC) method was also validated. The new UHPLC method uses a 1.7 µm particle size biphenyl column with aqueous ammonium formate and methanol as eluants. Analytes are identified by comparing their retention times and UV-visible spectra to those of reference standards. Calibration is performed in the presence of the dye matrix and analyte levels are determined from their peak areas. UHPLC and HPLC validation studies obtained linear calibration curves and excellent values for limits of detection, limits of quantitation, recovery, repeatability, and intermediate precision for all analytes. Survey analyses of 30 samples using the UHPLC and HPLC methods yielded results that were consistent within experimental error; however, the UHPLC method provided more accurate results for two analytes that coelute using the HPLC method. The UHPLC method satisfies the accuracy and precision requisites for routine certification of Y6. The UHPLC method is faster and generates less waste than the HPLC method while delivering better sensitivity, separation, and accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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

Carbon abundance and silicate mineralogy of anhydrous interplanetary dust particles

We have studied nineteen anhydrous chondritic interplanetary dust particles (IDPs) using analytical electron microscopy. We have determined a method for quantitative light element EDX analysis of small particles and have applied these techniques to a group of IDPs. Our results show that some IDPs have significantly higher bulk carbon abundances than do carbonaceous chondrites. We have also identified a relationship between carbon abundance and silicate mineralogy in our set of anhydrous IDPs. In general, these particles are dominated by pyroxene, olivine, or a subequal mixture of olivine and pyroxene. The pyroxene-dominated IDPs have a higher carbon abundance than those dominated by olivines. Members of the mixed mineralogy IDPs can be grouped with either the pyroxene- or olivine-dominated particles based on their carbon abundance. The high carbon, pyroxene-dominated particles have primitive mineralogies and bulk compositions which show strong similarities to cometary dust particles. We believe that the lower carbon, olivine-dominated IDPs are probably derived from asteroids. Based on carbon abundances, the mixed-mineralogy group represents particles derived from either comets or asteroids. We believe that the high carbon, pyroxene-rich anhydrous IDPs are the best candidates for cometary dust.

Thomas, Kathie L.↗

A Physical Model to Estimate Snowfall over Land using AMSU-B Observations

In this study, we present an improved physical model to retrieve snowfall rate over land using brightness temperature observations from the National Oceanic and Atmospheric Administration's (NOAA) Advanced Microwave Sounder Unit-B (AMSU-B) at 89 GHz, 150 GHz, 183.3 +/- 1 GHz, 183.3 +/- 3 GHz, and 183.3 +/- 7 GHz. The retrieval model is applied to the New England blizzard of March 5, 2001 which deposited about 75 cm of snow over much of Vermont, New Hampshire, and northern New York. In this improved physical model, prior retrieval assumptions about snowflake shape, particle size distributions, environmental conditions, and optimization methodology have been updated. Here, single scattering parameters for snow particles are calculated with the Discrete-Dipole Approximation (DDA) method instead of assuming spherical shapes. Five different snow particle models (hexagonal columns, hexagonal plates, and three different kinds of aggregates) are considered. Snow particle size distributions are assumed to vary with air temperature and to follow aircraft measurements described by previous studies. Brightness temperatures at AMSU-B frequencies for the New England blizzard are calculated using these DDA calculated single scattering parameters and particle size distributions. The vertical profiles of pressure, temperature, relative humidity and hydrometeors are provided by MM5 model simulations. These profiles are treated as the a priori data base in the Bayesian retrieval algorithm. In algorithm applications to the blizzard data, calculated brightness temperatures associated with selected database profiles agree with AMSU-B observations to within about +/- 5 K at all five frequencies. Retrieved snowfall rates compare favorably with the near-concurrent National Weather Service (NWS) radar reflectivity measurements. The relationships between the NWS radar measured reflectivities Z(sub e) and retrieved snowfall rate R for a given snow particle model are derived by a histogram matching technique. All of these Z(sub e)-R relationships fall in the range of previously established Z(sub e)-R relationships for snowfall. This suggests that the current physical model developed in this study can reliably estimate the snowfall rate over land using the AMSU-B measured brightness temperatures.

Kim, Min-Jeong↗

DUAL-PARTICLE DOSEMETER BASED ON ORGANIC SCINTILLATOR

Abstract Traditionally available handheld dosemeters are generally sensitive to only one type of radiation: neutrons or photons. Some dosemeters also rely on very specific attenuation correlations between response and dose, are not scalable in size and multiple dosemeters are required to characterise mixed-particle fields. The research presented here serves as a proof-of-concept for a method to simultaneously measure dose rates from neutrons and photons using a particle discriminating organic scintillation detector without the need for spectral deconvolution. The method was compared with traditional instruments and to simulation. Isotopic photon dose rates measured with this method were within 4% of simulated truth, whereas fission spectrum neutron dose rates were measured within 21%. Measurements of dose rates from both particles agree with simulated truth better than traditional instruments. This new method allows for measurement of dose equivalent from both neutrons and photons with a single instrument and no reliance on spectral deconvolution.

Miller, C. A.↗

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using the Hottest Particles in the Universe to Probe Icy Solar System Worlds

We present results of our Phase 1 NIAC Study to determine the feasibility of developing a competitive, low cost, low power, low mass passive instrument to measure ice depth on outer planet ice moons, such as Europa, Ganymede, Callisto, and Enceladus. Indirect measurements indicate that liquid water oceans are likely present beneath the icy shells of such moons (see e.g.,the JPL press release "The Solar System and Beyond is Awash in Water"), which has important astrobiological implications. Determining the thickness of these ice shells is challenging given spacecraft SWaP (Size, Weight and Power) resources. The current approach uses a suite of instruments, including a high power, massive ice penetrating radar. The instrument under study, called PRIDE (Passive Radio Ice Depth Experiment) exploits a remarkable confluence between methods from the high energy particle physics and the search for extraterrestrial life within the solar system. PRIDE is a passive receiver of a naturally occurring radio frequency (RF) signal generated by interactions of deep penetrating Extremely High Energy (> 10^18 eV) cosmic ray neutrinos. It could measure ice thickness directly, and at a significant savings to spacecraft resources. At RF frequencies the transparency of modeled Europan ice is up to many km, so an RF sensor in orbit can observe neutrino interactions to great depths, and thereby probe the thickness of the ice layer.

Exploration↗

Calorimetry with deep learning: particle simulation and reconstruction for collider physics

Using detailed simulations of calorimeter showers as training data, we investigate the use of deep learning algorithms for the simulation and reconstruction of single isolated particles produced in high-energy physics collisions. We train neural networks on single-particle shower data at the calorimeter-cell level, and show significant improvements for simulation and reconstruction when using these networks compared to methods which rely on currently-used state-of-the-art algorithms. We define two models: an end-to-end reconstruction network which performs simultaneous particle identification and energy regression of particles when given calorimeter shower data, and a generative network which can provide reasonable modeling of calorimeter showers for different particle types at specified angles and energies. We investigate the optimization of our models with hyperparameter scans. Furthermore, we demonstrate the applicability of the reconstruction model to shower inputs from other detector geometries, specifically ATLAS-like and CMS-like geometries. These networks can serve as fast and computationally light methods for particle shower simulation and reconstruction for current and future experiments at particle colliders.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Understanding Aitken Mode Aerosol Variability over the Southern Ocean and Antarctica: Insights from Cloud Condensation Nuclei Data

Aitken mode aerosol particles play an important role influencing cloud properties and sustenance, acting as a reservoir of potential cloud condensation nuclei against precipitation scavenging. However, there is limited data on Aitken mode aerosols. In this study, we develop a method to estimate Aitken mode aerosol concentrations and size distribution using cloud condensation nuclei measurements (CCN) and κ-Köhler theory. The performance of this method is evaluated using scanning mobility particle sizer (SMPS) data from recent field campaigns to demonstrate its skills and applicability. The method reasonably estimates Aitken- and accumulation-mode aerosol concentrations, achieving correlations of 0.7–0.9 with only modest biases (mean fractional bias within ±23% for Aitken mode and ±34% for accumulation-mode). This method is further applied to measurements collected over the Southern Ocean and Antarctica in recent years from multiple platforms, including ground sites, aircraft, and ships, to derive Aitken and accumulation-mode aerosol concentrations. Using the derived data, we examine the seasonal cycle, latitudinal variations, and vertical distribution of aerosols. Aitken mode aerosol concentrations are elevated over the Southern Ocean and Antarctica during the austral summer similar to the accumulation mode. In the austral summer, the free troposphere has more Aitken mode aerosols and fewer accumulation mode aerosols than the boundary layer, and thus likely serves as an important source of cloud-forming aerosol while also diluting the accumulation mode.

Kang, Litai [University of Washington] (ORCID:0000↗

Applications of Quantum Theory of Atomic and Molecular Scattering to Problems in Hypersonic Flow

The general status of a grant to investigate the applications of quantum theory in atomic and molecular scattering problems in hypersonic flow is summarized. Abstracts of five articles and eleven full-length articles published or submitted for publication are included as attachments. The following topics are addressed in these articles: fragmentation of heavy ions (HZE particles); parameterization of absorption cross sections; light ion transport; emission of light fragments as an indicator of equilibrated populations; quantum mechanical, optical model methods for calculating cross sections for particle fragmentation by hydrogen; evaluation of NUCFRG2, the semi-empirical nuclear fragmentation database; investigation of the single- and double-ionization of He by proton and anti-proton collisions; Bose-Einstein condensation of nuclei; and a liquid drop model in HZE particle fragmentation by hydrogen.

Malik, F. Bary↗