SEARCH · Engineering Papers

Results for “volumetric deformation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗

NASA NTRS1993

Strength and deformation properties of basaltic lava flows on planetary surfaces

Basaltic rocks are thought to constitute a volumetrically significant rock type on the Moon, Mercury, Mars, and Venus, in addition to the Earth. Spacecraft images of surfaces with known or suspected basaltic composition on these bodies, particularly on Venus, indicate that these rocks have been deformed in the brittle regime to form faults and perhaps dilatant cracks, in addition to folding and more distributed types of deformation. Predictions of brittle fracture or other types of deformation are made by comparing calculated stresses from a tectonic model to some criterion for rock strength. Common strength criteria used in the planetary science literature for near-surface deformation include a Griffith tensile-strength criterion for intact rock, a Mohr envelope for intact basalt, and a brittle strength envelope based on Byerlee's law of rock frictional resistance. However, planetary terrains of basaltic composition consist of much more than just intact basaltic rock. The aggregate basaltic material, termed the 'rock mass,' consists of both the intact rock and the associated fracture, faults, lithologic contacts, and other discontinuous surfaces. A basaltic rock mass is the relevant material for which strength properties must be defined and calculated model stresses must be compared to in order to more accurately predict brittle deformation. For example, the various strengths of a rock mass are less than that of intact material of the same composition. This means that tectonic models which compare stresses to intact failure strengths overestimate the stresses required for fracture and so underestimate the extent and magnitude of brittle deformation predicted in these models. On the other hand, rock mass shear strength can be greater than that predicted from Byerlee's law. The concept of rock mass strength is central to many engineering design studies in which calculated stresses are used to predict brittle fracture, and this experience indicates that brittle strength envelopes which assume properties for intact rock (Griffith parabolas) or sliding along a single, continuous surface (Byerlee's law) inadequately characterize the tensile, compressive, and shear strengths of rock masses. The criterion adopted here to relate stresses to rock mass fracture is based on a Griffith-type curve for tensile normal stress and a concave downward curve for compressive normal stress. It is the only available criterion that explicitly considers the weakening effects of discontinuities within the rock mass on the stress state required for fracture.

Schultz, Richard A.↗

DOE OSTI2024

Asymptotic quasisymmetric high-beta three-dimensional magnetohydrodynamic equilibria near axisymmetry

Quasisymmetry (QS), a hidden symmetry of the magnetic field strength, is known to support nested flux surfaces and provide superior particle confinement in stellarators. In this work, we study the ideal magnetohydrodynamic (MHD) equilibrium and stability of high-beta plasma in a large-aspect-ratio stellarator. In particular, we show that the lowest-order description of a near-axisymmetric equilibrium vastly simplifies the problem of three-dimensional quasisymmetric MHD equilibria, which can be reduced to a standard elliptic Grad–Shafranov equation for the flux function. We show that any large-aspect-ratio tokamak, deformed periodically in the vertical direction, is a stellarator with approximate volumetric QS. We discuss exact analytical solutions and numerical benchmarks. Finally, we discuss the ideal ballooning and interchange stability of some of our equilibrium configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

HFIR SiC-SiC Composite Clad Tube Bowing Test: Pre-Irradiation Characterization

This report describes the pre-irradiation characterization of the SiC-SiC composite clad tube bowing experiment that will be irradiated in the High Flux Isotope Reactor (HFIR). There are concerns that SiCSiC fuel cladding in light water reactors could undergo bowing because of the non-uniform fast neutron flux profiles. This experiment is combined with the SiC-SiC channel box experiment, which has similar concerns regarding control blade movements in boiling water reactors. The HFIR experiment aims to validate the predicted deformation and stresses in the composite clad specimen after the exposure of fast neutron flux gradients. Significant radial fast neutron flux gradients that exist in the permanent reflector of HFIR were thoroughly characterized by using detailed 3D neutronic calculations. The 3D displacement damage dose rate profile and the resulting volumetric swelling in SiC were used as inputs to structural analyses to determine the predicted deformation that will affect the specimens. The tube specimens were characterized by using traditional dimensional inspection and surface profilometry to provide detailed information regarding the pre-irradiation condition as the first step to evaluate the radiation-induced deformation. Furthermore, fine engraving markers were inscribed along all outer surfaces of the specimen and mapped by using a digital microscope and a multidimensional stage. This allows marker spacings to be accurately measured and compared with similar measurements that will be made post-irradiation to provide local radiation-induced strain mapping. The experiment successfully completed two HFIR cycle irradiations, cycles 492 and 493, which finished on July 25, 2021.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Assessing continuum plasticity postulates with grain stress and local strain measurements in triaxially compressed sand

Critical state and continuum plasticity theories have been used in research and engineering practice in soil and rock mechanics for decades. These theories rely on postulated relationships between material stresses and strains. Some classical postulates include coaxiality between stress and strain rates, stress–dilatancy relationships, and kinematic assumptions in shear bands. Although numerical and experimental data have quantified the strains and grain kinematics in such experiments, little data quantifying grain stresses are available. Here, we report the first-known grain stress and local strain measurements in triaxial compression tests on synthetic quartz sands using synchrotron X-ray tomography and 3D X-ray diffraction. We use these data to examine the micromechanics of shear banding, with a focus on coaxiality, stress-dilatancy, and kinematics within bands. Our results indicate the following: 1) elevated deviatoric stress, strain, and stress ratios in shear bands throughout experiments; 2) coaxial principal compressive stresses and strains throughout samples; 3) significant contraction along shear bands; 4) vanishing volumetric strain but nonvanishing stress fluctuations throughout samples at all stages of deformation. Our results provide some of the first-known in situ stress and strain measurements able to aid in critically evaluating postulates employed in continuum plasticity and strain localization theories for sands.

