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 181 records · Page 10

Lagrangian methods in plasma dynamics. II - Construction of Lagrangians for plasmas

Consideration of the construction of suitable Lagrangian functions for the dynamics of a cold plasma in such a way as to retain the relativistically covariant formalism. In one method, this is achieved by the introduction of a set of three variables which label the world lines of the particles. A second method results in a Clebsch-type representation. Sturrock's relativistic Lagrangian and Low's hot plasma Lagrangian are also briefly discussed in the context of the present work. The behavior of the canonical stress tensor is considered. The applicability of many of the general results in part I (Dougherty, 1970) is ensured by establishing the existence of the Lagrangian function.

Dougherty, J. P.↗

Evaluation of Aerosol Mixing State Classes in the GISS Modele-matrix Climate Model Using Single-particle Mass Spectrometry Measurements

Aerosol particles in the atmosphere are composed of multiple chemical species. The aerosol mixing state, which describes how chemical species are mixed at the single-particle level, provides critical information on microphysical characteristics that determine the interaction of aerosols with the climate system. The evaluation of mixing state has become the next challenge. This study uses aerosol time-of-flight mass spectrometry (ATOFMS) data and compares the results to those of the Goddard Institute for Space Studies modelE-MATRIX (Multiconfiguration Aerosol TRacker of mIXing state) model, a global climate model that includes a detailed aerosol microphysical scheme. We use data from field campaigns that examine a variety of air mass regimens (urban, rural, and maritime). At all locations, polluted areas in California (Riverside, La Jolla, and Long Beach), a remote location in the Sierra Nevada Mountains (Sugar Pine) and observations from Jeju (South Korea), the majority of aerosol species are internally mixed. Coarse aerosol particles, those above 1 micron, are typically aged, such as coated dust or reacted sea-salt particles. Particles below 1 micron contain large fractions of organic material, internally-mixed with sulfate and black carbon, and few external mixtures. We conclude that observations taken over multiple weeks characterize typical air mass types at a given location well; however, due to the instrumentation, we could not evaluate mass budgets. These results represent the first detailed comparison of single-particle mixing states in a global climate model with real-time single-particle mass spectrometry data, an important step in improving the representation of mixing state in global climate models.

climate models↗

Analysis of the Global Microwave Polarization Data of Clouds

Two NASA microwave radiometers, the satellite-borne GPM (Global Precipitation Measurement) Microwave Imager (GMI) and the aircraft-borne CoSMIR (the Conical Scanning Millimeter-wave Imaging Radiometer), measure vertically- and horizontally-polarized microwaves emitted by cloud particles and the Earth below, providing unique information on ice crystal properties in clouds. Their data reveal that non-spherical ice crystals are common and they fall in a preferred horizontally aligned orientation in convective and optically thick clouds especially near cloud top. A bin (particle-size-resolving) microphysical model with an ice crystal shape representation is developed to simulate the evolution of ice crystal properties (i.e., size, shape and orientation), where the radiation effect on microphysics (REM) is taken into account. Since REM represents the effect of all (e.g., both infrared and solar) radiation on ice crystal temperature, it relies upon the ice crystal proprieties that determine how an ice crystal receives radiation. Definitely, REM is different from the radiative effects that cause sensitivity at the microwave frequencies in the GPM and CoSMIR observations. Model results show that horizontally-oriented ice crystals grow faster than vertically-oriented (or spherical) ones due to REM. When both horizontally- and vertically-oriented ice crystals coexist in an air parcel, the model results show that the former grow by vapor deposition whereas the latter shrink by sublimation and disappear eventually. These modeling results are supported by the GMI data and the CoSMIR observations from MC3E (Midlatitude Continental Convective Clouds Experiment) and OLYMPEX (Olympic Mountains Experiment) on the prevalence of horizontally-oriented ice crystals. Moreover, the REM-induced precipitation explains the CloudSat observations of rare thin clouds in the tropical mid-troposphere as well as the common diamond dust in the high latitudes.

Xiping Zeng↗

An Evaluation of Biomass Burning Aerosol Mass, Extinction, and Size Distribution in Geos Using Observations From Camp2EX

