Results for “vulcan” — Page 9

NASA NTRS2019

The Effect of Turbulence Modeling on the Mixing Characteristics of Several Fuel Injectors at Hypervelocity Flow Conditions

CFD analysis is presented on the effects of turbulence modeling choices on the mixing characteristics and performance of three fuel injectors at hypervelocity flow conditions. The analyses were carried out with the VULCAN-CFD solver using Reynolds-Averaged Simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The three injectors are the baseline configurations used in the experiments and represent three categories of injectors typically considered individually or in combination for fueling high-speed propulsive devices. The current work discusses the impact of the turbulence model and the turbulent Schmidt number on the mixing flow field behavior and the mixing performance as described by the one-dimensional values of the Mach number, total pressure recovery, and the mixing efficiency. Because planar laser induced fluorescence (PLIF) images are available from the EIMP experiments, the sensitivity of the synthetic LIF signal to turbulence modeling choices is also examined to determine whether PLIF can be extended beyond its intended qualitative visualization purpose and used to guide CFD turbulence model and parameter selections. It is found that the mixing performance, as quantified using mixing efficiency, exhibits a strong sensitivity to both turbulence model choice and turbulent Schmidt number value. However, the synthetic LIF signal only demonstrates a modest level of sensitivity, which suggests that PLIF is of limited use for guiding CFD turbulence model and parameter selections.

Drozda, Tomasz G.↗

NASA NTRS2018

Comparison of Mixing Characteristics for Several Fuel Injectors at Mach 8, 12, and 15 Hypervelocity Flow Conditions

CFD analysis is presented of the mixing characteristics and performance of three fuel injectors at flight Mach numbers of 8, 12, and 15. The Reynolds-averaged simulations (RAS) were carried out using the VULCAN-CFD solver. The high Mach number flow conditions match those of the experiments conducted as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The EIMP aims to investigate scramjet fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies relevant to flight Mach numbers greater than 8. The injectors include a fuel placement device, a strut, and a fluidic vortical mixer, a ramp. These fuel injectors accomplish the necessary task of distributing and mixing fuel into the supersonic cross-flow, albeit via different strategies. For comparison, a flush-wall injector is also included. This type of injector generally represents the simplest method of introducing fuel into a scramjet combustor. The three injectors represent the baseline configurations of the EIMP experiments. The mixing parameters of interest, such as mixing efficiency and total pressure recovery, are computed from the RAS and compared for the three flight conditions and injector configurations. In addition to mixing efficiency and total pressure recovery, the combustion efficiency and thrust potential are also computed for the reacting simulations. Plotting the total pressure recovery and thrust potential as a function of mixing efficiency provides added insight into critical aspects of combustor performance as the flight condition and injector type are varied.

Shenoy, Rajiv R.↗

NASA NTRS2019

Efficient and Robust Weighted Least-Squares Cell-Average Gradient Construction Methods for the Simulation of Scramjet Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. The construction of cell-average gradients using a weighted linear least-squares method and the use of these gradients in the construction of the inviscid fluxes is the focus of this paper. A comparison of least-squares stencil construction methodologies is presented and approaches designed to minimize the number of cells used to augment/stabilize the least-squares stencil while preserving accuracy are explored. Due to our interest in hypersonic flow, a robust multidimensional cell-average gradient limiter procedure that is consistent with the stencil used to construct the cellaverage gradients is described. Canonical problems are computed to illustrate the challenges and investigate the accuracy, robustness and convergence behavior of the cell-average gradient methods on unstructured cell-centered finite-volume grids. Finally, thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet engine flowpath is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the preferred gradient method for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗

NASA NTRSDate not supplied

Automation of the Uncertainty Quantification Process Based on Probability Boxes with DAKOTA