42 ENGINEERING↗

DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

HFIR SiC Bowing Test Ready to Insert

This report describes the successful assembly of a High Flux Isotope Reactor (HFIR) irradiation experiment designed to assess radiation-induced lateral bowing of silicon carbide fiber–reinforced, silicon carbide matrix composite (SiC/SiC) components under a radial fast neutron flux gradient. Excessive bowing of a SiC/SiC channel box in a boiling water reactor could potentially interfere with control blade movements. Similar concerns exist for SiC/SiC fuel cladding in light water reactors. The experiment described herein will provide experimental validation of the structural response of a miniature SiC/SiC channel box and tube specimens with pressurized water reactor diameters during irradiation. The significant radial fast neutron flux gradients that exist in the permanent reflector of HFIR were characterized using detailed three-dimensional neutronic calculations. The three-dimensional displacement damage dose rate profile and the resulting volumetric swelling in SiC were used as inputs to structural analyses that predicted the deformation and stresses in the channel box specimen. The specimens were thoroughly characterized prior to irradiation using traditional dimensional inspection and surface profilometry so that these measurements can later be compared with similar measurements that will be made post-irradiation to determine radiation-induced deformations. Furthermore, fine engraving markers were inscribed along all outer surfaces of the specimen and mapped using a digital microscope and a three-dimensional stage. This technique allowed for accurate measurements of the marker spacings, which can be compared with similar measurements that will be made post-irradiation to provide local radiation-induced strain mapping. The experiment was successfully assembled and is scheduled for insertion during HFIR cycle 492, which is currently scheduled to run from May 25, 2021 to June 18, 2021.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Methodological guidelines on concrete degradation based on predictive models and the release of MOSAIC for industry use

