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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Improved Regional Moment Tensor Inversion for Moderately Large Earthquakes in the Western United States Using a 3D Earth Model Based on Full Waveform Tomography

The nature of seismic sources for moderately large (moment magnitude, M w 5.0–6.5) events are commonly characterized by their moment tensor (MT) solutions and obtained by inversion of regional distance (200–1600 km) long‐period (20–50 s) waveforms. Regional MT estimates are often calculated from average plane‐layered, one‐dimensional (1D) velocity models. However, 1D model calculations can produce misfits in the arrival times and waveform shapes that introduce errors, particularly at longer distances or for shorter periods, which are necessary for analyzing lower magnitude events. Approximate Earth models (e.g., 1D) representing broad areas may be inadequate, particularly in the crust and uppermost mantle of tectonically complex regions. In this study, we show how a three‐dimensional (3D) Earth model obtained from full waveform inversion tomography can improve waveform fits and decrease phase errors. We developed a platform and workflow to perform routine 3D MT inversions and inverted MTs for 25 earthquakes in the western United States and seven nuclear explosions using an average 1D and a recent 3D Earth model, WUS256 (Rodgers et al., 2022). Using the 3D model improves waveform fits (variance reduction and phase time shifts) compared with the 1D model, and the 3D MT solutions are stable across large distances. This study shows that 3D models obtained from full waveform tomography can improve MTs and source characterization especially at far regional distances (>800 km).

Geosciences↗
DOE OSTI2022

Computer Model Calibration with Time Series Data Using Deep Learning and Quantile Regression

Computer models play a key role in many scientific and engineering problems. One major source of uncertainty in computer model experiments is input parameter uncertainty. Computer model calibration is a formal statistical procedure to infer input parameters by combining information from model runs and observational data. The existing standard calibration framework suffers from inferential issues when the model output and observational data are high-dimensional dependent data, such as large time series, due to the difficulty in building an emulator and the nonidentifiability between effects from input parameters and data-model discrepancy. To overcome these challenges, we, in this work, propose a new calibration framework based on a deep neural network (DNN) with long short-term memory layers that directly emulates the inverse relationship between the model output and input parameters. Adopting the “learning with noise” idea, we train our DNN model to filter out the effects from data-model discrepancy on input parameter inference. We also formulate a new way to construct interval predictions for DNN using quantile regression to quantify the uncertainty in input parameter estimates. Through a simulation study and real data application with the Weather Research and Forecasting Model Hydrological modeling system (WRF-Hydro), we show our approach can yield accurate point estimates and well-calibrated interval estimates for input parameters.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Magneto-electroluminescence response in 2D and 3D hybrid organic–inorganic perovskite light emitting diodes

We have studied the magneto-electroluminescence (MEL) response in light emitting diodes based on 2D-(PEA) 2 PbI 4 and 3D-MAPbI 3 hybrid organic–inorganic perovskites at cryogenic temperatures. We found that the MEL is negative, i.e., the EL decreases with the applied field strength, B. In addition, the MEL(B) response has a Lorentzian line shape whose width depends on the perovskite used. We interpret the MEL(B) response using the “Δg mechanism” in which the spin of the injected electron–hole (e–h) pairs oscillates between singlet and triplet configurations due to different precession frequencies of the electron and hole constituents that originate from the difference, Δg, in the electron and hole gyromagnetic constants, g. In this model, the MEL(B) linewidth is inversely proportional to the spin lifetime and Δg. The model used is validated by directly measuring the spin lifetime of photogenerated e–h pairs using the circularly polarized pump–probe transmission technique with 100 fs resolution.

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

Permafrost, Peatland, and Cropland Regions Are Key to Reconciling North American Carbon Sink Estimates

