Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse methods”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Improving the MJO Forecast of S2S Operation Models by Correcting Their Biases in Linear Dynamics

The operational dynamic subseasonal to seasonal (S2S) models for Madden-Julian oscillation (MJO) forecasting mostly still suffer from systematic errors in capturing the MJO's key dynamic features, such as its growth rate and propagation speed. By deriving the linear dynamic operators using the linear inverse modeling (LIM) approach, we propose a method to partly correct the errors in MJO linear dynamic operators to improve the MJO predictions of three operational dynamic S2S models. Correcting the deficiencies of the too-fast decay rates and the unrealistic propagating phase speeds lead to MJO prediction skills being extended by approximately 2–4 days. The improvements are more significant for the models with larger biases in MJO amplitude and propagation. This approach in principle may be extendable to predictions of other types of climate variability such as ENSO on one hand, and possible inclusions of nonlinear dynamics effects on the other hand.

58 GEOSCIENCES↗
DOE OSTI2023

Applying Bayesian inference and deterministic anisotropy to retrieve the molecular structure ∣Ψ(R)∣2 distribution from gas-phase diffraction experiments

Abstract Currently, our general approach to retrieving molecular structures from ultrafast gas-phase diffraction heavily relies on complex ab initio electronic or vibrational excited state simulations to make conclusive interpretations. Without such simulations, inverting this measurement for the structural probability distribution is typically intractable. This creates a so-called inverse problem. Here we address this inverse problem by developing a broadly applicable method that approximates the molecular frame structure ∣Ψ( R , t )∣ 2 distribution independent of these complex simulations. We retrieve the vibronic ground state ∣Ψ( R )∣ 2 for both simulated stretched NO 2 and measured N 2 O. From measured N 2 O, we observe 40 mÅ coordinate-space resolution from 3.75 Å −1 reciprocal space range and poor signal-to-noise, a 50X improvement over traditional Fourier transform methods. In simulated NO 2 diffraction experiments, typical to high signal-to-noise levels predict 100–1000X resolution improvements, down to 0.1 mÅ. By directly measuring the width of ∣Ψ( R )∣ 2 , we open ultrafast gas-phase diffraction capabilities to measurements beyond current analysis approaches. This method has the potential to effectively turn gas-phase ultrafast diffraction into a discovery-oriented technique to probe systems that are prohibitively difficult to simulate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Effect of lithium incorporation on tweaking the electrocatalytic behavior of tantalum-based oxides

To study the effect of lithium incorporation on the electrochemical properties and catalytic behavior of tantalum-based oxides, a series of mesoporous lithium-contained tantalum oxides was developed adopting an inverse micelle-based sol gel synthesis method. Four various phases of amorphous Ta 2 O 5 , polycrystalline LiTaO 3 , Li 5.04 Ta 6.16 O 17.92 , and Li 3 TaO 4 were determined for samples synthesized at Li/Ta ratios of 0.5 (LTO-0.5), 1 (LTO-1), 2 (LTO-2), and 3 (LTO-3) respectively. The surface analysis revealed that LTO-2 sample is featured with lower amount of lithium atoms at the surface rather than the bulk of the sample. The electrochemical analysis of LTO samples showed that LTO-2 possesses lower charge transfer resistance extracted from impedance spectroscopy that consequently leads to a better electrocatalytic performance. The LTO-2 electrode also showed a low open circuit potential (OCP) combined with a relatively higher electroactive surface area, thus capable of running an electrooxidation reaction at a potential lower than all other LTO samples. The more negative OCP value of LTO-2 electrode in AA solution states a more facile electron transfer between the analyte and electrode compared to other synthesized samples. Mott–Schottky plots of LTO electrodes are characterized by positive slopes which indicate n-type semiconductor behavior. The charge carrier density (5.5 × 10 19 cm -3 ) in LTO-2 sample is the lowest among all other samples. The charge carrier density is dependent on the amount of applied lithium in the synthesis confirming various lithium diffusion behaviors in LTO samples. The lower carrier concentration in LTO-2 justifies its higher activity in electrooxidation providing less populated lower unoccupied molecular orbital (LUMO) for the accommodation of electrons extracted from the analyte higher occupied molecular orbital (HOMO). Therefore, LTO-2 showed an improved electrocatalytic performance due to lower charge transfer resistance (higher conductivity) offered by lower amount of lithium atom at its surface and higher oxygen vacancies in the bulk.

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