In light-water reactors (LWRs), the concrete biological shield (CBS) designates the concrete structure directly facing the reactor pressure vessel (RPV). The primary function of the CBS is to protect equipment and personnel from the neutron and gamma radiation exiting the RPV. Most CBSs in operation in the United States also provide structural support to the RPV. Based on a long-term structural analysis, the structural integrity of the CBS to transfer the in-service passive load to the foundation, dynamic loading during an earthquake, and thermal loading during a loss-of-coolant accident (LOCA) must be maintained. Based on test reactor data for accelerated conditions, the mechanical properties of irradiated concrete are affected when exposed to neutron fluence higher than ≈ 10 19 n.cm -2 (E > 0.1 MeV). At 80 years of operation, the fluence estimates at the surface of the CBS range between 1×10 19 and 7×10 19 n.cm -2 (E > 0.1 MeV): the fluence value depends mostly on the design. Although fluence attenuates through the CBS by virtue of the shielding properties of concrete, about 5% to 10% of the CBS wall depth is affected by high fluence. The main degradation mechanism is the so-called radiation-induced volumetric expansion (RIVE) caused by the amorphization of aggregate-forming minerals. RIVE causes large deformation of up to 18% in pure quartz, and it also reduces the structural properties of concrete. The structural effects of long-term exposure to fast neutron irradiation in LWRs was studied in the report entitled “Assessment of the Effect of the Irradiation-Induced Degradation on the In-Service Structural Performance of the Concrete Biological Shields” (M3LW-23OR0403013). This report illustrates that (1) irradiation-induced damage expands further than the region subject to a cumulative fluence of 10 19 n.cm -2 (E > 0.1 MeV) to reach the steel reinforcement location, and (2) the irradiation-induced damage region shows no substantial residual bearing capacity and should not be accounted for in the integrity assessment of the CBS.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1990

Seismotectonics of thin- and thick-skinned deformation in the Andean foreland from local network data - Evidence for a seismogenic lower crust

Local network data from San Juan, Argentina, provides new information about crustal seismicity in the Andean foreland above a horizontal segment of the subducted Nazca Plate. Two areas of foreland seismicity are found, one associated with the Sierras Pampeanas basement uplifts, and the other beneath, but not within, the Precordillera foreland fold-thrust belt. The Precordillera seismicity provides direct evidence for basement deformation beneath the sediments of the thrust belt and supports the idea that its eastern part is significantly modified by underlying basement deformation. In both areas, events are concentrated between 15 and 35 km depth and have volumetric, rather than planar, faultlike distributions. The depth distribution is unusually deep for intraplate earthquakes and suggests a brittle-ductile transition near 30-35 km.

Smalley, Robert, Jr.↗

DOE OSTI2021

Radiation effects on stress evolution and dimensional stability of large fusion energy structures

Here we assess the effects of neutron irradiation on the deformation and stress evolution of large-scale fusion energy structures. This is accomplished through non-linear finite element structural analysis of the coupled thermal and mechanical fields at the Beginning-Of-Life (BOL), at 45 dpa, and at 90 dpa. Radiation effects include volumetric swelling and the influence of radiation on the mechanical properties. The system studied here is a large section of a full inboard module of an integrated structure comprising the First Wall and Blanket (FW/B) of a Dual Cooled Lithium-Lead (DCLL) energy conversion unit in the Fusion Nuclear Science Facility (FNSF). The structural material is the ferritic/martensitic steel F82H. We analyze several radiation effects phenomena that can lead to significant impact on the mechanical design and lifetime of the structure. These include volumetric swelling effects on the assembly and disassembly of modules and the impact of spatial gradients in neutron damage and helium generation on deformation, stress evolution, and plastic strain accumulation. We show that the stress state is strongly influenced by the accumulation of swelling strains, radiation hardening and softening, and by spatial gradients in displacement damage and helium gas generation. Several key regions in the FW/B structure are identified where severe plastic strains accumulate and may be potential sites for failure, especially during cyclic reactor operations. We finally show that minimum gaps between Inboard (IB) blanket sectors must be present during assembly so as to avoid excessive stress and strain in FW/B modules critical regions and raising chances of failure. A range of structural displacements are shown, from as low as 4 mm at BOL for free side walls to as high as 46 mm at 90 dpa for fully constrained side walls.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Phase segregation and miscibility of TiO x nanocomposites in Gd-doped ceria solid electrolyte material

