Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large displacement”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Background Information on ARPA-E's REUSE Program (White Paper)

Under the REUSE Topic, ARPA-E seeks to fund the development of technologies to convert high-energy materials currently going to landfills to a high-energy content liquid product capable of displacing energy imports used for fuel or chemical production. The high-energy materials include plastics (#1-7 polymers, rubber, and composites) and paper. As discussed below, we estimate 30-34 MM ton plastic, 2-6 MM ton rubber, 18 MM ton of paper, and up to 0.1 MM ton composites are potentially available annually for this purpose. These numbers may be conservative based on changes in the plastic and paper export markets. ARPA-E anticipates deployment of multiple low-cost, simple, flexible, small-scale (100-500 ton per day) regional facilities using modular plants. This scale is consistent with the sources for high-energy materials, which include ~300 Material Recovery Facilities and industrial waste sources. The assumption is that such facilities can be more economical than the paradigm of large-scale facilities making purity products, due to cost for transporting and aggregating waste and the high operating costs (OPEX) and capital cost (CAPEX) for product purification. This document gives a brief technical review for multiple potential process technologies. The review is not intended to be comprehensive or limiting, only to provide an introduction to potential Applicants.

10 SYNTHETIC FUELS↗
DOE OSTI2020

White Paper: Background Information on ARPA-E's Reuse Program

This document provides supplemental information for ARPA-E’s exploratory research program, “Recycle Underutilized Solids to Energy” (REUSE). The goal is to provide additional technical information to prospective Applicants. Further information is available in a blog interview and webinar at https://arpa-e.energy.gov/?q=news-item/trash-treasure-reuse-creates-feedstock-plastic-waste. Award information, submission requirements, evaluation criteria, and other applicable information is provided in Funding Opportunity Announcements DE-FOA-0001953 and DE-FOA-0001954 (SBIR/STTR). REUSE-specific requirements are provided in Topic K of those FOAs, and is available at https://arpa-e-foa.energy.gov/Default.aspx?Search=0001953&SearchType=#FoaIde8647d89-1cac-4b58-8622-1b04de8958c4. Under the REUSE Topic, ARPA-E seeks to fund the development of technologies to convert high-energy materials currently going to landfills to a high-energy content liquid product capable of displacing energy imports used for fuel or chemical production. The high-energy materials include plastics (#1-7 polymers, rubber, and composites) and paper. As discussed below, we estimate 30-34 MM ton plastic, 2-6 MM ton rubber, 18 MM ton of paper, and up to 0.1 MM ton composites are potentially available annually for this purpose. These numbers may be conservative based on changes in the plastic and paper export markets. ARPA-E anticipates deployment of multiple low-cost, simple, flexible, small-scale (100-500 ton per day) regional facilities using modular plants. This scale is consistent with the sources for high-energy materials, which include ~300 Material Recovery Facilities and industrial waste sources. The assumption is that such facilities can be more economical than the paradigm of large-scale facilities making purity products, due to cost for transporting and aggregating waste and the high operating costs (OPEX) and capital cost (CAPEX) for product purification. This document gives a brief technical review for multiple potential process technologies. The review is not intended to be comprehensive or limiting, only to provide an introduction to potential Applicants.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Transformer coupled toroidal wave-heated remote plasma sources operating in Ar/NF 3 mixtures

Remote plasmas are used in semiconductor device manufacturing as sources of radicals for chamber cleaning and isotropic etching. In these applications, large fluxes of neutral radicals (e.g. F, O, Cl, H) are desired with there being negligible fluxes of potentially damaging ions and photons. One remote plasma source (RPS) design employs toroidal, transformer coupling using ferrite cores to dissociate high flows of moderately high pressure (up to several Torr) electronegative gases. In this paper, results are discussed from a computational investigation of moderate pressure, toroidal transformer coupled RPS sustained in Ar and Ar/NF 3 mixtures. Operation of the RPS in 1 Torr (133 Pa) of argon with a power of 1.0 kW at 0.5 MHz and a single core produces a continuous toroidal plasma loop with current continuity being maintained dominantly by conduction current. Operation with dual cores introduces azimuthal asymmetries with local maxima in plasma density. Current continuity is maintained by a mix of conduction and displacement current. Operation in NF 3 for the same conditions produces essentially complete NF 3 dissociation. Electron depletion as a result of dissociative attachment of NF 3 and NF x fragments significantly alters the discharge topology, confining the electron density to the downstream portion of the source where the NFx density has been lowered by this dissociation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Standardized Analysis Process Using Digital Image Correlation to Calculate In Situ Cladding Strain from Modified Burst Tests for Fuel Performance Code Validation

Historical data collection on nuclear fuel cladding materials has focused on generating a statistically significant amount of data to assess the material and its failure behavior. Furthermore, data generated to support material model and failure criteria development were previously posttest evaluations, so a large number of tests was required to gain new understanding. A way to expedite this process is to develop techniques capable of generating large, high-fidelity data sets from a single test with lower uncertainty or quantified uncertainty. One such example of this approach is Oak Ridge National Laboratory’s use of modified burst tests (MBTs) to analyze the mechanical behavior and failure conditions of cladding during a simulated reactivity-initiated accident (RIA). Each test incorporates digital image correlation (DIC) analysis techniques that are used to assess the accumulated strain in situ, as well as eventual cladding failure. This work has been fruitful in defining strain-to-failure conditions for materials like silicon carbide (SiC) fiber–reinforced/SiC matrix composite tubes (SiC/SiC), iron-chromium-aluminum (FeCrAl) alloy tubes, and chromium-coated Zircaloy-4 tubes. However, there are numerous DIC software available, including open-source and proprietary software. The different DIC software use various algorithms to process images and calculate displacement values. Using these different software and algorithms can lead to varying results, and perhaps larger-than-expected uncertainties. In the present study, previously published MBT data encompassing a variety of test conditions were reanalyzed with two different DIC software to assess the variance in the calculated strain results. The data consisted of SiC/SiC, FeCrAl, and chromium-coated Zircaloy-4 tubes. Plots of the calculated strains during the transient revealed good agreement between the two DIC software. The average root-mean-square errors between the two software was 0.20% strain, which is slightly larger than a previously reported error value for these tests. In conclusion, this variance in results is low enough that this analysis method can be used for code validation.

Reactivity-initiated accident↗
DOE OSTI2020

Robust paths to net greenhouse gas mitigation and negative emissions via advanced biofuels

Biofuel and bioenergy systems are integral to most climate stabilization scenarios for displacement of transport sector fossil fuel use and for producing negative emissions via carbon capture and storage (CCS). However, the net greenhouse gas mitigation benefit of such pathways is controversial due to concerns around ecosystem carbon losses from land use change and foregone sequestration benefits from alternative land uses. Here, we couple bottom-up ecosystem simulation with models of cellulosic biofuel production and CCS in order to track ecosystem and supply chain carbon flows for current and future biofuel systems, with comparison to competing land-based biological mitigation schemes. Analyzing three contrasting US case study sites, we show that on land transitioning out of crops or pasture, switchgrass cultivation for cellulosic ethanol production has per-hectare mitigation potential comparable to reforestation and severalfold greater than grassland restoration. In contrast, harvesting and converting existing secondary forest at those sites incurs large initial carbon debt requiring long payback periods. We also highlight how plausible future improvements in energy crop yields and biorefining technology together with CCS would achieve mitigation potential 4 and 15 times greater than forest and grassland restoration, respectively. Finally, we show that recent estimates of induced land use change are small relative to the opportunities for improving system performance that we quantify here. While climate and other ecosystem service benefits cannot be taken for granted from cellulosic biofuel deployment, our scenarios illustrate how conventional and carbon-negative biofuel systems could make a near-term, robust, and distinctive contribution to the climate challenge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Ordered hexagonal columnar liquid crystalline self-assembly of mesogen-free sulfonylated side-chain chiral polyethers and their high dielectric property

To achieve large spontaneous polarization for ferroelectric liquid crystalline polymers, we have designed a series of mesogen-free isotactic polyethers with sulfonylated side chains. Different from prior studies using linear alkyl side chains, in this work a methyl branch was implemented in the alkyl side chains of the chiral polyethers. As a result of the continuous twist from the methyl branches and strong dipole-dipole interactions among neighboring sulfonyl groups, a helical chain conformation with tilted side chains was achieved. When the branched side chains were short (e.g., C5) and double chirality existed in both the main chain and side chains, a crystalline structure was obtained. When the branched side chains were long (e.g., ≥ C8) and single chirality existed in the main chain, an ordered hexagonal columnar liquid crystalline phase was induced. Their high-field dielectric property was studied by electric displacement-electric field (D-E) loop tests; however, no ferroelectric behavior was observed up to 200 MV/m. It was considered that the dipole-dipole interaction among neighboring sulfonyl groups in the helical chains was too strong to allow easy ferroelectric flipping.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Subseasonal Variability of Elevated Dust Concentrations Over South Florida

Dust mass concentrations have been measured daily at Miami, Florida, in the summer months, and biweekly throughout the remaining months, since 1974. The 43-year record of dust mass concentrations indicate large daily, seasonal, and interannual variations, with most of the dust arriving within 5–8 episodes each summer. On average, dust arrives to Miami, Florida, 10 days after emission from North Africa, with measured concentrations depending on characteristics of the lower free-tropospheric winds due to the vast travel distance. Daily dust mass concentrations from July and August, the months that contribute the most to the annual mean, are used to characterize the synoptic conditions most favorable for dust transport. Two key regions are linked with the highest daily dust mass concentrations above Miami: (i) easterly winds, averaged over 850–500 mb, over the Tropical West Atlantic [15–25°N, 45–80°W], and (ii) southerly winds, similarly averaged, over the Florida Peninsula [20–30°N,75–80°W]. Winds within these two regions are enhanced when the North Atlantic subtropical high is displaced south and zonally elongated, relocating the western edge over Florida. A dust-transport-efficiency index, based on the maximum potential for dust to arrive above Miami with limited loss to deposition or mixing, identifies high-dust loading cases on the subseasonal scale. Monthly dust-transport-efficiency values agree well with the monthly dust trends over the 43-year time span. While seasonal dust loadings have been decreasing over Florida in the past decade, the transport efficiency has been increasing, possibly due to trends in the North Atlantic subtropical high.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Learning effective physical laws for generating cosmological hydrodynamics with Lagrangian deep learning

The goal of generative models is to learn the intricate relations between the data to create new simulated data, but current approaches fail in very high dimensions. When the true data-generating process is based on physical processes, these impose symmetries and constraints, and the generative model can be created by learning an effective description of the underlying physics, which enables scaling of the generative model to very high dimensions. In this work, we propose Lagrangian deep learning (LDL) for this purpose, applying it to learn outputs of cosmological hydrodynamical simulations. The model uses layers of Lagrangian displacements of particles describing the observables to learn the effective physical laws. The displacements are modeled as the gradient of an effective potential, which explicitly satisfies the translational and rotational invariance. The total number of learned parameters is only of order 10, and they can be viewed as effective theory parameters. Further, we combine N-body solver fast particle mesh (FastPM) with LDL and apply it to a wide range of cosmological outputs, from the dark matter to the stellar maps, gas density, and temperature. The computational cost of LDL is nearly four orders of magnitude lower than that of the full hydrodynamical simulations, yet it outperforms them at the same resolution. We achieve this with only of order 10 layers from the initial conditions to the final output, in contrast to typical cosmological simulations with thousands of time steps. This opens up the possibility of analyzing cosmological observations entirely within this framework, without the need for large dark-matter simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Combining different 3-D global and regional seismic wave propagation solvers towards box tomography in the deep Earth

SUMMARY In previous publications, we presented a general framework, which we called ‘box tomography’, that allows the coupling of any two different numerical seismic wave propagation solvers, respectively outside and inside a target region, or ‘box’. The goal of such hybrid wavefield computations is to reduce the cost of computations in the context of full-waveform inversion for structure within the target region, when sources and/or receivers are located at large distances from the box. Previously, we had demonstrated this approach with sources and receivers outside the target region in a 2-D acoustic spherical earth model, and demonstrated and applied this methodology in the 3-D spherical elastic Earth in a continental scale inversion in which all stations were inside the target region. Here we extend the implementation of the approach to the case of a 3-D global elastic earth model in the case where both sources and stations are outside the box. We couple a global 3-D solver, SPECFEM3D_GLOBE, for the computation of the wavefield and Green’s functions in a reference 3-D model, with a regional 3-D solver, RegSEM, for the computation of the wavefield within the box, by means of time-reversal mirrors. We briefly review key theoretical aspects, showing in particular how only the displacement is needed to be stored at the boundary of the box. We provide details of the practical implementation, including the geometrical design of the mirrors, how we deal with different sizes of meshes in the two solvers, and how we address memory-saving through the use of B-spline compression of the recorded wavefield on the mirror. The proposed approach is numerically efficient but also versatile, since adapting it to other solvers is straightforward and does not require any changes in the solver codes themselves, as long as the displacement can be recovered at any point in time and space. We present benchmarks of the hybrid computations against direct computations of the wavefield between a source and an array of stations in a realistic geometry centred in the Yellowstone region, with and without a hypothetical plume within the ‘box’, and with a 1-D or a 3-D background model, down to a period of 20 s. The ultimate goal of this development is for applications in the context of imaging of remote target regions in the deep mantle, such as, for example, Ultra Low Velocity Zones.

Geochemistry & Geophysics↗
DOE OSTI2020

In Situ Proteolysis Condition-Induced Crystallization of the XcpVWX Complex in Different Lattices

Although prevalent in the determination of protein structures; crystallography always has the bottleneck of obtaining high-quality protein crystals for characterizing a wide range of proteins; especially large protein complexes. Stable fragments or domains of proteins are more readily to crystallize; which prompts the use of in situ proteolysis to remove flexible or unstable structures for improving crystallization and crystal quality. In this work; we investigated the effects of in situ proteolysis by chymotrypsin on the crystallization of the XcpVWX complex from the Type II secretion system of Pseudomonas aeruginosa. Different proteolysis conditions were found to result in two distinct lattices in the same crystallization solution. With a shorter chymotrypsin digestion at a lower concentration; the crystals exhibited a P3 hexagonal lattice that accommodates three complex molecules in one asymmetric unit. By contrast; a longer digestion with chymotrypsin of a 10-fold higher concentration facilitated the formation of a compact P212121 orthorhombic lattice with only one complex molecule in each asymmetric unit. The molecules in the hexagonal lattice have shown high atomic displacement parameter values compared with the ones in the orthorhombic lattice. Taken together; our results clearly demonstrate that different proteolysis conditions can result in the generation of distinct lattices in the same crystallization solution; which can be exploited in order to obtain different crystal forms of a better quality

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

Effects of dilution in ionic liquid supercapacitors

Room-temperature ionic liquids (RTILs) are synthetic electrolytes that have a large electrochemical stability window, making them attractive candidates for electric double-layer capacitor (EDLC) applications. Due to their high viscosities and low ionic conductivities, RTILs are often diluted with organic solvent for practical use. Here, we study the effects of dilution on the performance of RTIL EDLCs using a simple mean-field model. We find that dilution diminishes the unfavorable hysteresis that results from a spontaneous surface charge separation (SSCS). As a result, the RTIL concentration can be used to modulate the proximity to the SSCS transition, and maximize capacitance. The interplay between the concentration and the correlation strength gives rise to complex zero-potential phase behavior, including a tricritical point and a λ-line, very similar to the Blume–Capel dilute Ising model. Additionally, electrodes that are solvophilic aid in the prevention of SSCS by drawing solvent molecules to the electrode and displacing ions. Solvophilic electrodes give rise to a phase transition at finite potential where the surface charge rapidly increases with a small increase in potential, leading to a substantial increase in capacitance and energy storage.

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

Discrete-Direct Model Calibration and Uncertainty Propagation Method Confirmed on Multi-Parameter Plasticity Model Calibrated to Sparse Random Field Data

A discrete direct (DD) model calibration and uncertainty propagation approach is explained and demonstrated on a 4-parameter Johnson-Cook (J-C) strain-rate dependent material strength model for an aluminum alloy. The methodology's performance is characterized in many trials involving four random realizations of strain-rate dependent material-test data curves per trial, drawn from a large synthetic population. The J-C model is calibrated to particular combinations of the data curves to obtain calibration parameter sets which are then propagated to “Can Crush” structural model predictions to produce samples of predicted response variability. These are processed with appropriate sparse-sample uncertainty quantification (UQ) methods to estimate various statistics of response with an appropriate level of conservatism. This is tested on 16 output quantities (von Mises stresses and equivalent plastic strains) and it is shown that important statistics of the true variabilities of the 16 quantities are bounded with a high success rate that is reasonably predictable and controllable. The DD approach has several advantages over other calibration-UQ approaches like Bayesian inference for capturing and utilizing the information obtained from typically small numbers of replicate experiments in model calibration situations—especially when sparse replicate functional data are involved like force–displacement curves from material tests. The DD methodology is straightforward and efficient for calibration and propagation problems involving aleatory and epistemic uncertainties in calibration experiments, models, and procedures.

42 ENGINEERING↗
DOE OSTI2026

Mesogenesis through the ephemeral dark decay of Beauty

Mesogenesis provides a path for generating the baryon asymmetry of the Universe, using only the CP violation furnished by the Standard Model in the decay of B mesons. While this is an intriguing possibility, it is largely constrained by the data on B meson branching fractions into baryons and missing energy carried into the dark sector. We point out that it is possible to make this branching fraction dominant only in the early Universe, through an ultralight scalar coupled to the dark sector and the Standard Model leptons. A scenario is examined where the thermal density of muons in the early Universe temporarily lowers the mass of a dark fermion, allowing for efficient B meson decays. This `dark’ decay channel is shut off later when the muon number density falls, making the scenario compatible with flavor data. Our model can be consistent with the LHC constraints on color-charged heavy bosons required to implement Mesogenesis; such states may be discovered in the future runs as their masses cannot be far above the current bounds. We also outline other possible signals, which can arise in future displaced vertex searches, long range force searches, and observations of neutron star binary mergers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Stereo-DIC Challenge 1.0 – Rigid Body Motion of a Complex Shape

Background Stereo-DIC is a widely used optical measurement technique that provides a dense full-field 3D measurement of the shape, displacement, and strain of a solid sample. When compared with 2D-DIC, Stereo-DIC provides greater flexibility and expands its use beyond flat, planar specimens. Furthermore, the widespread availability of commercial systems has led to the adoption of the technique throughout industry, academia, and government research labs. Objective Even though some research has been done to understand the effects of different experimental and stereo-DIC parameters, no reference is available to benchmark and compare the performance of current stereo-DIC algorithms to each other. Methods This paper provides the description and analysis of a carefully controlled 3D experiment and associated images used to compare the results from five subset based DIC software packages. Both the images and analysis codes used in this paper to compare the results are described here and are available for download and use for continued research. Results We show that over a very large range of motion, the 3D errors are very small, less than 80μm over a travel of ±20 mm out-of-plane and ±20 mm in-plane. While all codes performed similarly, there are important differences noted in the paper. Conclusion The image sets and results comparison software are hosted by the International DIC Society (www.iDICs.org) and are freely available for download and analysis for comparison with results in this paper. Furthermore, it is hoped that this set of images can be used for future research in improving stereo-DIC by future authors.

Algorithms comparison↗
DOE OSTI2024

Oxygen Storage Incorporated Into Net Power and the Allam–Fetvedt Oxy-Fuel sCO2 Power Cycle—Techno-Economic Analysis

Abstract With the planned future reliance on variable renewable energy, the ability to store energy for prolonged time periods will be required to reduce the disruption of market fluctuations. This paper presents a method to analyze a hybrid liquid-oxygen (LOx) storage/direct-fired supercritical carbon dioxide (sCO2) power cycle and optimize the economic performance over a diverse range of scenarios. The system utilizes a modified version of the NET Power process to produce energy when energy demand exceeds the supply while displacing much of the cost of the air separation unit (ASU) energy requirements through cryogenic storage of oxygen. The model uses marginal cost of energy data to determine the optimal times to charge and discharge the system over a given scenario. The model then applies ramp rates and other time-dependent factors to generate an economic model for the system without storage considerations. The size of the storage system is then applied to create a realistic model of the plant operation. From the real plant operation model, the amount of energy charged and discharged, the capital expenditures (CAPEX) of each system, energy costs and revenue and other parameters can be calculated. The economic parameters are then combined to calculate the net present value (NPV) of the system for the given scenario. The model was then run through the SMPSO genetic algorithm in Python for a variety of geographic regions and large-scale scenarios (high solar penetration) to maximize the NPV based on multiple parameters for each subsystem. The LOx storage requirements will also be discussed.

Engineering↗
DOE OSTI2024

Testing Protocol Development for the Fracture Toughness of Parts Built with Big Area Additive Manufacturing

The mechanical testing of additively manufactured parts has largely relied on the existing standards developed for traditional manufacturing. While this approach leverages the investment made in current standards development, it inaccurately assumes that the mechanical response of additive manufacturing (AM) parts is identical to that of parts manufactured through traditional processes. When considering thermoplastic, material extrusion AM, the differences in response can be attributed to an AM part’s inherent inhomogeneity caused by porosity, interlayer zones, and surface texture. Additionally, the interlayer bonding of parts printed with large-scale AM is difficult to adequately assess, as much testing is performed such that stress is distributed across many layer interfaces; therefore, the lack of AM-specific standards to assess interlayer bonding is a significant research gap. To quantify interlayer bonding via fracture toughness, double cantilever beam (DCB) testing has been used for some AM materials, and DCB has been generally used for a variety of materials including metal, wood, and laminates. Mode I DCB testing was performed on thermoplastic matrix composites printed with Big Area Additive Manufacturing (BAAM). Of particular interest was the notch shape and deflection speed during testing. The results examine the differences when using two notch types and three deflection speeds. The testing method introduced by the following paper differentiates itself from the ones described in the standards used by modernizing the methodology. This was conducted with the introduction of Digital Image Correlation (DIC) to gather displacement and load data simultaneously without human intervention.

Polymer Science↗
DOE OSTI2021

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty

Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. We develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)-based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one-dimensional problem and demonstrate it numerically for a three-dimensional aquifer model. Here we design a generalized preconditioned Crank-Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low-rank-based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two-dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.

54 ENVIRONMENTAL SCIENCES↗