Application of machine learning techniques for fast MeV x-ray spectra unfolding from filter stack spectrometer data

Recovery of MeV x-ray spectra from detector signals is difficult because the response matrix inversion is ill-conditioned and current methods are too slow for high-repetition-rate experiments. In this work, we make use of neural networks to unfold MeV x-ray spectra from measurements obtained with a filter stack spectrometer at rates of near 40 Hz. The neural network was trained on synthetic data and tested on both synthetic and experimental data, the latter obtained in two separate experiments performed at the Omega EP laser facility. We show here that this unfolding method has good performance on synthetic data and that it is a promising option for experimental data of up to 40 MeV. The accuracy on experimental data is verified by using a simple forward model to compare against measured values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Desmearing Bonse–Hart USANS data using Bayesian Gaussian process regression

Ultra-small-angle neutron scattering (USANS) enables access to micrometer-scale structures but is intrinsically affected by strong, anisotropic resolution smearing arising from slit-geometry optics. As a result, recovery of the intrinsic scattering intensity constitutes an ill-posed inverse problem, and commonly used iterative desmearing methods lack rigorous uncertainty quantification. We present a Bayesian desmearing framework for slit-geometry USANS based on Gaussian process regression. In this approach, the scattering intensity is modeled as a smooth random function, and the instrumental point spread function is incorporated explicitly as a forward operator. The resulting formulation yields a closed-form maximum a posteriori solution with well-defined credibility intervals. Computational benchmarks and experimental validation using combined USANS and small-angle neutron scattering (SANS) measurements demonstrate that the framework enables stable desmearing, suppresses experimental noise, and preserves physically meaningful structural features under realistic conditions.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗
DOE OSTI2025

Collective Modes in Multilayer Graphene/𝛼−RuCl 3 Heterostructures

Collective modes in multilayer graphene, such as plasmons and phonons, exhibit sensitivity to displacement fields and interlayer coupling, distinguishing them from their counterparts in single-layer graphene. Here, we engineer collective modes in charge-transfer heterostructures composed of multilayer graphene and 𝛼−RuCl 3 . In heterostructures with a single 𝛼−RuCl 3 interface, the charge transfer generates displacement fields up to 7 V/nm at the interface between 𝛼−RuCl 3 and the adjacent graphene layer—the highest value achieved through charge-transfer methods. As a result of the broken inversion symmetry, we discover enhanced nonlinear optical response and modified phonon selection rules. Conversely, we find that multilayer graphene sandwiched between two 𝛼−RuCl 3 flakes causes displacement fields to cancel. There, we achieve carrier densities as high as 8 ×10 13 cm −2 in multilayer graphene and restore the phonon selection rules to their unperturbed state. Meanwhile, we demonstrate that plasmonic properties derive from the depletion of multiple valence bands. As a result of the quasilinear band dispersion, these “Dirac multiband plasmons” are relatively unaffected by displacement fields. On the other hand, the inverted heterostructure sequence—two multilayer graphene sheets encapsulating 𝛼−RuCl 3 —activates significant alteration of the plasmons via interlayer plasmon-plasmon coupling. Hence, multilayer graphene and 𝛼−RuCl 3 heterostructures offer a gate-free platform for engineering collective modes derived from inversion symmetry and interlayer coupling.