To date, while the use of CFD is prevalent, very few efforts have been undertaken that truly attempt to document all (or even most) of the sources of uncertainty in the simulations. Instead, the current state-of-the-art relies heavily on the experience of the CFD practitioner to estimate the uncertainty associated with their simulations through simple sensitivity studies or subject matter expertise. This practice will have to be replaced with a formal uncertainty quantification (UQ) process if CFD is to play an expanded role in the design research and engineering community, test and evaluation community, and ultimately certification for flight. This is especially true for hypersonic air-breathing propulsion systems due to the environment, scale, and duration limitations of ground test facilities. Accounting for uncertainties in a formal manner is a tedious process. Moreover, the typical CFD practitioner is not likely to be familiar with formal UQ methods. Hence, a major obstacle that has prevented the adoption of UQ methods for engineering design and development work is the lack of a tool set to automate most (if not all) of the UQ workflow. Towards this end, the SANDIA package DAKOTA (which has been developed to drive both UQ and optimization processes) will be tightly wrapped around the VULCAN-CFD code to automate the uncertainty quantification process. The automated process will be applied to an isolator turbulence model validation exercise that has previously been documented using a manual approach to the UQ process. Hence, the focus of this paper will be documenting the level to which automation can hide the UQ process details from the CFD practitioner rather than the UQ method itself.

CFD↗

NASA NTRSDate not supplied

CFD Validation Study of a Hypersonic Cone-Slice-Flap Configuration

Model validation is the process of determining the degree of accuracy between the real world and the model. The result of model validation can either be used to improve the model through calibration or quantify the model-form uncertainty. This work focuses on providing the model-form uncertainty through an area metric for a hypersonic cone-slice-flap geometry configuration given uncertainty in both the simulation and experimental data. A procedure will be presented that can give an accurate representation of the model-form uncertainty using a small number of high-fidelity runs. The work will also assess the impact between using different high-fidelity solvers and different turbulence models.The goal of performing this comparison is to provide a quantifiable measurement for the accuracy of each solver and turbulence model for this type of design. A low-fidelity analysis will also be performed to get a model-form uncertainty for the lower fidelity analysis. The two high-fidelty CFD solvers used are VULCAN-CFD and FUN3D and the low-fidelity results will come from Cart3D. The experimental data came from the 20-inch Mach 6Tunnel located at NASA Langley Research Center.

Laura M. White↗

NASA NTRSDate not supplied

Comparisons Between NO PLIF Imaging and CFD Simulations of Mixing Flowfields for High-Speed Fuel Injectors

The current work compares experimentally and computationally obtained nitric oxide (NO) planar laser-induced fluorescence (PLIF) images of the mixing flowfields for three types of high-speed fuel injectors: a strut, a ramp, and a rectangular flushwall. These injection devices, which exhibited promising mixing performance at lower flight Mach numbers, are currently being studied as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The EIMP aims to investigate scramjet fuel injection and mixing physics, and improve the understanding of underlying physical processes relevant to flight Mach numbers greater than eight. In the experiments conducted in the NASA Langley Arc-Heated Scramjet Test Facility (AHSTF), the injectors are placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number air flow at the entrance of a scramjet combustor. Helium is used as an inert substitute for hydrogen fuel. Both schlieren and PLIF techniques are applied to obtain mixing flowfield flow visualizations. The experimental PLIF is obtained by using a UV laser sheet to interrogate a plane of the flow by exciting fluorescence from the NO molecules that are present in the AHSTF air. Consequently, the absence of signal in the resulting PLIF images is an indication of pure helium (fuel). The computational PLIF is obtained by applying a fluorescence model for NO to the results of the Reynolds-averaged simulations (RAS) of the mixing flowfield carried out using the VULCAN-CFD solver. This approach is required because the PLIF signal is a nonlinear function of not only NO concentration, but also pressure, temperature, and the flow velocity. This complexity allows additional flow features to be identified and compared with those obtained from the computational fluid dynamics (CFD) simulations, however, such comparisons are only semiquantitative. Three-dimensional image reconstruction, similar to that used in magnetic resonance imaging, is also used to obtain images in the streamwise and spanwise planes from select cross-stream PLIF plane data. Synthetic schlieren is also computed from the RAS data. Good agreement between the experimental and computational results provides increased confidence in the CFD simulations for investigations of injector performance.

