Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

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 OSTI2023

Unsupervised Learning Based Interaction Force Model for Nonspherical Particles in Incompressible Flows

This project provides a neural network-based interaction force model for gas-solid flows from low to intermediate Reynolds numbers and concentration, which can be linked to MFiX-DEM. We have constructed a database of the interaction force between the irregular-shaped particles using a spherical harmonic method and the fluid phase based on the particle-resolved direct numerical simulation (PR-DNS) with immersed boundary-based gas kinetic scheme. Unsupervised learning method, i.e., variational auto-encoder (VAE) has been applied to extract the primitive shape factors determining the drag force, lifting forces, and torque. The interaction force model has been trained and validated with a simple but effective multi-layer feed-forward neural network: multi-layer perceptron (MLP), which will be concatenated after the encoder of the previously trained VAE for geometry feature extraction for single, irregular particles. We have trained transpose convolutional neural networks with the PR-DNS data to predict the velocity and pressure gradient of the single particle systems and utilized them to calculate drag force of multi-particle systems. This model can provide high computational efficiency because it does not require collecting multiparticle system data from PR-DNS.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Seismic anisotropy, dominant slip systems and phase transitions in the lowermost mantle

SUMMARY The presence of seismic anisotropy at the base of the Earth's mantle is well established, but there is no consensus on the deformation mechanisms in lower mantle minerals that could explain it. Strong anisotropy in magnesium post-perovskite (pPv) has been invoked, but different studies disagree on the dominant slip systems at play. Here, we aim to further constrain this by implementing the most recent results from atomistic models and high-pressure deformation experiments, coupled with a realistic composition and a 3-D geodynamic model, to compare the resulting deformation-induced anisotropy with seismic observations of the lowermost mantle. We account for forward and reverse phase transitions from bridgmanite (Pv) to pPv. We find that pPv with either dominant (001) or (010) slip can both explain the seismically observed anisotropy in colder regions where downwellings turn to horizontal flow, but only a model with dominant (001) slip matches seismic observations at the root of hotter large-scale upwellings. Allowing for partial melt does not change these conclusions, while it significantly increases the strength of anisotropy and reduces shear and compressional velocities at the base of upwellings.

58 GEOSCIENCES↗
DOE OSTI2025

PREFACE: A Search for Long-Lived Particles at the Large Hadron Collider

The Standard Model (SM) fails to explain many problems (neutrino masses, dark matter, and matter–antimatter asymmetry, among others) that may be resolved with new particles beyond the SM. No observation of such new particles may be explained either by their exceptionally high mass or by considerably small coupling to SM particles. The latter case implies relatively long lifetimes. Such long-lived particles (LLPs) then to have signatures different from those of SM particles. Searches in the “central region” are covered by the LHC general purpose experiments. The forward small angle region far from the interaction point (IP) is unexplored. Such particles are expected to have the energy as large as E = O(1 TeV) and Lorentz time dilation factor γ = E/m ≈ 10 2 –10 3 (with m the particle mass) hence long enough decay distances. A new class of specialized LHC detectors dedicated to LLP searches has been proposed for the forward regions. Among these experiments, FASER is already operational, and FACET is under consideration at a location 100 m from the LHC IP5 (the CMS detector intersection). However, some features of FACET require a specially enlarged beam pipe, which cannot be implemented for LHC Run 4. In this study, we explore a simplified version of the proposed detector PREFACE compatible with the standard LHC beam pipe in the HL-LHC Run 4. Realistic G EANT 4 simulations are performed and the background is evaluated. An initial analysis of the physics potential with the PREFACE geometry indicates that several significant channels could be accessible with sensitivities comparable to FACET and other LLP searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

INN Airfoil Tools

This software package contains code to perform airfoil inverse shape design using a sparsified invertible neural network model. This is a specialized deep learning architecture that enables both network weights to be shared to evaluate forward or inverse passes through the network. As a design tool, this network allows practitioners to specify desired aerodynamic and structural properties and generate shapes that satisfy those properties. The invertible neural network is accompanied with tuned parameters as well as software tools to train new models as needed.