Moore, Samuel L. [Columbia University, New York, N↗
DOE OSTI2021

pnnl/neuromancer

Dynamics-based deep learning methods to modernize current scientific computing methods. Neuromancer is currently capable of solving inverse problems for a system of ordinary differential equations. The functionality includes system identification and constrained optimal control of unknown or partially known ODEs.

Skomski, Elliott↗
DOE OSTI2025

pnnl/MCRASTA

McRasta (Markov Chain Rate and State Analysis) was developed to estimate parameter uncertainty in constitutive friction models via Bayesian inverse and Markov Chain Monte Carlo (MCMC) methods.

Fichera, Marissa [Pacific Northwest National Labor↗
DOE OSTI2024

Seismic Monitoring near Ithaca, New York, Reveals Nonuniform Distribution of Microseismicity in an Intraplate Region

Abstract Cornell University intends to use a deep direct-use geothermal system to heat its Ithaca, New York, campus. In preparation for this project, the Cornell Seismic Network has been monitoring the background seismicity in this intraplate region since 2019. From January 2020 to June 2023, 95 events were detected within 20 km of the proposed geothermal well site, with local magnitudes ranging from −1.02 to 0.56. None of these events appear in regional or national catalogs. Events locate in a narrow geographic band, with one-fourth exhibiting multimodal hypocentral probability peaks both near the surface and at 1–4 km depth. We relocate events with a joint hypocenter and 1D velocity model inversion, in addition to a fully nonlinear method, and then compare observations with synthetic waveforms. Together, these approaches provide strong evidence for >95% of events locating at the surface or within the 3-km-thick sedimentary sequence. We explore how anthropogenic activity and regional topographic stress may contribute to frequent surficial events. This information is critical for characterizing the background microseismicity for comparison during future geothermal operations. Ithaca’s geology of Paleozoic sediments overlying Precambrian crystalline basement is typical of many continental interiors, so these results also provide insight into intraplate microseismicity patterns.

Geochemistry & Geophysics↗
DOE OSTI2021

Nuclear Data to Reduce Uncertainties in Reactor Antineutrino Measurements

The large quantities of antineutrinos produced through the decay of fission fragments in nuclear reactors provide an opportunity to study the properties of these particles and investigate their use in reactor monitoring. The reactor antineutrino spectra are measured using specialized, large area detectors that detect antineutrinos through inverse beta decay, electron elastic scattering, or coherent elastic neutrino nucleus scattering; although, inverse beta decay is the only demonstrated method so far. Reactor monitoring takes advantage of the differences in the antineutrino yield and spectra resulting from uranium and plutonium fission providing an opportunity to estimate the fissile material composition in the reactor. Recent experiments reveal a deviation between the measured and calculated antineutrino flux and spectra (the reactor anomaly) indicating either the existence of yet undiscovered neutrino physics, uncertainties in the reactor source term calculation, incorrect nuclear data, or a combination of all three. To address the nuclear data that impact the antineutrino spectrum calculations and measurements, an international group of over 180 experts in antineutrino physics, reactor analysis, detector development, and nuclear data came together during the Workshop on Nuclear Data for Reactor Antineutrino Measurements (WoNDRAM) to discuss nuclear data needs and achieve concordance on a set of recommended priorities for nuclear data improvements. Three topical sessions focused on the reactor source term, the antineutrino spectrum, and the detector response, provided a forum to gain consensus amongst the participants on the most important data improvements to address two goals: 1) understand the reactor anomaly and 2) improve the ability to monitor reactors using antineutrinos. This report summarizes the outcomes of the workshop discussions and the recommendations for nuclear data efforts that reduce reactor antineutrino measurement uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Advancing Development of Emissions Detection (Final Report)

