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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Accelerating full-waveform inversion using source stacking: synthetic experiments at the global scale in a realistic 3-D earth model

SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.

Geochemistry & Geophysics↗
DOE OSTI2025

Scattering-based structural inversion of soft materials via Kolmogorov–Arnold networks

Small-angle scattering techniques are indispensable tools for probing the structure of soft materials. However, traditional analytical models often face limitations in structural inversion for complex systems, primarily due to the absence of closed-form expressions of scattering functions. To address these challenges, we present a machine learning framework based on the Kolmogorov–Arnold Network (KAN) for directly extracting real-space structural information from scattering spectra in reciprocal space. This model-independent, data-driven approach provides a versatile solution for analyzing intricate configurations in soft matter. By applying the KAN to lyotropic lamellar phases and colloidal suspensions—two representative soft matter systems—we demonstrate its ability to accurately and efficiently resolve structural collectivity and complexity. Here, our findings highlight the transformative potential of machine learning in enhancing the quantitative analysis of soft materials, paving the way for robust structural inversion across diverse systems.

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

Probabilistic inversion of circular phase spectra: application to two-station phase-velocity dispersion estimation in western Canada

SUMMARY Periodic directional and temporal measurements are common in seismology, and necessitate specific statistical analyses that are appropriate for circular quantities. In this work, we explore the use of a von Mises distribution as a representation of errors on circular seismological observations. Specifically, we automate the estimation of surface-wave phase-velocity dispersion for the teleseismic two-station method, which generally suffers from a 2π phase ambiguity. The use of Bayesian inverse techniques, which aim to rigorously quantify model parameter uncertainty, have become widespread throughout seismology over the last decade. Here, we apply Bayesian inversion to measurements of surface-wave phase spectra in order to estimate 1-D, path-averaged Earth structure between station pairs. The dispersion curve and associated uncertainties are additional results of the inversion, which can then be used as input for subsequent analyses (e.g. tomography). We demonstrate this technique through application to surface-wave recordings from long-running seismic stations throughout western Canada. Our results for over 10 000 station pairs reveal first-order tectonic features consistent with previous studies, which provides confidence in our approach as well as an initial step towards resolving a full 3-D seismic velocity model for the region. This work also presents a foundation for the inversion of surface-wave phase spectra to estimate 3-D Earth structure directly. Finally, the ideas presented in this work are not limited to the inversion of surface-wave phase spectra, but can also be considered for Bayesian geophysical inversion of any circular quantities.

Gosselin, Jeremy M. (ORCID:0000000203754102)↗
DOE OSTI2021

Dairy Methane Emissions in California's San Joaquin Valley Inferred With Ground‐Based Remote Sensing Observations in the Summer and Winter

Abstract The dairy industry in the San Joaquin Valley (SJV) is one of California’s largest methane (CH 4 ) sources. Reducing dairy emissions is a priority for the state’s climate change plans. Observations of current dairy CH 4 emissions are key to monitoring actions taken toward this goal. To help support this, we present new ground‐based measurements of atmospheric column‐averaged CH 4 mixing ratio (XCH 4 ) gradients across a group of 600 dairies in the central SJV using EM27/SUN solar spectrometers. We used measurements from the 2019 summer and 2020 winter seasons for a top‐down emission inversion based on the WRF‐STILT model. Our top‐down estimates of the region’s dairy emissions range from 90% to 183% of the current CALGEM inventory’s emissions of 277 Gg/yr. In contrast to the strong temperature dependence found by earlier dairy CH 4 emission studies, we also find that our top‐down emissions during the winter measurement days are comparable to the summer measurement days, possibly due to seasonal changes in dairy management practices and meteorological conditions. Furthermore, we find significant interday variability in our measurements and find that our emission estimates overlap with earlier top‐down studies and bottom‐up inventories in this region. Our study demonstrates how analysis of ground‐based remotely sensed CH 4 gradient observations can help improve our understanding of CH 4 sources at scales relevant to mitigation policy. It also reflects the need for long‐term monitoring of CH 4 emissions in the region and at individual facilities to better understand their emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

WUS324: Multiscale Full Waveform Inversion Approaching Convergence Improves Waveform Fits While Imaging Seismic Structure of the Western United States

