Search NASA⌕ Search

SEARCH · Search NASA

Results for “inverse analysis”

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

Combining translational and rotational seismic motions to invert local-scale seismic data for time-variable moment tensors: do rotational motions help for high-frequency seismic data produced by underground explosions?

SUMMARY We present an analysis of combining translational and rotational seismic data in an inversion for the time-variable source time functions corresponding to the components of the seismic moment tensor. We conduct a series of numerical experiments where the data are simulated by a combination of an underground explosion and a co-located double couple shear source and recorded on surface-mounted seismometers within 1–2 km of the source. The experiments are designed to mimic explosion seismology experiments, and thus the data are in the 1–10 Hz frequency range and contain very few surface waves. We use a Monte Carlo method to propagate Earth model uncertainty into the estimates of seismic source parameters. In our experiments, we find that the uncertainty of the estimated seismic source parameters increases when we add rotational seismic motions to the inversion when using a constant number of data channels. In this case, the increased degree of uncertainty in the final results is most likely due to the near-surface Earth model uncertainty that we introduce in our simulations. However, for a fixed number of seismic stations, adding rotational seismic motions to the inversion acts to decrease the uncertainty of the estimated seismic source parameters, most likely due to the increase in the number of data channels used in the inversion.

Poppeliers, Christian (ORCID:0000000159526849)↗
DOE OSTI2020

A Generalized Stress Inversion Approach With Application to Residual Stress Estimation

In this work develop a generalized stress inversion technique (or the generalized inversion method) capable of recovering stresses in linear elastic bodies subjected to arbitrary cuts. Specifically, given a set of displacement measurements found experimentally from digital image correlation (DIC), we formulate a stress estimation inverse problem as a partial differential equation-constrained optimization problem. We use gradient-based optimization methods, and we accordingly derive the necessary gradient and Hessian information in a matrix-free form to allow for parallel, large-scale operations. By using a combination of finite elements, DIC, and a matrix-free optimization framework, the generalized inversion method can be used on any arbitrary geometry, provided that the DIC camera can view a sufficient part of the surface. We present numerical simulations and experiments, and we demonstrate that the generalized inversion method can be applied to estimate residual stress.

42 ENGINEERING↗
DOE OSTI2025

The weak land carbon sink hypothesis

Over the past three decades, assessments of the contemporary global carbon budget consistently report a strong net land carbon sink. Here, we review evidence supporting this paradigm and quantify the differences in global and Northern Hemisphere estimates of the net land sink derived from atmospheric inversion and satellite-derived vegetation biomass time series. Our analysis, combined with additional synthesis, supports a hypothesis that the net land sink is substantially weaker than commonly reported. At a global scale, our estimate of the net land carbon sink is 0.8 ± 0.7 petagrams of carbon per year from 2000 through 2019, nearly a factor of two lower than the Global Carbon Project estimate. With concurrent adjustments to ocean (+8%) and fossil fuel (−6%) fluxes, we develop a budget that partially reconciles key constraints provided by vegetation carbon, the north-south CO 2 gradient, and O 2 trends. We further outline potential modifications to models to improve agreement with a weaker land sink and describe several approaches for testing the hypothesis.