This document is the final report to the U.S. Department of Energy (DOE for contract DE-FE0031873) awarded to Colorado State University (CSU). CSU and partners at Harrisburg University of Science and Technology, University of Texas Arlington, and University of Texas at Austin organized several testing rounds to provide knowledge focused on advancing the detection capabilities of emissions monitoring devices. With this funding opportunity the research group began by establishing the Advancing the Development of Emission Detection (ADED) program, with the goal focused on enhancing the accuracy, reliability, and field applicability of methane detection technologies. This effort aimed to address the critical challenges of identifying and quantifying methane emissions while enabling industry stakeholders to meet regulatory compliance and environmental sustainability goals. The program engaged with industry, government, and technology stakeholders to promote adoption and consensus on testing techniques for methane detection solutions. Methane, a potent greenhouse gas, contributes significantly to global warming, and the oil and natural gas (O&G) sector is a primary source of methane emissions. Regulatory measures such as leak detection and repair (LDAR) programs have been implemented to address emissions. However, traditional LDAR approaches, reliant on handheld and component-level measurements, are resource-intensive. To address these limitations and move with evolving regulations, advanced methane technologies are emerging. These solutions include ground-based sensors, mobile systems (e.g., drones, vehicles, and aircraft), and satellite-based platforms. They offer innovative capabilities for autonomous monitoring, larger spatial coverage, and emission quantification using methods such as tracer gas techniques and inverse modeling with Gaussian plume analysis. The ADED program began with creating protocols for methane controlled release (CR) testing these continuous monitoring (CM) and survey technologies that detect and monitor methane emissions at O&G facilities. The protocols were then implemented throughout testing of CM and survey devices at CSU’s Methane Emissions Technology Evaluation Center (METEC) facility from 2021 through 2024. As apart of the protocol, solutions that tested under the ADED program installed their solutions at METEC, documented their system under test, and provided detection reports to the METEC team for analysis. The METEC team would provide the solutions with analyzed reports of their emissions and ground truth data of the releases conducted during their testing session. Under the ADED program, CMs were also tested at O&G facilities for a six week test run of challenge release (ChR) releases. The findings from the ADED program underscore the critical role of collaborative research and innovation in tackling methane emissions, offering a pathway for the oil and gas sector to achieve significant environmental and economic benefits. Results from METEC testing saw improvement of performance and accuracy across all solutions over the extent of the ADED experiments. The results also showed a variance in CM solution performance between CRs and ChRs. That variance pushed the team to further analyze the differences between CR testing environments and field conditions. With the drive from regulations and that variance in field conditions, the ADED team began designing a new CR testing protocol and additions to the METEC testing facility. The METEC team is furthering the progress made through the ADED program with awarded funding from DE-FE0032276. This funding pushes the development of METEC’s addition with new equipment, allowing for an updated facility layout. METEC still facilitates for traditional facilities, with a legacy pad, while expanding an new design based on how O&G infrastructure has Final Report - Contract Number: DE-FE0031873 been changed over the last decade. Throughout the ADED program the team has also been working with international partners to ensure staying in the trend globally. International partners have been essential in moving the new protocol forward to implement into CR testing at the METEC facility in Spring 2025.

42 ENGINEERING↗
DOE OSTI2024

Early Inference of Nuclear Technology-Directed Research Activities of Authors from Scientific Publications

Nuclear research articles can provide information about early nuclear proliferation indicators such as influential research entities and technology capability levels of a country, but detection of nuclear activities typically occurs after they have started. We investigate the extent to which nuclear research articles can be used to infer whether a research entity will acquire or develop a nuclear technology before it happens. Early detection of nuclear proliferation or technology development indicators from data is challenging due to partial observability, sparse and unlabeled information, and confounding signals from multiple concurrent activities. This paper presents the early detection problem as a sequential decision-making, goal inference problem, where the objective is to characterize and predict an individual’s, organization’s, or a country’s intent (unobserved goal-directed behavior) towards developing a nuclear capability from partially observed sequences of their research publications, using inverse reinforcement learning and Bayesian goal inference methods. A computational framework is presented, and its application demonstrated using 29,196 Scopus records for a case study related to a civil nuclear capability. The case study results serve as a proof-of-concept demonstration for inference of technology-directed research activity of authors who publish in the nuclear domain. The inference method, combined with advanced computing, may be used to assess and monitor activities pertaining to early developmental stages of a nuclear technology or capability, which in turn can help to identify and prioritize activities with nuclear proliferation potential for further investigation.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

thornado-transport: Anderson- and GPU-accelerated nonlinear solvers for neutrino-matter coupling

