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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Shock and statistical acceleration of energetic particles in the interplanetary medium

Definite evidence for particle acceleration in the solar wind came around a decade ago. Two likely sources are known to exist: particles may be accelerated by the turbulence resulting from the superposition of Alfven and Magnetosonic waves (Statistical Acceleration) or they may be accelerated directly at shock fronts formed by the interaction of fast and slow solar wind (CIR's) or by traveling shocks due to sporadic coronal mass ejections. Naurally both mechanisms may be operative. In this work the acceleration problem was tackled numerically using Helios 1 and 2 data to create a realistic representation of the Heliospheric plasma. Two 24 hour samples were used: one where there are only wave like fluctuations of the field (Day 90 Helios 1) and another with a shock present in it (Day 92 of Helios 2) both in 1976 during the STIP 2 interval. Transport coefficients in energy space have been calculated for particles injected in each sample and the effect of the shock studied in detail.

Valdes-Galicia, J. F.↗
DOE OSTI2021

Predicting Frigid Mixed-Phase Clouds for Pristine Coastal Antarctica

Supercooled water is common in the clouds near coastal Antarctica and occasionally occurs at temperatures at or below -30°C. Yet the ice physics in most regional and global numerical models will glaciate out these clouds. This presents a challenge for the simulation of highly supercooled clouds that were observed at McMurdo, Antarctica during the Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE) project during 2015–2017. The polar optimized version of the Weather Research and Forecasting model (Polar WRF) with the recently developed two-moment P3 microphysics scheme was used to simulate observed supercooled liquid water cases during March and November 2016. Nudging of the simulations to observed rawinsonde profiles and Antarctic automatic weather station observations provided increased realism and much greater cloud water amounts. Furthermore, sensitivity tests that adjust the ice physics for extremely low ice nucleating particle (INP) concentrations decrease cloud ice and increases the cloud liquid water closer to observed amounts. In these tests, a liquid layer near cloud top is simulated, in agreement with observations. Accurate representation of INP concentrations appears to be critical for the simulation of coastal Antarctic clouds.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

Research in computer science

Synopses are given for NASA supported work in computer science at the University of Virginia. Some areas of research include: error seeding as a testing method; knowledge representation for engineering design; analysis of faults in a multi-version software experiment; implementation of a parallel programming environment; two computer graphics systems for visualization of pressure distribution and convective density particles; task decomposition for multiple robot arms; vectorized incomplete conjugate gradient; and iterative methods for solving linear equations on the Flex/32.

Ortega, J. M.↗
NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

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.

Discrete Element Method↗
DOE OSTI2026

AI-assisted object condensation clustering for calorimeter shower reconstruction at CLAS12

Several nuclear physics studies using the CLAS12 detector rely on the accurate reconstruction of neutrons and photons from its forward angle calorimeter system. These studies often place restrictive cuts when measuring neutral particles due to an overabundance of false clusters created by the existing calorimeter reconstruction software. In this work, we present a new AI approach to clustering CLAS12 calorimeter hits based on the object condensation framework. The model learns a latent representation of the full detector topology using GravNet layers, serving as the positional encoding for an event’s calorimeter hits which are processed by a Transformer encoder. This unique structure allows the model to contextualize local and long range information, improving its performance. Evaluated on one million simulated $e^-$ $+$ $p$ collision events, our method significantly improves cluster trustworthiness: the fraction of reliable neutron clusters, increasing from 8.88% to 30.73%, and photon clusters, increasing from 51.07% to 64.73%. In conclusion, our study also marks the first application of AI clustering techniques for hodoscopic detectors, showing potential for usage in many other experiments.

Calorimeters↗
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↗
DOE OSTI2020

Derivation of K-matrix reaction theory in a discrete basis formalism

