Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse 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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

On the Feasibility of Geophysical Methods for CO 2 Monitoring in the North Dakota CarbonSAFE Project

Conference presentation at American Institute of Chemical Engineers (AIChE) Annual Meeting, Phoenix, AZ, November 13–18, 2022. An analysis of rock physics crossplots and diagnostic models suggests that the baseline seismic data may be sensitive to changes in lithology and porosity of the Broom Creek Formation and that these changes may be detectable by prestack seismic inversion.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Particle-in-cell simulations of expanding high energy density plasmas with laser ray tracing

The design and analysis of high energy density (HED) laser experiments typically rely on radiation hydrodynamics simulations. However, some laser–plasma interaction regimes are not collisional and cannot be adequately modeled with hydrodynamics. For example, strongly driven magnetic reconnection and magnetized collisionless shock experiments possess extended hydrodynamic or even kinetic properties, necessitating first-principles kinetic simulations. In this paper, we present the benchmarking and first results obtained with a laser-ray-tracing and inverse bremsstrahlung absorption module implemented in the particle-in-cell code PSC. The simulation results are compared to radiation hydrodynamic simulations using the FLASH code as well as analytical estimates. We successfully benchmark the energy deposition model and overall hydrodynamic evolution of the systems. We also consider possible kinetic effects that may be expected from laser-target ablation in the HED regime, including non-local transport and two-temperature effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Topography inversion in scanning tunneling microscopy of single-atom-thick materials from penetrating substrate states

Scanning tunneling microscopy (STM) is one of the indispensable tools to characterize surface structures, but the distinction between atomic geometry and electronic effects based on the measured tunneling current is not always straightforward. In particular, for single-atomic-thick materials (graphene or boron nitride) on metallic substrates, counterintuitive phenomena such as a larger tunneling current for insulators than for metal and a topography opposite to the atomic geometry are reported. Using first-principles density functional theory calculations combined with analytical modeling, we reveal the critical role of penetrating states of metallic substrates that surpass 2D material states, hindering the measurement of intrinsic 2D materials states and leading to topography inversion. Our finding should be instrumental in the interpretation of STM topographies of atomic-thick materials and in the development of 2D material for (opto)electronic and various quantum applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗
DOE OSTI2026

Machine learning inversion of interatomic force constants from single-crystal inelastic neutron scattering

Atomic vibrations govern many macroscopic properties of materials, but experiments to comprehensively probe them remain challenging. Inelastic neutron scattering (INS) is a powerful technique to map phonon dispersions in crystals, especially when leveraging modern time-of-flight (ToF) spectrometers with large detectors. However, efficiently and robustly extracting interatomic force constants (FCs) parameterizing phonon dynamics from experimental spectra remains a bottleneck due to the complexity and high dimensionality of ToF INS datasets. Here, we present a machine learning approach for the direct inversion of FCs from single-crystal INS measurements. The framework leverages synthetic training data generated using universal machine-learned force fields and an efficient physics-based forward model. We benchmark two neural architectures–one emphasizing structured latent representation learning and the other direct, supervised spectral regression–across simulated datasets for two materials under idealized and noisy conditions. The latent-representation model is subsequently applied to experimental single-crystal INS data on germanium. The model is shown to reproduce FCs derived from both first-principles simulations and from iterative optimization, and furthermore achieves reliable inference even from sparse, single-orientation measurements representing short data acquisitions. Analysis of the learned latent space reveals semantically continuous and physically interpretable encodings that support strong cross-domain generalization. By bridging theoretical and experimental domains, we establish a path toward rapid inversion of experimental spectra and data-driven interpretation of temperature-dependent lattice dynamics.

42 ENGINEERING↗
DOE OSTI2022

Evaluation of the Fitch Wind-Farm Wake Parameterization with Large-Eddy Simulations of Wakes Using the Weather Research and Forecasting Model