Randerson, James T. [Univ. of California, Irvine, ↗
DOE OSTI2023

Light Elements $R$-matrix Analyses with the SAMMY code towards the Foundation of Charged-particle Nuclear Data Libraries [Slides]

This presentation covers newly developed SAMMY module for inverse channel transformation. Additionally covered is the R-matrix analysis of 7 Be compound nucleus and the R-­matrix analysis of 17 O compound nucleus. Further touched on is the evaluated Nuclear Data File generation and processing with the AMPX code. The presentation concludes with talks on future evaluation work and tests on light nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Variational encoder geostatistical analysis (VEGAS) with an application to large scale riverine bathymetry

Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2025

Quasinormal Coupled-Mode Analysis of Dynamic Gain in Exceptional-Point Lasers

One of the key features of lasers operating near exceptional points (EPs) is that the gain medium can support an oscillating population inversion above a pump threshold, leading to self-modulated laser dynamics. This unusual behavior opens up new possibilities for frequency comb generation and temporal modulation. However, the dynamic population inversion couples signals with different frequencies and is difficult to capture using conventional temporal coupled-mode theory (TCMT) based on stationary saturable gain. In this paper, we develop a perturbative coupled-mode analysis framework to capture the spatial-temporal dynamics of near-EP lasers. By decomposing discrete frequency generation into multiple excitations of resonant modes, our analysis establishes a minimal physical model that translates the local distribution of dynamic population-inversion into a resonant modal interpretation of laser gain. Furthermore, this work enables the exploration of unique properties in this self-time-modulated systems, such as time-varying scattering and nonreciprocal transmission.

Absorption↗
DOE OSTI2025

Machine learning inversion from small-angle scattering for charged polymers

We develop Monte Carlo simulations for uniformly charged polymers and a machine learning algorithm to interpret the intra-polymer structure factor of the charged polymer system, which can be obtained from small-angle scattering experiments. The polymer is modeled as a chain of fixed-length bonds, where the connected bonds are subject to bending energy, and there is also a screened Coulomb potential for charge interaction between all joints. The bending energy is determined by the intrinsic bending stiffness, and the charge interaction depends on the interaction strength and screening length. All three contribute to the stiffness of the polymer chain and lead to longer and larger polymer conformations. The screening length also introduces a second length scale for the polymer besides the bending persistence length. To obtain the inverse mapping from the structure factor to these polymer conformation and energy-related parameters, we generate a large data set of structure factors by running simulations for a wide range of polymer energy parameters. We use principal component analysis to investigate the intra-polymer structure factors and determine the feasibility of the inversion using the nearest neighbor distance. We employ Gaussian process regression to achieve the inverse mapping and extract the characteristic parameters of polymers from the structure factor with low relative error.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Measurement of mix at the fuel–ablator interface in indirectly driven capsule implosions on the National Ignition Facility

The interface between the capsule ablator and fuel ice layer is susceptible to hydrodynamic instabilities. The subsequent mixing of hot ablator material into the ice reduces fuel compression at stagnation and is a candidate for reduced capsule performance. The ability to diagnose ice–ablator mix is critical to understanding and improving stability at this interface. Combining the crystal backlighter imager with the single line of sight camera on the National Ignition Facility (NIF) allows direct measurement of ice–ablator mix by providing multiple quasi-monochromatic radiographs of layered capsule implosions per experiment with high spatial (∼12 μm) and temporal (∼35 ps) resolution. The narrow bandwidth of this diagnostic platform allows radiography of the inner edge of the capsule limb close to stagnation without capsule self-emission contaminating the data and removes opacity uncertainties typically associated with the spectral content of the radiograph. Analysis of radiographic data via a parameterized forward-fitting Abel inversion technique provides measurements of the distribution of mix mass inwards from the ice–ablator interface. The sensitivity of this mix measurement technique was demonstrated by applying it to layered experiments in which the stability of the ice–ablator interface was expected to vary significantly. Additional experiments suggest that high-density carbon capsules that employ a buried-layer dopant profile suffer from mixing at the innermost doped–undoped interface. Data from these experiments suggest that opacity models used in hydrodynamic simulations of NIF experiments can potentially over-predict the opacity of doped capsules. LLNL-JRNL-850535-DRAFT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Can a transport model predict inverse signatures in Lithium metal batteries without modifying kinetics?

In this study, a one-dimensional transport model is developed and analyzed that predicts the inverse signatures obtained in lithium metal batteries. This simple approach’s ability to predict inverse signatures stems from the competing interplay between moving boundary rates and mass transfer limitations to yield inverse signatures without modifying the kinetics. The numerical scheme used for the present model simulations is presented in detail which has been further used to study the effect of design parameters on the prevalence and strength of inverse signatures. It was found that the proposed model and the analysis is more pertinent to thick symmetric cells, commonly used for in-depth fundamental studies.

model, transport, dendrite, lithium metal battery,↗
DOE OSTI2021

Overview of Hydraulic Fracturing Test Site 2 in the Permian Delaware Basin (HFTS-2)

Here, the Hydraulic Fracturing Test Site 2 (HFTS-2) is a large collaborative field-based R&D program in the Permian Delaware Basin, funded by the US Department of Energy through the National Energy Technology Laboratory (NETL) and the E&P industry, with support from academia. The projects' main objective is to improve the understating of the hydraulic fracturing process through utilization of advanced diagnostics and collection of through-fracture cores to provide undisputable evidence and attributes of the created hydraulic fractures. At the HFTS-2, in excess of $30 million was used to perform hydraulic fracturing research focusing on the Wolfcamp formation at a field site hosted and operated by Occidental. In addition to the research data collected by the project, Occidental provided a significant amount of background data for about a dozen existing wells in the test area as well as access to previously collected core. Additional technical and laboratory support was provided by the program members. Building on learnings and unanswered questions from HFTS-1 in the Permian Midland basin, the HFTS-2 used eight new producing wells and two existing (parent) wells to perform hydraulic fracturing research. Multiple science wells were drilled to sample and characterize the subsurface, including the collection of 540 feet of core in a vertical pilot hole and 948 feet of high-angle through-fracture core. The project installed permanent fiber optic cables in 3 wells to monitor near wellbore signals during fracturing and to collect cross-well strain measurements. Additional advanced diagnostics included a significant formation evaluation program on the vertical whole core, multi array moment tensor inversion capable microseismic survey, multi-well time-lapse geochemistry analysis, analysis of proppant distribution in producing child and slant core well, and others. We will provide an overview of the HFTS-2 project, including list of the consortium members, details of the test site, experiments performed, and technologies tested.

58 GEOSCIENCES↗
DOE OSTI2022

Accurate and Rapid Forecasts for Geologic Carbon Storage via Learning-Based Inversion-Free Prediction

Carbon capture and storage (CCS) is one approach being studied by the U.S. Department of Energy to help mitigate global warming. The process involves capturing CO 2 emissions from industrial sources and permanently storing them in deep geologic formations (storage reservoirs). However, CCS projects generally target “green field sites,” where there is often little characterization data and therefore large uncertainty about the petrophysical properties and other geologic attributes of the storage reservoir. Consequently, ensemble-based approaches are often used to forecast multiple realizations prior to CO 2 injection to visualize a range of potential outcomes. In addition, monitoring data during injection operations are used to update the pre-injection forecasts and thereby improve agreement between forecasted and observed behavior. Thus, a system for generating accurate, timely forecasts of pressure buildup and CO 2 movement and distribution within the storage reservoir and for updating those forecasts via monitoring measurements becomes crucial. This study proposes a learning-based prediction method that can accurately and rapidly forecast spatial distribution of CO 2 concentration and pressure with uncertainty quantification without relying on traditional inverse modeling. The machine learning techniques include dimension reduction, multivariate data analysis, and Bayesian learning. The outcome is expected to provide CO 2 storage site operators with an effective tool for timely and informative decision making based on limited simulation and monitoring data.

58 GEOSCIENCES↗
DOE OSTI2022

A Statistical Interpolation Code for Ocean Analysis and Forecasting

Abstract We present a data assimilation package for use with ocean circulation models in analysis, forecasting, and system evaluation applications. The basic functionality of the package is centered on a multivariate linear statistical estimation for a given predicted/background ocean state, observations, and error statistics. Novel features of the package include support for multiple covariance models, and the solution of the least squares normal equations either using the covariance matrix or its inverse—the information matrix. The main focus of this paper, however, is on the solution of the analysis equations using the information matrix, which offers several advantages for solving large problems efficiently. Details of the parameterization of the inverse covariance using Markov random fields are provided and its relationship to finite-difference discretizations of diffusion equations are pointed out. The package can assimilate a variety of observation types from both remote sensing and in situ platforms. The performance of the data assimilation methodology implemented in the package is demonstrated with a yearlong global ocean hindcast with a 1/4° ocean model. The code is implemented in modern Fortran, supports distributed memory, shared memory, multicore architectures, and uses climate and forecasts compliant Network Common Data Form for input/output. The package is freely available with an open source license from www.tendral.com/tsis/ .

Srinivasan, Ashwanth↗
DOE OSTI2021

LL18-Signal Analysis-NDD2Ad (Final Report)

Under favorable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially-distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection, and scaling with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning, and possible application to studies of repeating events in induced, volcanic, and glacial seismicity. Our results place limits on the widely-held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗
DOE OSTI2025

Dynamical phase transition in a strongly hybridized phonon-triplon chain

In this article, we study a dimerized spin-1/2 chain, such as CuGeO 3 , hosting triplon excitations coupled to optical phonons under weak terahertz laser driving. Both phonons and triplons weakly lose energy into the surrounding baths, forming a non-equilibrium steady state. In the strong phonon-triplon coupling regime, phonons near the two-triplon continuum hybridize strongly with triplons. Using mean-field Lindblad dynamics, we show that strong hybridization induces sharp first-order phase transitions—either single or simultaneous double—in the emission spectrum, mainly due to nonlinearities in the chain and the instantaneous breaking of spatial inversion symmetry, as captured by the phonon displacement. Using mean-field Floquet analysis of harmonic modes in both sectors, we analytically confirm the existence of these phase transitions. Furthermore, we map the complete steady-state phase diagram by varying key control parameters and provide experimentally relevant parameters for observing these transitions in laser-driven CuGeO 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Source-Resolved Inversion of Elemental Carbon Emissions in California Using Log-Space Bayesian Inference

Elemental carbon (EC), operationally quantified by thermal-optical analysis, is widely used as a proxy for black carbon (BC) relevant to short-term climate forcing and public health. Current EC emission inventories remain highly uncertain, with persistent discrepancies between bottom-up and top-down estimates. In this study, we develop a source-resolved, log-space Bayesian inversion framework applied to estimate California’s statewide EC emissions in 2019. By integrating surface EC measurements from the EPA’s Air Quality System network with high-resolution source contributions simulated by a chemical transport model, we identify a one-third underestimation in the existing statewide EC inventory, requiring an increase of the total from a prior of 8.58 [5.49–13.75] Gg year–1 to a posterior estimate of 12.78 [10.71–15.37] Gg year–1. This discrepancy is primarily driven by substantial underestimations in the power and industrial and off-road mobile sectors. Furthermore, population-weighted exposure analysis reveals a marked sectoral divergence between emission mass and health burden: off-road mobile sources dominate both emissions and exposure, accounting for 31% of statewide exposure, while residential wood combustion contributes 26% of total exposure despite comprising only 19% of total emissions, due to its source proximity to population. These findings underscore the need to update sector-specific EC speciation profiles and demonstrate that mitigation strategies targeting off-road mobile sources and residential wood combustion are critical for reducing EC-related health impacts in California.

Zhang, Jie↗
DOE OSTI2021

Tuning valley degeneracy with band inversion

Valley degeneracy is a key feature of the electronic structure that benefits the thermoelectric performance of a material. Despite recent studies which claim that high valley degeneracy can be achieved with inverted bands, our analysis of rock-salt IV–VI compounds using first-principles calculations and k•p perturbation theory demonstrates that mere band inversion is an insufficient condition for high valley degeneracy; rather, there is a critical degree to which the bands must be inverted to induce multiple carrier pockets. The so-called “band inversion parameter” is formalized as a chemically-tunable property, offering a design route to achieving high valley degeneracy in compounds with inverted bands. In conclusion, we predict that the valley degeneracy of rock-salt IV–VI compounds can be increased from N V = 4 to N V = 24, which could result in a corresponding increase in the thermoelectric figure of merit zT .

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

Surface structures of magnetostrictive D0 3 -Fe 3 Ga(0 0 1)

First-principles total energy calculations and scanning tunneling microscopy experiments were performed to study the surface reconstruction of the magnetostrictive Fe 3 Ga alloy. The inverse magnetostrictive behavior was evaluated in the bulk by compressing and stretching its lattice parameter, showing an increase in magnetic moments as strain increases. Surface analysis demonstrates two thermodynamically stable surfaces, the (1 × 1) and (3 × 1). The (1 × 1) is an ideal FeGa terminated surface, whereas the (3 × 1) is also FeGa terminated but it has a first-layer Fe atom substituted by a Ga atom every three unit-cells, forming a row-like surface structure. Tersoff–Hamann scanning tunneling microscopy simulations were obtained and compared with experimental results. We found good agreement between theory and experiment, in which the distance between rows is ~12.3 Å. Theoretical findings suggest that the substrate-induced strain may increase the stability of the (3 × 1) reconstruction. Analysis of the magnetic moments in the reconstructions showed that their behavior is affected by a surface effect, as well as by the inverse magnetostriction of the structure. A good understanding of the FeGa surface reconstructions is an important step towards further improvements in magnetic storage devices and sensors.

36 MATERIALS SCIENCE↗