Electro-chemo-mechanical (ECM) coupling refers to mechanical deformation due to electrochemically driven compositional change in a solid. An ECM actuator producing micrometre-size displacements and long-term stability at room temperature was recently reported, comprising a 20 mol% Gd-doped ceria (20GDC), a solid electrolyte membrane, placed between two working bodies made of TiO x /20GDC (Ti-GDC) nanocomposites with Ti concentration of 38 mol%. The volumetric changes originating from oxidation or reduction in the local TiO x units are hypothesized to be the origin of mechanical deformation in the ECM actuator. Studying the Ti concentration-dependent structural changes in the Ti-GDC nanocomposites is therefore required for (i) understanding the mechanism of dimensional changes in the ECM actuator and (ii) maximizing the ECM response. Here, the systematic investigation of the local structure of the Ti and Ce ions in Ti-GDC over a broad range of Ti concentrations using synchrotron X-ray absorption spectroscopy and X-ray diffraction is reported. The main finding is that, depending on the Ti concentration, Ti atoms either form a cerium titanate or segregate into a TiO 2 anatase-like phase. The transition region between these two regimes with Ti(IV) concentration between 19% and 57% contained strongly disordered TiO x units dispersed in 20GDC containing Ce(III) and Ce(IV) and hence rich with oxygen vacancies. As a result, this transition region is proposed to be the most advantageous for developing ECM-active materials.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Analysis of time-dependent strain heterogeneity in shales using X-ray microscopy and digital volume correlation

In this paper, we report the results of a 4-dimensional (3 spatial dimensions plus time) in-situ creep experiment on a milli-core shale plug from Barnett formation in Texas. An outstanding challenge in shale mechanics is connecting observable characteristics to predictable mechanical and hydraulic behavior. We employ novel imaging and analysis techniques to statistically link sample mineralogy and structure to observed deformation. In particular, we focus on mechanisms of porosity reduction during transient response to step changes in stress, analogous to the aseismic deformation that occurs after an episode of fracturing. Our main goal in this study is to take the first steps to visualize and capture the small-scale deformation mechanisms involved in creep of shales, using imaging techniques. To do so, we conducted a uniaxial creep experiment inside an Xradia Versa 520 micro-CT system for 12 h at room temperature under a constant uniaxial stress of 25 MPa. We employed a Digital Volume Correlation (DVC) method to measure micro-scale deforspmation of the sample by comparing the resulting volumetric images with a high-resolution scan of the sample obtained prior to the creep experiment. The results indicate that the time-dependent deformation in the compliant clay plus kerogen matrix constituents facilitates rigid grain rotation and compaction of intergranular porosity at the boundary between rigid grains and matrix.

58 GEOSCIENCES↗

DOE OSTI2025

Localized orientation gradients in additively manufactured stainless steel 316H structures

The high solidification rates during additive manufacturing cause highly localized thermal and strain gradients. The effect of these gradients on the evolution of local orientation misorientations within a grain is not well understood. In this study, stainless steel 316H parts were fabricated via laser powder bed fusion using three different energy densities: 43, 71, and 135 J/mm 3 . Electron backscatter diffraction showed that the maximum misorientations of the grains can be up to 25° along the build direction. Misorientation gradients (RM g ) within grains are process-dependent and can change from 0.036°/μm to 0.015°/μm with increased volumetric energy densities. In conclusion, the characterized misorientation gradients are an indication of the level of dislocations and, to an extent, the plastic deformation resulting from the rapid solidification during laser powder bed fusion.

316H↗

NASA NTRS2001

High-Frequency Focused Water-Coupled Ultrasound Used for Three-Dimensional Surface Depression Profiling