Abstract Wind-farm parameterizations in weather models can be used to predict both the power output and farm effects on the flow; however, their correctness has not been thoroughly assessed. We evaluate the wind-farm parameterization of the Weather Research and Forecasting Model with large-eddy simulations (LES) of the wake performed with the same model. We study the impact on the velocity and turbulence kinetic energy (TKE) of inflow velocity, roughness, resolution, number of turbines (one or two), and inversion height and strength. We compare the mesoscale with the LES by spatially averaging the LES within areas correspondent to the mesoscale horizontal spacing: one covering the turbine area and two downwind. We find an excellent agreement of the velocity within the turbine area between the two types of simulations. However, within the same area, we find the largest TKE discrepancies because in mesoscale simulations, the turbine-added TKE has to be highest at the turbine position to be advected downwind. Within the downwind areas, differences between velocities increase as the wake recovers faster in the LES, whereas for the TKE both types of simulations show similar levels. From the various configurations, the impact of inversion height and strength is small for these heights and inversion levels. The highest impact for the one-turbine simulations appears under the low-speed case due to the higher thrust, whereas the impact of resolution is low for the large-eddy simulations but high for the mesoscale simulations. Our findings demonstrate that higher-fidelity simulations are needed to validate wind-farm parameterizations.

17 WIND ENERGY↗
DOE OSTI2020

Search for TiO and Optical Nightside Emission from the Exoplanet WASP-33b

With a temperature akin to an M dwarf, WASP-33b is among the hottest Jupiters known, making it an ideal target for high-resolution optical spectroscopy. By analyzing both transmission and emission spectra, we aim to substantiate previous reports of atmospheric TiO and a thermal inversion within the planet’s atmosphere. We observed two transits and six arcs of the phase curve with the Echelle SpectroPolarimetric Device for the Observation of Stars (ESPaDOnS) on the Canada–France–Hawaii Telescope and High Resolution Echelle Spectrometer (HIRES) on the Keck telescope, which provide high spectral resolution and ample wavelength coverage. We employ the Doppler cross-correlation technique to search for the molecular signatures of TiO and H{sub 2}O in these spectra, using models based on the TiO line list of Plez. Though we cannot exclude line-list-dependent effects, our data do not corroborate previous indications of a thermal inversion. Instead we place a 3σ upper limit of 10{sup −9} on the volume mixing ratio of TiO for the T–P profile we consider. While we are unable to constrain the volume mixing ratio of water, our strongest constraint on TiO comes from dayside emission spectra. This apparent absence of a stratosphere sits in stark contrast to previous observations of WASP-33b as well as theoretical predictions for the atmospheres of highly irradiated planets. The discrepancy could be due to variances between line lists, and we stress that detection limits are only as good as the line list employed, and are only valid for the specific T–P profile considered due to the strong degeneracy between lapse rate (dT/dlogP) and molecular abundance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Physics Informed Neural Networks as Computational Physics Emulators

This report is a brief overview and evaluation of Physics Informed Neural Networks (PINNs). Karniadakis and co-workers, e.g., Karniadakis et al. (2021) assert that the PINNs approach integrates seamlessly both data and mathematical physics models, even in partially understood, uncertain and high-dimensional contexts. They further claim that PINNs are effective and efficient for ill-posed and inverse problems, and when combined with domain decomposition, are scalable to large problems and a tool to discover hidden physics. While they demonstrate the capabilities in specific academic instances, their overarching claims about PINNs seem to be an overstatement, at least at the current time. We have briefly considered a few of the limitations of PINNs in this investigation. It is not clear to us if the PINNs approach can ever be competitive with approaches that use specialized algorithms to achieve high-accuracy solutions of governing equations and other techniques that can combine observational data with such solutions. As an example of the latter, consistent with the principles of Bayesian inference, data assimilation, or more generally data-model fusion, is a process that fuses observational data typically with a computational model that respects certain constraints such as conservation laws. For example, improvements in observational network combined with data assimilation have been key in improving weather predictions over the past four decades Kalnay (2003).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Fossil fuel CO 2 emissions over metropolitan areas from space: A multi-model analysis of OCO-2 data over Lahore, Pakistan

