Engineering Papers⌕ Search

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.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Theory of molecular rate processes in the presence of intense laser radiation

The present paper deals with the influence of intense laser radiation on gas-phase molecular rate processes. Representations of the radiation field, the particle system, and the interaction involving these two entities are discussed from a general rather than abstract point of view. The theoretical methods applied are outlined, and the formalism employed is illustrated by application to a variety of specific processes. Quantum mechanical and semiclassical treatments of representative atom-atom and atom-diatom collision processes in the presence of a field are examined, and examples of bound-continuum processes and heterogeneous catalysis are discussed within the framework of both quantum-mechanical and semiclassical theories.

George, T. F.↗
NASA NTRS2000

How Well Can Aerosol Measurements from the Terra Morning Polar Orbiting Satellite Represent the Daily Aerosol Abundance and Properties?

The Terra mission, launched at the dawn of 1999, and Aqua mission to be launched soon, will possess innovative measurements of the aerosol daily spatial distribution, distinguish between dust, smoke and regional pollution and measure aerosol radiative forcing of climate. Their polar orbit gives daily global coverage, however measurements are acquired at specific time of the day. To what degree can present measurements from Terra taken between 10:00 and 11:30 AM local time, represent the daily average aerosol forcing of climate? Here we answer this question using 7 years of data from the distributed ground based 50-70 instrument Aerosol Robotic Network (AERONET) This (AERONET) half a million measurement data set shows that Terra aerosol measurements represent the daily average values within 5%. The excellent representation is found for large dust particles or small aerosol particles from Fires or regional pollution and for any range of the optical thickness, a measure of the amount of aerosol in the atmosphere.

Kaufman, Y. J.↗
NASA NTRS2007

Validation and Determination of Ice Water Content - Radar Reflectivity Relationships during CRYSTAL-FACE: Flight Requirements for Future Comparisons

In order for clouds to be more accurately represented in global circulation models (GCM), there is need for improved understanding of the properties of ice such as the total water in ice clouds, called ice water content (IWC), ice particle sizes and their shapes. Improved representation of clouds in models will enable GCMs to better predict for example, how changes in emissions of pollutants affect cloud formation and evolution, upper tropospheric water vapor, and the radiative budget of the atmosphere that is crucial for climate change studies. An extensive cloud measurement campaign called CRYSTAL-FACE was conducted during Summer 2002 using instrumented aircraft and a variety of instruments to measure properties of ice clouds. This paper deals with the measurement of IWC using the Harvard water vapor and total water instruments on the NASA WB-57 high-altitude aircraft. The IWC is measured directly by these instruments at the altitude of the WB-57, and it is compared with remote measurements from the Goddard Cloud Radar System (CRS) on the NASA ER-2. CRS measures vertical profiles of radar reflectivity from which IWC can be estimated at the WB-57 altitude. The IWC measurements obtained from the Harvard instruments and CRS were found to be within 20-30% of each other. Part of this difference was attributed to errors associated with comparing two measurements that are not collocated in time an space since both aircraft were not in identical locations. This study provides some credibility to the Harvard and CRS-derived IWC measurements that are in general difficult to validate except through consistency checks using different measurement approaches.

Sayres, D. S.↗
NASA NTRSDate not supplied

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗
DOE OSTI2021

Monte Carlo Hauser-Feshbach computer code system to model nuclear reactions: YAHFC

A computer program framework, YAHFC, to model low-energy nuclear reactions is presented. The framework allows for reactions with incident particles ranging from protons/neutrons to alphas and is designed to address reactions that ultimately lead to the formation of compound nuclear systems that then decay statistically as outlined in concepts of Hauser and Feshbach. Additionally, instead of a reaction, it is also possible to model the decay of a nuclear system with an initial excitation and population. The code models nuclear decays with a Monte Carlo process that tracks the decay of each state. This allows for an exact representation of the spectra for all emitted particles in each of the final exit channels and the possibility of generating reaction data for simulation purposes. The program is interfaced with the optical model code system FRESCOX to calculate transmission coefficients as well as the effects of coupled channels and other direct excitations via the distorted wave Born approximation (DWBA). Modules are included to account for nuclear processes such as width corrections, pre-equilibrium emission, and fission. The program is controlled by a series of input commands and while a set of input parameters exists for each projectile and target, the input commands allow for complete control over each input parameter. Extensive data files are produced and a program is provided that converts YAHFC data files into nuclear data library entries in the generalized nuclear data structure (GNDS).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Automated characterization of spatial and dynamical heterogeneity in supercooled liquids via implementation of machine learning