Biomass burning aerosol impacts aspects of the atmosphere and Earth system through radiative forcing, serving as cloud condensation nuclei, and air quality. Despite its importance, the representation of biomass burning aerosol is not always accurate in numerical weather prediction and climate models or reanalysis products. Using observations collected as part of the Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex) in August through October of 2019, aerosol concentration and optical properties are evaluated within the Goddard Earth Observing System (GEOS) and its underlying aerosol module, GOCART. In the operational configuration, GEOS assimilates aerosol optical depth observations at 550 nm to constrain aerosol fields. Particularly for biomass burning aerosol, without the assimilation of aerosol optical depth, aerosol extinction is underestimated compared to observations collected in the Philippines region during the CAMP2Ex campaign. The assimilation process adds excessive amounts of carbon to account for the underestimated extinction, resulting in positive biases in the mass of black and organic carbon, especially within the boundary layer, relative to in situ observations from the Langley Aerosol Research Group Experiment. Counteracting this, GEOS is deficient in sulphate and nitrate aerosol just above the boundary layer. Aside from aerosol mass, extinction within GEOS is a function of ambient relative humidity and an assumed particle size distribution. The relationship between dry and ambient extinction in GEOS reveals that hygroscopic growth is too aggressive within the model for biomass burning aerosol. An additional concern lies in the assumed particle size distribution for GEOS, which has a mode radius that is too small for organic carbon. Variability in the observed particle size distribution for biomass burning aerosol within a single flight also illuminates the fact that a single assumed particle size distribution is not sufficient and that for a proper representation, a more advanced aerosol module with GEOS may be necessary.

Aerosol Mass↗

Evaluation of GEOS Aerosols Using CAMP2Ex Observations

Biomass burning aerosol impacts aspects of the atmosphere and Earth system through direct and semi-direct effects, as well as influencing air quality. Despite its importance, the representation of biomass burning aerosol is not always accurate in numerical weather prediction and climate models or reanalysis products. Using observations collected as part of the Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex) in August through October of 2019, aerosol concentration and optical properties are evaluated within the Goddard Earth Observing System (GEOS) and its underlying aerosol module, GOCART. In the operational configuration, GEOS assimilates aerosol optical depth observations at 550 nm from AERONET and MODIS to constrain aerosol fields. Particularly for biomass burning aerosol, without the assimilation of aerosol optical depth, aerosol extinction is underestimated compared to observations collected in the Philippines region during the CAMP2Ex campaign. The assimilation process adds excessive amounts of carbon to account for the underestimated extinction, resulting in positive biases in the mass of black and organic carbon, especially within the boundary layer, relative to in situ observations from the Langley Aerosol Research Group Experiment. Counteracting this, GEOS is deficient in sulfate and nitrate aerosol just above the boundary layer. Aerosol extinction within GEOS is a function of the mass of different aerosol species, the ambient relative humidity, the assumed spectral optical properties, and particle size distribution per species. The relationship between dry and ambient extinction in GEOS reveals that hygroscopic growth is too high within the model for biomass burning aerosol. An additional concern lies in the assumed particle size distribution for GEOS, which has a single mode radius that is too small for organic carbon. Variability in the observed particle size distribution for biomass burning aerosol within a single flight also illuminates the fact that a single assumed particle size distribution is not sufficient and that for a proper representation, a more advanced aerosol module within GEOS may be necessary.

Allison Collow↗

A General Materials Data Science Framework for Quantitative 2D Analysis of Particle Growth from Image Sequences

Abstract Phase transformations are a challenging problem in materials science, which lead to changes in properties and may impact performance of material systems in various applications. We introduce a general framework for the analysis of particle growth kinetics by utilizing concepts from machine learning and graph theory. As a model system, we use image sequences of atomic force microscopy showing the crystallization of an amorphous fluoroelastomer film. To identify crystalline particles in an amorphous matrix and track the temporal evolution of the particle dispersion, we have developed quantitative methods of 2D analysis. 700 image sequences were analyzed using a neural network architecture, achieving 0.97 pixel-wise classification accuracy as a measure of the correctly classified pixels. The growth kinetics of isolated and impinged particles were tracked throughout time using these image sequences. The relationship between image sequences and spatiotemporal graph representations was explored to identify the proximity of crystallites from each other. The framework enables the analysis of all image sequences without the requirement of sampling for specific particles or timesteps for various materials systems.

36 MATERIALS SCIENCE↗

Models Transport Saharan Dust too Low in the Atmosphere: a Comparison of the MetUM and CAMS Forecasts with Observations