Persistent discrepancies between bottom-up, terrestrial biosphere models (TBMs), and top-down, atmospheric inversions, have made it difficult to quantify the magnitude of the North American terrestrial carbon sink. Previous studies have compared aggregated continent-scale estimates of carbon fluxes from TBMs and inversions for all of North America, but this provides limited insights into finer-scale mismatches that contribute to the overall discrepancies. Here we evaluate agreement between TBM and inversion carbon flux estimates at 1° × 1° resolution to provide more direct insights into where models disagree and what underlying factors drive discrepancies. We find that the additional carbon uptake estimated by inversions, in just 16% of the area of North America, is large enough to account for the discrepancy between TBMs and inversions across the whole continent. The majority of these differences occur in permafrost, peatland, and cropland regions. In these regions, we find a higher likelihood of potential biases in the weaker sink estimates from TBMs, suggesting that the stronger sink implied by inversions is more likely to be realistic. However, the current observational coverage is insufficient for fully assessing the causes of discrepancies or the magnitude of biases in either approach. Encouragingly, improved representation of agricultural processes in a TBM led to better agreement with inversions in croplands. Efforts to accurately model cropland dynamics will help improve agreement between TBMs and inversions. Overall, this work presents a clear path for reconciling the discrepancies between inversion and TBM estimates of the North American carbon sink that have persisted for two decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

ML-based Dimension Reduction Strategies

Deep learning (DL)--based surrogate models have achieved success in various applications in carbon capture and storage (CCS). However, the model training on high-dimensional spaces is computationally expensive and impractical for large-scale and complex geological models, because the models usually contain hundreds of thousands to millions of grid cells, each with a set of parameters. Furthermore, the high cost of generating training data with sufficient variation is another limitation of model training on high-dimensional spaces, which may result in overfitting and reduce the model efficiency and prediction performance. We proposed the workflow incorporating dimension reduction methods and deep learning models, which aim to extract the latent variables of input parameters and output state variables, and then build the mapping function at the latent spaces. The proposed workflow can significantly reduce the computational complexity in solving both forward and inverse problems compared to models trained on high-dimensional spaces. Dimensionality reduction models showed great potential in workflows for fast reservoir simulation, history matching, prior model generation, visualization, and more, ultimately enhancing DL model performance in related SMART Work Packages.

Hosseini, Seyyed↗
DOE OSTI2022

A combined experimental and numerical approach that eliminates the non-uniqueness associated with the Johnson-Cook parameters obtained using inverse methods

Abstract Johnson-Cook constitutive model is a commonly used material model for machining simulations. The model includes five parameters that capture the initial yield stress, strain-hardening, strain-rate hardening, and thermal softening behavior of the material. These parameters are difficult to determine using experiments since the conditions observed during machining (such as high strain-rates of the order of $$10^5$$ 10 5 /sec - $$10^6$$ 10 6 /sec) are challenging to recreate in the laboratory. To address this problem, several researchers have recently proposed inverse approaches where a combination of experiments and analytical models are used to predict the Johnson-Cook parameters. The errors between the measured cutting forces, chip thicknesses and temperatures and those predicted by analytical models are minimized and the parameters are determined. In this work, it is shown that only two of the five Johnson-Cook parameters can be determined uniquely using inverse approaches. Two different algorithms, namely, Adaptive Memory Programming for Global Optimization (AMPGO) and Particle Swarm Optimization (PSO), are used for this purpose. The extended Oxley’s model is used as the analytical tool for optimization. For determining a parameter’s value, a large range for each parameter is provided as an input to the algorithms. The algorithms converge to several different sets of values for the five Johnson-Cook parameters when all the five parameters are considered as unknown in the optimization algorithm. All of these sets, however, yield the same chip shape and cutting forces in FEM simulations. Further analyses show that only the strain-rate and thermal softening parameters can be determined uniquely and the three parameters present in the strain-hardening term of the Johnson-Cook model cannot be determined uniquely using the inverse method. A combined experimental and numerical approach is proposed to eliminate this determine all parameters uniquely.

42 ENGINEERING↗
DOE OSTI2025

Quantifying uncertainty in regional-scale seismic moment tensors