Glaws, Andrew↗
DOE OSTI2020

Computational modeling and neutron imaging to understand interface shape and solute segregation during the vertical gradient freeze growth of BaBrCl:Eu

In this work we apply continuum models to analyze phase change, heat transfer, fluid flow, solute transport, and segregation in order to understand prior neutron imaging observations of the vertical gradient freeze growth of Eu-doped BaBrCl. The models provide a rigorous framework in which to understand the mechanisms that are responsible for the complicated evolution of interface shape and dopant distribution in the growth experiment. We explain how a transition in the solid/liquid interface shape from concave to convex is driven by changes in radial heat transfer caused by furnace design. We also provide a mechanistic explanation of how dynamic growth conditions and changes of the flow structure in the melt result in complicated segregation patterns in this system. A growth pause caused by controller lock-up is shown to result in a band of solute depletion in accordance with classical theory. However, changing flow patterns during growth result in a non-monotonic axial distribution of solute that cannot be explained by simple application of classical segregation models. We assert that the approach presented here, namely the use of rigorous models in conjunction advanced diagnostics, such as neutron imaging, provides an exciting path forward for process optimization and control, accelerating the incremental advances that have, in the past, typically relied on empiricism, experience, and intuition.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Increasing the Reproducibility and Replicability of Supervised AI/ML in the Earth Systems Science by Leveraging Social Science Methods

Artificial intelligence (AI) and machine learning (ML) pose a challenge for achieving science that is both reproducible and replicable. The challenge is compounded in supervised models that depend on manually labeled training data, as they introduce additional decision-making and processes that require thorough documentation and reporting. We address these limitations by providing an approach to hand labeling training data for supervised ML that integrates quantitative content analysis (QCA)—a method from social science research. The QCA approach provides a rigorous and well-documented hand labeling procedure to improve the replicability and reproducibility of supervised ML applications in Earth systems science (ESS), as well as the ability to evaluate them. Specifically, the approach requires (a) the articulation and documentation of the exact decision-making process used for assigning hand labels in a “codebook” and (b) an empirical evaluation of the reliability” of the hand labelers. In this paper, we outline the contributions of QCA to the field, along with an overview of the general approach. We then provide a case study to further demonstrate how this framework has and can be applied when developing supervised ML models for applications in ESS. With this approach, we provide an actionable path forward for addressing ethical considerations and goals outlined by recent AGU work on ML ethics in ESS.

58 GEOSCIENCES↗
DOE OSTI2023

Phase coexistence implications of violating Newton’s third law

Newton’s third law, action = reaction, is a foundational statement of classical mechanics. However, in natural and living systems, this law appears to be routinely violated for constituents interacting in a nonequilibrium environment. Here, we use computer simulations to explore the macroscopic phase behavior implications of breaking microscopic interaction reciprocity for a simple model system. We consider a binary mixture of attractive particles and introduce a parameter that is a continuous measure of the degree to which interaction reciprocity is broken. In the reciprocal limit, the species are indistinguishable, and the system phase separates into domains with distinct densities and identical compositions. Furthermore, increasing nonreciprocity is found to drive the system to explore a rich assortment of phases, including phases with strong composition asymmetries and three-phase coexistence. Many of the states induced by these forces, including traveling crystals and liquids, have no equilibrium analogs. By mapping the complete phase diagram for this model system and characterizing these unique phases, our findings offer a concrete path forward toward understanding how nonreciprocity shapes the structures found in living systems and how this might be leveraged in the design of synthetic materials.

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

Non-Gaussianity constraints with anisotropic μ distortion measurements from Planck