Algorithms for neutrino-matter coupling in core-collapse supernovae (CCSNe) are investigated in the context of a spectral two-moment model, which is discretized in space with the discontinuous Galerkin method, integrated in time with implicit-explicit (IMEX) methods, and implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado). The model considers electron neutrinos and antineutrinos and tabulated opacities from Bruenn (1985), which includes neutrino-electron scattering and pair processes. The nonlinear system arising from implicit time discretization of the equations governing neutrino-matter coupling is iterated to convergence using Anderson-accelerated fixed-point methods, which avoid formation of Jacobians and inversion of dense linear systems. Numerical experiments show that, for a given tolerance, a nested iteration scheme which aims to reduce opacity evaluations can lower the computational cost. Our initial port to GPUs, using both OpenMP and OpenACC, shows an overall speedup of up to ~ 100× when compared to results using a single CPU core. These results indicate that the algorithms implemented in thornado are well-suited to GPU acceleration.

Laiu, Paul↗
DOE OSTI2020

Coupled Time-Lapse Full-Waveform Inversion for Subsurface Flow Problems Using Intrusive Automatic Differentiation

We describe a novel framework for estimating subsurface properties, such as rock permeability and porosity, from time-lapse observed seismic data by coupling full-waveform inversion (FWI), subsurface flow processes, and rock physics models. For the inverse modeling, we handle the back propagation of gradients by an intrusive automatic differentiation strategy that offers three levels of user control: (1) At the wave physics level, we adopted the discrete adjoint method in order to use our existing high-performance FWI code; (2) at the rock physics level, we used built-in automatic differentiation operators from the TensorFlow backend; (3) at the flow physics level, we implemented customized partial differential equation (PDE) operators for the multiphase flow equations. The three-level coupled inversion strategy strikes a good balance between computational efficiency and programming efforts, and when the gradients are chained together, it constitutes a coupled inverse system. Our numerical experiments demonstrate that the three-level coupled inverse problem is superior in terms of accuracy to a traditional decoupled inversion strategy. Additionally, our method is able to simultaneously invert for parameters in empirical relationships such as the rock physics models. Our proposed inverted model can be used for reservoir performance prediction and reservoir management/optimization purposes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Two-dimensional plasma density evolution local to the inversion layer during sawtooth crash events using Beam Emission Spectroscopy

Here, we present methods for analyzing Beam Emission Spectroscopy (BES) data to obtain the plasma density evolution associated with rapid sawtooth crash events at the DIII-D tokamak. BES allows coverage over a 2D spatial plane, inherently local measurements, with fast time responses, and, therefore, provides a valuable new channel for data during sawtooth events. A method is developed to remove sawtooth-induced edge-light pulses contained in the BES data. The edge light pulses appear to be from the D α emission produced by edge recycling during sawtooth events, and are large enough that traditional spectroscopic filtering and data analysis techniques are insufficient to deduce physically meaningful quantities. A cross-calibration of 64 BES channels is performed by using a novel method to ensure accurate measurements. For the large-amplitude density oscillations observed, we discuss and use the non-linear relationship between the BES signal δI/ I 0 and the plasma density variation δn e / n e0 . The 2D BES images cover an 8 × 20 cm 2 region around the sawtooth inversion layer and show large-amplitude density oscillations, with additional significant spatial variations across the inversion layer that grows and peaks near the time of the temperature crash. The edge light removal technique and method of converting large-amplitude δI/ I 0 to δn e / n e0 presented here may help analyze other impulsive MHD phenomena in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fast and Scalable FFT-Based GPU-Accelerated Algorithms for Block-Triangular Toeplitz Matrices with Application to Linear Inverse Problems Governed by Autonomous Dynamical Systems