PLIF↗

NASA NTRS2020

An Initial Concept for a JUMP Mating Mechanism

Lunar landers currently under development for crew and cargo missions to the Moon are limited by the available boosters to launch them from Earth. Most landers are being developed for commercial launch vehicles (CLVs) in the Falcon Heavy / Vulcan / New Glenn / Atlas V class, limiting their performance with many lunar surface studies consequently restricted to 5-10 metric ton payloads. The JUMP Lander (Joinable Undercarriage to Maximize Payload) attempts to overcome this limitation by performing an in-space mating, linking multiple landers together. This research provides an initial proof of concept demonstration that two or more CLV-launched lunar landers can be mated in-space to produce a larger lander with greater payload capacity. Multiple rockets have been connected hundreds if not thousands of times in spaceflight history. An example is the space shuttle. The External Tank-Orbiter interface is an example where the 100-ton Orbiter is bolted to the side of the External Tank. This example serves as a concept metaphor for JUMP where a system of bolts is used as the mating system. This paper performs a high-level sizing analysis of the bolt, including the structural connection to each lander and the physical dimensions of the bolt, given the projected lander and payload masses. It sizes the necessary motors to perform the in-space mating and describes the motor to bolt integration. It describes the mounting and translation for the bolt-motor assembly, including structural supports, physical stops, and range of travel. It discusses active-active and active-passive implementations for the mating system and post-mating options to jettison no-longer-needed components of the assembly, including a discussion of benefits of jettisoning or not jettisoning. Finally, forward work for the mating mechanism and overall JUMP Lander concept is described.

Lunar Lander↗

NASA NTRSDate not supplied

Adaptation of Metal Additive Manufacturing Processes for the International Space Station

The In-Space Manufacturing (ISM) project at NASA Marshall Space Flight Center, in a partnership with the company, Made in Space, has previously investigated 3D printing of polymer materials on-orbit. In recent years, the project has begun exploring the potential for metal additive manufacturing (AM) on future space missions to reduce logistics and enable point-of-use manufacturing for sparing and repair. This presentation will provide an overview of constraints for demonstrating a manufacturing process on the International Space Station (ISS) as well as trades of available metal AM processes and their potential for in-space use. There are currently two processes in development as payloads for an ISS technology demonstration: wire+arc additive manufacturing (the Vulcan payload from Made in Space, Inc.) and bound metal deposition (the Fabrication Laboratory from Techshot, Inc). An update on both of these systems, results to date, and future development efforts will be presented. Relevant modeling work to evaluate operation of certain aspects of the processes in a microgravity environment will also be summarized.

manufacturing↗

NASA NTRSDate not supplied

CFD Validation Study of a Hypersonic Cone-Slice-Flap Configuration

Model validation is the process of determining the degree of accuracy between the real world and the model. The result of model validation can be used to either improve the model through calibration or quantify the model-form uncertainty. This work focuses on estimating the model-form uncertainty through an area metric for a hypersonic cone-slice-flap geometry configuration given the uncertainty in both the simulation and experimental data. A procedure that can give an accurate representation of the model-form uncertainty using a small number of high-fidelity runs is outlined. The work also assesses the impact of using different high-fidelity solvers and different turbulence models. The goal of performing this comparison is to provide a quantifiable measurement of the accuracy of each solver and turbulence model for this type of design. A low-fidelity analysis is also performed to get a model-form uncertainty for this type of analysis. The two high-fidelity CFD solvers used here are VULCAN-CFD and FUN3D, and the low-fidelity results comes from Cart3D.The experimental data comes from the 20-inch Mach 6 Tunnel located at NASA Langley Research Center. The work here shows that the models used tend to under-predict the drag and moment aerodynamic coefficients for this type of design leading to a model-form uncertainty estimate which can account for up to an 86% of the total uncertainty in the model predictions.