ABSTRACT Primordial non-Gaussianity can source μ-distortion anisotropies that are correlated with the large-scale temperature and polarization signals of the cosmic microwave background (CMB). A measurement of μT and μE correlations can therefore be used to constrain it on wavelengths of perturbations not directly probed by the standard CMB anisotropies. We carry out a first rigorous search for μ-distortion anisotropies with Planck data, applying the well-tested constrained ILC component-separation method combined with the needlet framework. We correlate the reconstructed μ map with the CMB anisotropies to derive constraints on the amplitude fNL of the local form bispectrum, specifically on the squeezed configurations with effective wavenumbers $k_s \simeq {740}\, \mathrm{Mpc^{-1}}$ and $k_L \simeq {0.05}\, \mathrm{Mpc^{-1}}$, improving previously estimated constraints by more than an order of magnitude. This enhancement is owing to the fact that we are able to use the full multipole information by carefully controlling biases and systematic effects in the analysis. We also for the first time incorporate constraints from measurements of μE correlations, which further tighten the limits. A combination of the derived Planck μT and μE power spectra yields |fNL| ≲ 6800 (95 per cent c.l.) on this highly squeezed bispectrum. This is only ≃3 times weaker than the anticipated constraint from Litebird. Furthermore we show that a combination of Litebird with Planck can improve the expected future constraint by ${\simeq}20{{\%}}$. These limits can be used to constrain multifield inflation models and primordial black hole formation scenarios, thus providing a promising novel avenue forward in CMB cosmology.

Rotti, Aditya (ORCID:0000000341245226)↗
DOE OSTI2022

Development and Commercialization of an IDAES-Based Power Plant Performance Monitoring and Optimization System

DOE and NETL have created an advanced, open-source computational platform through the Institute for the Design of Advanced Energy Systems (IDAES). The IDAES platform is a very extensive modeling environment that can be used for a broad range of power plant and process applications. MapEx Software is developing and commercializing a software application that makes it easier to set-up and run IDAES-based analyses. The MapEx-developed software application will replace the need for custom Python language coding with a user-friendly graphical user interface where the user can construct a flowsheet diagram of the IDAES model by inserting icons representing the plant equipment onto the screen. This application will make the implementation of modeling and optimization of existing fossil-fired power plants more straight-forward and less time-consuming. The effort focuses on performance monitoring and optimization of plant operations for the existing coal-fired power plant fleet but is built on a structure that allows expansion into the broad range of applications where IDAES methods may be applied.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Photoproduction of K0K0-bar with the GlueX Experiment

In this dissertation, we study photoproduction of K0 ?K 0 through the KSKLp and KSKSp ?nal states with the GlueX Phase-I (GlueX-I) data set. We measure the spin-density matrix elements (SDMEs) and di?erential cross section of ?(1020) ! KSKL to better understand photoproduction of light vector mesons. A mass independent Partial Wave Analysis (PWA) of the KSKL and KSKS meson spectrum below 2 GeV is also undertaken to study the meson spectrum with JPC = even++ and odd??. The ?(1020) di?erential cross section is measured at E = 8:2 ? 8:8 GeV and ?t = 0:15 ? 1:00 GeV2. The di?erential cross section is well described by an exponential decay with slope 4:44?0:01 GeV?2 and the integrated cross section is determined to be 295:7?0:4 nb, only statistical uncertainties are quoted. Both measurements are consistent and far more precise than the previous measurement by Ballam et al. [1]. The ?(1020) SDMEs were measured in nine bins of ?t in the same range. At low ?t, we ?nd the data were consistent with s-channel helicity conservation, SCHC, i. e. the only non-zero SDMEs were ?1 1?1 = ?Im(?2 1?1) = 1=2. At higher ?t, the SDMEs deviate from SCHC and the measurements indicate this is due to natural parity exchange since we observe that the SDMEs are consistent with zero unnatural exchange. We also ?find that the contribution from helicity double-flip amplitudes is consistent with zero. The measured SDMEs are in poor agreement with theoretical predictions put forward by JPAC [35]. These measurements will serve as input to re?ne models of production processes, which will be essential for the interpretation of possible signals of exotic mesons in GlueX. Analysis of the spectrum above the ?(1020) indicates the presence of at least three particles at ? 1:50, ? 1:75, and 2:20 GeV. We employ simple parametrizations of the resonance line shape based on relativistic Breit-Wigner functions considering cases with and without interference. The resonance at ? 1:50 GeV is not identifi?ed as any specifi?c resonance, but could be due to interference between the ?(1450) and !(1420). Our models favor the resonance at ~1:75 GeV as the X(1750) rather than the ?(1680). We ?find that introducing a third Breit-Winger into the models improves the ?t quality with the ?2/ndf going from 1.75 to 1.59 and 1.41 to 1.18 for models without and with interference respectively. A PWA of the KSKL spectrum indicates that the spectrum is consistent with exclusively spin-1 contribution up to ? 1:6 GeV. Above ? 1:6 GeV we ?find evidence for a small spin-3 contribution. We also ?nd that the spectrum predominantly positive reflectivity up to ? 1:6 GeV, after which the spectrum becomes a nearly equal mix of both. The KSKS system shows a rich spectrum with multiple resonance structures but is more statistically limited than the KSKL system. Modelling the line shape we ?find evidence at greater than 10? level in favor of a resonance at ? 1:75 GeV. The model parameters suggest this resonance is the f0(1710). The PWA of the KSKS system was inconclusive. However, using a small set of amplitudes suggests that the spin-2 contribution to the spectrum is primarily in the range 1:2?1:6 GeV, where f2(1270), a2(1320) and f02 (1525) are expected.

