Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse prediction”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Metal dissolution from first principles: Potential-dependent kinetics and charge transfer

Metal dissolution and its inverse process are integral to both corrosion and electrodeposition; however, many mechanistic details regarding the dissolution process are challenging to decipher. These include how ion dissolution kinetics and charge transfer are influenced by the competition between metal and solvent interactions under an electrode potential. Here, in this work, we introduce a computational framework based on density functional theory with grand-canonical treatment of electrons to directly predict the potential energy landscape for metal dissolution at a constant potential. Using aluminum as an example, we demonstrate that dissolution kinetics is governed by competing kinetics between two physical processes associated with metal–metal bond breaking and ion-migration within the electrical double layer, respectively. We identify a kinetic transition between regimes dominated by each of these processes and show that this transition depends on the operating electrode potential, among other key factors. It is further found that kinetics and thermodynamics of these processes can be described with a simple, one-parameter Marcus-theory-type model. Beyond offering new understanding of charge transfer during dissolution, our simulation protocol provides a recipe for directly predicting other important quantities in electrochemical reactions from first principles that are difficult to measure, such as the symmetry factor.

36 MATERIALS SCIENCE↗
NASA NTRS1990

The analysis and simulation of compressible turbulence

Compressible turbulent flows at low turbulent Mach numbers are considered. Contrary to the general belief that such flows are almost incompressible, (i.e., the divergence of the velocity field remains small for all times), it is shown that even if the divergence of the initial velocity field is negligibly small, it can grow rapidly on a non-dimensional time scale which is the inverse of the fluctuating Mach number. An asymptotic theory which enables one to obtain a description of the flow in terms of its divergence-free and vorticity-free components has been developed to solve the initial-value problem. As a result, the various types of low Mach number turbulent regimes have been classified with respect to the initial conditions. Formulae are derived that accurately predict the level of compressibility after the initial transients have disappeared. These results are verified by extensive direct numerical simulations of isotropic turbulence.

Erlebacher, Gordon↗
NASA NTRS1990

The analysis and simulation of compressible turbulence

Compressible turbulent flows at low turbulent Mach numbers are considered. Contrary to the general belief that such flows are almost incompressible (i.e., the divergence of the velocity field remains small for all times), it is shown that even if the divergence of the initial velocity field is negligibly small, it can grow rapidly on a nondimensional time scale which is the inverse of the fluctuating Mach number. An asymptotic theory which enables one to obtain a description of the flow in terms of its divergence-free and vorticity-free components has been developed to solve the initial-value problem. As a result, the various types of low Mach number turbulent regimes have been classified with respect to the initial conditions. Formulae are derived that accurately predict the level of compressibility after the initial transients have disappeared. These results are verified by extensive direct numerical simulations of isotropic turbulence.

Erlebacher, Gordon↗
NASA NTRS2010

Ares I Reaction Control System Propellant Feedline Decontamination Modeling

The objective of the work presented here is to quantify the effects of purge gas temperature, pressure, and mass flow rate on Hydrazine (Hz) decontamination rates of the Ares I Roll Control System and Reaction Control System. A survey of experts in this field revealed the absence of any decontamination rate prediction models. Three basic decontamination methods were identified for analysis and modeling. These include low pressure eduction, high flow rate purge, and pulse purge. For each method, an approach to predict the Hz mass transfer rate, as a function of system pressure, temperature, and purge gas mass flow rate, is developed based on the applicable physics. The models show that low pressure eduction is two orders of magnitude more effective than the high velocity purge, which in turn is two orders of magnitude more effective than the pure diffusion component of pulse purging of deadheads. Eduction subjects the system to low pressure conditions that promote the extraction of Hz vapors. At 120 F, Hz is saturated at approximately 1 psia. At lower pressures and 120 F, Hz will boil, which is an extremely efficient means to remove liquid Hz. The Hz boiling rate is predicted by equating the rate at which energy is added to the saturated liquid Hz through heaters at the tube outer wall with the energy removed from the liquid through evaporation. Boil-off fluxes were predicted by iterating through the range of local pressures with limits set by the minimum allowed pressure of 0.2 psia and maximum allowed wall temperature of 120 F established by the heaters, which gives a saturation pressure of approximately 1.0 psia. Figure 1 shows the resulting boil-off fluxes as a function of local eduction pressure. As depicted in figure 1, the flux is a strong inverse function of eduction pressure, and that minimizing the eduction pressure maximizes the boil-off flux. Also, higher outer wall temperatures lead to higher boil-off fluxes and allow for boil-off over a greater range of eduction pressures.