Abstract We report a new model of radially anisotropic crustal and upper mantle structure of the western United States (WUS324) obtained from full waveform inversion of earthquake data. We ran three multiscale inversion stages beyond model WUS256 (Rodgers et al., 2022, https://doi.org/10.1029/2022jb024549 ) allowing them to approach convergence to fit a larger data set to a shorter minimum period of 16 s. WUS324 is based on 324 total iterations from its starting model, significantly more (16 times) than previous studies. Waveform misfit reductions are 66%–70% for both the inversion data and an independent validation data set providing confidence in the predictive power of the model. WUS324 provides much better fits and reveals shear wavespeed, v S , structure of this large region with more detail than previous waveform tomography models. We show representative images demonstrating the resolution of diverse seismic structure across this highly heterogeneous region including oceanic lithosphere, subducting slabs and continental magmatism.

58 GEOSCIENCES↗
DOE OSTI2026

Wake-Resolving Acoustic Tomography: Advances through Numerical Covariance Methods

Acoustic tomography offers path-integrated measurements of atmospheric velocity and temperature fluctuations with high spatial resolution. Classical implementations of time-dependent stochastic inversion rely on homogeneous, isotropic covariance models that are poorly suited to the anisotropic structure of wind turbine wakes. By directly estimating heterogeneous covariances from large-eddy simulations (LESs) into the time-dependent stochastic inversion operator, we relax implicit assumptions in the analytical models used historically. Retrievals using these LES-informed models improve agreement with true fields in variance, turbulent kinetic energy, and spectral content compared to analytical and precursor-based covariance models. The results indicate that LES-informed covariance models can enhance the accuracy of acoustic tomography retrievals in complex, anisotropic flows such as wind turbine wakes in some cases and highlight instances where analytical models still offer competitive performance, despite their simplifying assumptions.

17 WIND ENERGY↗
DOE OSTI2024

Inverse design of a pyrochlore lattice of DNA origami through model-driven experiments

Sophisticated statistical mechanics approaches and human intuition have demonstrated the possibility of self-assembling complex lattices or finite-size constructs. However, attempts so far have mostly only been successful in silico and often fail in experiment because of unpredicted traps associated with kinetic slowing down (gelation, glass transition) and competing ordered structures. Theoretical predictions also face the difficulty of encoding the desired interparticle interaction potential with the experimentally available nano- and micrometer-sized particles. To overcome these issues, we combine SAT assembly (a patchy-particle interaction design algorithm based on constrained optimization) with coarse-grained simulations of DNA nanotechnology to experimentally realize trap-free self-assembly pathways. In this paper, we use this approach to assemble a pyrochlore three-dimensional lattice, coveted for its promise in the construction of optical metamaterials, and characterize it with small-angle x-ray scattering and scanning electron microscopy visualization.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
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 OSTI2026

Inverse Lieb materials: altermagnetism and more

The Lieb lattice, originally proposed for cuprate superconductors, has gained new attention in the emerging field of altermagnetism as a minimal analytical model for the latter. While the inverse Lieb lattice (ILL) was once considered purely theoretical, several materials with this motif have recently been discovered. Its unique geometry supports complex magnetic orders driven by geometric frustration, offering high tunability. In this work, we provide comprehensive insights into ILL magnetic phases and establish guidelines for identifying altermagnets. Using a Heisenberg model, we first construct phase diagrams to elucidate the mechanisms underlying experimental magnetic phases. We then bridge theory and experiment via density functional theory (DFT) calculations on existing ILL compounds, finding results consistent with experimental data. Notably, we identified a trend linking d-shell filling to magnetic order, where d x (x≤5) configurations show a propensity for altermagnetism. We also highlight Sr 2 CrO 2 Cr 2 OAs 2 as a promising metallic altermagnet with highly anisotropic J 2 exchange and a high Néel temperature (~600 K). Finally, our magnon spectra calculations confirm that chiral splittings correlate directly with anisotropies between inequivalent J 2 interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Connecting Hindered Transport in Porous Media across Length Scales: From Single-Pore to Macroscopic

Hindered mass transport is widely observed in various porous media, including filtration membranes and porous rocks. However, there is no universal model capable of predicting transport in porous media due to the heterogeneity of porous structures and the complexity of the underlying microscopic mechanisms, including electrostatic and hydrodynamic effects. Here, we used a highly ordered porous medium (i.e. a silica inverse opal film) as a model system to directly explore the effects of geometric parameters (i.e. pore size, pore throat size and tracer particle size) and microscopic interaction parameters (e.g. controlled by ionic strength) on nanoparticle transport in porous environments using three-dimensional single-particle tracking. We found a linear scaling relation between the macroscopic diffusion coefficient and microscopic diffusion behaviors and a combination of parameters associated with microscopic features and phenomena, including both geometric effects and particle-level interactions. The proportionality coefficient relating micro and macro behaviors was complex and related to the connectivity of the matrix and pore-size variation, which could lead to tortuous diffusion pathways, hindering macroscopic transport. Furthermore, these findings provide important insights into nanoparticles and macromolecules transport in various porous media.

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

Probing the Fine Symmetry Breaking in High-Temperature Superconductor Bi 2 Sr 2 CaCu 2 O 8+δ with Angle-Resolved Nonreciprocal Transport

After decades of research, symmetry breaking in high-temperature cuprate superconductors remains a key issue to resolve and is relevant to understanding their exotic quantum phases. In the prototypical cuprate superconductor, Bi 2 Sr 2 CaCu 2 O 8+δ (Bi2212), the possible symmetry breaking has been mostly examined microscopically with scanning tunneling microscopy and photoemission spectroscopy. However, macroscopic evidence and the direct implications for electronic transport have remained elusive. Using superconductivity-enhanced nonreciprocal transport, we report macroscopic evidence of inversion symmetry breaking in Bi2212. While the inversion symmetry breaking is subtle, its effect on nonreciprocal transport is significantly enhanced by the vortex motion during the superconducting transition, leading to a robust manifestation of inversion symmetry breaking in macroscopic transport. Combining angle-resolved nonreciprocal transport and 3D tight-binding model calculations, we derive that the inversion symmetry breaking is due to subtle crystal distortions that give rise to both in-plane and out-of-plane polar axes. Our work not only establishes nonreciprocal transport as a sensitive macroscopic probe to provide electrical transport-based evidence of fine symmetry breaking in Bi2212, but also paves the way for novel device applications such as high-temperature superconducting diodes and superconducting spintronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Simulating Heisenberg interactions in the Ising model with strong drive fields

The time evolution of an Ising model with large driving fields over discrete time intervals is shown to be reproduced by an effective XXZ-Heisenberg model at leading order in the inverse field strength. For specific orientations of the drive field, the dynamics of the XXX-Heisenberg model is reproduced. Finally, these approximate equivalences, valid above a critical driving field strength set by dynamical phase transitions in the Ising model, are expected to enable quantum devices that natively evolve qubits according to the Ising model to simulate more complex systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Observation and Simulation of Methane Plumes During the Morning Boundary Layer Transition

Abstract Methane (CH 4 ) contributes significantly to global warming. However, accurate identification of CH 4 sources for reducing CH 4 emissions is often hampered by inadequate accuracy and spatiotemporal coverage of CH 4 detection, and lack of accurate CH 4 forward modeling used in top‐down inversion systems. In this study, a field experiment was conducted in Pampa, Texas using two CH 4 sensors (LI‐COR and OGI camera) to detect CH 4 releases. We investigated whether high‐resolution simulations using the Weather Research and Forecasting (WRF) model with greenhouse gases (WRF‐GHG) could accurately simulate the CH 4 plumes in the presence of evolving atmospheric boundary layer from sunrise to noon. CH 4 plumes showed substantial variation in time. At a release rate of ∼17.5 kg hr −1 , the maximum enhancement of CH 4 measured by LI‐COR was 2.6 ppm at sunrise (7:36 a.m.), 250 m from the release location. Within half an hour after sunrise, this enhancement decreased to 0.3–0.4 ppm. The enhancement was 0.2 ppm by 10:00 a.m. and further dropped to less than 0.1 ppm after 11:30 a.m. Due to the low temperature at sunrise, the OGI camera failed to detect the CH 4 plume. The WRF‐GHG large‐eddy simulation (LES) with 32 m grid spacing successfully reproduced these CH 4 enhancements. In situ measurements together with numerical simulations illustrate the impact of the transition from a stable boundary layer in the early morning to a convective boundary layer at noon on the dispersion of CH 4 plumes. Additionally, CH 4 plumes from a cattle farm in Oklahoma are briefly examined using the same modeling approach.

Hu, Xiao‐Ming [Center for Analysis and Prediction ↗
DOE OSTI2024

Deep Learning for Full Waveform Inversion of Elastic Active-Source Seismic Data to Estimate P-Wave Velocity Models

Seismic imaging methods are critical for Global Security and Energy & Homeland Security missions and activities that rely on subsurface characterization, but traditional methods remain computationally expensive and require significant labor hours and expertise to execute. Within the past few years, machine learning (ML), namely deep learning (DL), has been used to develop data-driven end-to-end full waveform inversion (FWI) methods to estimate 2D P-wave velocity (Vp) models in a fraction of the time as conventional FWI. These methods, however, are trained on simplistic acoustic wave seismic data and Vp models that are not realistic nor representative of real-world observations, leaving a large gap between the state-of-the-art and deployable, feasible, and practical DL FWI methods. Here, we generate a synthetic active-source, 3D, elastic wave seismic data set and a variety of Vp models with realistic geologic structure for training DL FWI methods. We evaluate six different methods that have performed well for acoustic DL FWI or medical imaging tasks using our more realistic dataset. We find that these six trained models do not match the performance of published acoustic end-to-end DL FWI methods, indicating more training data may be needed, physics may need to be incorporated to achieve good accuracy at the sacrifice of the end-to-end advantage, and/or novel methods need to be developed to enable end-to-end DL FWI methods to perform well for real-world seismic data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

In-depth analysis of solar models with high-metallicity abundances and updated opacity tables

Context. As a result of the high-quality constraints available for the Sun, we are able to carry out detailed combined analyses using neutrino, spectroscopic, and helioseismic observations. These studies lay the ground for future improvements of the key physical components of solar and stellar models because ingredients such as the equation of state, the radiative opacities, or the prescriptions for macroscopic transport processes of chemicals are then used to study other stars in the Universe. Aims. We study the existing degeneracies in solar models using the recent high-metallicity spectroscopic abundances by comparing them to helioseismic and neutrino data and discuss the effect on their properties of changes in the micro and macro physical ingredients. Methods. We carried out a detailed study of solar models computed with a high-metallicity composition from the literature based on averaged 3D models that were claimed to resolve the solar modelling problem. We compared these models to helioseismic and neutrino constraints. Results. The properties of the solar models are significantly affected by the use of the recent OPLIB opacity tables and the inclusion of macroscopic transport. The properties of the standard solar models computed using the OPAL opacities are similar to those for which the OP opacities were used. We show that a modification of the temperature gradient just below the base of the convective zone is required to remove the discrepancies in solar models, particularly in the presence of macroscopic mixing. This can be simulated by a localised increase in the opacity of a few percent. Conclusions. We conclude that the existing degeneracies and issues in solar modelling are not removed by using an increase in the solar metallicity, in contradiction to what has been suggested in the recent literature. Therefore, standard solar models cannot be used as an argument for a high-metallicity composition. While further work is required to improve solar models, we note that direct helioseismic inversions indicate a low metallicity in the convective envelope, in agreement with spectroscopic analyses based on full 3D models.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Robust Optimal Experimental Design of Infinite-Dimensional Bayesian Nonlinear Inverse Problems

Abstract. We consider robust optimal experimental design (ROED) for nonlinear Bayesian inverse problems governed by partial differential equations (PDEs). An optimal design is one that maximizes some utility quantifying the quality of the solution of an inverse problem. However, the optimal design is dependent on elements of the inverse problem such as the simulation model, the prior, or the measurement error model. ROED aims to produce an optimal design that is aware of the additional uncertainties encoded in the inverse problem and remains optimal even after variations in them. We follow a worst-case scenario approach to develop a new framework for robust optimal design of nonlinear Bayesian inverse problems. The proposed framework (a) is scalable and designed for infinite-dimensional Bayesian nonlinear inverse problems constrained by PDEs; (b) develops efficient approximations of the utility, namely the expected information gain; (c) employs eigenvalue sensitivity techniques to develop analytical forms and efficient evaluation methods of the gradient of the utility with respect to the uncertainties against which we wish to be robust; and (d) employs a probabilistic optimization paradigm that properly defines and efficiently solves the resulting combinatorial max-min optimization problem. The effectiveness of the proposed approach is illustrated for optimal sensor placement problem in an inverse problem governed by an elliptic PDE.

Chowdhary, Abhijit↗
DOE OSTI2023

Machine learning-based inversion for acoustic impedance with large synthetic training data: Workflow and data characterization

Where wells are sparse or training data are difficult to label with high-quality wireline-derived impedance logs, machine learning (ML)-based inversion of acoustic impedance typically depends on small training data sets, leading to biased prediction. We have advanced a novel workflow that applies large synthetic seismic training data to reduce facies-related bias. Using a geologically realistic model as the truth model, we randomly select sparse seed wells to perform sequential Gaussian simulation (SGS) for impedance models of the same geometry and simulate facies variability. We implement random forest regression on 30 features extracted from the synthetic volume. We observe that more seed wells tend to reduce facies-induced bias by sampling more types of facies, resulting in a better prediction. We then focus on the responses of SGS models to facies changes, the number of seed wells necessary for a useful synthetic model, and how much a synthetic model can help ML-based inversion. Here, we observe that the SGS synthetic training model outperforms well-direct training in general. For modeled clastic shore-zone systems in Miocene Gulf of Mexico, two or more seed wells are necessary for a significant reduction of root-mean-square error and outliners, and improvement of facies imaging. In a field-data test, we apply a similar workflow to quantitatively predict acoustic impedance, which is then converted to a sand-volume map at a high-frequency sequence (10–100 m), revealing detailed facies and sandstone patterns. Such results are valuable in many geologic and engineering applications, such as hydrocarbon and CO 2 reservoir prospecting, reserve estimation, simulation, etc.

3D seismic↗