The usual derivations of the S and K matrices for two-particle reactions proceed through the Lippmann–Schwinger equation with formal definitions of the incoming and outgoing scattering states. Here we present a simpler alternative derivation that is carried out completely in the Hamiltonian representation, using a discrete basis of configurations for the scattering channels as well as the quasi-bound configurations of the combined fragments. We use matrix algebra to derive an explicit expression for the K matrix in terms of the Hamiltonian of the internal states of the compound system and the coupling between the channels and the internal states. The formula for the K matrix includes explicitly a real dispersive shift matrix to the internal Hamiltonian that is easily computed in the formalism. That expression is applied to derive the usual form of the S matrix as a sum over poles in the complex energy plane. Some extensions and limitations of the discrete-basis Hamiltonian formalism are discussed in the concluding remarks and in the Appendix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Blueprints for Training Information Bottlenecks for Collider Analyses

Dimensionality reduction is a crucial aspect of data analysis in high energy physics, even if accompanied by information loss. Several methods, including histogram- and kernel-based analyses, are only computationally feasible for low-dimensional data. Furthermore, simulation models used in HEP can often only be validated for low-dimensional data. We provide several blueprints for using machine learning to create low-dimensional data representations (continuous event variables and discrete classification labels) for use in signal discovery and parameter estimation tasks. We also describe how to design the learned representation to facilitate a) searches with unknown model parameters and b) validation of simulation models in data control regions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

An alternative representation of the ice canopy for calculating microwave brightness temperatures over a thunderstorm

Passive microwave brightness temperatures (TBs) at 92 and 183 GHz from an aircraft thunderstorm overflight are compared with values calculated from radar-derived hydrometeor profiles and a modified proximity sounding. Two methods for modeling particles in the ice canopy are contrasted. The first is a 'traditional' approach employing Marshall-Palmer ice spheres. The second, or 'alternative', method partitions 20 percent of the ice water content into a Marshall-Palmer component for graupel and hail, and 80 percent into a modified gamma spherical particle size distribution function representing ice crystals. Results from the alternative approach are superior to those from the traditional method in the anvil and mature convective core. In the decaying convective region, the traditional approach yields better agreement with observed magnitudes. Neither method, however, matches the geometry of the observed TB depression associated with the decaying convective core. This is likely due to the presence of graupel, which is not detected as a special signature in radar reflectivity, but does diminish TBs through scattering. Brightness temperatures at the relatively high microwave frequencies considered are shown to be very sensitive to the ice-particle size distribution.

Muller, Bradley M.↗
NASA NTRS1995

Bounce-Averaged Hamiltonian for Charged Particles in an Axisymmetric but Nondipolar Model Magnetosphere

In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M, J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. The L value of a field line is (by definition) inversely proportional to the flux enclosed by the corresponding magnetic shell of equatorial radius r(sub 0), and the L value at the neutral line (r(sub 0) = b) is denoted L*. The azimuthal coordinate phi measures magnetic local time. The best functional representation found for the normalized difference (L(exp 3)a(exp 3)/mu(sub E))(B(sub m) - B(sub 0)) between mirror-point field B(sub m) and equatorial field B(sub 0) along any field line is a 5-term expansion in powers (2/3 through 6/3) of the quantity X equivalent to (La/mu(sub E))(exp 1/2)K, where K equivalent to (J(exp 2)/8m(sub 0)M)(exp 1/2) is an adiabatically conserved quantity independent of particle energy, m(sub 0) is the rest mass of the particle, and a is the radius of the Earth. This functional form is motivated by results for limiting cases in which particles mirror very near and very far from the magnetic equator. Expansion coefficients corresponding to various powers of X are obtained from least squares fits to numerically computed results for X as a function of L and B(sub m). These are accurately expressible as fourth-order polynomials in (r(sub 0)/b)(exp 3), hence indirectly as functions of L/L* = 3La/2b. This representation, which leads (except for a manageably small region of parameter space) to better than 1% accuracy in the specification of B(sub m) as a function of K and L, allows bounce-averaged guiding center simulations to be performed without actually tracing the bounce motions of individual particles. Bounce-averaged drifts L' (meridional) and phi' (azimuthal) are proportional to derivatives of the Hamiltonian H (sum of kinetic and potential energies) with respect to phi and L, respectively. Our formulation thus provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Schulz, Michael↗
DOE OSTI2026