We examine the ability of three different inversion methods: 1) first motion (FM) inversions, 2) amplitude inversions, and 3) full waveform (FW) time variable moment tensor (TVMT) inversions, to recover an accurate source mechanism for simulated data for an earthquake as well as an explosion. The ability of inversion methods, such as the ones described above, to recover accurate models representative of the data depends on both a priori information as well as the quality of data being inverted. Therefore, we examine the effect that station geometry, geologic model, and noise level has on the inversion and estimated source mechanism. We find that FM data can provide more robust solutions than amplitude data and are not as sensitive to inaccurate earth models, especially in low-noise cases, and overall FW inversions are the most accurate out of the methods examined, but can still be biased by inaccurate earth models.

58 GEOSCIENCES↗
DOE OSTI2024

Bayesian calibration and uncertainty quantification of a rate-dependent cohesive zone model for polymer interfaces

In this work we present a rate-dependent cohesive zone model for the fracture of polymeric interfaces and performs a Bayesian calibration, an uncertainty quantification, and a sensitivity analysis for the model. The proposed cohesive zone model accounts for both reversible elastic and irreversible rate-dependent separation sliding deformation at the interface. The viscous dissipation due to the irreversible opening at the interface is modeled using elastic-viscoplastic kinematics that incorporates the effects of strain rate. Inverse calibration of parameters for such complex models through trial and error is challenging due to the large number of parameters of the model. Moreover, the calibrated parameter values are often non-unique and uncertain when the available experimental data is limited. To tackle this challenge, we employ a Bayesian calibration approach to identify parameters from experimental data, the resulting parameters significantly enhance the accuracy of the model. To quantify the uncertainty associated with the inverse parameter estimation, a modular Bayesian approach is employed to calibrate the unknown model parameters, accounting for the parameter uncertainty of the cohesive zone model. The advantages of the Bayesian calibration over a deterministic parameter fit are demonstrated. Further, to quantify the model uncertainties, such as incorrect assumptions or missing physics, a discrepancy function is introduced, which significantly improves the model’s prediction. Finally, the total uncertainty of the model is quantified in a predictive setting. A sensitivity analysis is performed to assess how changes in the input variables of the model affect the peak load, facilitating the identification of a concise set of highly influential parameters. The present approach can be used for calibration and uncertainty quantification for other complex computational mechanics models. It should also facilitate the designing of interface materials under uncertainty.

42 ENGINEERING↗
DOE OSTI2023

Compensating attenuation effects in full-waveform inversion with dissipation-dispersion decoupling

Seismic attenuation poses challenges to velocity model building from acoustic/elastic full-waveform inversion (FWI). In particular, when constructing the FWI gradient for velocity inversion, recent studies have indicated that the high-attenuation structure could distort the gradient by damping the amplitudes and shifting the kinematic phases, resulting in an imbalanced update and thus unreliable velocity model. These Q effects are particularly significant in the reflection acquisition geometry due to the “double-damping” issue. Here, we develop a Q-compensated FWI algorithm for constructing a Q-free FWI gradient. By using a recently developed viscoacoustic wave propagator, this compensation can be done conveniently by keeping the dispersion term and flipping the dissipation term in the wave equation when we simulate the forward and time-reversed adjoint wavefields. The resultant gradient obtained by interacting these two wavefields has correct kinematics and Q-free amplitudes. This Q-compensated FWI can balance the update between Q- and no- Q-areas, which we determine using synthetic examples. In addition, we illustrate how to take advantage of the dissipation-dispersion decoupling to determine the anomalous Q value in the Q model building workflow via Q-compensated reverse time migration.

Geochemistry & Geophysics↗
DOE OSTI2022

Proxy-Based Bayesian Inversion Of Poroelastic Simulations To Interpret Strain Tensor Data Measured During Well Testing