Pasch, James J.↗
DOE OSTI2022

Scaled Vecchia Approximation for Fast Computer-Model Emulation

Many scientific phenomena are studied using computer experiments consisting of multiple runs of a computer model while varying the input settings. Gaussian processes (GPs) are a popular tool for the analysis of computer experiments, enabling interpolation between input settings, but direct GP inference is computationally infeasible for large datasets. We adapt and extend a powerful class of GP methods from spatial statistics to enable the scalable analysis and emulation of large computer experiments. Specifically, we apply Vecchia’s ordered conditional approximation in a transformed input space, with each input scaled according to how strongly it relates to the computer-model response. The scaling is learned from the data by estimating parameters in the GP covariance function using Fisher scoring. Our methods are highly scalable, enabling estimation, joint prediction, and simulation in near-linear time in the number of model runs. In several numerical examples, our approach substantially outperformed existing methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Diffusion Model-Guided Inverse Design of Bimetallic Catalysts for Ammonia Decomposition

In the past decade, artificial intelligence and deep learning have played increasingly prominent roles in materials design and discovery. Among these, generative AI models, known for their ability to create unique and complex structures, have emerged as state-of-the-art tools for materials screening due to their high efficiency and low computational cost. In catalysis, one of the major challenges is identifying promising material candidates within an immense chemical space. This challenge can be addressed using generative approaches, such as diffusion-based inverse design models. In this study, we present a machine learning-guided workflow that employed a diffusion model for the inverse design of bimetallic alloy catalysts for low-carbon ammonia decomposition, a key reaction for ammonia emission control and sustainable hydrogen production. Catalyst candidates were evaluated using nitrogen adsorption energy as the key descriptor, inspired by multiscale modeling. The proposed workflow identified low-cost, environmentally friendly catalysts with excellent catalytic performance, which have been validated theoretically and experimentally. Our framework decoupled the generative and property-prediction components, enhancing both flexibility and accuracy in the catalytic material design process.

Adsorption↗
DOE OSTI2021

Emergence of spin singlets with inhomogeneous gaps in the kagome lattice Heisenberg antiferromagnets Zn-barlowite and herbertsmithite

The kagome Heisenberg antiferromagnet formed by frustrated spins arranged in a lattice of corner-sharing triangles is a prime candidate for hosting a quantum spin liquid (QSL) ground state consisting of entangled spin singlets. However, the existence of various competing states makes a convincing theoretical prediction of the QSL ground state difficult, calling for experimental clues from model materials. The kagome lattice materials Zn-barlowite (ZnCu 3 (OD) 6 FBr) and herbertsmithite (ZnCu 3 (OD) 6 Cl 2 ) do not exhibit long-range order and are considered the best realizations of the kagome Heisenberg antiferromagnet known so far. In this study, we use 63 Cu nuclear quadrupole resonance combined with the inverse Laplace transform to locally probe the inhomogeneity of delicate quantum ground states affected by disorder. We present direct evidence for the gradual emergence of spin singlets with spatially varying excitation gaps, but even at temperatures far below the super-exchange energy scale their fraction is limited to ~60% of the total spins. Theoretical models need to incorporate the role of disorder to account for the observed inhomogeneously gapped behaviour.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Machine learning and LHC event generation

First-principle simulations are at the heart of the high-energy physics research program. They link the vast data output of multi-purpose detectors with fundamental theory predictions and interpretation. This review illustrates a wide range of applications of modern machine learning to event generation and simulation-based inference, including conceptional developments driven by the specific requirements of particle physics. New ideas and tools developed at the interface of particle physics and machine learning will improve the speed and precision of forward simulations, handle the complexity of collision data, and enhance inference as an inverse simulation problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2000