In this work, we present an efficient and scalable algorithm for performing matrix-vector multiplications (matvecs) for block Toeplitz matrices. Such matrices, which are shift-invariant with respect to their blocks, arise in the context of solving inverse problems governed by autonomous systems, and time-invariant systems in particular. In this article, we consider inverse problems that infer unknown parameters from observational data of a linear time-invariant dynamical system given in the form of partial differential equations (PDEs). Matrix-free Newton-conjugate-gradient methods are often the gold standard for solving these inverse problems, but they require numerous actions of the Hessian on a vector. Matrix-free adjoint-based Hessian matvecs require solution of a pair of linearized forward/adjoint PDE solves per Hessian action, which may be prohibitive for large-scale inverse problems. Time invariance of the forward PDE problem leads to a block Toeplitz structure of the discretized parameter-to-observable (p2o) map defining the mapping from inputs (parameters) to outputs (observables) of the PDEs. This block Toeplitz structure enables us to exploit two key properties: (1) compact storage of the p2o map and its adjoint, and (2) efficient fast Fourier transform–based Hessian matvecs. The proposed algorithm is mapped onto large multi-GPU clusters and achieves more than 80% of peak bandwidth on NVIDIA A100 GPUs. Excellent weak scaling is shown for up to 48 A100 GPUs. For the targeted problems, the implementation executes Hessian matvecs within fractions of a second, which is orders of magnitude faster than can be achieved by conventional matrix-free Hessian matvecs via forward/adjoint PDE solves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Physics-informed Karhunen-Loeve and Neural Network Approximations for Solving Inverse Differential Equation Problems

Here we present the PI-CKL-NN method for parameter estimation in differential equation (DE) models given sparse measurements of the parameters and states. In the proposed approach, the space- or time-dependent parameters are approximated by Karhunen-Loeve (KL) expansions that are conditioned on the parameters’ measurements, and the states are approximated by deep neural networks (DNNs). The unknown weights in the KL expansions and DNNs are found my minimizing the cost function that enforces the measurements of the states the DE constraint. Regularization is achieved by adding the l2 norm of the conditional KL coefficients into the loss function. Our approach assumes that the parameter fields are correlated in space or time and enforces the statistical knowledge (the mean and the covariance function) in addition to the DE constraints and measurements as opposed to the physics-informed neural network (PINN) and other similar physics-informed machine learning methods where only DE constraints and data are used for parameter estimation. We use the PI-CKL-NN method for parameter estimation in an ordinary differential equation with an unknown time-dependent parameter and the one- and two-dimensional partial differential diffusion equations with unknown space-dependent diffusion coefficients. We also demonstrate that PI-CKL-NN is more accurate than the PINN method, especially when the observations of the parameters are very sparse

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Estimating Watershed Subsurface Permeability From Stream Discharge Data Using Deep Neural Networks

Subsurface permeability is a key parameter in watershed models that controls the contribution from the subsurface flow to stream flows. Since the permeability is difficult and expensive to measure directly at the spatial extent and resolution required by fully distributed watershed models, estimation through inverse modeling has had a long history in subsurface hydrology. The wide availability of stream surface flow data, compared to groundwater monitoring data, provides a new data source to infer soil and geologic properties using integrated surface and subsurface hydrologic models. As most of the existing methods have shown difficulty in dealing with highly nonlinear inverse problems, we explore the use of deep neural networks for inversion owing to their successes in mapping complex, highly nonlinear relationships. We train various deep neural network (DNN) models with different architectures to predict subsurface permeability from stream discharge hydrograph at the watershed outlet. The training data are obtained from ensemble simulations of hydrographs corresponding to an permeability ensemble using a fully-distributed, integrated surface-subsurface hydrologic model. The trained model is then applied to estimate the permeability of the real watershed using its observed hydrograph at the outlet. Our study demonstrates that the permeabilities of the soil and geologic facies that make significant contributions to the outlet discharge can be more accurately estimated from the discharge data. Their estimations are also more robust with observation errors. Compared to the traditional ensemble smoother method, DNNs show stronger performance in capturing the nonlinear relationship between permeability and stream hydrograph to accurately estimate permeability. Our study sheds new light on the value of the emerging deep learning methods in assisting integrated watershed modeling by improving parameter estimation, which will eventually reduce the uncertainty in predictive watershed models.

54 ENVIRONMENTAL SCIENCES↗