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 91 records · Page 5

A full year of aerosol size distribution data from the central Arctic under an extreme positive Arctic Oscillation: insights from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition

Abstract. The Arctic environment is rapidly changing due to accelerated warming in the region. The warming trend is driving a decline in sea ice extent, which thereby enhances feedback loops in the surface energy budget in the Arctic. Arctic aerosols play an important role in the radiative balance and hence the climate response in the region, yet direct observations of aerosols over the Arctic Ocean are limited. In this study, we investigate the annual cycle in the aerosol particle number size distribution (PNSD), particle number concentration (PNC), and black carbon (BC) mass concentration in the central Arctic during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. This is the first continuous, year-long data set of aerosol PNSD ever collected over the sea ice in the central Arctic Ocean. We use a k-means cluster analysis, FLEXPART simulations, and inverse modeling to evaluate seasonal patterns and the influence of different source regions on the Arctic aerosol population. Furthermore, we compare the aerosol observations to land-based sites across the Arctic, using both long-term measurements and observations during the year of the MOSAiC expedition (2019–2020), to investigate interannual variability and to give context to the aerosol characteristics from within the central Arctic. Our analysis identifies that, overall, the central Arctic exhibits typical seasonal patterns of aerosols, including anthropogenic influence from Arctic haze in winter and secondary aerosol processes in summer. The seasonal pattern corresponds to the global radiation, surface air temperature, and timing of sea ice melting/freezing, which drive changes in transport patterns and secondary aerosol processes. In winter, the Norilsk region in Russia/Siberia was the dominant source of Arctic haze signals in the PNSD and BC observations, which contributed to higher accumulation-mode PNC and BC mass concentrations in the central Arctic than at land-based observatories. We also show that the wintertime Arctic Oscillation (AO) phenomenon, which was reported to achieve a record-breaking positive phase during January–March 2020, explains the unusual timing and magnitude of Arctic haze across the Arctic region compared to longer-term observations. In summer, the aerosol PNCs of the nucleation and Aitken modes are enhanced; however, concentrations were notably lower in the central Arctic over the ice pack than at land-based sites further south. The analysis presented herein provides a current snapshot of Arctic aerosol processes in an environment that is characterized by rapid changes, which will be crucial for improving climate model predictions, understanding linkages between different environmental processes, and investigating the impacts of climate change in future Arctic aerosol studies.

54 ENVIRONMENTAL SCIENCES↗

Nonlinear spectral analysis of ion acoustic solitons arising from a streaming charged object using the numerical inverse scattering transform

Observations of nonlinear coherent plasma phenomena in spacecraft and terrestrial experiments often rely on visual identification of solitary modes because nonlinear coherent modes are broadband in a Fourier transform based analysis. We implement an alternative spectral decomposition known as the inverse scattering transform and demonstrate its ability to successfully isolate nonlinear modes on the homogeneous and forced Korteweg–De Vries equation which models these nonlinear modes. We also demonstrate for the first time that this decomposition is useful when forcing is applied. This is because the stable modes generated by localized forcing are similar to the homogeneous solutions where the inverse scattering transform is a rigorous decomposition. This spectral technique is then applied to simulations of ion acoustic waves generated by a 10 mm spherical debris object interacting with the ionospheric plasma orbiting at 2000 km altitude. The algorithm is found to successfully detect the resultant solitons. This demonstrates the feasibility of using this spectral technique as a real time analysis tool for screening spacecraft data for nonlinear solitary modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aquifer Monitoring Using Ambient Seismic Noise Recorded With Distributed Acoustic Sensing (DAS) Deployed on Dark Fiber

