Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium reconstruction”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ringing Artifact Reduction Method for Ultrasound Reconstruction Using Multi-Agent Consensus Equilibrium

Non-destructive characterization of multi-layered structures that can be accessed from only a single side is important for applications such as well-bore integrity inspection. Existing methods related to Synthetic Aperture Focusing Technique (SAFT) rapidly produce acceptable results but with significant artifacts. Recently, ultrasound model-based iterative reconstruction (UMBIR) approaches have shown significant improvements over SAFT. However, even these methods produce ringing artifacts due to the high fractional-bandwidth of the excitation signal.In this paper, we propose a ringing artifact reduction method for ultrasound image reconstruction that uses a multi-agent consensus equilibrium (RARE-MACE) framework. Our approach integrates a physics-based forward model that accounts for the propagation of a collimated ultrasonic beam in multi-layered media, a spatially varying image prior, and a denoiser designed to suppress the ringing artifacts that are characteristic of reconstructions from high-fractional bandwidth ultrasound sensor data. We test our method on simulated and experimental measurements and show substantial improvements in image quality compared to SAFT and UMBIR.

Alanazi, Abdulrahman M↗
DOE OSTI2024

Reconstruction of magnetic island electron temperature in mixed second and third harmonic electron cyclotron emission conditions

Here we develop a method to use the mixed third and second harmonic electron cyclotron emission (ECE) signal in the DIII-D tokamak to reconstruct the electron temperature profile of a rotating magnetic island. The third harmonic ECE is removed by extracting the rotating-island-associated fluctuations in the mixed signal, and the extracted fluctuation is combined with the equilibrium temperature obtained from other diagnostics after correcting for the third harmonic reabsorption. The accuracy of the reconstruction is studied by considering a DIII-D shot where an unmixed signal from an island is available on the low field side of the plasma and a mixed signal from the same island is available from the high field side. It is found that the reconstruction method successfully reproduces the island shape and temperature perturbation magnitude without the distortion caused by third harmonic ECE mixing. However, the radial location of the reconstructed island is somewhat displaced relative to the location of the q = 2 surface in the axisymmetric equilibrium reconstruction, resulting in a corresponding inaccuracy in the absolute temperature of the island. It is conjectured that this may arise from an inaccuracy of the reconstructed axisymmetric equilibrium in this region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Climate reconstructions for the Last Glacial Maximum from a simple cirque glacier in Fiordland, New Zealand

Mountain glacier records offer important constraints on the timing and magnitude of climate variability during the last glacial cycle. Existing moraine chronologies from the central Southern Alps indicate maximum ice extent was achieved during marine isotope stages 3–4, followed by repeated advances of similar, but gradually declining extent during marine isotope stage 2, until onset of the glacial termination. Questions remain over the precise role of climate in driving these changes, as most existing moraine chronologies come from large, complex former valley glacier systems, where non-climatic influences such as changing bed topography and proglacial lake formation may have influenced glacier length changes. In this work, we address this problem via a new cosmogenic 10Be chronology and equilibrium line altitude reconstruction from a cirque glacier situated in Fiordland, New Zealand. Our chronology shows moraine deposition at 32 ± 11 ka, 18.7 ± 0.2 ka, 18.1 ± 0.1 ka, and c. 17.2 ± 0.3 ka. The simple geometry of the former glacier supports the role of climate in driving a net decline of regional ice volume during marine isotope stages 3–2. Close spacing and good preservation of the 19–17 ka moraines permits 2D glacier reconstruction which suggests the equilibrium line altitude remained depressed by c.1130 m (equivalent to 5.8 ± 0.6 °C colder than present) during this interval. Onset of warming after 17.2 ± 0.2 ka is consistent with climate proxy evidence for a sustained southward shift in the southern westerly winds, which may have promoted deglaciation via shifting ocean currents and promoting increases in atmospheric carbon dioxide.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Neural net modeling of equilibria in NSTX-U