To interface with other solids, many surfaces are engineered via methods such as plating, coating, and machining to produce a functional surface ensuring successful end products. In addition, subsurface properties such as hardness, residual stress, deformation, chemical composition, and microstructure are often linked to surface characteristics. Surface topography, therefore, contains the signatures of the surface and possibly links to volumetric properties, and as a result serves as a vital link between surface design, manufacturing, and performance. Hence, surface topography can be used to diagnose, monitor, and control fabrication methods. At the NASA Glenn Research Center, the measurement of surface topography is important in developing high-temperature structural materials and for profiling the surface changes of materials during microgravity combustion experiments. A prior study demonstrated that focused air-coupled ultrasound at 1 MHz could profile surfaces with a 25-m depth resolution and a 400-m lateral resolution over a 1.4-mm depth range. In this work, we address the question of whether higher frequency focused water-coupled ultrasound can improve on these specifications. To this end, we employed 10- and 25-MHz focused ultrasonic transducers in the water-coupled mode. The surface profile results seen in this investigation for 25-MHz water-coupled ultrasound, in comparison to those for 1-MHz air-coupled ultrasound, represent an 8 times improvement in depth resolution (3 vs. 25 m seen in practice), an improvement of at least 2 times in lateral resolution (180 vs. 400 m calculated and observed in practice), and an improvement in vertical depth range of 4 times (calculated).

Roth, Don J.↗

NASA NTRS2012

Demonstration of Repeat-Pass POLINSAR Using UAVSAR: The RMOG Model

In this paper we show our first POLINSAR results using the Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) developed by the Jet Propulsion Laboratory (JPL). UAVSAR is a L-band repeat-pass polarimetric and interferometric system designed for measuring vegetation structure and monitoring crustal deformations. In order to extract canopy height from POLINSAR data and account for temporal decorrelation, we formulate a physical model of the temporal-volumetric coherence, random motion over ground (RMOG) model. Canopy height extracted from single-baseline UAVSAR data using the RMOG model is shown to be in agreement with canopy height measured by the Land, Vegetation, and Ice Sensor (LVIS) lidar.

polarimetry↗

DOE OSTI2026

Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems

The Material Point Method (MPM) has been increasingly used to simulate large strain deformations. Linear interpolation functions are commonly used to perform the spatial integration. It is well-known that the discontinuities in the interpolation function derivatives induce shock-like artifacts known as ‘cell-crossing’ error. These errors compound with volumetric locking errors when used with hydromechanical formulations for porous media, where different velocity fields are used for each phase. The capabilities of higher-order MPM frameworks have not been explored for real-scale geotechnical problems. As such, this paper aims to assess, validate, and further discuss a higher-order B-spline MPM (BS-MPM) framework. First, the BS-MPM framework is verified against the large-strain oedometer consolidation problem. Second, the framework is validated against a real-scale slope failure experiment triggered by pore water pressure recharge. Landslide features that are captured using the higher-order framework are specifically highlighted, and results (e.g., pore water pressure and deformation) are validated with field measurements. A generally convergent numerical solution is observed when using cubic interpolation functions. Third, a footing penetration problem is simulated using the multi-patch BS-MPM. Trends are examined with respect to penetration velocity and variation in hydraulic conductivity. The BS-MPM framework ultimately presents a stabilized numerical solution that captures plausible hydromechanical interaction trends important in geotechnical engineering applications.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Fourier-based three-dimensional multistage transformer for aberration correction in multicellular specimens

High-resolution tissue imaging is often compromised by sample-induced optical aberrations that degrade resolution and contrast. Although wavefront sensor-based adaptive optics (AO) can measure these aberrations, such hardware solutions are typically complex, expensive to implement and slow when serially mapping spatially varying aberrations across large fields of view. Here we introduce AOViFT (adaptive optical vision Fourier transformer)—a machine learning-based aberration sensing framework built around a three-dimensional multistage vision transformer that operates on Fourier domain embeddings. AOViFT infers aberrations and restores diffraction-limited performance in puncta-labeled specimens with substantially reduced computational cost, training time and memory footprint compared to conventional architectures or real-space networks. We validated AOViFT on live gene-edited zebrafish embryos, demonstrating its ability to correct spatially varying aberrations using either a deformable mirror or postacquisition deconvolution. By eliminating the need for the guide star and wavefront sensing hardware and simplifying the experimental workflow, AOViFT lowers technical barriers for high-resolution volumetric microscopy across diverse biological samples.

Alshaabi, Thayer [Howard Hughes Medical Institute,↗

NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