Abstract A computational approach by an implementation of the principle component analysis (PCA) with K -means and Gaussian mixture (GM) clustering methods from machine learning algorithms to identify structural and dynamical heterogeneities of supercooled liquids is developed. In this method, a collection of the average weighted coordination numbers ( W C N s ‾ ) of particles calculated from particles’ positions are used as an order parameter to build a low-dimensional representation of feature (structural) space for K -means clustering to sort the particles in the system into few meso-states using PCA. Nano-domains or aggregated clusters are also formed in configurational (real) space from a direct mapping using associated meso-states’ particle identities with some misclassified interfacial particles. These classification uncertainties can be improved by a co-learning strategy which utilizes the probabilistic GM clustering and the information transfer between the structural space and configurational space iteratively until convergence. A final classification of meso-states in structural space and domains in configurational space are stable over long times and measured to have dynamical heterogeneities. Armed with such a classification protocol, various studies over the thermodynamic and dynamical properties of these domains indicate that the observed heterogeneity is the result of liquid–liquid phase separation after quenching to a supercooled state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Fracture mechanics approach to TRISO fuel particle failure analysis

Weibull stress-based methods for failure probability assessment have been developed and analyzed to assess the integrity of tristructural isotropic (TRISO) fuel particles during fuel life cycles and accident operating conditions. While simple, these methods entail a number of drawbacks when stress concentrates near crack tips, including finite element mesh size dependency when the Weibull stress is averaged over the finite element domain. Fracture mechanics approaches involving the use of interaction integrals eliminate this lack of mesh convergence and produce consistent fracture predictions. In this work, we use an interaction integral approach to computing stress intensity factors for a crack in the inner pyrolytic carbon layer perpendicular to the silicon carbide layer, which is simplified representation of a failure mode in TRISO particles. Further, the interface between these two TRISO layers has been shown to become porous, which we simulate by considering a transition of mechanical properties over such porous length, i.e. the layers are modeled as a functionally graded material. Aspects such as porosity and thermal and irradiation eigenstrains are considered in computing the stress intensity factor from a fracture mechanics approach and compared with the known Weibull stress failure approach. The methodology introduced in this paper enables a more general fracture probability assessment in TRISO particles and eliminates the need to identify best suited parameters when using local or averaged stress-based failure criterion. Finally, the numerical sensitivity studies show how parameters such as the porous transition zone length, the material stiffness, and creep affect the probability of TRISO fuel particle failure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Helicity is a topological invariant of massless particles: C = - 2 h

There is an elementary but indispensable relationship between the topology and geometry of massive particles. The geometric spin s is related to the topological dimension of the internal space V by dim⁡ V = 2⁢s + 1. This breaks down for massless particles, which are geometrically characterized by their helicity ℎ, but all have 1D internal spaces. We show that a subtler relation exists between the topology and geometry of massless particles. Wave functions of massless particles are sections of nontrivial line bundles over the light cone whose topology is completely characterized by their first Chern number C. We prove that in general C = -2⁢ℎ. In doing so, we also exhibit a method of generating all massless bundle representations via an Abelian group structure of massless particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Zwanziger’s pairwise little group on the celestial sphere

We generalize Zwanziger’s pairwise little group to include a boost subgroup. We do so by working in the celestial sphere representation of scattering amplitudes. We propose that due to late time soft photon and graviton exchanges, matter particles in the asymptotic states in massless QED and gravity transform under the Poincaré group with an additional pair of pairwise celestial representations for each pair of matter particles. We demonstrate that the massless abelian and gravitational exponentiation theorems are consistent with the proposed pairwise Poincaré transformation properties. For massless QED we demonstrate that our results are consistent with the effects of the Faddeev-Kulish dressing and the abelian exponentiation theorem for celestial amplitudes found in arXiv:2012.04208. We discuss electric and magnetic charges simultaneously as it is especially natural to do so in this formalism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2013

Automatic Parameter Tuning for the Morpheus Vehicle Using Particle Swarm Optimization

A high fidelity simulation using a PC based Trick framework has been developed for Johnson Space Center's Morpheus test bed flight vehicle. There is an iterative development loop of refining and testing the hardware, refining the software, comparing the software simulation to hardware performance and adjusting either or both the hardware and the simulation to extract the best performance from the hardware as well as the most realistic representation of the hardware from the software. A Particle Swarm Optimization (PSO) based technique has been developed that increases speed and accuracy of the iterative development cycle. Parameters in software can be automatically tuned to make the simulation match real world subsystem data from test flights. Special considerations for scale, linearity, discontinuities, can be all but ignored with this technique, allowing fast turnaround both for simulation tune up to match hardware changes as well as during the test and validation phase to help identify hardware issues. Software models with insufficient control authority to match hardware test data can be immediately identified and using this technique requires very little to no specialized knowledge of optimization, freeing model developers to concentrate on spacecraft engineering. Integration of the PSO into the Morpheus development cycle will be discussed as well as a case study highlighting the tool's effectiveness.

Birge, B.↗
NASA NTRS2005

Representation of Vegetation and Other Nonerodible Elements in Aeolian Shear Stress Partitioning Models for Predicting Transport Threshold