Neural networks (NNs) offer a path towards synthesizing and interpreting data on faster timescales than traditional physics-informed computational models. In this work we develop two NNs relevant to equilibrium and shape control modeling, which are part of a suite of tools being developed for the National Spherical Torus Experiment-Upgrade for fast prediction, optimization, and visualization of plasma scenarios. The networks include Eqnet, a free-boundary equilibrium solver trained on the EFIT01 (Equilibrium FITtting 01) reconstruction algorithm, and Pertnet, which is trained on the Gspert code and predicts the non-rigid plasma response, a nonlinear term that arises in shape control modeling. The NNs are trained with different combinations of inputs and outputs in order to offer flexibility in use cases. In particular, Eqnet can use magnetic diagnostics as inputs and act as an EFIT-like reconstruction algorithm, or, by using pressure and current profile information the NN can act as a forward Grad–Shafranov equilibrium solver. This forward-mode version is envisioned to be implemented in the suite of tools for simulation of plasma scenarios. The reconstruction-mode version gives some performance improvements compared to the online reconstruction code real-time EFIT, especially when vessel eddy currents are significant. Here, we report strong performance for all NNs indicating that the models could reliably be used within closed-loop simulations or other applications. Some limitations are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2015

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗
DOE OSTI2025

Unraveling the Surface Termination and Evolution of Surface States for Electrocatalyst PtSn 4 in Alkaline HER

Semimetal PtSn 4 has been experimentally demonstrated as a promising topological electrocatalyst for the hydrogen evolution reaction (HER) under both acidic and alkaline conditions. While two possible mechanisms have been proposed to explain its activity, the role of its surface states in HER remains unclear. It is indeed in question how the surface states of this alloy evolve as HER proceeds. In this study, we investigate the surface termination that sustains conducting surface states on PtSn 4 , and we track their evolution during HER catalysis. We show that a reconstructed surface with a Sn-poor termination reproduces the scanning tunneling microscopy pattern observed in experiments and sustains a conducting surface. Through phase diagram and geometric structure analysis, we outline the HER profile following the Volmer–Heyrovsky mechanism. As hydrogen atoms adsorb onto the surface, the structure undergoes further reconstruction to an equilibrium phase with a coverage of two hydrides per unit cell. Meanwhile, the surface electronic bands evolve in response to interactions with the adsorbed hydrogen atoms. A hybridization diagram is further proposed for understanding the surface state evolution based on wave function and chemical bonding analyses. While the Pt atoms serve as conventional sites for hydrogen binding, the surface states of PtSn 4 are essential for stabilizing the hydrogen antibonding states via in-phase electronic interactions with the Sn components. This stabilization results in frontier surface bands that are responsible for driving the HER catalysis. Here, our findings provide a detailed description for the direct involvement of surface states on PtSn4 when employed as a catalyst for HER.

catalysts↗
DOE OSTI2023

An optical-input Maximum Likelihood Estimation feedback system demonstrated on tokamak horizontal equilibrium control

A readily parallelized Maximum Likelihood Estimation (MLE) algorithm with linear computational complexity is demonstrated in real time using only measurements from an extreme ultraviolet (EUV) diagnostic to control the horizontal position of a tokamak plasma. A set of trial emissivity profiles are parameterized by the control quantity of interest (R m ), and the MLE is identified from the profile which minimizes the signal reconstruction residual. The algorithm depends on an empirically determined likelihood function with exponential form. EUV emission (λ ≈ 15eV-1keV) is captured in a poloidal plane by four 16-channel AXUV diodes mounted at different poloidal angles with radial and angular resolution sufficient to discern plasma equilibrium evolution in HBT-EP. Calculations of the plasma major radius by the system are consistent within diagnostic uncertainty for the majority of the discharge with those of: a weighted average of vertical soft X-ray or EUV chords, magnetic sensors, and an equilibrium reconstruction. The feedback system corrects for a horizontal displacement of the major radius equal to 20% of the plasma minor radius by adjusting the vertical field produced from 40 in-vessel control coils in real time. The MLE calculation is performed on a GPU in a 15 μs cycle, with similar performance in this application to a simple weighted average of vertical chords. Finally, results demonstrate horizontal position control using magnetic actuators and an optical observer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Physics‐based iterative reconstruction for dual‐source and flying focal spot computed tomography

