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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A new technique to retrieve aerosol vertical profiles using micropulse lidar and ground-based aerosol measurements

Accurately characterizing the vertical distribution of aerosols and their cloud-forming properties is crucial for understanding aerosol-cloud interactions and their impact on climate. This study presents a novel technique for retrieving vertical profiles of aerosols, cloud condensation nuclei (CCN), and ice nucleating particles (INP) by combining micropulse lidar, radiosonde, and ground-based aerosol measurements. Herein, the technique is applied to data collected by our team at Texas A&M University during the Tracking Aerosol Convection Interactions ExpeRiment (TRACER) campaign. Ground-based aerosol size distribution and CCN counter data are used to estimate the value of the aerosol hygroscopicity parameter, κ. The derived κ, together with Mie scattering theory and the relative humidity profile from the radiosonde, is used to estimate aerosol size growth and the associated increase in backscatter at each altitude. We then correct the lidar backscatter to dry conditions to produce the dry aerosol backscatter coefficient profile. The dry aerosol backscatter coefficient profile is linearly scaled to collocated surface measurements of aerosols, CCN, and INP to produce corresponding vertical profiles. Combining lidar backscatter profiles with aerosol and cloud nucleation measurements leads to a more realistic representation of vertical distributions of aerosol properties. The method could be readily applied to lidar measurements in future field campaigns.