Abstract Groundwater is a critical resource for human activities worldwide, and a vital component of many natural ecosystems. However, the state and dynamics of water‐bearing aquifers remain uncertain, mostly due to the paucity of subsurface data at high spatial and temporal resolution. Here, we show that analysis of infrastructure‐generated ambient seismic noise acquired on distributed acoustic sensing (DAS) arrays has potential as a tool to track variations in seismic velocities (d v / v ) caused by groundwater level fluctuations. We analyze 5 months of ambient noise acquired along an unused, 23 km‐long telecommunication fiber‐optic cable in the Sacramento Valley, CA, a so‐called “dark fiber." Three array subsections, ∼6 km apart, are processed and the stretching technique is applied to retrieve daily d v / v beneath each location. Near the Sacramento river, d v / v variations in the order of 2%–3% correlate with precipitation events and fluctuations in river stage of ∼1.5 m. In contrast, regions away (2.5 km) from the river do not experience large d v / v variations. These observations reveal short‐scale spatial variability in aquifer dynamics captured by this approach. Dispersion analysis and surface wave inversion of noise gathers reveal that seismic velocity perturbations occur at depths of 10–30 m. Rock physics modeling confirms that observed d v / v are linked to pore pressure changes at these depths, caused by groundwater table fluctuations. Our results suggest that DAS combined with ambient noise interferometry provides a means of tracking aquifer dynamics at high spatial and temporal resolutions at local to regional scales, relevant for effective groundwater resource management.

58 GEOSCIENCES↗

Advances in soft X-ray RIXS for studying redox reaction states in batteries

Redox (reduction and oxidation) chemistry provides the fundamental basis for numerous energy-related electrochemical devices. Detecting the electrochemical redox chemistry is pivotal but challenging because it requires independent probes of the cationic and anionic redox states at different electrochemical states. Here, the synchrotron-based soft X-ray mapping of resonant inelastic X-ray scattering (mRIXS) has recently emerged as a powerful tool for exploring such states in electrochemical devices, especially batteries. High-efficiency mRIXS covers the energy range of the absorption edge with the extra dimension of information on the emitted photon energies. In this frontier article, we review recent representative demonstrations of utilizing soft X-ray mRIXS for detecting the novel chemical state during electrochemical operation and for quantifying the cationic redox reactions through inverse partial fluorescence yield analysis (mRIXS-iPFY). More importantly, the non-divalent states of oxygen in electrodes involving oxygen redox reactions could be reliably captured by mRIXS, with its reversibility quantified by the intensity variation of the characteristic mRIXS feature through a super-partial fluorescence yield analysis (mRIXS-sPFY). These recent demonstrations inspire future perspectives on using mRIXS for studying the complex phenomena in energy materials, with both technical and scientific challenges in RIXS theory, in situ/operando experiments, and spatially resolved RIXS imaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Emulation of the calculations of final r -process abundance patterns with a neural network

This work explores the construction of a fast emulator for the calculation of the final pattern of nucleosynthesis in the rapid neutron capture process (the r-process). An emulator is built using a feed-forward artificial neural network (ANN). We train the ANN with nuclear data and relative abundance patterns. We take as input the β-decay half-lives and the one-neutron separation energy of the nuclei in the rare-earth region. The output is the final isotopic abundance pattern. In this work, we focus on the nuclear data and abundance patterns in the rare-earth region to reduce the dimension of the input and output space. We show that the ANN can capture the effect of the changes in the nuclear physics inputs on the final r-process abundance pattern in the adopted astrophysical conditions. We employ the deep ensemble method to quantify the prediction uncertainty of the neural network emulator. The emulator achieves a speed-up by a factor of about 20 000 in obtaining a final abundance pattern in the rare-earth region. The emulator may be utilized in statistical analyses such as uncertainty quantification, inverse problems, and sensitivity analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Predictive Inverse Model for Advective Heat Transfer in a Short–Circuited Fracture: Dimensional Analysis, Machine Learning, and Field Demonstration

Identifying fluid flow maldistribution in planar geometries is a well–established problem in subsurface science/engineering. Of particular importance to the thermal performance of enhanced (or “engineered”) geothermal systems is identifying the existence of nonuniform (i.e., heterogeneous) permeability and subsequently predicting advective heat transfer. Here, machine learning via a genetic algorithm (GA) identifies the spatial distribution of an unknown permeability field in a two–dimensional Hele–Shaw geometry (i.e., parallel plates). The inverse problem is solved by minimizing the L2 norm between simulated residence time distribution (RTD) and measurements of an inert tracer breakthrough curve (BTC) (C–Dot nanoparticle). Principal component analysis (PCA) of spatially correlated permeability fields enabled reduction of the parameter space by more than a factor of 10 and restricted the inverse search to reservoir–scale permeability variations. Thermal experiments and tracer tests conducted at the mesoscale Altona Field Laboratory (AFL) demonstrate that the method accurately predicts the effects of extreme flow channeling on heat transfer in a single bedding–plane rock fracture. However, this is only true when the permeability distributions provide adequate matches to both tracer RTD and frictional pressure loss. Without good agreement to frictional pressure loss, it is still possible to match a simulated RTD to measurements, but subsequent predictions of heat transfer are grossly inaccurate. Here, the results of this study suggest that it is possible to anticipate the thermal effects of flow maldistribution, but only if both simulated RTDs and frictional pressure loss between fluid inlets and outlets are in good agreement with measurements.