Urban areas, where gathering more than 55% of the global population, alone contributed to more than 70% of anthropogenic fossil fuel carbon dioxide (CO 2ff ) emissions. Accurate quantification of CO 2ff emissions from urban areas is of great importance to the formulation of global warming mitigation policies to achieve carbon neutrality by 2050. Satellite-based inversion techniques are unique among “top-down” approaches, potentially allowing us to track CO 2ff emission changes over cities globally. However, its accuracy is still limited by incomplete background information, cloud blockages, aerosol contaminations, and uncertainties in models and emission inventories used as prior. To evaluate the current potential of space-based quantification techniques, we present the first attempt to monitor long-term changes in CO 2ff emissions based on the OCO-2 satellite measurements of column-averaged dry-air mole fractions of CO 2 (X CO2 ) over a fast-growing Asian metropolitan area: Lahore, Pakistan. We first examined the OCO-2 data availability at global scale. About 17% of OCO-2 soundings are marked as high-quality soundings by quality flags over the global 70 most populated cities over the period 2014-2019. Cloud blockage and aerosol contamination are the two main causes of data loss. As an attempt to recover additional soundings, we evaluated the effectiveness of OCO-2 quality flags at the city level by comparing three flux quantification methods (WRF-Chem, X-STILT, and flux cross-sectional integration method), all based on the Open-Data Inventory for Anthropogenic Carbon dioxide (ODIAC) product. The satellite/bottom-up emissions (OCO-2/ODIAC) ratios of the high-quality tracks better converged across the three methods compared to the all-data tracks with reduced uncertainties in emissions. Thus, OCO-2 quality flags are useful filters of low-quality OCO-2 retrievals at local scales, although originally designed for global-scale studies. All three methods consistently suggested that the ratio medians are greater than 1, which implies that the ODIAC slightly underestimated the CO2ff emissions over Lahore. Additionally, our estimation of the a posteriori CO2ff emission trend was about 734 kt C/year (i.e., an annual 6.7% increase). 10,000 Monte Carlo simulations of the Mann-Kendall upward trend test showed that less than 10% prior uncertainty for 8 tracks (or less than 20% prior uncertainty for 25 tracks) is required to achieve a greater-than-50% trend significant possibility at a 95% confidence level. It implies that the trend is driven by the prior and not due to the assimilation of OCO-2 retrievals. The key to improving the role of satellite data in CO 2 emission trend detection lies in collecting more frequent high-quality tracks near metropolitan areas to achieve significant constraints from X CO2 retrievals.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Black carbon aerosol number and mass concentration measurements by picosecond short-range elastic backscatter lidar

Abstract Black carbon aerosol emissions are recognized as contributors to global warming and air pollution. There remains, however, a lack of techniques to remotely measure black carbon aerosol particles with high range and time resolution. This article presents a direct and contact-free remote technique to estimate the black carbon aerosol number and mass concentration at a few meters from the emission source. This is done using the Colibri instrument based on a novel technique, referred to here as Picosecond Short-Range Elastic Backscatter Lidar (PSR-EBL). To address the complexity of retrieving lidar products at short measurement ranges, we apply a forward inversion method featuring radiometric lidar calibration. Our method is based on an extension of a well-established light-scattering model, the Rayleigh–Debye–Gans for Fractal-Aggregates (RDG-FA) theory, which computes an analytical expression of lidar parameters. These parameters are the backscattering cross-sections and the lidar ratio for black carbon fractal aggregates. Using a small-scale Jet A-1 kerosene pool fire, we demonstrate the ability of the technique to quantify the aerosol number and mass concentration with centimetre range-resolution and millisecond time-resolution.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Framework for simulating gauge theories with dipolar spin systems