Laura M White↗

NASA NTRSDate not supplied

Comparisons Between NO PLIF Imaging and CFD Simulations of Mixing Flowfields for High-Speed Fuel Injectors

The current work compares experimentally and computationally obtained nitric oxide (NO) planar laser-induced fluorescence (PLIF) images of the mixing flowfields for three types of high-speed fuel injectors: a strut, a ramp, and a rectangular flushwall. These injection devices, which exhibited promising mixing performance at lower flight Mach numbers, are currently being studied as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The EIMP aims to investigate scramjet fuel injection and mixing physics, and improve the understanding of underlying physical processes relevant to flight Mach numbers greater than eight. In the experiments conducted in the NASA Langley Arc-Heated Scramjet Test Facility (AHSTF), the injectors are placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number air flow at the entrance of a scramjet combustor. Helium is used as an inert substitute for hydrogen fuel. Both schlieren and PLIF techniques are applied to obtain mixing flowfield flow visualizations. The experimental PLIF is obtained by using a UV laser sheet to interrogate a plane of the flow by exciting fluorescence from the NO molecules that are present in the AHSTF air. Consequently, the absence of signal in the resulting PLIF images is an indication of pure helium (fuel). The computational PLIF is obtained by applying a fluorescence model for NO to the results of the Reynolds-averaged simulations (RAS) of the mixing flowfield carried out using the VULCAN-CFD solver. This approach is required because the PLIF signal is a nonlinear function of not only NO concentration, but also pressure, temperature, and the flow velocity. This complexity allows additional flow features to be identified and compared with those obtained from the computational fluid dynamics (CFD) simulations, however, such comparisons are only semiquantitative. Three-dimensional image reconstruction, similar to that used in magnetic resonance imaging, is also used to obtain images in the streamwise and spanwise planes from select cross-stream PLIF plane data. Synthetic schlieren is also computed from the RAS data. Good agreement between the experimental and computational results provides increased confidence in the CFD simulations for investigations of injector performance.

PLIF↗

NASA NTRSDate not supplied

Adaptation of Metal Additive Manufacturing Processes for the International Space Station

The In-Space Manufacturing (ISM) project at NASA Marshall Space Flight Center, in a partnership with the company, Made in Space, has previously investigated 3D printing of polymer materials on-orbit. In recent years, the project has begun exploring the potential for metal additive manufacturing (AM) on future space missions to reduce logistics and enable point-of-use manufacturing for sparing and repair. This paper provides an overview of constraints for demonstrating a manufacturing process on the International Space Station (ISS) as well as information on previous trades of available metal AM processes and their potential for in-space use. There are currently two processes in development as payloads for an ISS technology demonstration: wire+arc additive manufacturing (the Vulcan payload from Made in Space, Inc.) and bound metal additive manufacturing (the Fabrication Laboratory from Techshot, Inc). An update on both of these systems, key results to date, and future development efforts will be presented. Relevant modeling work, performed by NASA Ames Research Center, to evaluate operation of certain aspects of the bound metal AM process in a microgravity environment will also be summarized.

in-space manufacturing↗

NASA NTRS2022

S-band Network Analysis and Strategies for LEO Multi-CubeSat Science Missions