Nonlinear Shock Acceleration and Photon Emission in Supernova Remnants

We have extended a simple model of nonlinear diffusive shock acceleration (Berezhko & Ellison 1999: Ellison &, Berezhko 1999a) to include the injection and acceleration of electrons and the production of photons from bremsstrahlung, synchrotron, inverse Compton, and pion-decay processes. We argue that, the results of this model, which is simpler to use than more elaborate ones, offer a significant improvement, over test-particle, power-law spectra which are often used in astrophysical applications of diffusive shock acceleration. With an evolutionary supernova remnant (SNR) model to obtain shock parameters as functions of ambient interstellar medium parameters and time, we predict broad-band continuum photon emission from supernova remnants in general, and SN1006 in particular, showing that our results compare well with the more complete time-dependent and spherically symmetric nonlinear model of Berezhko, Ksenofontov, & Petukhov (1999a). We discuss the implications nonlinear shock acceleration has for X-ray line emission, and use our model to describe how ambient conditions determine the TeV/radio flux ratio, an important parameter for gamma-ray observations of radio SNRs.

Ellison, Donald C.↗
DOE OSTI2022

Latitude, Elevation, and Mean Annual Temperature Predict Peat Organic Matter Chemistry at a Global Scale

Peatlands contain a significant fraction of global soil carbon, but how these reservoirs will respond to the changing climate is still relatively unknown. A global picture of the variations in peat organic matter chemistry will aid our ability to gauge peatland soil response to climate. Here, the goal of this research is to test the hypotheses that (a) peat carbohydrate content, an indicator of soil organic matter reactivity, will increase with latitude and decrease with mean annual temperatures, (b) while peat aromatic content, an indicator of recalcitrance, will vary inversely, and (c) elevation will have a similar effect to latitude. We used Fourier Transform Infrared Spectroscopy to examine variations in the organic matter functional groups of 1034 peat samples collected from 10 to 20, 30–40, and 60–70 cm depths at 165 individual sites across a latitudinal gradient of 79°N–65°S and from elevations of 0–4,773 m. Carbohydrate contents of high latitude peat were significantly greater than peat originating near the equator, while aromatic content showed the opposite trend. For peat from similar latitudes but different elevations, the carbohydrate content was greater and aromatic content was lower at higher elevations. Higher carbohydrate content at higher latitudes indicates a greater potential for mineralization, whereas the chemical composition of low latitude peat is consistent with their apparent relative stability in the face of warmer temperatures. The combination of low carbohydrates and high aromatics at warmer locations near the equator suggests the mineralization of high latitude peat until reaching recalcitrance under a new temperature regime.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Constraining baryon annihilation in the hadronic phase of heavy-ion collisions via event-by-event fluctuations

We point out that the variance of net-baryon distribution normalized by the Skellam distribution baseline, κ 2 [$B-\overline{B}$]/$\langle B+\overline{B} \rangle$, is sensitive to the possible modification of (anti)baryon yields due to $B\overline{B}$ annihilation in the hadronic phase. The corresponding measurements can thus place stringent limits on the magnitude of the $B\overline{B}$ annihilation and its inverse reaction. We perform Monte Carlo simulations of the hadronic phase in Pb-Pb collisions at the LHC via the recently developed subensemble sampler + UrQMD afterburner and show that the effect survives in net-proton fluctuations, which are directly accessible experimentally. The available experimental data of the ALICE Collaboration on net-proton fluctuations disfavors a notable suppression of (anti)baryon yields in $B\overline{B}$ annihilations predicted by the present version of UrQMD if only global baryon conservation is incorporated. On the other hand, the annihilations improve the data description when local baryon conservation is imposed. The two effects can be disentangled by measuring κ 2 [$B-\overline{B}$]/$\langle B+\overline{B} \rangle$, which at the LHC is notably suppressed by annihilations but virtually unaffected by baryon number conservation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior

The traditional goal of inverse self-assembly is to design interactions that drive particles toward a desired target structure. However, achieving successful self-assembly also requires tuning the thermodynamic conditions under which the structure is stable. In this work, we extend the inverse design paradigm to explicitly address this challenge by developing a framework for inverse thermodynamics, i.e., the design of interaction potentials that realize specific thermodynamic behavior. As a step in this direction, using patchy particle mixtures as a model system, we demonstrate how precise control over both bonding topology and bond energetics enables the programming of targeted phase behavior. In particular, we establish design principles for azeotropic demixing and show how to create mixtures that exhibit azeotropy at any prescribed composition. Our predictions are validated through Gibbs-ensemble simulations [Panagiotopoulos, Mol. Phys. 1987, 61, 813−826]. These results highlight the necessity of coupling structural design with thermodynamic engineering, and provide a blueprint for controlling complex phase behavior in multicomponent systems.

Azeotropes↗
DOE OSTI2020

Broadband Tunable Mid-infrared Plasmon Resonances in Cadmium Oxide Nanocrystals Induced by Size-Dependent Nonstoichiometry

A central theme of nanocrystal (NC) research involves synthesis of dimension-controlled NCs and study of size-dependent scaling laws governing their optical, electrical, magnetic, and thermodynamic properties. In this paper, we describe the synthesis of monodisperse CdO NCs that exhibit high quality-factor (up to 5.5) mid-infrared (MIR) localized surface plasmon resonances (LSPR) and elucidate the inverse scaling relationship between carrier concentration and NC size. The LSPR wavelength is readily tunable between 2.4 and similar to 6.0 μm by controlling the size of CdO NCs. Structural and spectroscopic characterization provide strong evidence that free electrons primarily originate from self-doping due to NC surface-induced nonstoichiometry. The ability to probe and to control NC stoichiometry and intrinsic defects will pave the way toward predictive synthesis of doped NCs with desirable LSPR characteristics.

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

Using Machine Learning to Predict Cloud Turbulent Entrainment–Mixing Processes

Different turbulent entrainment–mixing mechanisms between clouds and environment are essential to cloud–related processes; however, accurate representation of entrainment–mixing in weather/climate models still poses a challenge. This study exploits the use of machine learning (ML) to address this challenge. Four ML (Light Gradient Boosting Machine [LGB], eXtreme Gradient Boosting, Random Forest, and Support Vector Regression) are examined and compared. It is found that LGB performs best, and thus is selected to understand the impact of entrainment–mixing on microphysics using simulation data from Explicit Mixing Parcel Model. Compared with traditional parameterizations, the trained LGB provides more accurate microphysical properties (number concentration and cloud droplet spectral dispersion). The partial dependences of predicted microphysics on features exhibit a strong alignment with physical mechanisms and expectations, as determined by the interpreting method, thus overcoming the limitations of the “black box” scheme. The underlying mechanisms are that the smaller number concentration and larger spectral dispersion correspond to more inhomogeneous entrainment–mixing. Specifically, number concentration after entrainment–mixing is positively correlated with adiabatic number concentration and liquid water content affected by entrainment–mixing, and inversely correlated with adiabatic volume mean radius. Spectral dispersion after entrainment–mixing is negatively correlated with liquid water content affected by entrainment–mixing, turbulent dissipation rate and relative humidity of entrained air. Sensitivity analysis further suggests that number concentration is mainly determined by cloud microphysical properties whereas spectral dispersion is influenced by both cloud microphysical properties and environmental variables. The results indicate that the LGB scheme has the potential to enhance the representation of entrainment–mixing in weather/climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Electron transfer in nonpolar media

The Marcus theory of electron transfer views fluctuating orientations of permanent dipoles as the nuclear mode bringing the donor and acceptor into the tunneling resonance. Electronic polarization of the solvent is excluded as the fast mode adiabatically following the electronic density. This view, valid for solids, does not apply to molecular liquids where molecular translations (density fluctuations) modulate the induction interaction of the donor-acceptor complex with the solvent. This mechanism of promoting radiation-less electronic transitions is considered here in the framework of the perturbation liquid-state theory. The reorganization energy of electron transfer in nonpolar solvents is nonzero and reaches the values of 0.1–0.3 eV for typical molecular sizes and solvents used in applications. The reorganization energy scales quadratically with the molecular polarizability of the solvent and decays as the inverse fifth power with the size of the donor and acceptor. The combination of the entropic character of the density fluctuations, driven by re-packing of molecular cores, with the short range of induction solute–solvent interactions leads to the violation of the fluctuation–dissipation theorem for the variance of the donor-acceptor energy gap. An explicit, approximately hyperbolic, dependence of the reorganization energy on temperature is predicted. It leads to a non-Arrhenius kinetic law for the rate constant of electron transfer in nonpolar liquid solvents.