The long runtimes of 3D poroelastic numerical simulators makes it impractical to interpret deformation datasets using many inversion schemes. Recent advances in instrumentation have made it possible to measure the strain tensor during well testing, but the lack of robust inversion methods is limiting the ability to interpret these data. We have developed an inversion workflow that reduces the number of computations required to complete a Bayesian inversion using DREAMzs. The workflow trains a KNN model using output from the poroelastic simulator, and then uses the KNN model as a proxy for the simulator during inversion. The workflow also includes a strategy for ensuring the results from the proxy model converge to the results from the simulator, ensuring the accuracy of the final results. An idealized example configured to represent a well test in a deep aquifer is used to verify that the workflow correctly identifies parameters and characterizes noise. Field data measured using strainmeters during an injection test at an oil reservoir in Oklahoma are used to evaluate performance with a real dataset. The workflow identified 265 history matching solutions out of 1240 total simulation runs (21% acceptance ratio), and the results are used to characterize posterior parameter distribution and evaluate the prediction uncertainty. This approach makes it feasible to invert strain data measured during well testing and this has the potential to improve the characterization of aquifers and reservoirs.

Roudini, Soheil↗
DOE OSTI2021

Spatial distributions of X CO 2 seasonal cycle amplitude and phase over northern high-latitude regions

Satellite-based observations of atmospheric carbon dioxide (CO 2 ) provide measurements in remote regions, such as the biologically sensitive but undersampled northern high latitudes, and are progressing toward true global data coverage. Recent improvements in satellite retrievals of total column-averaged dry air mole fractions of CO 2 (X CO 2 ) from the NASA Orbiting Carbon Observatory 2 (OCO-2) have allowed for unprecedented data coverage of northern high-latitude regions, while maintaining acceptable accuracy and consistency relative to ground-based observations, and finally providing sufficient data in spring and autumn for analysis of satellite-observed X CO 2 seasonal cycles across a majority of terrestrial northern high-latitude regions. Here, we present an analysis of X CO 2 seasonal cycles calculated from OCO-2 data for temperate, boreal, and tundra regions, subdivided into 5° latitude by 20° longitude zones. We quantify the seasonal cycle amplitudes (SCAs) and the annual half drawdown day (HDD). OCO-2 SCAs are in good agreement with ground-based observations at five high-latitude sites, and OCO-2 SCAs show very close agreement with SCAs calculated for model estimates of X CO 2 from the Copernicus Atmosphere Monitoring Services (CAMS) global inversion-optimized greenhouse gas flux model v19r1 and the CarbonTracker2019 model (CT2019B). Model estimates of X CO 2 from the GEOS-Chem CO 2 simulation version 12.7.2 with underlying biospheric fluxes from CarbonTracker2019 (GC-CT2019) yield SCAs of larger magnitude and spread over a larger range than those from CAMS, CT2019B, or OCO-2; however, GC-CT2019 SCAs still exhibit a very similar spatial distribution across northern high-latitude regions to that from CAMS, CT2019B, and OCO-2. Zones in the Asian boreal forest were found to have exceptionally large SCA and early HDD, and both OCO-2 data and model estimates yield a distinct longitudinal gradient of increasing SCA from west to east across the Eurasian continent. In northern high-latitude regions, spanning latitudes from 47 to 72° N, longitudinal gradients in both SCA and HDD are at least as pronounced as latitudinal gradients, suggesting a role for global atmospheric transport patterns in defining spatial distributions of X CO2 seasonality across these regions. GEOS-Chem surface contact tracers show that the largest X CO 2 SCAs occur in areas with the greatest contact with land surfaces, integrated over 15–30d. The correlation of X CO 2 SCA with these land surface contact tracers is stronger than the correlation of X CO 2 SCA with the SCA of CO 2 fluxes or the total annual CO 2 flux within each 5° latitude by 20° longitude zone. This indicates that accumulation of terrestrial CO 2 flux during atmospheric transport is a major driver of regional variations in X CO 2 SCA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

BayesMT: A Probabilistic Bayesian Framework for the Seismic Moment Tensor

Moment tensors (MTs) have long been used in earthquake and explosion source analysis, and there has been a renewed interest in how they can inform us about the seismic source, particularly in the geophysical monitoring community due to its application in event identification and yield analysis. However, parameter uncertainties in seismic MT inversion are rarely available. The inverse procedure often does not quantify MT model errors such as event location, data noise and Earth model that are essential for estimating solution robustness. To address this need, we propose to adopt the Bayesian probabilistic framework to incorporate uncertainties in MT inversions. In this study, we present the theoretical background of a probabilistic Bayesian framework for MT inversion accounting for model and measurements errors and illustrate the implementation of the method using a synthetic example.