This paper presents the network architecture utilized at NASA’s Goddard Space Flight Center (GSFC) to support a mission set of five 6U CubeSats. These five CubeSats employ a multi-waveform Software Defined Radio (SDR) platform developed by Vulcan Wireless for use with NASA’s Space Relay (SR) and Direct to Earth (DTE) networks. The performance testing of the SDR is discussed via a comprehensive S-band and communication link analysis. The tested capabilities of the SDR and antenna components are reviewed in terms of the mission requirements for each CubeSat. Size Weight and Power (SWaP), required availability, and access times are discussed. The measured data from experimental compatibility testing is incorporated into detailed simulations of the CubeSat mission set to verify desired performance over the mission lifetime. The model is also used to investigate the impact of potential adverse effects on the communication links such as interference and weather conditions. This paper also reviews potential improvements from future technological advances and commercial partnerships. A collaborative investigation between GSFC and Oklahoma State University is presented in which a qualitative analysis of Hybrid RF/Optical communication strategies is performed. The data rate improvements of optical communication techniques are weighed against attitude control and science mission requirements for CubeSats, and network architectures/switching strategies are discussed. A separate analysis reviews the potential benefits of ground station partnerships, which aligns with NASA’s future goal to include commercial partners in its Earth and space network architectures, e.g. LunaNet.

Space networks↗

NASA NTRSDate not supplied

F-ANG+: A 3-D Augmented-Stencil Face-Averaged Nodal-Gradient Cell-Centered Finite-Volume Method for Hypersonic Flows

We describe the extension of a 2-D simplified face-averaged nodal-gradient (F-ANG) method to 3-D and demonstrate that the 3-D simplified F-ANG method is accomplished by augmenting the nodecentered gradient least squares stencil. This augmented stencil F-ANG method is shown to result in advection and diffusion schemes that are stable for hexahedral, prismatic, pyramidal and tetrahedral cells without having to resort to cell-averaged nodal gradients. In addition, we describe the modifications to the augmented stencil required to support the use of wall function boundary conditions. Finally we describe a consistent, face-stencil based multi-dimensional limiter procedure (MLP), and show it to be fully consistent and compatible with the linearity-preserving unstructured- MUSCL (LP-U-MUSCL) scheme for all values of kappa. These methods and schema are implemented in the cell-centered finite-volume code VULCAN-CFD, which is then used to investigate whether the robustness improvements demonstrated in 2-D carry over to 3-D by computing hypersonic flows using mixed-element grids as well as highly adapted tetrahedral grids.

Weighted Least-Squares↗

NASA NTRSDate not supplied

Deep Neural Network Based Unsteady Flamelet Progress Variable Approach in a Supersonic Combustor

Higher dimensional flamelet manifolds are essential in capturing the coupled effects of pressure gradients and unsteady chemical kinetics observed in supersonic combustion applications. Previous studies have validated the feasibility of using deep neural networks as an alternative to computation-ally intensive multidimensional flamelet table storage and lookup. This approach has demonstrated a significant reduction in memory footprint and enabled the use of larger dimensional tabulated manifolds for supersonic combustion in canonical problems. In this study, the Unsteady Flamelet Progress Variable (UFPV)-ANN model implemented in the VULCAN-CFD code is validated by the Burrows-Kurkov supersonic mixing/combustion configuration. The well characterized experimental problem consists of hydrogen injection into a supersonic vitiated crossflow that results in a lifted flame structure. The initial model consists of a 4-dimensional table where the independent variables Z, C, Xst, P are tabulated using an unsteady flamelet code with boundary conditions corresponding to the vitiated air conditions. The results show the development of a lifted flame structure and over-all acceptable agreement with finite-rate chemistry (FRC) simulation and the experimental data. Moreover, direct mapping between the independent variables and the flamelet table is replaced by a deep neural network for significant memory reduction. The results indicate that the UFPV-ANN approach can retrieve the same solution as the memory intensive lookup table approach.

Flamelet↗

NASA NTRSDate not supplied

Flight and Direct to Earth/Space Relay Communication System Architecture for GSFC CubeSat Missions