We investigate the dust forecasts from two operational global atmospheric models in comparison with in situ and remote sensing measurements obtained during the AERosol properties – Dust (AER-D) field campaign. Airborne elastic backscatter lidar measurements were performed on board the Facility for Airborne Atmospheric Measurements during August 2015 over the eastern Atlantic, and they permitted us to characterise the dust vertical distribution in detail, offering insights on transport from the Sahara. They were complemented with airborne in situ measurements of dust size distribution and optical properties, as well as datasets from the Cloud–Aerosol Transport System (CATS) spaceborne lidar and the Moderate Resolution Imaging Spectroradiometer (MODIS). We compare the airborne and spaceborne datasets to operational predictions obtained from the Met Office Unified Model (MetUM) and the Copernicus Atmosphere Monitoring Service (CAMS). The dust aerosol optical depth predictions from the models are generally in agreement with the observations but display a low bias. However, the predicted vertical distribution places the dust lower in the atmosphere than highlighted in our observations. This is particularly noticeable for the MetUM, which does not transport coarse dust high enough in the atmosphere or far enough away from the source.We also found that both model forecasts underpredict coarse-mode dust and at times overpredict fine-mode dust, but as they are fine-tuned to represent the observed optical depth, the fine mode is set to compensate for the underestimation of the coarse mode. As aerosol–cloud interactions are dependent on particle numbers rather than on the optical properties, this behaviour is likely to affect their correct representation. This leads us to propose an augmentation of the set of aerosol observations available on a global scale for constraining models, with a better focus on the vertical distribution and on the particle size distribution. Mineral dust is a major component of the climate system; therefore, it is important to work towards improving how models reproduce its properties and transport mechanisms.

Debbie OSullivan↗

Graph neural network for neutrino physics event reconstruction

Liquid argon time projection chamber (LArTPC) detector technology offers a wealth of high-resolution information on particle interactions, and leveraging that information to its full potential requires sophisticated automated reconstruction techniques. Here, this article describes NUGRAPH 2, a graph neural network for low-level reconstruction of simulated neutrino interactions in a LArTPC detector. Simulated neutrino interactions in the MicroBooNE detector geometry are described as heterogeneous graphs, with energy depositions on each detector plane forming nodes on planar subgraphs. The network utilizes a multihead attention message-passing mechanism to perform background filtering and semantic labeling on these graph nodes, identifying those associated with the primary physics interaction with 98.0% efficiency and labeling them according to particle type with 94.9% efficiency. The network operates directly on detector observables across multiple two-dimensional representations but utilizes a three-dimensional-context-aware mechanism to encourage consistency between these representations. Model inference takes 0.12 s / event on a CPU and 0.005 s / event batched on a GPU. This architecture is designed to be a general-purpose solution for particle reconstruction in neutrino physics, with the potential for deployment across a broad range of detector technologies, and offers a core convolution engine that can be leveraged for a variety of tasks beyond the two described in this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Gas-Granular Simulation Framework for Spacecraft Landing Plume-Surface Interaction and Debris Transport Analysis

The Gas-Granular Flow Solver (GGFS) multi-phase flow computational framework has been developed to enable simulations of particle flows complex extra-terrestrial regolith materials. Particle flows of interest include the damage of unprepared landing sites from rocket plume impingement on Moon, Mars, and asteroids. The flow solver implements an Eulerian-Eulerian two-fluid model with fluid representation of the gas phase and granular phase to avoid the need to model billions of particle interactions. The granular phase is modeled as an Eulerian fluid with constituent physics closure models derived from first-principle Discrete Element Model (DEM) particle interaction simulations that capture the complex, non-linear granular particle interaction effects. Granular phase constituent models have been developed and integrated that address the complex, non-linear granular material mechanics complexities resulting from both: the irregular, jagged particle shapes and poly-disperse mixture effects encountered in extra-terrestrial regolith, with lunar regolith as the extreme. The GGFS capabilities are being integrated into a proven NASA plume-surface interaction and debris transport simulation framework featuring the Loci/CHEM CFD program and Debris Transport Analysis (DTA) post-processing tools for applications in robotic and human Moon and Mars lander development. Integration of the three simulation tool components. Loci/CHEM, GGFS, and DTA, into a coordinated simulation framework will enable time-accurate spacecraft landing simulations that account for the alteration of the landing surface through plume-induced cratering and the resulting redirection of plume impingement flow and debris transport. Initial implementation of this simulation framework and application examples will be presented.

Liever, Peter A.↗

Masked Particle Modeling on Sets: Towards Self-Supervised High Energy Physics Foundation Models