Purpose For single‐source helical Computed Tomography (CT), both Filtered‐Back Projection (FBP) and statistical iterative reconstruction have been investigated. However, for dual‐source CT with flying focal spot (DS‐FFS CT), a statistical iterative reconstruction that accurately models the scanner geometry and acquisition physics remains unknown to researchers. Therefore, our purpose is to present a novel physics‐based iterative reconstruction method for DS‐FFS CT and assess its image quality. Methods Our algorithm uses precise physics models to reconstruct from the native cone‐beam geometry and interleaved dual‐source helical trajectory of a DS‐FFS CT. To do so, we construct a noise physics model to represent data acquisition noise and a prior image model to represent image noise and texture. In addition, we design forward system models to compute the locations of deflected focal spots, the dimension, and sensitivity of voxels and detector units, as well as the length of intersection between x‐rays and voxels. The forward system models further represent the coordinated movement between the dual sources by computing their x‐ray coverage gaps and overlaps at an arbitrary helical pitch. With the above models, we reconstruct images by an advanced Consensus Equilibrium (CE) numerical method to compute the maximum a posteriori estimate to a joint optimization problem that simultaneously fits all models. Results We compared our reconstruction with Siemens ADMIRE, which is the clinical standard hybrid iterative reconstruction (IR) method for DS‐FFS CT, in terms of spatial resolution, noise profile, and image artifacts through both phantoms and clinical scan datasets. Experiments show that our reconstruction has a higher spatial resolution, with a Task‐Based Modulation Transfer Function (MTF task ) consistently higher than the clinical standard hybrid IR. In addition, our reconstruction shows a reduced magnitude of image undersampling artifacts than the clinical standard. Conclusions By modeling a precise geometry and avoiding data rebinning or interpolation, our physics‐based reconstruction achieves a higher spatial resolution and fewer image artifacts with smaller magnitude than the clinical standard hybrid IR.

Wang, Xiao↗
DOE OSTI2024

Validation of NSFsim as a Grad-Shafranov equilibrium solver at DIII-D

Plasma shape is a significant factor that must be considered for any Fusion Pilot Plant (FPP) as it has significant consequences for plasma stability and core confinement. A new simulator, NSFsim, has been developed based on a historically successful code, DINA [1], offering tools to simulate both transport and plasma shape. Specifically, NSFsim is a free boundary equilibrium and transport solver and has been configured to match the properties of the DIII-D tokamak. This paper is focused on validating the Grad-Shafranov (GS) solver of NSFsim by analyzing its ability to recreate the plasma shape, the poloidal flux distribution, and the measurements of the simulated diagnostic signals originating from flux loops and magnetic probes in DIII-D. Five different plasma shapes are simulated to show the robustness of NSFsim to different plasma conditions; these shapes are Lower Single Null (LSN), Upper Single Null (USN), Double Null (DN), Inner Wall Limited (IWL), and Negative Triangularity (NT). The NSFsim results are compared against real measured signals, magnetic profile fits from EFIT [2], and another plasma equilibrium simulator, GSevolve [3]. EFIT reconstructions of shots are readily available at DIII-D, but GSevolve was manually ran by us to provide simulation data to compare against.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Factors controlling the oxygen isotopic composition of lacustrine authigenic carbonates in Western China: implications for paleoclimate reconstructions