58 GEOSCIENCES↗
DOE OSTI2024

The key role of forest disturbance in reconciling estimates of the northern carbon sink

Northern forests are an important carbon sink, but our understanding of the driving factors is limited due to discrepancies between dynamic global vegetation models (DGVMs) and atmospheric inversions. We show that DGVMs simulate a 50% lower sink (1.1 ± 0.5 PgC yr –1 over 2001–2021) across North America, Europe, Russia, and China compared to atmospheric inversions (2.2 ± 0.6 PgC yr –1 ). We explain why DGVMs underestimate the carbon sink by considering how they represent disturbance processes, specifically the overestimation of fire emissions, and the lack of robust forest demography resulting in lower forest regrowth rates than observed. We reconcile net sink estimates by using alternative disturbance-related fluxes. We estimate carbon uptake through forest regrowth by combining satellite-derived forest age and biomass maps. We calculate a regrowth flux of 1.1 ± 0.1 PgC yr –1 , and combine this with satellite-derived estimates of fire emissions (0.4 ± 0.1 PgC yr –1 ), land-use change emissions from bookkeeping models (0.9 ± 0.2 PgC yr –1 ), and the DGVM-estimated sink from CO 2 fertilisation, nitrogen deposition, and climate change (2.2 ± 0.9 PgC yr –1 ). The resulting ‘bottom-up’ net flux of 2.1 ± 0.9 PgC yr –1 agrees with atmospheric inversions. The reconciliation holds at regional scales, increasing confidence in our results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Manifold Learning-Based Polynomial Chaos Expansions for High-Dimensional Surrogate Models

In this work we introduce a manifold learning-based method for uncertainty quantification (UQ) in systems describing complex spatiotemporal processes. Our first objective is to identify the embedding of a set of high-dimensional data representing quantities of interest of the computational or analytical model. For this purpose, we employ Grassmannian diffusion maps, a two-step nonlinear dimension reduction technique which allows us to reduce the dimensionality of the data and identify meaningful geometric descriptions in a parsimonious and inexpensive manner. Polynomial chaos expansion is then used to construct a mapping between the stochastic input parameters and the diffusion coordinates of the reduced space. An adaptive clustering technique is proposed to identify an optimal number of clusters of points in the latent space. The similarity of points allows us to construct a number of geometric harmonic emulators which are finally utilized as a set of inexpensive pretrained models to perform an inverse map of realizations of latent features to the ambient space and thus perform accurate out-of-sample predictions. Thus, the proposed method acts as an encoder-decoder system which is able to automatically handle very high-dimensional data while simultaneously operating successfully in the small-data regime. The method is demonstrated on two benchmark problems and on a system of advection-diffusion-reaction equations which model a first-order chemical reaction between two species. In all test cases, the proposed method is able to achieve highly accurate approximations which ultimately lead to the significant acceleration of UQ tasks.

42 ENGINEERING↗
DOE OSTI2021

Enhanced mechanical property evaluation using innovative data analytics capability

This report focuses on efforts to improve Multiphysics Object-Oriented Simulation Environ-ment (MOOSE) for mechanical property evaluation using data analytics. These efforts involveimprovements to the stochastic tools module (STM) for stochastic simulations of MOOSE mul-tiphysics model and the development of inverse optimization capabilities. The report gives anoverview of the STM and describes recent updates to its core capabilities and theory on itsreduced-order model (ROM) schemes. Examples are also provided showing the impact of theseupdates and exhibits the usefulness of ROMs. An overview of the gradient based inverse opti-mziation algorithms are given along with examples of their application to source identification.Inverse optimization will provide a new methodology in STM for fitting model parameters toexperimental data.

97 MATHEMATICS AND COMPUTING↗