Gauge theories appear broadly in physics, ranging from the standard model of particle physics to long-wavelength descriptions of topological systems in condensed matter. However, systems with sign problems are largely inaccessible to classical computations and also beyond the current limitations of digital quantum hardware. In this work, we develop an analog approach to simulating gauge theories with an experimental setup that employs dipolar spins (molecules or Rydberg atoms). We consider molecules fixed in space and interacting through dipole-dipole interactions, avoiding the need for itinerant degrees of freedom. Each molecule represents either a site or gauge degree of freedom, and Gauss's law is preserved by a direct and programmatic tuning of positions and internal state energies. This approach can be regarded as a form of analog systems programming and charts a path forward for near-term quantum simulation. As a first step, we numerically validate this scheme in a small-system study of U(1) quantum link models in (1+1) dimensions with link spin S = 1/2 and S =1 and illustrate how dynamical phenomena such as string inversion and string breaking could be observed in near-term experiments. Our work brings together methods from atomic and molecular physics, condensed matter physics, high-energy physics, and quantum information science for the study of nonperturbative processes in gauge theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Calculation of Critical Heat Flux Using an Inverse Heat Transfer Method to Support TREAT Experiment Analysis

Heat transfer between cladding and coolant during transient scenarios remains a critical area of uncertainty in understanding nuclear reactor safety. To advance the understanding of transient and accident scenarios involving critical heat flux (CHF), an in-pile experiment for the Transient Reactor Test (TREAT) facility at Idaho National Laboratory was developed. The experiment named CHF-Static Environment Rodlet Transient Test Apparatus (CHF-SERTTA), consists of a hollow borated stainless steel heater rod submerged in a static water pool which is heated via the (n, a) reaction in boron-10. This paper presents a novel inverse heat transfer method to determine CHF by using the optimization and uncertainty software DAKOTA to calibrate a RELAP5-3D model of CHF-SERTTA to temperature measurements obtained from a thermocouple welded to the surface of the rod.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Development of the detector array for photons, protons, and exotic residues

Accurate neutron capture predictions are needed in fundamental and applied science. Empirical models are constrained largely for stable nuclei, but predictions far from stability are needed. Neutron transfer in (d,p) reactions in inverse kinematics with rare isotope beams offer a means to access key aspects of the capture process. To this end, the detector array for photons, protons, and exotic residues (DAPPER) has been designed, built, tested, and used for measurement at Texas A&M University.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Machine learning framework for quantum sampling of highly constrained, continuous optimization problems

In recent years, there is growing interest in using quantum computers for solving combinatorial optimization problems. In this work, we developed a generic, machine learning-based framework for mapping continuous-space inverse design problems into surrogate quadratic unconstrained binary optimization (QUBO) problems by employing a binary variational autoencoder and a factorization machine. The factorization machine is trained as a low-dimensional, binary surrogate model for the continuous design space and sampled using various QUBO samplers. Using the D-Wave Advantage hybrid sampler and simulated annealing, we demonstrate that by repeated resampling and retraining of the factorization machine, our framework finds designs that exhibit figures of merit exceeding those of its training set. We showcase the framework’s performance on two inverse design problems by optimizing (i) thermal emitter topologies for thermophotovoltaic applications and (ii) diffractive meta-gratings for highly efficient beam steering. This technique can be further scaled to leverage future developments in quantum optimization to solve advanced inverse design problems for science and engineering applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Spillway preheating

In traditional models only an order one fraction of energy is transferred from the inflaton to radiation through nonperturbative resonance production in preheating immediately after inflation, due to backreaction effects. We propose a particle production mechanism that could improve the depletion of the inflaton energy density by up to four orders of magnitude. The improvement comes from the fast perturbative decays of resonantly produced daughter particles. They act as a “spillway” to drain these daughter particles, reducing their backreaction on the inflaton and keeping the resonant production effective for a longer period. Thus we dub the scenario “spillway preheating”. We also show that the fraction of energy density remaining in the inflaton has a simple inverse power-law scaling in the scenario. In general, spillway preheating is a much more efficient energy dissipation mechanism, which may have other applications in model building for particle physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Underlying simplicity of 5f unoccupied electronic structure