Chen, Bo [Texas A & M Univ., College Station, TX (↗
DOE OSTI2023

A guide to diagnosing colored resonances at hadron colliders

We present a comprehensive study on how to distinguish the properties of heavy dijet resonances at hadron colliders. A variety of spins, chiral couplings, charges, and QCD color representations are considered. Distinguishing the different color representations is particularly difficult at hadron colliders. To determine the QCD color structure, we consider a third jet radiated in a resonant dijet event. We show that the relative rates of three-jet versus two-jet processes are sensitive to the color representation of the resonance. We also show analytically that the antennae radiation pattern of soft radiation depends on the color structure of dijet events and develops an observable that is sensitive to the antennae patterns. Finally, we exploit a Convolutional Neural Network with Machine Learning techniques to differentiate the radiation patterns from different colored resonances and find encouraging results to discriminate them. We demonstrate our results numerically at a 14 TeV LHC, and the methodology presented here should be applicable to other future hadron colliders.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Optimization performance, fidelity, and cost: SIAM VQE

This dataset contains files storing results from classically-simulated quantum subroutines within a dynamical mean-field theory workflow, and jupyter notebooks processing the data in these files to generate plots. The files store: (1) Results from variational quantum eigensolver (VQE) simulations searching for optimal parameters allowing parametrized quantum circuits to prepare approximations to ground states of different Anderson impurity models (AIMs) (2) Results from simulations of a quantum Lanczos algorithm (QLA) estimating the Lanczos coefficients defining the continued-fraction representation of an (AIM) Green’s function Description: Any file named vqe_gs_results* stores approximations to the ground state and energy of a given AIM estimated using three different methods: (1) Numerical diagonalization (2) Ideal VQE simulation (3) VQE simulation with sampling noise For each VQE simulations metadata about the optimization (optimization results plus number of quantum circuits that would have been executed on real hardware) is also stored. Any file named qla_dos_results* estimations for the Lanczos coefficients defining the Green’s function of an AIM. The stored estimations are achieved using different methods: (1) Numerical Lanczos algorithm from initial states obtained from numerical diagonalization (2) Simulated quantum Lanczos algorithm from initial states prepared from parametrized quantum circuits yielded by corresponding ideal and noisy VQE subroutines. The dataset is used and described in M. Karabin et al., "Quantum solver for single-impurity Anderson models with particle-hole symmetry", Phys. Rev. Research 8, 033066 (2026). DOI: https://doi.org/10.1103/7ys3-tl4l

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2011

Micro-Scale Distributions of Major and Trace Elements in Chondrites

The Hayabusa spacecraft has successfully returned to Earth after two touchdowns on the surface of Asteroid 25143 Itokawa. This asteroid is classified as an S-type and inferred to consist of materials similar to ordinary chondrites or primitive achondrites [1]. More than 1500 particles have been identified consisting of olivine, pyroxene, plagioclase, Fe sulfide and Fe metal, with compositions consistent with being of LL origin. While the chondritic components are familiar to us, the level of detail to which the Itokawa samples will be exposed to will be unprecedented given that the samples are reasonably large and accessible to a wide variety of techniques. In many ways, we expect that our knowledge base of the comparator chondrites will be found to be wanting. Chondrites are the building blocks of the solar system. However, these rocks are essentially breccias and they are quite variable in bulk element compositions as well as compositions of the individual components. We have initiated a program of analysis for chondrites focusing on major and trace element distributions between the mineral components and the matrix. The issues to be addressed include the homogeneity of matrix and chondrule components and the representivity of any given sample to the bulk meteorite. This may be particularly important given the limited numbers of Itokawa grains that may be available for a specific analysis. As an initial study, we have taken thin sections of carbonaceous chondrites to study the representivity of the matrix compositions. Spot locations were constrained to limited regions of the sections so as to assess the variability of a local scale. Further work will be required to assess variability over a centimeter scale.

Ireland, T. R.↗
DOE OSTI2024

Modeling of Macroscopic Dynamics in Three-Dimensional Plasma Configurations (Final Technical Report)

A numerical simulation-based study to investigate the stability properties of inherently three-dimensional plasma magnetic configurations is presented. These configurations, called stellarators and heliotrons, use external helical coils to provide the twisting magnetic field that is required to confine ionized plasma particles. In contrast to the tokamak configuration, which nominally has the property of axisymmetry, electrical current does not need to be induced in stellarator and heliotron plasma. Experiments have shown that these three-dimensional configurations are remarkably robust to macroscopic instability, even when driven to relatively high levels of plasma pressure. Our understanding of how the plasma is able to withstand macroscopic instability will be enhanced by the work proposed here. The effort will generalize existing numerical modeling capability, and it will apply the new capability to model macroscopic dynamics in these three-dimensional configurations. The development work has been applied to the NIMROD simulation code (https://nimrodteam.org) that was originally based on axisymmetric problem domains. It has been generalized to accommodate the three-dimensional shaping of the plasma volume in stellarators and heliotrons. Its use of equilibrium data has also been generalized, and the implementation allows different numerical representation of the magnetic field. Analysis has been performed to inform the selection of numerical methods, and a method to accelerate the iterative solution of the necessary algebraic systems has been developed and implemented. Verification of the new NIMSTELL code on benchmark problems has also been an important part of the effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Application of Three-Component PIV to a Hovering Rotor Wake

The key to accurate predictions of rotorcraft aerodynamics, acoustics, and dynamics lies in the accurate representation of the rotor wake. The vortical wake computed by rotorcraft CFD analyses typically suffer from numerical dissipation before the first blade passage. With some a priori knowledge of the wake trajectory, grid points can be concentrated along the trajectory to minimize the dissipation. Comprehensive rotorcraft analyses based on lifting-line theory rely on classical vortex models and/or semi-empirical information about the tip vortex structure. Until the location, size, and strength of the trailed tip vortex can be measured over a range of wake ages, the analyses will continue to be adjusted on a trial and error basis in order to correctly predict blade airloads, acoustics, dynamics, and performance. Using the laser light sheet technique, tip vortex location can be acquired in a straightforward manner. Measuring wake velocities and vortex core size, however, has been difficult and tedious using point-measurement techniques such as laser velocimetry. Recently, the Particle Image Velocimetry (PIV) technique has proven to be an efficient method for acquiring velocity measurements over relatively large areas and volumes of a rotor wake. The work reported to date, however, has been restricted to 2-component velocity measurements of the rotor wake. Three-component velocity measurements of a hovering rotor wake were acquired at NASA Ames Research Center in May 1999. This experiment represents a major step toward understanding the detailed structure of a rotor wake. This paper will focus primarily on the experimental technique used in acquiring this data. The accuracy and limitations of the current technique will also be discussed. Representative velocity field measurements will be included.

Yamauchi, Gloria K.↗
DOE OSTI2020

High energy modular bootstrap, global symmetries and defects

We derive Cardy-like formulas for the growth of operators in different sectors of unitary 2 dimensional CFT in the presence of topological defect lines by putting an upper and lower bound on the number of states with scaling dimension in the interval [Δ – δ, Δ + δ] for large Δ at fixed δ. Consequently we prove that given any unitary modular invariant 2D CFT symmetric under finite global symmetry G (acting faithfully), all the irreducible representations of G appear in the spectra of the untwisted sector; the growth of states is Cardy like and proportional to the “square” of the dimension of the irrep. In the Schwarzian limit, the result matches onto that of JT gravity with a bulk gauge theory. If the symmetry is non-anomalous, the result applies to any sector twisted by a group element. For c > 1, the statements are true for Virasoro primaries. Furthermore, the results are applicable to large c CFTs. We also extend our results for the continuous U(1) group.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗
DOE OSTI2023

Direct prediction of inelastic neutron scattering spectra from the crystal structure*

Abstract Inelastic neutron scattering (INS) is a powerful technique to study vibrational dynamics of materials with several unique advantages. However, analysis and interpretation of INS spectra often require advanced modeling that needs specialized computing resources and relevant expertise. This difficulty is compounded by the limited experimental resources available to perform INS measurements. In this work, we develop a machine-learning based predictive framework which is capable of directly predicting both one-dimensional INS spectra and two-dimensional INS spectra with additional momentum resolution. By integrating symmetry-aware neural networks with autoencoders, and using a large scale synthetic INS database, high-dimensional spectral data are compressed into a latent-space representation, and a high-quality spectra prediction is achieved by using only atomic coordinates as input. Our work offers an efficient approach to predict complex multi-dimensional neutron spectra directly from simple input; it allows for improved efficiency in using the limited INS measurement resources, and sheds light on building structure-property relationships in a variety of on-the-fly experimental data analysis scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Improving GCM Representation of Convective Cloud Microphysics by Using ARM Raman Lidar and Cloud Radar Observations

The overall objective of this Department of Energy (DOE) Atmospheric System Research (ASR) funded project is to improve the representation of convective cloud microphysics in global climate models (GCMs) and check it by comparing the model simulations with observations. We have derived the cloud ice water content by synthesizing ARM Raman lidar (RL) and cloud radar observations at the Atmospheric Radiation Measurement (ARM) sites. Noting that the simulated anvil clouds in terms of their macro- and micro-physical properties are sensitive to the parameterization of convective microphysical processes, observed ice water content in anvil clouds provides a useful constraint on these parameterizations. We have improved the convective microphysics parameterization scheme by (1) considering sedimentation for cloud ice crystals that do not fall in the original scheme, (2) applying a new terminal velocity parameterization that depends on the environmental conditions for convective snow, (3) adding a new hydrometeor category, “rimed ice,” to the original four-class (cloud liquid, cloud ice, rain, and snow) scheme, and (4) allowing convective clouds to detrain snow particles into stratiform clouds. We have examined the impact of improved convective cloud microphysics parameterization on the simulated global climate from GCM simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Navigating collinear superspace

We introduce a new set of effective field theory rules for constructing Lagrangians with \( \mathcal{N} \) = 1 supersymmetry in collinear superspace. In the standard superspace treatment, superfields are functions of the coordinates \( \left({x}^{\mu },{\theta}^{\alpha },{\theta}^{\dagger \overset{\cdot }{\alpha }}\right) \) , and supersymmetry preservation is manifest at the Lagrangian level in part due to the inclusion of auxiliary F - and D -term components. By contrast, collinear superspace depends on a smaller set of coordinates ( x μ , η , η † ), where η is a complex Grassmann number without a spinor index. This provides a formulation of supersymmetric theories that depends exclusively on propagating degrees of freedom, at the expense of obscuring Lorentz invariance and introducing inverse momentum scales. After establishing the general framework, we construct collinear superspace Lagrangians for free chiral matter and non-Abelian gauge fields. For the latter construction, an important ingredient is a superfield representation that is simultaneously chiral, anti-chiral, and real; this novel object encodes residual gauge transformations on the light cone. Additionally, we discuss a fundamental obstruction to constructing inter- acting theories with chiral matter; overcoming these issues is the subject of our companion paper, where we introduce a larger set of superfields to realize the full range of interactions compatible with \( \mathcal{N} \) = 1. Along the way, we provide a novel framing of reparametrization invariance using a spinor decomposition, which provides insight into this important light-cone symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A new large-scale suspended sediment model and its application over the United States

Abstract. Suspended sediment plays a vital role in the regional and global cycling of carbon and nutrients by carrying carbon and nutrients from headwaters into lowland rivers and the oceans. Sediment transport through river systems is often fundamentally modified by human activities such as reservoir management. However, a physically based representation of sediment transport is still missing in most existing Earth system models (ESMs), which are essential tools for modeling and predicting Earth system changes. Here, we introduce a multi-process river sediment module for ESMs, which includes (1) hillslope soil erosion and sediment discharge into streams, (2) sediment transport processes through river networks, (3) reservoir operation based on the inflows from upstream areas and water demand from downstream areas, and (4) sediment trapping by reservoirs. All model parameters are estimated a priori without calibration. We apply this new sediment modeling framework to the contiguous United States and validate it against historical observations of monthly streamflow and sediment discharges at 35 river gauges. The model reasonably well captures the long-term balance and seasonal variations of suspended sediment in large river systems. Furthermore, our model results show that suspended sediment discharge in managed rivers is affected more by reservoirs' direct trapping of sediment particles than by their flow regulation. This new sediment module enables future modeling of the transportation and transformation of carbon and nutrients carried by the fine sediment along the river–ocean continuum to close the global carbon and nutrient cycles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of 208 Pb with Minimal Modeling Assumptions

The symmetry energy and its density dependence are crucial inputs for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of 0.283 ± 0.071 fm for the neutron-skin thickness of 208 Pb, implying a slope parameter L = 106 ± 37 MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy S 0 , L, and $R^{^{208}Pb}_{skin}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER clearly favor smaller values of the neutron skin and L, as well as negative symmetry incompressibilities. Finally, combining astrophysical data with PREX-II and chiral effective field theory constraints yields S 0 = $33.0^{+2.0}_{-1.8}$ MeV , L = $53^{+14}_{-15}$ MeV, and $R^{^{208}Pb}_{skin} = 0.17^{+0.04}_{-0.04}$ fm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Physics-constrained machine learning for electrodynamics without gauge ambiguity based on Fourier transformed Maxwell’s equations

We utilize a Fourier transformation-based representation of Maxwell’s equations to develop physics-constrained neural networks for electrodynamics without gauge ambiguity, which we label the Fourier–Helmholtz–Maxwell neural operator method. In this approach, both of Gauss’s laws and Faraday’s law are built in as hard constraints, as well as the longitudinal component of Ampère–Maxwell in Fourier space, assuming the continuity equation. An encoder–decoder network acts as a solution operator for the transverse components of the Fourier transformed vector potential, $\hat{A}_⟂(k,t)$, whose two degrees of freedom are used to predict the electromagnetic fields. This method was tested on two electron beam simulations. Among the models investigated, it was found that a U-Net architecture exhibited the best performance as it trained quicker, was more accurate and generalized better than the other architectures examined. We demonstrate that our approach is useful for solving Maxwell’s equations for the electromagnetic fields generated by intense relativistic charged particle beams and that it generalizes well to unseen test data, while being orders of magnitude quicker than conventional simulations. We show that the model can be re-trained to make highly accurate predictions in as few as 20 epochs on a previously unseen data set.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2019

Micro-scale Artificial Weave Generation Capabilities for Thermal Protection System Material Modeling

Thermal Protection System (TPS) modeling requires accurate representation and prediction of the thermomechanical behavior of ablative materials. State-of-the-art TPS materials such as Phenolic Impregnated Carbon Ablator (PICA) have a proven flight record and demonstrate exceptional capabilities for handling extreme aerothermal heating conditions. The constant push for lightweight materials that are flexible in their design and performance, and hence allow for a wide range of mission profiles, has led NASA over the past years to develop its Heatshield for Extreme Entry Environment Technology (HEEET). HEEET is based primarily on a dual layer woven carbon fiber architecture and the technology has successfully been tested in arc-jet facilities. These recent developments have sparked interest in the accurate micro-scale modeling of composite weave architectures, to predict the structural response of macro-scale heatshields upon atmospheric entry. This effort can be extended to incorporate in-depth failure mechanics analyses as a result of local thermal gradients or high-velocity particle impact.

Visser, Sander J.↗
DOE OSTI2021

Flipped SU(5) GUT phenomenology: proton decay and $\mathbf {g_\mu - 2}$

We consider proton decay and $\mathbf {g_\mu - 2}$ in flipped SU(5) GUT models. We first study scenarios in which the soft supersymmetry-breaking parameters are constrained to be universal at some high scale M in above the standard GUT scale where the QCD and electroweak SU(2) couplings unify. In this case the proton lifetime is typically > rsim10 36 years, too long to be detected in the foreseeable future, and the supersymmetric contribution to $\mathbf {g_\mu - 2}$ is too small to contribute significantly to resolving the discrepancy between the experimental measurement and data-driven calculations within the Standard Model. However, we identify a region of the constrained flipped SU(5) parameter space with large couplings between the 10- and 5-dimensional GUT Higgs representations where p → e + π 0 decay may be detectable in the Hyper-Kamiokande experiment now under construction, though the contribution to $\mathbf {g_\mu - 2}$ is still small. A substantial contribution to $\mathbf {g_\mu - 2}$ is possible, however, if the universality constraints on the soft supersymmetry-breaking masses are relaxed. We find a ‘quadrifecta’ region where observable proton decay co-exists with a (partial) supersymmetric resolution of the $\mathbf {g_\mu - 2}$ discrepancy and acceptable values of m h and the relic LSP density.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Line operators in Chern-Simons-Matter theories and Bosonization in Three Dimensions II: Perturbative analysis and all-loop resummation

We study mesonic line operators in Chern-Simons theories with bosonic or fermionic matter in the fundamental representation. In this paper, we elaborate on the classification and properties of these operators using all loop resummation of large N perturbation theory. We show that these theories possess two conformal line operators in the fundamental representation. One is a stable renormalization group fixed point, while the other is unstable. They satisfy first-order chiral evolution equations, in which a smooth variation of the path is given by a factorized product of two mesonic line operators. The boundary operators on which the lines can end are classified by their conformal dimension and transverse spin, which we compute explicitly at finite ’t Hooft coupling. We match the operators in the bosonic and fermionic theories. Finally, we extend our findings to the mass deformed theories and discover that the duality still holds true.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1993

Energy storage and dissipation in the magnetotail during substorms. I - Particle simulations. II - MHD simulations

2D electromagnetic particle simulations are used to investigate the dynamics of the tail during development of substorms under the influence of the pressure in the magnetospheric boundary layer and the dawn-to-dusk electric field. It is shown that pressure pulses result in thinning of the tail current sheet as the magnetic field becomes pinched near the region where the pressure pulse is applied. The pinching leads to the tailward flow of the current sheet plasma and the eventual formation and injection of a plasmoid. Surges in the dawn-to-dusk electric field cause plasma on the flanks to convect into the center of the current sheet, thereby thinning the current sheet. The pressure in the magnetospheric boundary laser is coupled to the dawn-to-dusk electric field through the conductivity of the tail. Changes in the predicted evolution of the magnetosphere during substorms due to changes in the resistivity are investigated under the assumption that MHD theory provides a suitable representation of the global or large-scale evolution of the magnetotail to changes in the solar wind and to reconnection at the dayside magnetopause. It is shown that the overall evolution of the magnetosphere is about the same for three different resistivity distributions with plasmoid formation and ejection in each case.

Winglee, R. M.↗