Linera, Gabriel Rodriguez↗
DOE OSTI2024

Measurement of the double-differential cross section of muon-neutrino charged-current interactions with low hadronic energy in the NOvA Near Detector

The NOvA collaboration reports cross-section measurements for $\nu_{\mu}$ charged-current interactions with low hadronic energy (maximum kinetic energy of 250 MeV for protons and 175 MeV for pions) in the NOvA Near Detector. The results are presented as a double-differential cross section as a function of the direct observables of the final-state muon kinematics. Results are also presented as a single-differential cross section as a function of the derived square of the four-momentum transfer, $Q^{2}$, and as a function of the derived neutrino energy. The data correspond to an accumulated 8.09$\times10^{20}$ protons-on-target (POT) in the neutrino mode of the NuMI beam, with a narrow band of neutrino energies peaked at 1.8 GeV. The analysis provides a sample of neutrino-nucleus interactions with an enhanced fraction of quasi-elastic and two-particle-two-hole (2p2h) interactions. This enhancement allows quantitative comparisons with various nuclear models. We find strong disagreement between data and theory-based models in various regions of the muon kinematic phase space, especially in the forward muon direction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Hadronic interaction model sibyll 2.3d and extensive air showers

We present a new version of the hadron interaction event generator SIBYLL. While the core ideas of the model have been preserved, the new version handles the production of baryon pairs and leading particles in a new way. In addition, production of charmed hadrons is included. Updates to the model are informed by high-precision measurements of the total and inelastic cross sections with the forward detectors at the LHC that constrain the extrapolation to ultrahigh energy. Minimum-bias measurements of particle spectra and multiplicities support the tuning of fragmentation parameters. This paper demonstrates the impact of these changes on air-shower observables such as X max and N μ , drawing comparisons with other contemporary cosmic-ray interaction models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Machine Learning Self-Diffusion Prediction for Lennard-Jones Fluids in Pores

Predicting the diffusion coefficient of fluids under nanoconfinement is important for many applications including the extraction of shale gas from kerogen and product turnover in porous catalysts. Due to the large number of important variables, including pore shape and size, fluid temperature and density, and the fluid–wall interaction strength, simulating diffusion coefficients using molecular dynamics (MD) in a systematic study could prove to be prohibitively expensive. Here, we use machine learning models trained on a subset of MD data to predict the self-diffusion coefficients of Lennard-Jones fluids in pores. Our MD data set contains 2280 simulations of ideal slit pore, cylindrical pore, and hexagonal pore geometries. We use the forward feature selection method to determine the most useful features (i.e., descriptors) for developing an artificial neutral network (ANN) model with an emphasis on easily acquired features. Our model shows good predictive ability with a coefficient of determination (i.e., R 2 ) of ~0.99 and a mean squared error of ~2.9 × 10 –5 . Finally, we propose an alteration to our feature set that will allow the ANN model to be applied to nonideal pore geometries.

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