42 ENGINEERING↗

CO 2 zonal injection rate allocation and plume extent evaluation through wellbore temperature analysis

Temperature analysis during a pause in injection operations, known as warmback analysis, has been used in the petroleum industry for evaluating the injection conformance and estimating the location of the flooded front in applications, such as waterflooding oil reservoirs. Here in this work, methods are introduced to extend the application of temperature warmback analysis to estimate the zonal CO 2 injection rate and zonal CO 2 plume extent during geologic CO 2 storage in a saline aquifer. First, novel analytical solutions are developed to model transient temperature in the aquifer during the injection and subsequent shut-in periods considering two-phase flow (gaseous CO 2 and aqueous brine) conditions in the aquifer. The solution involves a discretization of the aquifer into regions; the energy and mass conservation equations for the regions are solved simultaneously considering appropriate boundary conditions at the interfaces. Two solutions techniques are presented: multi-region and three-region solutions. Inverse models are developed accordingly to evaluate the injection profile and estimate the extent of the plume front in the reservoir during the injection period. The multi-region solution results in an inversion approach that requires regression analysis. However, the three-region formulation results in a simple graphical technique for inverse modeling. The analytical solutions are validated against a thermally coupled reservoir simulation tool using different synthetic cases for CO 2 injection in deep saline aquifers. The results of the developed solutions provide a good match with numerical results during forward and inverse modeling.

54 ENVIRONMENTAL SCIENCES↗

Utah FORGE: Southwestern Utah Magnetotelluric (MT) Data

This comprehensive magnetotellurics (MT) dataset, which covers southwestern Utah, integrates 600 sites from various surveys, including those from the Utah FORGE, SubTER, and Play Fairway projects, all of which are linked below. The core of this dataset is the use of a 3D finite element (FE) algorithm for inversion imaging, allowing for detailed analysis of the MT data. Provided here are all .edi files, the 3D finite element mesh (with a readme file), graphic slices, node graphics, resistivity graphics, and reports on the acquisition areas and inversion.

15 GEOTHERMAL ENERGY↗

Effect of solvothermal synthesis parameters on the crystallite size and atomic structure of cobalt iron oxide nanoparticles

We here investigate how the synthesis method affects the crystallite size and atomic structure of cobalt iron oxide nanoparticles. By using a simple solvothermal method, we first synthesized cobalt ferrite nanoparticles of ca. 2 and 7 nm, characterized by Transmission Electron Microscopy (TEM), Small Angle X-ray scattering (SAXS), X-ray and neutron total scattering. The smallest particle size corresponds to only a few spinel unit cells. Nevertheless, Pair Distribution Function (PDF) analysis of X-ray and neutron total scattering data shows that the atomic structure, even in the smallest nanoparticles, is well described by the spinel structure, although with significant disorder and a contraction of the unit cell parameter. These effects can be explained by the surface oxidation of the small nanoparticles, which is confirmed by X-ray near edge absorption spectroscopy (XANES). Neutron total scattering data and PDF analysis reveal a higher degree of inversion in the spinel structure of the smallest nanoparticles. Neutron total scattering data also allow magnetic PDF (mPDF) analysis, which shows that the ferrimagnetic domains correspond to ca. 80% of the crystallite size in the larger particles. A similar but less well-defined magnetic ordering was observed for the smallest nanoparticles. Finally, we used a co-precipitation synthesis method at room temperature to synthesize ferrite nanoparticles similar in size to the smallest crystallites synthesized by the solvothermal method. Structural analysis with PDF demonstrates that the ferrite nanoparticles synthesized via this method exhibit a significantly more defective structure compared to those synthesized via a solvothermal method.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A Butterfly-Accelerated Volume Integral Equation Solver for Broad Permittivity and Large-Scale Electromagnetic Analysis