Abstract In the carbonate-water system, at equilibrium, the oxygen isotopic composition of carbonate is dependent not only on the temperature but also on the isotopic composition of host water in which the carbonate is formed. In this study, lake surface sediment and water samples were collected from 33 terminal lakes in Western China to evaluate controls on the oxygen isotopic composition of lacustrine authigenic carbonates (δ 18 O carb ) and its spatial distribution. Our results show that water oxygen isotopic composition (δ 18 O water ) rather than lake summer water temperature (T water ), is the main determinant of δ 18 O carb , irrespective of whether oxygen isotope equilibrium is achieved. There are significant linear correlations between δ 18 O carb and elevation, as well as that between δ 18 O carb and latitude for lakes located on the Tibetan Plateau. In Western China, the spatial distribution of δ 18 O carb is consistent with that of δ 18 O water , and is ultimately controlled by the isotopic composition of local precipitation (δ 18 O precipitation ) that depends on the source of water vapor. Therefore, changes in δ 18 O carb can be predominantly interpreted as variations of δ 18 O water , which in turn represent changes in δ 18 O precipitation for paleoclimate reconstructions in this region, and may be relevant to studies of other areas.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1973

Equilibrium condensation in a solar nebula

In attempts to reconstruct the environment of condensation of solar system materials, particularly exemplified by certain meteorite components, the relative temperatures of the gas and the solid are of critical importance. The relationships that determine the heat balance in a circumsolar grain-gas system are examined. Fundamental considerations show that regardless of opacity or gas density, the gas will always be at a higher temperature than the solid in such regions of the system where condensation is possible. Implications of the characteristic temperature differential between the gas and the condensing solid are discussed.

Arrhenius, G.↗
DOE OSTI2022

Cosmogenic 10Be constraints on deglacial snowline rise in the Southern Alps, New Zealand