The CubeSat platform is finding increasing use in space science applications due to its low cost and comparative ease of launch. It is becoming a key scientific discovery tool in low Earth orbit (LEO) and beyond, including geosynchronous equatorial orbit (GEO), the Lagrange Points, Lunar missions, and more. The increasing complexity of these missions and their scientific goals must be supported by equal advancements in communications technology. Higher data rates and greater reliability are required every year. However, the reduced Size, Weight, and Power (SWaP) constraints of CubeSat platforms introduce unique challenges in the area of satellite communications. There is currently a lack of communication equipment tailored specifically to the CubeSat platform. This lack of standardized, tested equipment extends development time and reduces mission confidence. Furthermore, missions utilizing the CubeSat platform are often subject to more difficult design constraints. Antenna placement, size, and pointing are often subordinate to the requirements of the payload instruments and mission goals. Traditional link margin estimation techniques are insufficient in these cases, as they emphasize worst case scenarios. In reality the actual link parameters may vary widely even during a single pass. This presents new challenges in predicting communications performance and scheduling ground station contacts, but also new opportunities for improving efficiency. This paper presents the integration, testing, and validation process for a new software defined radio (SDR) designed for the CubeSat platform in conjunction with Vulcan Wireless, Inc. The SDR is planned for use on 5 upcoming CubeSat missions at NASAs Goddard Space Flight Center (GSFC) including a Geosynchronous Transfer Orbit (GTO) mission and it may also serve as a standard and well-tested option for future missions by enabling a standardized, rapid and low cost CubeSat communication system network integration process. Detailed simulations have been developed to estimate the communication performance of these missions, taking the unique antenna placements and attitude behavior of each satellite into account. These simulations allow a much more accurate analysis of the expected link margin, which varies considerably during each pass for the NASA Space Relay (SR) and Direct to Earth (DTE) network. The modelling procedures are outlined, and the results are used to predict communications performance of the missions.

Space Networks↗

NASA NTRSDate not supplied

Beyond Gravity Single Launch Payload Accommodations Shock Testing

The Vulcan Centaur launch vehicle is purpose-built to meet all the U.S. government’s National Security Space Launch (NSSL) needs. A shock test for the Centaur V’s payload accommodations (PLA), single launch (SL) configuration was successfully conducted at NASA’s Marshall Space Flight Center (MSFC) in Huntsville, Alabama, in May 2021. This test was a joint effort between MSFC, Beyond Gravity and United Launch Alliance. The primary purpose of the shock test was to characterize the shock environment induced by the separation event on components attached to the PLA and its interfaces. Additional test requirements included: (1) verifying the payload fairing (PLF) did not contact the flight hardware during the mechanical separation of the PLF horizontal separation system (HSS) and (2) evaluating the performance of the HSS spring packages during jettison. The data acquired during the PLF SL shock test will be used to validate the shock predictions from pre-test analysis of the launch vehicle and the PLF separation analysis. This paper provides an overview of the PLA SL shock test planning and execution as well as the required catch system, instrumentation and results.

separation↗

NASA NTRSDate not supplied

Beyond Gravity Single Launch Payload Accommodations Shock Testing

The Vulcan Centaur launch vehicle is purpose-built to meet all the U.S. government’s National Security Space Launch (NSSL) needs. A shock test for the Centaur V’s payload accommodations (PLA), single launch (SL) configuration was successfully conducted at NASA’s Marshall Space Flight Center (MSFC) in Huntsville, Alabama, in May 2021. This test was a joint effort between MSFC, Beyond Gravity and United Launch Alliance. The primary purpose of the shock test was to characterize the shock environment induced by the separation event on components attached to the PLA and its interfaces. Additional test requirements included: (1) verifying the payload fairing (PLF) did not contact the flight hardware during the mechanical separation of the PLF horizontal separation system (HSS) and (2) evaluating the performance of the HSS spring packages during jettison. The data acquired during the PLF SL shock test will be used to validate the shock predictions from pre-test analysis of the launch vehicle and the PLF separation analysis. This paper provides an overview of the PLA SL shock test planning and execution as well as the required catch system, instrumentation and results.

separation↗

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↗