The presence of nonerodible elements is well understood to be a reducing factor for soil erosion by wind, but the limits of its protection of the surface and erosion threshold prediction are complicated by the varying geometry, spatial organization, and density of the elements. The predictive capabilities of the most recent models for estimating wind driven particle fluxes are reduced because of the poor representation of the effectiveness of vegetation to reduce wind erosion. Two approaches have been taken to account for roughness effects on sediment transport thresholds. Marticorena and Bergametti (1995) in their dust emission model parameterize the effect of roughness on threshold with the assumption that there is a relationship between roughness density and the aerodynamic roughness length of a surface. Raupach et al. (1993) offer a different approach based on physical modeling of wake development behind individual roughness elements and the partition of the surface stress and the total stress over a roughened surface. A comparison between the models shows the partitioning approach to be a good framework to explain the effect of roughness on entrainment of sediment by wind. Both models provided very good agreement for wind tunnel experiments using solid objects on a nonerodible surface. However, the Marticorena and Bergametti (1995) approach displays a scaling dependency when the difference between the roughness length of the surface and the overall roughness length is too great, while the Raupach et al. (1993) model's predictions perform better owing to the incorporation of the roughness geometry and the alterations to the flow they can cause.

King, James↗
DOE OSTI2021

Assessment of filtered cameras for quantitative 2D analysis of divertor conditions during detachment in JET L-mode plasmas

We report that estimates for 2D distributions of electron temperature, T e , electron density, n e , and atomic deuterium density, n o , in the JET divertor volume have been inferred from deuterium Balmer line intensity ratios obtained from tomographic reconstructions of divertor camera measurements. This enables also investigation of ionization, S ion , and recombination, S rec , rates. The analysis shows a decrease of T e to 0.5–1.0 eV throughout the outer divertor during detachment in low-confinement (L-mode) plasmas. Simultaneously, the high-n e region and the n 0 distribution in the outer divertor are observed to elongate and shift from the outer strike point towards the X-point. The observations are in qualitative agreement and follow the same sequence with modelling predictions of EDGE2D-EIRENE simulations of a density scan. While the method was found to provide good representation of the evolution of volumetric recombination during detachment, in agreement with the simulations, the movement of the ionization front upstream could not be followed due to lack of spatial overlap between the ionization region and the necessary emission distributions. Consequently, the representation of the ionization conditions and the particle balance in the detached outer divertor are compromised.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Neutrino many-body flavor evolution: The full Hamiltonian

We study neutrino flavor evolution in the quantum many-body approach using the full neutrino-neutrino Hamiltonian, including the usually neglected terms that mediate nonforward scattering processes. Working in the occupation number representation with plane waves as single-particle states, we explore the time evolution of simple initial states with up to N = 10 neutrinos. We discuss the time evolution of the Loschmidt echo, one body flavor and kinetic observables, and the one-body entanglement entropy. For the small systems considered, we observe “thermalization” of both flavor and momentum degrees of freedom on comparable time scales, with results converging towards expectation values computed within a microcanonical ensemble. We also observe that the inclusion of nonforward processes generates a faster flavor evolution compared to the one induced by the truncated (forward) Hamiltonian. Published by the American Physical Society 2024

Cirigliano, Vincenzo (ORCID:000000029056754X)↗
DOE OSTI2020

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2014

Investigating the Sensitivity of Surface-Level Nitrate Seasonality in Antarctica to Primary Sources Using a Global Model

Determining the sources of total nitrate (TNIT NO-3 + HNO3) reaching Antarctica is a long-standing challenge. Here we analyze the monthly sensitivity of surface-level TNIT in Antarctica to primary sources using a global 3-D chemical transport model, GEOS-Chem, and its adjoint. Modeled seasonal variation of TNIT concentrations shows good agreement with several measurement studies, given that the lack of post-depositional processing in the model leads to an expected underestimate of maximum values in November through January. Remote NOx sources have the greatest impact May-July, during when the model background concentrations are sensitive to NOx emissions from fossil fuel combustion, soil, and lightning originating from 25S to 65S. In this season, NOx is transported to Antarctica as TNIT, which is formed above continental source regions at an altitude of 5-11 km. In other seasons, more NOx is transported as a reservoir species (e.g., peroxyacetyl nitrate, PAN) through the free troposphere, transforming into TNIT within a cone of influence that extends to 35S and above 4 km altitude. Photolysis of PAN over Antarctica is the main driver of modeled NOx seasonality. Stratospheric production and loss of tracers are relatively unimportant in monthly sensitivities in GEOS-Chem, driving only a few percent of surface level variability of TNIT. A small peak concentration in August is captured by the model, although some measured values in August fall outside the range of simulated concentrations. Modifications to the model to represent sedimentation of polar stratospheric clouds (PSCs) lead to increased surface level August TNIT concentrations. However, this simple representation does not explicitly account for PSC particle deposition or disappearance of the tropopause in the middle of winter, and thus the influence of stratospheric nitrate sources estimated in this study is likely a lower bound.

stratosphere↗
NASA NTRS2019

Electromagnetic Scattering by Discrete Random Media. III: The Vector Radiative Transfer Equation

A vector radiative transfer equation with an additional source term typical of dense media is obtained. The analysis includes (i) the derivation of an integral equation for the correlation matrix of the exciting field coefficients accounting for the correlation between the particles, (ii) the derivation of an integral representation for the specific coherency dyadicin terms of this matrix, and (iii) the simplification of the integral equation for the correlation matrix and of the integral representation for the specific coherency dyadic by employing a series of approximations which are characteristic of sparse media.

Adrian Doicu↗