Prognostic simulations of mixed-phase clouds with model AC-1D v1.0: the impact of aerosol types and freezing parameterizations on ice crystal budgets

Mixed-phase clouds at high latitudes contribute to the uncertainty in predicting cloud feedbacks and climate sensitivity, mainly due to the complexity of microphysical processes that influence the partitioning between the supercooled liquid and ice phases, and hence, cloud radiative effects on regional scales. Particularly in Arctic mixed-phase clouds, the activation of ice-nucleating particles (INPs) from various aerosol populations remains a leading source of uncertainty. We developed an aerosol-cloud one-dimensional (AC-1D) model, which provides a novel framework to prognostically treat INP and ice crystal budgets while explicitly accounting for polydisperse and multicomponent aerosol that activate INPs following different freezing parameterizations. The AC-1D model is informed by large-eddy simulations to probe the impact of INP representation on predicted ice crystal number concentrations (N i ) and ice crystal budgets in mixed-phase Arctic stratus. We apply three immersion freezing (IMF) parameterizations, two time-independent (singular) and one time-dependent (classical nucleation theory), to predict the evolution of the INP reservoir and resulting ice crystal budget from polydisperse mineral dust, organic (humic-like substances), and sea spray aerosol particle size distributions. Our analysis focuses on how variations in aerosol number concentration and cloud system parameters such as cloud cooling rate, cloud-top entrainment rate, and ice crystal fall speed influence the INP reservoir and ice crystal budgets. Furthermore, this study investigates the competitive ice nucleation dynamics in mixed aerosol environments and provides a process-level quantification of the INP budget terms, which directly controls ice crystal budgets. For all studied case scenarios, the aerosol types and associated particle size distributions significantly impact INP and N i , and the choice between a time-dependent and a singular freezing description yields orders-of-magnitude differences in the predicted INP and N i over the 10 h simulation time, reflecting typical cloud lifetimes. Our results show that the influence of cloud cooling, INP entrainment, and sedimentation varies significantly depending on the chosen freezing parameterization. These findings underscore the critical need for robust IMF parameterizations and precise cloud system observations to enhance the accuracy of models in predicting mixed-phase cloud structure and evolution.

Arctic clouds↗
DOE OSTI2022

Nondestructive neutron imaging diagnosis of acidic gas reduction catalyst after 400-Hour operation in natural gas furnace

Residential natural gas furnaces are the primary space-heating devices in US homes, leading to substantial environmental impact caused by the acidic components in the furnace combustion gases. Here, to experimentally demonstrate acidic gas reduction in a furnace, a monolithic catalyst was fabricated and was called the AGR. A commercially-available condensing furnace was retrofitted with the AGR, and a 400-hour reliability and durability test was conducted. The results showed that the AGR significantly reduced acidic gases in the flue gas and produced condensate with neutral pH. Challenges were also revealed: inappropriate condensate drainage caused incomplete combustion and amorphous carbon deposits. To nondestructively survey the internal state of the AGR, neutron computed tomography (NCT) was employed to produce spatially resolved 2D and 3D representations of the 2-L, aged AGR component. The NCT results confirm the integrity of the AGR component, which consists of two blocks, without deformation or damage to AGR channels. The distribution of the particle accumulation in the middle of the top block was visibly heavier than the entrance and exit of the top block. The representative cross-section views revealed significant aggregation in the central region but not at the rim. The spatially resolved details provide in-depth diagnosis, evaluation, and understanding of the AGR. The insights can enable new AGR designs that realize a uniform and self-cleaning flow pattern, alleviate significant particle aggregation, and thus enhance AGR-enabled furnace performance. The neutron imaging method demonstrates the good potential that can diagnose faults and improve the design, optimization, and production processes of novel catalysts and other components or systems with heavy metal shell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Three-dimensional plasmoid-mediated reconnection and turbulence in Hall magnetohydrodynamics