Abstract We propose masked particle modeling (MPM) as a self-supervised method for learning generic, transferable, and reusable representations on unordered sets of inputs for use in high energy physics (HEP) scientific data. This work provides a novel scheme to perform masked modeling based pre-training to learn permutation invariant functions on sets. More generally, this work provides a step towards building large foundation models for HEP that can be generically pre-trained with self-supervised learning and later fine-tuned for a variety of down-stream tasks. In MPM, particles in a set are masked and the training objective is to recover their identity, as defined by a discretized token representation of a pre-trained vector quantized variational autoencoder. We study the efficacy of the method in samples of high energy jets at collider physics experiments, including studies on the impact of discretization, permutation invariance, and ordering. We also study the fine-tuning capability of the model, showing that it can be adapted to tasks such as supervised and weakly supervised jet classification, and that the model can transfer efficiently with small fine-tuning data sets to new classes and new data domains.

Heinrich, Lukas (ORCID:0000000240487584)↗

Microscopic Scattering Approach to In-Gap States

We develop a microscopic scattering formalism to describe Yu-Shiba-Rusinov (YSR) states due to a single Cr adatom on the Bi-terminated surface of beta Bi2Pd, by combining ab initio Wannier functions with a real-space Green's function approach in the Bogoliubov-de Gennes formalism[1]. Our framework reproduces key scanning tunneling spectroscopy features, including a single particle-hole asymmetric YSR peak and isotropic dIdV maps around the impurity. Decomposing the YSR states reveals contributions from four nearly degenerate C4v representations, with energy broadening masking their individual signatures. Spin-orbit coupling induces partial spin polarization, while the spatial asymmetry between particle and hole components arises from Cr d-Bi p hybridization. These results highlight the importance of realistic band structures and microscopic modeling for interpreting STM data for magnetic in-gap states on superconductors. Further advances examining layered 2D material surfaces, such as NbSe2, will be described[2]. For this system the superconducting properties are obtained from a full anisotropic Eliashberg calculation of the superconducting order parameter along with the charge density wave gap. Additional features associated with proposals to measure the dynamics of these individual YSR states will be presented. [1] arXiv:2507.08740 [2] arXiv:2507.11856

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Deep learning-based segmentation of lithium-ion battery microstructures enhanced by artificially generated electrodes

Accurate 3D representations of lithium-ion battery electrodes, in which the active particles, binder and pore phases are distinguished and labeled, can assist in understanding and ultimately improving battery performance. Here, we demonstrate a methodology for using deep-learning tools to achieve reliable segmentations of volumetric images of electrodes on which standard segmentation approaches fail due to insufficient contrast. We implement the 3D U-Net architecture for segmentation, and, to overcome the limitations of training data obtained experimentally through imaging, we show how synthetic learning data, consisting of realistic artificial electrode structures and their tomographic reconstructions, can be generated and used to enhance network performance. We apply our method to segment x-ray tomographic microscopy images of graphite-silicon composite electrodes and show it is accurate across standard metrics. We then apply it to obtain a statistically meaningful analysis of the microstructural evolution of the carbon-black and binder domain during battery operation.

25 ENERGY STORAGE↗

Exact matrix product state representation and convergence of a fully correlated electronic wavefunction in the infinite-basis limit

Here In this paper we present the exact representation of a fully correlated electronic wavefunction as the single-particle basis approaches completeness. It consists of a half-infinite chain of matrices of exponentially increasing size. The complete basis limit is illustrated numerically using the density-matrix renormalization-group method by computing the core-valence entanglement in the C 2 ground state in increasing subsets of cc-pVTZ and pVQZ bases until convergence is reached.

36 MATERIALS SCIENCE↗

Advanced Closures for Nonlinear Fluid Simulations of Plasmas (Final Report)

Advanced Closures for Nonlinear Fluid Simulations of Plasmas E. Held, Utah State University (Principal Investigator) J.-Y. Ji, Utah State University (Co-Investigator) Utah State University (USU) continued to develop and apply closures for hybrid fluid/kinetic simulations of magnetized fusion plasmas. The closures developed by USU's fusion theory and computation group use the drift ordering and are referred to as "continuum" and "general moment" approaches. They differ in their representation for the distribution functions which probabilistically describe the density of particles in 5 dimensional phase space. The two chosen velocity variables, needed in studies of kinetic physics for magnetized plasmas, are the speed and the component of a particle's velocity along the magnetic field. In magnetized plasmas, the charged particle trajectories are tied closely to the magnetic field lines. Both approaches account for binary interactions between particles using the Coulomb collision operator. Theoretical and computational development along both lines provided avenues for verification between the two approaches. This funding allowed USU to further develop physically accurate and numerically efficient models that advance our understanding of kinetic physics in experiments like the International Thermonuclear Experimental Reactor (ITER). The ITER device will likely prove the physical practicality of thermonuclear fusion by producing 500 megawatts of clean, abundant fusion power. This work helped to realize this goal by advancing our simulation capability using hybrid fluid/kinetic models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Difference in Fractional Occurrences of Precipitation Categories in Terms of Cloud Properties