In this work, a butterfly-accelerated volume integral equation (VIE) solver is proposed for fast and accurate electromagnetic (EM) analysis of scattering from heterogeneous objects. The proposed solver leverages the hierarchical off-diagonal butterfly (HOD-BF) scheme to construct the system matrix and obtain its approximate inverse, used as a preconditioner. Complexity analysis and numerical experiments validate the O(N log 2 N) construction cost of the HOD-BF-compressed system matrix and O(N log 1.5 N) inversion cost for the preconditioner, where N is the number of unknowns in the high-frequency EM scattering problem. For many practical scenarios, the proposed VIE solver requires less memory and computational time to construct the system matrix and obtain its approximate inverse compared to a H matrix-accelerated VIE solver. The accuracy and efficiency of the proposed solver have been demonstrated via its application to the EM analysis of large-scale canonical and real-world structures comprising of broad permittivity values and involving millions of unknowns.

42 ENGINEERING↗

Use of Gas Adsorption and Inversion Methods for Shale Pore Structure Characterization

The analysis of porosity and pore structure of shale rocks has received special attention in the last decades as unconventional reservoir hydrocarbons have become a larger parcel of the oil and gas market. A variety of techniques are available to provide a satisfactory description of these porous media. Some techniques are based on saturating the porous rock with a fluid to probe the pore structure. In this sense, gases have played an important role in porosity and pore structure characterization, particularly for the analysis of pore size and shapes and storage or intake capacity. In this review, we discuss the use of various gases, with emphasis on N2 and CO2, for characterization of shale pore architecture. We describe the state of the art on the related inversion methods for processing the corresponding isotherms and the procedure to obtain surface area and pore-size distribution. The state of the art is based on the collation of publications in the last 10 years. Limitations of the gas adsorption technique and the associated inversion methods as well as the most suitable scenario for its application are presented in this review. Finally, we discuss the future of gas adsorption for shale characterization, which we believe will rely on hybridization with other techniques to overcome some of the limitations.

04 OIL SHALES AND TAR SANDS↗

Chemically Recyclable Analogs of Styrene–Butadiene Copolymers Enabling Perfectly Linear Ethylene–Styrene Materials with Random Phenyl Distribution

Copolymerization of cyclopentene (CP) and 4-phenylcyclopentene (4PCP) at a full range of comonomer feed ratios is reported using Ru-based ring-opening metathesis polymerization (ROMP) yielding homogeneous copolymers analogous to poly(styrene-ran-1,4-butadiene) and poly(ethylene-ran-styrene) copolymers following hydrogenation under mild conditions. In all cases, total monomer conversions of 86%–92% yielded copolymers with compositions within 4% of monomer feeds. Analysis of equilibrium copolymerization thermodynamics, rarely performed on two cycloolefin monomers with low ring strain energies, provides rational design strategies for negotiating two monomers with different equilibrium monomer concentrations. Inverse-gated decoupled 13 C NMR analysis of dyad sequences on the resulting copolymer microstructures concludes a near-random distribution of comonomer units. The copolymers produced from ROMP have number-average molar masses up to 60 kg mol –1 , moderate dispersities (1.5 ± 0.1), and high trans olefin content (86% ± 2%) while glass-transitions temperatures follow the Fox equation and span the full range between homopolymer extremities of PCP (−96 °C) and P4PCP (17 °C). Unlike most prevulcanized elastomers, these materials undergo facile chemical recycling to monomer, producing complete ring-closing metathesis depolymerization (RCMD) of the polymer back to the CP comonomers. Quantitative olefin hydrogenation produced perfectly linear polyethylene with 4%–16% of the backbone carbons containing a phenyl pendant, analogous to ES copolymers with up to 71.5% w/w styrene units but with random distribution of the aromatic pendants. Thermal properties of these materials are discussed, which span from semicrystalline to amorphous, and with T g values notably less than the reported ES copolymer analogs at similar compositions.

animal feed↗

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↗

Forward and inverse modeling of fault transmissibility in subsurface flows

Characterizing physical properties of faults, such as their transmissibility, is crucial for performing predictive numerical simulation of subsurface flows, such as those encountered in petroleum engineering and remediation of subsurface contamination. Here, this paper provides a complete investigation of the inverse problem for fault transmissibility in subsurface flow models, under appropriate assumptions on fault structure. In particular, the following aspects are considered: 1) fault modeling and well-posedness of the forward problem; 2) finite element (FEM) discretizations of the forward problem and their rigorous a priori convergence analysis; 3) Well-posedness of the Bayesian inverse problem, FEM discretization of the infinite dimensional Bayesian inverse formulation, and its rigorous a priori analysis. Moreover, computation of the maximum a posteriori (MAP) point via fast inexact Newton-conjugate gradient optimization and a Laplace approximation of the Bayesian posterior are also presented. Numerical results illustrate the use of the proposed fault model in forward and inverse problems for subsurface flows in two dimensional domains with multiple faults.