Plasmoid instability accelerates reconnection in collisional plasmas by transforming a laminar reconnection layer into numerous plasmoids connected by secondary current sheets in two dimensions (2D) and by fostering self-generated turbulent reconnection in three dimensions (3D). In large-scale astrophysical and space systems, plasmoid instability likely initiates in the collisional regime but may transition into the collisionless regime as the fragmentation of the current sheet progresses toward kinetic scales. Hall magnetohydrodynamics (MHD) models are widely regarded as a simplified yet effective representation of the transition from collisional to collisionless reconnection. However, plasmoid instability in 2D Hall MHD simulations often leads to a single-X-line reconnection configuration, which significantly differs from fully kinetic particle-in-cell simulation results. This study shows that single-X-line reconnection is less likely to occur in 3D compared to 2D. Moreover, depending on the Lundquist number and the ratio between the system size and the kinetic scale, Hall MHD can also realize 3D self-generated turbulent reconnection. We analyze the features of the self-generated turbulent state, including the energy power spectra and the scale dependence of turbulent eddy anisotropy.

Physics↗
DOE OSTI2024

Direct solution of Minkowski-space Bethe-Salpeter equation in the massive Wick-Cutkosky model

Here in order to solve the Bethe-Salpeter equation (BSE) in the Minkowski space, we first introduce the Nakanishi integral representations of the Bethe-Salpeter amplitude and the Bethe-Salpeter wave function. We then derive the explicit integral equations for the corresponding spectral functions from the BSE for states of two scalar particles bound by a scalar-particle exchange interaction, where the propagators of constituents are allowed to be fully dressed. These integral equations are subsequently solved numerically in the variation of the Wick-Cutkosky model with massive exchange particles, where an algorithm of adaptive mash grid is proposed. The equations and algorithm we develop here serve as the foundation of Minkowski-space formulation of BSE for bound states of fermions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Rapid SACR Observations of Convection at Bankhead National Forest (RAPID) Field Campaign Report

Improving our representation of convective cell processes requires better quantification of convective clouds throughout their entire life cycle. This includes gaining a clearer understanding of the controls on key convective cloud properties, such as updraft intensity, particle size distributions, rainfall rates, and hydrometeor species. Our inability to improve convective cloud process modeling stems, in part, from a limited understanding of convective cell properties, particularly given how rapidly these storms evolve. This lack of detailed observations in the most intense and organized convective storms is especially significant, as large errors remain in representing these clouds, which are critical for severe weather prediction and Earth system model performance. Cloud and precipitation radars are essential tools for studying cloud microphysics and dynamics, particularly in deeper convective clouds. The recent U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s third Mobile Facility (AMF3) Bankhead National Forest (BNF) deployment provides a unique opportunity to investigate important land–atmosphere interactions, as well as the environmental controls on deep convective cloud processes, in a location favorable for frequent convection. We operate the X/Ka-band Scanning ARM Cloud Radar (X/Ka SACR) using a scan strategy optimized to capture these rapidly evolving clouds and their properties, thereby improving studies of deep convective cloud processes. This effort is strengthened by a complementary and coordinated partnership with ongoing university and multi-agency radar activities collocated in north Alabama—a unique opportunity to examine clouds and precipitation from a lifetime-centric perspective.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Denoising diffusion models with geometry adaptation for high fidelity calorimeter simulation

Simulation is crucial for all aspects of collider data analysis, but the available computing budget in the High Luminosity LHC era will be severely constrained. Generative machine learning models may act as surrogates to replace physics-based full simulation of particle detectors, and diffusion models have recently emerged as the state of the art for other generative tasks. Here, we introduce CaloDiffusion, a denoising diffusion model trained on the public CaloChallenge datasets to generate calorimeter showers. Our algorithm employs 3D cylindrical convolutions, which take advantage of symmetries of the underlying data representation. To handle irregular detector geometries, we augment the diffusion model with a new geometry latent mapping (GLaM) layer to learn forward and reverse transformations to a regular geometry that is suitable for cylindrical convolutions. The showers generated by our approach are nearly indistinguishable from the full simulation, as measured by several different metrics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