The relationship between the precipitation stage and water cloud properties was studied in the mid-latitudes using space-borne radar and radiometer data, with a focus on a comparison between continental (over China) and oceanic (over the northwest Pacific) clouds. More specifically, fractional occurrences of the precipitation categories were investigated in terms of two-dimensional representations of cloud parameters, such as combinations of optical depth - effective particle radius, and liquid water path - droplet number density. As a result, systematic changes in the transition regarding pairs of the cloud parameters were well captured for all precipitation categories. The transition pattern was similar for both land and oceanic clouds. Notably, the ranges of all values of cloud parameters (except cloud droplet number density) were considerably larger for oceanic clouds for the precipitating category.

drizzle↗

Interaction networks for the identification of boosted H → $b\overline b$ decays

We develop an algorithm based on an interaction network to identify high-transverse-momentum Higgs bosons decaying to bottom quark-antiquark pairs and distinguish them from ordinary jets that reflect the configurations of quarks and gluons at short distances. The algorithm's inputs are features of the reconstructed charged particles in a jet and the secondary vertices associated with them. Describing the jet shower as a combination of particle-to-particle and particle-to-vertex interactions, the model is trained to learn a jet representation on which the classification problem is optimized. The algorithm is trained on simulated samples of realistic LHC collisions, released by the CMS Collaboration on the CERN Open Data Portal. The interaction network achieves a drastic improvement in the identification performance with respect to state-of-the-art algorithms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evaluation of two-particle properties within finite-temperature self-consistent one-particle Green’s function methods: Theory and application to GW and GF2

One-particle Green’s function methods can model molecular and solid spectra at zero or non-zero temperatures. One-particle Green’s functions directly provide electronic energies and one-particle properties, such as dipole moment. However, the evaluation of two-particle properties, such as $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$, can be challenging because they require a solution of the computationally expensive Bethe–Salpeter equation to find two-particle Green’s functions. We demonstrate that the solution of the Bethe–Salpeter equation can be completely avoided. Applying the thermodynamic Hellmann–Feynman theorem to self-consistent one-particle Green’s function methods, we derive expressions for two-particle density matrices in a general case and provide explicit expressions for GF2 and GW methods. Such density matrices can be decomposed into an antisymmetrized product of correlated one-electron density matrices and the two-particle electronic cumulant of the density matrix. Cumulant expressions reveal a deviation from ensemble representability for GW, explaining its known deficiencies. We analyze the temperature dependence of $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$ for a set of small closed-shell systems. Interestingly, both GF2 and GW show a non-zero spin contamination and a non-zero fluctuation of the number of particles for closed-shell systems at the zero-temperature limit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Convergence study of wakefield simulations with GdfidL and ECHO3D

The interaction of charged particle beams with vacuum chamber components gives rise to electromagnetic wakefields, whose frequency-domain representation is known as beam coupling impedance. Geometric impedance arising from discontinuities and transitions in the vacuum chamber is the focus of this study. Minimizing this impedance is essential to mitigate adverse collective effects in modern storage rings operating with high-intensity particle beams. Accurate and reliable impedance simulations is a key factor of the vacuum chamber design. This paper presents the results of a convergence study of two widely used electromagnetic solvers, GdfidL and ECHO3D, applied to key vacuum-chamber components of the National Synchrotron Light Source II (NSLS-II) storage ring. Detailed comparisons are performed for several geometries, including flange absorbers, RF bellows, button-type beam position monitors, and an in-vacuum undulator (IVU). The results show notable differences in convergence and computational efficiency between the two codes. While GdfidL provides highly resolved results and serves as a common benchmark tool, ECHO3D yields consistent results with coarser meshes, significantly reducing simulation time and memory demands. Simulations with a full-geometry IVU model demonstrate that simplified taper-transition models can miss important impedance contributions. In conclusion, these findings provide practical guidelines for efficient and accurate impedance modeling to optimize design of vacuum chamber components for accelerators.

36 MATERIALS SCIENCE↗