We report geochronological dating of glacial landforms, such as terminal and lateral moraines, is useful for determining the extent and timing of past glaciation and for reconstructing the magnitude and rate of past climate changes. In the Southern Alps of New Zealand, well-dated glacial geomorphological records constrain the last glacial cycle across much of the Waitaki River basin (e.g. Ōhau, Pukaki, Tekapo) but its southern sector such as the Ahuriri River valley remains comparatively unconstrained. Recently, there has been debate on the scale and rapidity of mountain glacier retreat during the last glacial termination, particularly the 20–17 ka period in New Zealand. Missing from this debate is well-constrained equilibrium-line altitude (ELA) and associated temperature reconstructions, particularly over the period around 17 ka, which can help us to develop a more complete picture of how past temperature changes drove glacier retreat. Here we report the first glacial chronology dataset from the Last Glacial Maximum (LGM) and subsequent deglaciation from the Ahuriri River valley, Southern Alps, New Zealand (44°23'54"S, 169°39'48"E) based on 38 beryllium-10 ( 10 Be) surface-exposure ages from terminal moraine systems and glaciated bedrock situated at the lower and middle sections of the valley. Our results show that the former Ahuriri Glacier reached its maximum extent at 19.8 ± 0.3 ka, which coincides with the global Last Glacial Maximum. By 16.7 ± 0.3 ka, the glacier had retreat ~18 km up-valley suggesting at least ~43% glacier-length loss relative to its full LGM extent. This deglaciation was accompanied by the formation of a shallow proglacial lake. Using the accumulation area ratio (AAR) method, we estimate that the ELA was lower than present by ~880 m (~1120 m a.s.l.) at 19.8 ± 0.3 ka, and ~770 m lower (~1230 m a.s.l.) at 16.7 ± 0.3 ka. Applying an estimate for temperature lapse rate, this ELA anomaly implies that local air temperature was 5 ± 1 °C colder than present (1981–2010) at 19.8 ± 0.3 ka, while it was 4.4 ± 0.9 °C colder at 16.7 ± 0.3 ka, assuming no change in precipitation. The substantial glacier retreat in response to a relatively small accompanying increases in ELA (110 m) and temperature (0.6 °C) may have been a result of the high glacier-length sensitivity of this glacier system due to its low gradient of former ice surface. Our low warming estimate differs markedly from other deglaciation studies, specifically from Rakaia River valley, which reports a much larger temperature increase at the onset of the last deglaciation. This precisely-dated moraine record along with reconstructed ELA as proxies for atmospheric conditions, provides new insight into post LGM glacier behaviour and climate conditions in New Zealand.

58 GEOSCIENCES↗
DOE OSTI2024

Iterative Stress Reconstruction Algorithm to Estimate Three-Dimensional Residual Stress Fields in Manufactured Components

Residual stress (RS) significantly impacts the mechanical performance of components. Measurement of RS often provides incomplete data in terms of components of stress and spatial density. Employing such fields in finite element simulations results in significant modification of the field to achieve equilibrium and compatibility among strains. To overcome this, an iterative stress reconstruction algorithm (ISRA) is developed to estimate 3D RS fields that satisfy equilibrium, are stress component-wise complete, and represent the characterized data sampled. An Al 7075-T651 plate and an additively manufactured (AM) A36 steel wall are considered for RS reconstruction using measurement data from the literature. A maximum variation of ~2.5 MPa in the Al plate, and ~10 MPa in the steel wall are observed between the reconstructed and measured stresses. Furthermore, unknown stress components emerge and reach significant magnitudes (upto ~2.3 MPa in the Al plate and ~45 MPa in the AM wall) during ISRA. Indeed, it is found that minor errors in measurement or data processing are eliminated through the physical requirements during ISRA. Employing a reconstructed RS field is hence not just more accurate given its compatibility, but it additionally corrects for minor errors in measurement. Furthermore, it is found that spatially dense measurement data result in convergence with fewer iterations. Finally, although ISRA yields a nonunique solution dependent on boundary conditions, measurement errors, fitting errors, and mesh density, it accommodates for uncertainties and inaccuracies in measurement, as opposed to failing to reach a physically realistic converged solution.

42 ENGINEERING↗
NASA NTRS2004

High-Order Simulation of Non-Linear Oscillations and Shocks in the Solar Atmosphere

The solar atmosphere presents a rich source of highly non-linear magneto-hydrodynamic phenomena: strong gradients and forcing terms result in both large shocks and oscillations. The additional requirements of energy balance and initialization in hydrostatic equilibrium compound the challenge of this problem. A wealth of observational data allows us to check the results of our simulations. The problem of simulating the solar atmosphere provides, in addition to an interesting system in its own right, a challenging testbed for high-order shock-capturing methods. We discuss the challenge of simulating solar atmospheric phenomena, concentrating on various high-order central methods ranging from second to fourth order. Our method is based on the central-upwind scheme of Kurganov, Noelle and Petrova, which we extend to high order via various interpolants. We investigate various initial data for our simulations, corresponding to observed conditions in different regions of the solar surface: the normal quiet sun and sunspots. When non-oscillatory using second- and third-order methods, we are able to reproduce non-trivial observational results. In particular we find a correlation between initial data and both the shock speeds and particle oscillation spectra that match observations in the corresponding regions. When using fourth-order WENO interpolants, we find that while the individual shock profiles at any given time appear non-oscillatory, spurious oscillations appear in the fields after long time integrations. The issue of initialization in hydrostatic equilibrium raises difficult issues. Careful treatment of the gravitational source term can reduce violations of hydrostatic equilibrium, but difficulties remain, primarily due to discontinuities in the piecewise-polynomial reconstructions. In the case of the Euler equations in gravity, the use of high-order methods reduces the violation of hydrostatic equilibrium to a sufficiently low level for the above described results to be obtained. Including various energy flux terms in the equations, however, leads to severe instability when hydrostatic balance is violated. We discuss various strategies to address this problem, including adaptive grid methods.

Bryson, S.↗
NASA NTRSDate not supplied

Using Multiplicative Algebraic Reconstruction Techniques (MART) to Derive Instrument Requirements for Computed Tomography Imaging Spectrographs

Computed Tomography Imaging Spectrographs (CTISs), which are generally slitless or large-aperture spectrographs that observe a wide, dispersed field of view in multiple dispersion angles or diffraction orders, provide a unique opportunity to capture spectral information over an extended source, such as the Sun, but require software reconstruction techniques to be fully utilized. A useful data product that can be recovered from CTIS observations is intensity maps of solar features in single spectral lines. These intensity maps can then be used to determine temperature, density, abundance, and equilibrium properties of the emitting plasma. We apply the multiplicative algebraic reconstruction technique (MART) to example data from a variety of CTIS configurations to determine the capability and limitations of the method to return spectrally-pure maps. By completing this study, we aim to establish a path to derive requirements for CTIS instruments that rely on reconstruction techniques to meet their science objectives.

Amy Winebarger↗
DOE OSTI2025

Transferable dispersion-aware machine learning interatomic potentials for multilayer transition metal dichalcogenide heterostructures

Stacking atomically thin transition metal dichalcogenides (TMDs) into heterostructures enables exploration of exotic quantum phases, particularly through twist-angle-controlled moiré superlattices. These structures exhibit novel electronic and optical behaviors driven by atomic-scale structural reconstruction. However, studying such systems with DFT is computationally demanding due to their large unit cells and van der Waals (vdW) interactions between layers. To address this, we develop a transferable neural network potential (NNP) that includes long-range vdW corrections up to 12Å with minimal overhead. Trained on vdW-corrected DFT data for Mo- and W-based TMDs with S, Se, and Te, the NNP accurately models monolayers, bilayers, heterostructures, and their interaction with h-BN substrates. It reproduces equilibrium structures, energy landscapes, phonon dispersions, and matches experimental atomic reconstructions in twisted WS2 and MoS2/WSe2 systems. We demonstrate that our NNP achieves DFT-level accuracy and high computational efficiency, enabling large-scale simulations of TMD-based moiré superlattices both with and without substrates.

materials for energy and catalysis↗
DOE OSTI2023

Projected Multi-Agent Consensus Equilibrium (PMACE) With Application to Ptychography

Multi-Agent Consensus Equilibrium (MACE) formulates an inverse imaging problem as a balance among multiple update agents such as data-fitting terms and denoisers. However, each such agent operates on a separate copy of the full image, leading to redundant memory use and slow convergence when each agent affects only a small subset of the full image. In this article, we extend MACE to Projected Multi-Agent Consensus Equilibrium (PMACE), in which each agent updates only a projected component of the full image, thus greatly reducing memory use for some applications. We describe PMACE in terms of an equilibrium problem and an equivalent fixed-point problem and show that in most cases the PMACE equilibrium is not the solution of an optimization problem. To demonstrate the value of PMACE, we apply it to the problem of ptychography, in which a sample is reconstructed from the diffraction patterns resulting from coherent X-ray illumination at multiple overlapping spots. In our PMACE formulation, each spot corresponds to a separate data-fitting agent, with the final solution found as an equilibrium among all the agents. In conclusion, our results demonstrate that the PMACE reconstruction algorithm generates more accurate reconstructions at a lower computational cost than existing ptychography algorithms when the spots are sparsely sampled.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Consensus Equilibrium for Subsurface Delineation

Heterogeneity and insufficient site characterization limit our knowledge of the subsurface. Inversion techniques, which minimize the mismatch between observations and model predictions, have become an essential tool of subsurface characterization. Most optimization-based approaches fail to incorporate various implicit priors and capture the geological complexity. We overcome these limitations by deploying the plug-and-play and consensus equilibrium (CE) strategies, which provide a flexible framework for image reconstruction. Our CE methodology for spatial delineation of geologic formations consists of an image denoiser and a variational auto-encoder (deep learning-based emulator). The former ameliorates the reconstruction noise, yielding well-defined geological structures; its mathematical equivalence with the proximal operator allows the deployment of advanced denoisers (e.g., CNN-based denoiser) that do not correspond to a regularization objective. The latter defines a geology prior that imposes a geological constraint, for example, continuity and shape of geological features, onto the reconstructed image. Here, we conduct a series of numerical experiments dealing with transient two-dimensional flow driven by a pumping well and natural hydraulic head gradient. They demonstrate the CE framework's ability to delineate, both probabilistically and deterministically, complex subsurface environments with sufficient quality.

42 ENGINEERING↗