14 SOLAR ENERGY↗
DOE OSTI2022

Clump survival and migration in VDI galaxies: an analytical model versus simulations and observations

ABSTRACT We address the nature of the giant clumps in high-z galaxies that undergo violent disc instability, distinguishing between long-lived and short-lived clumps. We study the evolution of long-lived clumps during migration through the disc via an analytical model tested by simulations and confront theory with CANDELS-HST observations. The clump ‘bathtub’ model, which considers gas and stellar gain and loss, involves four parameters: the accretion efficiency α, the star formation rate (SFR) efficiency ϵd, and the outflow mass-loading factors for gas and stars, η and ηs. The corresponding time-scales are comparable to the migration time, two-three orbital times. The accretion-rate dependence on clump mass, gas, and stars, allows an analytical solution involving exponential growing and decaying modes. For the fiducial parameter values there is a main evolution phase where the SFR and gas mass are constant and the stellar mass is rising linearly with time. This makes the inverse specific SFR an observable proxy for clump age. When η or ϵd are high, or α is low, the decaying mode induces a decline of SFR and gas mass till the migration ends. Later, the masses and SFR approach an hypothetical exponential growth with a constant specific SFR. The model matches simulations with different, moderate feedbacks, both in isolated and cosmological settings. The observed clumps agree with our predictions, indicating that the massive clumps are long-lived and migrating. A challenge is to model feedback that is non-disruptive in massive clumps but suppresses SFR to match the galactic stellar-to-halo mass ratio.

Dekel, Avishai↗
DOE OSTI2022

Parameter estimation for X-ray scattering analysis with Hamiltonian Markov Chain Monte Carlo

Bayesian-inference-based approaches, in particular the random-walk Markov Chain Monte Carlo (MCMC) method, have received much attention recently for X-ray scattering analysis. Hamiltonian MCMC, a state-of-the-art development in the field of MCMC, has become popular in recent years. It utilizes Hamiltonian dynamics for indirect but much more efficient drawings of the model parameters. We described the principle of the Hamiltonian MCMC for inversion problems in X-ray scattering analysis by estimating high-dimensional models for several motivating scenarios in small-angle X-ray scattering, reflectivity, and X-ray fluorescence holography. Hamiltonian MCMC with appropriate preconditioning can deliver superior performance over the random-walk MCMC, and thus can be used as an efficient tool for the statistical analysis of the parameter distributions, as well as model predictions and confidence analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

On the Properties of Cosmological Ionization Fronts

We investigate the properties of cosmological ionization fronts during the epoch of reionization using the Cosmic Reionization on Computers simulations. By analyzing reionization timing maps, we characterize ionization front velocities and curvatures and their dependence on the density structure of the intergalactic medium (IGM). The velocity distribution of ionization fronts in the simulations indicates that while the barrier-crossing analytical model captures the overall shape in high-velocity regions, it fails to reproduce the low-velocity tail, highlighting the non-Gaussian nature of the IGM’s density field. Ionization front velocities are inversely correlated with local density, propagating faster in underdense regions and more slowly in overdense environments. Faster ionization fronts also lead to higher postionization temperatures, reaching a plateau at ∼2 × 10$^{4}$ K for velocities exceeding 3000 km s$^{−1}$. Examining curvature statistics further establishes a connection between ionization front structure and the normalized density contrast ν, with trends in overdense regions aligning well with barrier-crossing model predictions, while deviations appear in underdense environments due to model limitations. These results provide a detailed characterization of ionization front dynamics and their interaction with the underlying density field, bridging small-scale reionization physics with large-scale observables such as the 21 cm signal and the IGM’s thermal history.

79 ASTRONOMY AND ASTROPHYSICS↗