97 MATHEMATICS AND COMPUTING↗

A novel closed-form inversion of the convection–diffusion equation for rapid convection, diffusion, and source profile estimation

To simplify and routinize particle transport analysis in fusion devices, a novel closed form linear inversion of the 1-D convection diffusion equation to estimate diffusion and convection profiles D(r ⃗ ), v(r ⃗ ) and source distribution s(r ⃗ ), of a single species from measured data is derived and demonstrated on synthetic data. Profile estimates of D(r ⃗ ), v(r ⃗ ), s(r ⃗ ) and their uncertainties are given as a matrix expression constructed directly from the incoming density data of the transported species in space and time, as well as physics assumptions such as particle conservation and experimental geometry. The derived matrix expression can be applied to a pumped or non-pumped recycling species, or a non-recycling species that is effectively “pumped” by plasma-facing surfaces.

Hinson, Edward [ORNL] (ORCID:000000019713140X)↗

Inverse design for waveguide dispersion with a differentiable mode solver

Inverse design of optical components based on adjoint sensitivity analysis has the potential to address the most challenging photonic engineering problems. However, existing inverse design tools based on finite-difference-time-domain (FDTD) models are poorly suited for optimizing waveguide modes for adiabatic transformation or perturbative coupling, which lies at the heart of many important photonic devices. Among these, dispersion engineering of optical waveguides is especially challenging in ultrafast and nonlinear optical applications involving broad optical bandwidths and frequency-dependent anisotropic dielectric material response. In this work, we develop gradient back-propagation through a general-purpose electromagnetic eigenmode solver and use it to demonstrate waveguide dispersion optimization for second harmonic generation with maximized phase-matching bandwidth. This optimization of three design parameters converges in eight steps, reducing the computational cost of optimization by ∼100x compared to exhaustive search and identifying new designs for broadband optical frequency doubling of laser sources in the 1.3–1.4 µm wavelength range. Furthermore, we demonstrate that the computational cost of gradient back-propagation is independent of the number of parameters, as required for optimization of complex geometries. This technique enables practical inverse design for a broad range of previously intractable photonic devices.

Gray, Dodd (ORCID:000000030469599X)↗

Efficient 3-D velocity model building using joint inline and crossline plane-wave wave-equation migration velocity analyses

SUMMARY Wave-equation migration velocity analysis (WEMVA) is an image-domain inversion method for velocity model building. Automatic plane-wave WEMVA (PWEMVA) calculates the moveouts of plane-wave common-image gathers (CIGs) by searching a best-fitting parabola with semblance analysis and backprojects residual CIG moveouts into wavefield wave paths with a reflection tomographic kernel. However, 3-D PWEMVA is very computationally expensive because 3-D reflection tomographic inversion requires at least five 3-D reverse-time migrations per iteration and stores two types of source wavefields at model boundaries. We develop a joint inline and crossline PWEMVA method for efficient 3-D velocity model building. We alternatively implement the inline and crossline PWEMVAs with a constraint for each other, in which we iteratively construct the 3-D velocity model update through 1-D spline interpolation of 2-D gradients. The inline and crossline joint inversion is practical since PWEMVA only inverts for low-wavenumber velocity perturbations along wave paths, and the method can take less than 1 per cent of the computational cost of full 3-D PWEMVA. To construct unaliased plane waves for our joint inline and crossline PWEMVA, we develop a 3-D data interpolation method in the frequency–wavenumber (FK) domain to recover regularly and randomly missing traces. The method minimizes the misfit on sufficiently localized data subsets with iterative optimal step lengths and a gradient preconditioner that iteratively selects dominant dips along different azimuths. In numerical experiments, we use a 3-D synthetic seismic data set and a land 3-D field seismic data set acquired at the Farnsworth CO2-EOR (enhanced oil recovery) field to demonstrate the efficacy of our velocity model building and data interpolation methods.

Liu, Xuejian↗