Parameter estimation with maximal updated densities

A recently developed measure-theoretic framework solves a stochastic inverse problem (SIP) for models where uncertainties in model output data are predominantly due to aleatoric (i.e., irreducible) uncertainties in model inputs (i.e., parameters). The subsequent inferential target is a distribution on parameters. Another type of inverse problem is to quantify uncertainties in estimates of “true” parameter values under the assumption that such uncertainties should be reduced as more data are incorporated into the problem, i.e., the uncertainty is considered epistemic. A major contribution of this work is the formulation and solution of such a parameter identification problem (PIP) within the measure-theoretic framework developed for the SIP. The approach is novel in that it utilizes a solution to a stochastic forward problem (SFP) to update an initial density only in the parameter directions informed by the model output data. In other words, this method performs “selective regularization” only in the parameter directions not informed by data. Further, the solution is defined by a maximal updated density (MUD) point where the updated density defines the measure-theoretic solution to the PIP. Another significant contribution of this work is the full theory of existence and uniqueness of MUD points for linear maps with Gaussian distributions. Data-constructed Quantity of Interest (QoI) maps are also presented and analyzed for solving the PIP within this measure-theoretic framework as a means of reducing uncertainties in the MUD estimate. We conclude with a demonstration of the general applicability of the method on two problems involving either spatial or temporal data for estimating uncertain model parameters. The first problem utilizes spatial data from a stationary partial differential equation to produce a MUD estimate of an uncertain boundary condition. The second problem utilizes temporal data obtained from the state-of-the-art ADvanced CIRCulation (ADCIRC) model to obtain a MUD estimate of uncertain wind drag coefficients for a simulated extreme weather event near the Shinnecock Inlet located in the Outer Barrier of Long Island, NY, USA.

42 ENGINEERING↗
DOE OSTI2021

Comparison of the Forward and Reverse Photocycle Dynamics of Two Highly Similar Canonical Red/Green Cyanobacteriochromes Reveals Unexpected Differences

Cyanobacteriochromes (CBCRs) are cyanobacterial photoreceptors that exhibit photochromism between two states: a thermally stable dark-adapted state and a metastable light-adapted state with bound linear tetrapyrrole (bilin) chromophores possessing 15Z and 15E configurations, respectively. The photodynamics of canonical red/green CBCRs have been extensively studied; however, the time scales of their excited-state lifetimes and subsequent ground-state evolution rates widely differ and, at present, remain difficult to predict. Here, we compare the photodynamics of two closely related red/green CBCRs that have substantial sequence identity (~68%) and similar chromophore environments: AnPixJg2 from Anabaena sp. PCC 7120 and NpR6012g4 from Nostoc punctiforme. Using broadband transient absorption spectroscopy on the primary (125 fs to 7 ns) and secondary (7 ns to 10 ms) time scales together with global analysis modeling, our studies revealed that AnPixJg2 and NpR6012g4 have comparable quantum yields for initiating the forward ( 15Z P r → 15E P g ) and reverse ( 15E P g → 15Z P r ) reactions, which proceed through monotonic and nonmonotonic mechanisms, respectively. In addition to small discrepancies in the kinetics, the secondary reverse dynamics resolved unique features for each domain: intermediate shunts in NpR6012g4 and a Meta-Gf intermediate red-shifted from the 15Z P r photoproduct in AnPixJg2. Altogether, this study supports the conclusion that sequence similarity is a useful criterion for predicting pathways of the light-induced evolution and quantum yield of generating primary intermediate Φ p within subfamilies of CBCRs, but more studies are still needed to develop a comprehensive molecular level understanding of these processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Opening the Door to Grid-Interactive Efficient Buildings with Energy Codes