Using a simple empirical model based upon the bremsstrahlung isochromat spectroscopy of elemental Th, it is possible to explain the recent high energy resolution fluorescence detection measurements of UF 4 (n = 2) and UCd 11 (n = 3) as well as the new inverse photoelectron spectroscopy of Pu 2 O 3 (n = 5), where n is the 5f occupation number. Furthermore, a critical issue in this analysis is the assumption that the Th 5f states are essentially empty, which will be confirmed both experimentally and computationally. Thus, for 5f systems, this simple model provides a unified and consistent picture of 5f unoccupied density of states in simple, localized systems, as the 5f occupation varies in the early part of the series, for n = 0, 2, 3, and 5.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Learning generative neural networks with physics knowledge

Deep generative neural networks have enabled modeling complex distributions, but incorporating physics knowledge into the neural networks is still challenging and is at the core of current physics-based machine learning research. To this end, we propose a physics generative neural network (PhysGNN), a new class of generative neural networks for learning unknown distributions in a physical system described by partial differential equations (PDE). PhysGNN couples PDE systems with generative neural networks. It is a fully differentiable model that allows back-propagation of gradients through both numerical PDE solvers and generative neural networks, and is trained by minimizing the discrete Wasserstein distance between generated and observed probability distributions of the PDE outputs using the stochastic gradient descent method. Moreover, PhysGNN does not require adversarial training like standard generative neural networks, which offers better stability than adversarial training. We show that PhysGNN can learn complex distributions in stochastic inverse problems, where conventional methods such as maximum likelihood estimation and momentum matching methods may be inapplicable when little knowledge is known about the form of unknown distributions or the physical model is too complex. Furthermore, our method allows physics-based generative neural network training for learning complex distributions in the context of differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Direct CP violation and the ΔI=1/2 rule in K→ππ decay from the standard model

We present a lattice QCD calculation of the Δ I = 1 / 2 , K → π π decay amplitude A 0 and ϵ ' , the measure of direct C P violation in K → π π decay, improving our 2015 calculation [1] of these quantities. Both calculations were performed with physical kinematics on a 32 3 × 64 lattice with an inverse lattice spacing of a - 1 = 1.3784 ( 68 ) GeV . However, the current calculation includes nearly 4 times the statistics and numerous technical improvements allowing us to more reliably isolate the π π ground state and more accurately relate the lattice operators to those defined in the standard model. We find Re ( A 0 ) = 2.99 ( 0.32 ) ( 0.59 ) × 10 - 7 GeV and Im ( A 0 ) = - 6.98 ( 0.62 ) ( 1.44 ) × 10 - 11 GeV , where the errors are statistical and systematic, respectively. The former agrees well with the experimental result Re ( A 0 ) = 3.3201 ( 18 ) × 10 - 7 GeV . These results for A 0 can be combined with our earlier lattice calculation of A 2 [2] to obtain Re ( ϵ ' / ϵ ) = 21.7 ( 2.6 ) ( 6.2 ) ( 5.0 ) × 10 - 4 , where the third error represents omitted isospin breaking effects, and Re ( A 0 ) / Re ( A 2 ) = 19.9 ( 2.3 ) ( 4.4 ) . The first agrees well with the experimental result of Re ( ϵ ' / ϵ ) = 16.6 ( 2.3 ) × 10 - 4 . A comparison of the second with the observed ratio Re ( A 0 ) / Re ( A 2 ) = 22.45 ( 6 ) , demonstrates the standard model origin of this “ Δ I = 1 / 2 rule” enhancement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