Building codes represent standard design practice in the construction industry and must continually evolve to account for advancing technologies and innovative practices. Energy codes have historically focused on energy efficiency within buildings and across their systems. However, many of tomorrow's technologies go beyond efficiency and target increased flexibility. These demand flexibility (DF) measures can postpone or reduce building electric load based on price or other grid signals, and include smart appliances, connected lighting, and connected mechanical systems. Expanding codes to enable such grid-interactive efficient building (GEB) has the potential to influence their realization at scale, which will support renewable-energy grid-integration and building-sector decarbonization. However, considering DF measures in code development creates a new set of challenges for codes and the practices contained therein. This paper considers the role of building energy codes in enabling GEB. Specifically, it reviews the status of demand flexibility (DF) measures in current commercial building energy code. It examines the national model code development process and identifies components barring the consideration and inclusion of DF measures - including code scope, characterization and analysis of proposed new prescriptive measures, and accounting for time-of-day and geographic differences in their benefits. The paper presents findings from code development analyses that indicate the cost benefit of DF measures and the limitations associated with current code development conventions. To encourage building flexibility and improved energy resilience moving forward, recommendations are made for removing code development barriers and sanctioning DF measure consideration in future model codes.

Building Energy Codes and Standards, Grid-Interact↗
DOE OSTI2023

β + Gamow-Teller Strengths from Unstable 14 O via the (d, 2 He) Reaction in Inverse Kinematics

For the first time, the (d, 2 He) reaction was successfully used in inverse kinematics to extract the Gamow-Teller transition strength in the β + direction from an unstable nucleus. The new technique was made possible by the use of an active-target time-projection chamber and a magnetic spectrometer, and opens a path to addressing a range of scientific challenges, including in astrophysics and neutrino physics. Here, in this Letter, the nucleus studied was 14 O, and the Gamow-Teller transition strength to 14 N was extracted up to an excitation energy of 22 MeV. The data were compared to shell-model and state-of-the-art coupled-cluster calculations. Shell-model calculations reproduce the measured Gamow-Teller strength distribution up to about 15 MeV reasonably well, after the application of a phenomenological quenching factor. In a significant step forward to better understand this quenching, the coupled-cluster calculation reproduces the full strength distribution well without such quenching, owing to the large model space, the inclusion of strong correlations, and the coupling of the weak interaction to two nucleons through two-body currents.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Many-body expansion based machine learning models for octahedral transition metal complexes

Abstract Graph-based machine learning (ML) models for material properties show great potential to accelerate virtual high-throughput screening of large chemical spaces. However, in their simplest forms, graph-based models do not include any 3D information and are unable to distinguish stereoisomers such as those arising from different orderings of ligands around a metal center in coordination complexes. In this work we present a modification to revised autocorrelation descriptors, a molecular graph featurization method, for predicting spin state dependent properties of octahedral transition metal complexes (TMCs). Inspired by analytical semi-empirical models for TMCs, the new modeling strategy is based on the many-body expansion (MBE) and allows one to tune the captured stereoisomer information by changing the truncation order of the MBE. We present the necessary modifications to include this approach in two commonly used ML methods, kernel ridge regression and feed-forward neural networks. On a test set composed of all possible isomers of binary TMCs, the best MBE models achieve mean absolute errors (MAEs) of 2.75 kcal mol −1 on spin-splitting energies and 0.26 eV on frontier orbital energy gaps, a 30%–40% reduction in error compared to models based on our previous approach. We also observe improved generalization to previously unseen ligands where the best-performing models exhibit MAEs of 4.00 kcal mol −1 (i.e. a 0.73 kcal mol −1 reduction) on the spin-splitting energies and 0.53 eV (i.e. a 0.10 eV reduction) on the frontier orbital energy gaps. Because the new approach incorporates insights from electronic structure theory, such as ligand additivity relationships, these models exhibit systematic generalization from homoleptic to heteroleptic complexes, allowing for efficient screening of TMC search spaces.

Meyer, Ralf (ORCID:0000000322360261)↗