Search NASASearch

SEARCH · Search NASA

Results for “Topological Phases of Matter”

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.

1,514 records · Page 2

Diffusion Quantum Monte Carlo Calculation of the Austenite and Martensite Phases of NiTi

NiTi is a promising material for smart and active technologies due to its exhibition of the shape memory effect, superelasticity, and biocompatibility. The shape memory effect is tied to the reversible transition between the austenite and martensite phases. A major research direction is to alloy NiTi with Zr, Hf, Pd, Pt, etc., in order to tune the martensitic transition temperature (MTT). Modeling the MTT from first principles is challenging because the lattice dynamics is complicated by anharmonicity and various low-energy structures. Using density functional theory, the energy difference between the austenite and martensite phases of NiTi varies by up to 100 meV/atom depending on the choice of density functional, which is of the same order of the energy difference itself. Consequently, free energy calculations with different functionals can result in estimates of the MTT that vary by several hundred K. Using diffusion quantum Monte Carlo, we calculated the energy difference between the B2 and B19' structures of NiTi to be 70.9 +- 2.5 meV/atom.

Kevin K Ly

NIAC Phase I Final Report Lunar South Pole Oxygen Pipeline (LSPOP)

The Lunar South Pole Oxygen Pipeline (LSPOP) Phase I NIAC, is to analyze the feasibility of constructing a pipeline at the Moon’s South Pole for transporting gaseous oxygen from point of generation to point of use. A lunar pipeline has never been pursued and will revolutionize lunar surface operations for the Artemis program and reduce cost and risk. Our starting concept is for a 5 km pipeline to transport oxygen gas from an oxygen production source, for example from our molten regolith electrolysis (MRE, currently at TRL 5, see references[1],[2],[3],[4],[5],[6] for technical discussions of the MRE process) extraction site, or any other source, to an oxygen storage/liquification plant near a lunar base. The pipeline is designed to: 1) be constructed robotically from regolith-derived metals with minimal material transferred from Earth, 2) be repairable robotically, 3) have an oxygen flow rate of ~2 kg/hour, which is commensurate with the NASA initial projected need of 10,000 kg/year, 4) operate with minimal power over the lifetime of the pipeline, with a goal of high operational reliability and the ability to survive in the lunar environment for > 10 years. We compare this concept to a pipeline concept with pipe produced on Earth and then transported to and assembled on the Moon. Both approaches will use earth-based compressors, valves, and fittings which are integrated into the pipeline on the Moon.

NIAC Phase I

A Battery Charger and State of Charge Indicator

A battery charger which has a full wave rectifier in series with a transformer isolated 20 kHz dc-dc converter with high frequency switches, which are programmed to actively shape the input dc line current to be a mirror image of the ac line voltage is discussed. The power circuit operates at 2 kW peak and 1 kW average power. The BC/SCI has two major subsystems: (1) the battery charger power electronics with its controls; and (2) a microcomputer subsystem which is used to acquire battery terminal data and exercise the state of charge software programs. The state of charge definition employed is the energy remaining in the battery when extracted at a 10 kW rate divided by the energy capacity of a fully charged new battery. The battery charger circuit is an isolated boost converter operating at an internal frequency of 20 kHz. The switches selected for the battery charger are the single most important item in determining its efficiency. The combination of voltage and current requirements dictate the use of high power NPN Darlington switching transistors. The power circuit topology is a three switch design which utilizes a power FET on the center tap of the isolation transformer and the power Darlingtons on each of the two ends. An analog control system is employed to accomplish active input current waveshaping as well as the necessary regulation.

Latos, T. S.

Characterization of a New Phase and Its Effect on the Work Characteristics of a Near-Stoichiometric Ni30Pt20Ti50 High-Temperature Shape Memory Alloy (HTSMA)

A new phase observed in a nominal Ni 30 Pt 20 Ti 50 (at.%) high temperature shape memory alloy has been characterized using transmission electron microscopy and 3-D atom probe tomography. This phase forms homogenously in the B2 austenite matrix by a nucleation and growth mechanism and results in a concomitant increase in the martensitic transformation temperature of the base alloy. Although the structure of this phase typically contains a high density of faults making characterization difficult, it appears to be trigonal (-3m point group) with a o ~1.28 nm and c o ~1.4 nm. Precipitation of this phase increases the microhardness of the alloy substantially over that of the solution treated and quenched single-phase material. The effect of precipitation strengthening on the work characteristics of the alloy has been explored through load-biased strain-temperature testing in the solution-treated condition and after aging at 500 ºC for times ranging from 1 to 256 hours. Work output was found to increase in the aged alloy as a result of an increase in transformation strain, but was not very sensitive to aging time. The amount of permanent deformation that occurred during thermal cycling under load was small but increased with increasing aging time and stress. Nevertheless, the dimensional stability of the alloy at short aging times (1-4 hours) was still very good making it a potentially useful material for high-temperature actuator applications.

A Garg

Wave Rotor Enhanced Nuclear (WREN) Propulsion: NASA Innovative Advanced Concepts (NIAC) - Phase I Final Report

Nuclear Thermal Propulsion (NTP) is identified as one of the preferred propulsion technologies for manned missions throughout the solar system (NASA MSFC).[1, 2] The state-ofthe-art NTP cycle is based on a solid core Nuclear Engine for Rocket Vehicle Application (NERVA)[3] class technology (Fig. 1) that is envisioned to provide a specific impulse of 900 seconds doubling chemical rocket performance (450 seconds). Even with this impressive increase, the NTP NERVA designs still have issues providing adequate initial to final mass fractions for high ΔV missions.[4] Nuclear Electric Propulsion (NEP) can provide extremely high Isp (2,000 to over 10,000 seconds) but with only low thrust and limits on mass to power ratios. The need for an electric power source also adds the issue of heat rejection in space where thermal energy conversion is at best 30-40% under ideal conditions. NASA Space Technology Mission Directorate (STMD) has recently expressed interest in finding advanced nuclear propulsion technology through the NASA Go:Thrust RFI.[5, 6] A novel Wave Rotor (WR) topping cycle has been proposed for our NASA NIAC concept. It promises to deliver similar thrust as NERVA class NTP propulsion, but with Isp in the 1,200-2,000 second range. Coupled with an NEP cycle, the duty cycle Isp can further be increased (1,800-4,000 seconds) with minimal addition of dry mass. This bimodal design enables fast transit trajectories for manned missions to Mars and revolutionizes the deep space exploration of our solar system.

Nuclear Thermal Propulsion

First-Principles Investigation of Phase Transitions in RuNb, RuTa, and Ru-Nb-Ta Ultrahigh Temperature Shape Memory Alloys

Ultrahigh temperature shape memory alloys (UHT-SMAs) have transition temperatures above 600 C, and they have found applications for sensing and actuating devices in aerospace industry. Among very few such UHT-SMAs currently known are Ru-based alloys such as RuNb and RuTa, whose martensite structures and phase transitions are totally different from those of NiTi-based SMAs and were poorly understood. In this work, we carried out a systematical study of RuNb using first-principles total energy calculations and molecular dynamics (MD) simulations. The transition paths and mechanisms in cubic → tetragonal →monoclinic transitions are revealed. and the transition sequence and martensitic transition temperatures are determined (MTTs) by evaluating the Gibbs free energies using thermodynamic integration. The calculated MTTs are in very good agreement with the experimental data. We found that the monoclinic phase at the second transition has the 𝑃21/m symmetry instead of experimentally identified 𝑃2/m. Our calculations demonstrate that RuTa has very similar phase transitions to those of RuTa. Furthermore, we studied the Ru0.5Nb0.25Ta0.25 ternary. Our results can explain the measured significant drop of MTT (~ 200 C) in the second transition for Ru0.5Nb0.25Ta0.25 compared with those of RuNb and RuTa, while in the first transition its MTT is between MTTs of RuNb and RuTa. The insights gained by this study and the verified ab initio methods for accurate MTT calculations can be applied to fast screen and quantitatively design novel UHT-SMAs having similar properties with desirable MTTs and much reduced cost.

Zhigang Wu

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG

Lyotropic Liquid Crystalline Phase Behavior of Boron Nitride Nanotube Aqueous Dispersions

Boron nitride nanotubes (BNNTs) are gaining significant interest due to their outstanding mechanical and thermal properties, as well as their potential to serve as a model nanorod system. Processing BNNT liquid crystalline (LC) dispersions enables precise control over BNNT orientation in macroscopic assemblies, while their low absorption in the visible spectrum facilitates studying LCs at exceptionally high concentrations. Here, we investigate the behavior of BNNTs in aqueous solutions stabilized by the surfactant sodium deoxycholate (SDC), examining the effect of BNNT purity and BNNT−SDC concentrations on lyotropic LC formation. We disperse up to 15wt% BNNT in SDC solutions and use polarized light microscopy to detail the transition from an isotropic state to a biphasic regime, where isotropic and nematic domains coexist due to phase separation, to a single fully nematic phase. Cryogenic electron microscopy provides direct evidence of BNNT alignment within nematic domains. Our results show that enhanced depletion-induced attractions, driven by increased surfactant concentration, lower the threshold concentration of BNNT required to form nematic domains. In contrast, low surfactant concentrations relative to BNNT result in insufficiently coated nanotube surfaces, leading to poor dispersions and BNNT aggregation. Additionally, we fabricate well-aligned BNNT films from aqueous LC dispersions. Our findings advance the understanding of BNNT LCs, offering insight into controlling their orientation and highlighting their potential for high-performance materials.

Joe Khoury

The Effect of Temperature and Pressure on the Distribution of Iron Group Elements Between Metal and Olivine Phases in the Process of Differentiation of Protoplanetary Material

The distribution patterns of Ni, Co, Mn, and Cr were studied in olivines of various origins: from meteorites (chondrites, achondrites, pallasites), which are likely analogs of the protoplanetary material, to peridotite inclusions in kimberlite pipes, which are analogs of mantle material. According to X-ray microanalysis data, each genetic group of olivines is characterized by a specific concentration of these elements. Nickel is concentrated (up to 0.34 percent) in peridotite olivines, while manganese is concentrated in meteoritic olivines. The maximum chromium content (0.2 percent) was found in ureilites, which were formed under reducing conditions. Experiments at pressures of 20 to 70 kbar and temperatures of 1100 to 2000°C have shown that in a mixture of olivine and Ni metal or NiO nickel enters the silicate phase (up to 4 percent), displacing Fe into the metallic phase. Equilibrium temperatures were estimated from the Fe, Ni distribution coefficients between the metal and olivine: 1500 K for pallasites, 1600 K for olivine-bronzite H6 chondrites, 1200 K for olivine-hypersthene L6, 900 K for LL6, and 1900 K for ureilites (at P = 1 atm). The equilibrium conditions of peridotites are close to T = 1800 K and P over 100 kbar. The distribution patterns of the transition elements are explained on the basis of physical-chemical properties. It is concluded that there is a sharp difference between the conditions of differentiation of the protoplanetary material at the time meteorites were formed and the conditions of differentiation of the planets into concentric layers.

A P Vinogradov

NASA Hybrid Thermally Efficient Core (HyTEC) Phase 2 Q3’GFY25 to Q2’GFY26 Unlimited Rights Annual Report for CLIN 1 and 2

The HyTEC Phase 2 Project focuses on the development and demonstration of high-power density, small engine core gas turbine engine technologies. The focus of HyTEC on these small engine core technologies will provide direct benefits to the next single-aisle class aircraft in terms of thermal efficiency, as well as integrate with other technologies, through increased hybridization, as they mature to provide substantial continuous fuel burn reductions during the aircraft lifecycle. These advances will strengthen the US position in the commercial aviation engine market and enable cost benefits in commercial aviation. In Government fiscal year 26 the HyTEC Project was closed and contract moved to the Subsonic Vehicle Technologies and Tools Project. GE Aerospace is executing two major work plans as part of the HyTEC Phase 2 contract. The first is Contract Line-Item Number (CLIN) 001 and its objective is to mature High Pressure Turbine (HPT) aerodynamics technology to TRL 5 through multiple rigs so that it may be incorporated into a complete TRL 6 evaluation and unlock learnings for the engine core demonstration as applicable. CLIN 002, the second work plan, is to design, procure and test a cost shared portion of technologies on the TRL6 Compact Core test vehicle. The compact core demonstrated through HyTEC Phase 2 along with the HPT rig TRL advancement directly matures the technologies being developed as part of the CFM RISE demonstrator program which is targeting a more than 20% fuel burn reduction at the engine level versus today’s state of the art. The RISE program’s objective is to advance both the novel Open Fan architecture, compact core technologies and advanced systems such as hybrid electric needed to achieve the desired fuel burn reduction to TRL6 ahead of a new product introduction in the single-aisle class aircraft in the 2030s timeframe.

gas turbine engine

Cryogenic Flow Boiling in Microgravity: Effects of Reduced Gravity on Two-Phase Fluid Physics and Heat Transfer

With the growing interest in space exploration, cryogenic technologies involving two-phase flow and heat transfer are in high demand to successfully procure advanced space applications such as fuel depots and nuclear thermal propulsion (NTP) systems for deep space missions. However, the unique and extreme thermal properties of cryogenic fluids introduce distinct flow boiling fluid physics and energy transport phenomena, which differ significantly from those observed with conventional fluids. Understanding the unique two-phase physics in cryogenic flow boiling remains an ongoing challenge. Furthermore, the lack of readily available microgravity cryogenic steady-state heat transfer data hinders the assessment of gravitational effects on cryogenic flow boiling. This study aims to elucidate the gravitational effects on two-phase fluid physics and heat transfer by conducting the first-ever experimental measurement of cryogenic flow boiling performance using a steady-state heated method in a reduced gravity environment. Parabolic flight experiments were performed to acquire both heat transfer measurements and high-speed video of interfacial behaviors, under varying gravity levels (microgravity, hypergravity, Lunar gravity, and Martian gravity). The experiments involved flow boiling of liquid nitrogen (LN 2 ) with a near-saturated inlet along a circular heated tube of dimensions 8.5-mm inner diameter and 680-mm heated length. The operating parameters varied are mass velocity of 398.3 - 1342.8 kg/m2s, inlet quality of -0.08 to -0.01, and inlet pressure of 413.68 - 689.48 kPa. Captured microgravity flow patterns range from bubbly to annular, all having vapor structures that are larger than those under higher gravity levels. Under microgravity, absence of buoyancy yields symmetrical vapor structures without flow stratification, laying a physical foundation for the distinct two-phase heat transfer trends during LN 2 flow boiling in microgravity. Transient data collected during the flight parabolas exhibited decreasing heated wall temperature as the aircraft transitioned from hypergravity to microgravity phases. The temperature variation indicated an enhancement in flow boiling heat transfer with decreasing gravity levels and a reduction with increasing gravity levels. The effect of reduced gravity on cryogenic flow boiling heat transfer coefficient (HTC) is discussed based on steady state heat transfer analysis. Seminal HTC correlations are evaluated against the measured microgravity HTC data, of which one is identified for superior accuracy in predicting microgravity data. Finally, a new HTC correlation is proposed to improve accuracy of microgravity predictions, yet there still exists room for further improvement with future terrestrial flow boiling experiments at different flow orientations relative to Earth gravity.

Microgravity

Assessment of Aircraft Electric Engine Regulatory Special Conditions

This report provides a regulatory gap assessment of the Federal Aviation Administration (FAA) and the European Union Aviation Safety Agency (EASA) Special Conditions (SC) for electric engines. This information is presented to support ongoing international harmonization efforts for electric engine certification, to establish a level of safety consistent across regulatory bodies and equivalent to that of existing combustion engines. HS Advanced Concepts LLC assembled a team of subject-matter experts to develop a comprehensive comparison of the regulatory text of EASA SC E-19 and the FAA special conditions for the certification of electric engines. The team included Mr. Chris Baczynski of BAC Aerospace Inc. and Mr. Zach Kaufman, formerly with magniX.

aircraft certification

Photophoretic Propulsion Enabling Mesosphere Exploration NIAC Phase I Final Report

This Phase I report presents a comprehensive study on photophoretic flyers—innovative, ultralight, solar-powered vehicles that harness photophoretic forces generated via Knudsen pumping to achieve sustained flight in the mesosphere (50–80 km altitude). By integrating advanced materials such as nanocardboard— characterized by its extremely low areal density (~1 g/m²) and high bending stiffness—with ultrathin light-absorbing coatings, the project converts incident solar radiation directly into a directed thrust. Extensive experimental investigations, coupled with high-fidelity computational fluid dynamics (CFD) simulations using ANSYS Fluent, validate the concept across various three-dimensional geometries, including spherical, conical, and rocket-shaped configurations. These simulations bridge the gap between free-molecular and continuum flow regimes, demonstrating that optimized designs can generate lift forces sufficient to support kilogram-scale payloads even in low-pressure environments. At the heart of this innovation is the use of Knudsen pumping, where temperature gradients across porous surfaces induce directional gas flow, creating a modest overpressure that provides lift. The report introduces an analytical framework that interpolates between the well-known low-Reynolds number drag regime and the high-Reynolds number momentum theory. This model accurately predicts lift based on design parameters such as microchannel dimensions, porous wall geometry, areal density, and nozzle exit area. For instance, simulations indicate that 10-meter-scale structures with carefully engineered porous walls can achieve the necessary pressure differential to support scientifically significant payloads (~1 kg). The study also explores a hybrid propulsion approach that combines solar buoyancy with photophoretic lift. Initially, solar heating creates a buoyant force that elevates the flyer to mesospheric altitudes. Once in the optimal pressure range, the photophoretic mechanism—powered by Knudsen pumping—takes over as the primary source of lift, ensuring stable, long-duration flight. This dual-mode operation not only facilitates the deployment of photophoretic flyers but also broadens the potential applications for mesospheric exploration. In addition to propulsion, the report investigates the integration of photophoretic thrusters for trajectory control of existing research balloons in the upper stratosphere. By dynamically adjusting the nozzle orientation and controlling the flow-through velocity, these thrusters provide precise maneuverability, enabling the flyers to counteract atmospheric disturbances and adjust their flight paths in real time. For example, a photophoretic thruster approximately 7.5 by 7.5 meters in size could be unfolded below a payload gondola of a 60 million-cubic-foot zero-pressure balloon. Such a thruster can provide horizontal speed control of approximately 1 m/s using only sunlight and no moving parts (except those needed to track the Sun and control the jet direction). Importantly, photophoretic thrusters operate more efficiently at higher altitudes, which is complementary to known trajectory control techniques, such as propellers and tethered wings, which are more effective at lower altitudes. Finally, the report identifies three scientific research thrusts where mesospheric aircraft technology can have a profound impact: atmospheric tides, characterization of gravity waves, and investigation of mesospheric instabilities. Overall, the findings of this Phase I project represent a significant advancement in photophoretic propulsion technology. By demonstrating that large-scale, ultralight structures can be powered solely by solar radiation—via carefully engineered Knudsen pumping—this work lays a robust foundation for scalable, near-space flight architectures. Future refinements in material fabrication, structural optimization, and integrated trajectory control are expected to further enhance performance, paving the way for operational demonstrations that could revolutionize atmospheric science, remote sensing, and communication networks.

Knudsen Pump

Thermal Cycling and Isothermal Deformation Response of Polycrystalline NiTi: Simulations vs. Experiment

A recent microstructure-based FEM model that couples crystal-based plasticity, the B2<-> MB190 phase transformation and anisotropic elasticity at the grain scale is calibrated to recent data for polycrystalline NiTi (49.9 at.% Ni). Inputs include anisotropic elastic properties, texture and differential scanning calorimetry data, as well as a subset of recent isothermal deformation and load-biased thermal cycling data. The model is assessed against additional experimental data. Several experimental trends are captured - in particular, the transformation strain during thermal cycling monotonically increases and reaches a peak with increasing bias stress. This is achieved, in part, by modifying the martensite hardening matrix proposed by Patoor et al. [Patoor E, Eberhardt A, Berveiller M. J Phys IV 1996;6:277]. Some experimental trends are underestimated - in particular, the ratcheting of macrostrain during thermal cycling. This may reflect a model limitation that transformation-plasticity coupling is captured on a coarse (grain) scale but not on a fine (martensitic plate) scale.

Phase Transformations

Ejecta Management in a Safe Lithium Ion Battery Design

As lithium-ion battery energy densities continue to rise, managing the heat and pressure generated during failure events has become increasingly critical. Safety systems must effectively relieve pressure without releasing sparks, flames, or particulate matter, requiring robust filtration solutions. The challenge is compounded by the reduced free volume available for gas expansion in high-density designs, which increases the demands on these filters. While significant progress has been made through experimental studies and modeling efforts to understand the behavior of ejecta during battery failures, there remains a pressing need for practical, rule-of-thumb sizing parameters. These parameters would help correlate high-energy waste streams with appropriate filter design, ensuring reliable containment and safety. This talk will explore recent experimental findings in this area and discuss the potential pathways for developing these essential sizing guidelines.

Ejecta Management

Rewet Temperature Correlations for Liquid-Nitrogen Boiling Pipe Flows Across Varying Flow Conditions and Orientations

In many convective liquid-vapor phase change heat transfer engineering applications, cryogenic fluids are widely used in industrial processes, spacecraft and cryosurgery systems, and so on. For example, cryogens are usually used as liquid fuels such as liquid hydrogen and oxygen in the rocket industry, liquid nitrogen (LN2) and helium are frequently used to cool superconducting magnetic device for medical applications. In these systems, proper transport, handling, and storage of cryogenic fluids are of extreme importance. Among all the cryogenic transport processes performed in room temperatures, quenching, also termed chilldown, is a unavoidable initial, transient phase-change heat transfer process that brings the system down to the cryogenic condition. The Leidenfrost temperature or rewet temperature that signals the end of film boiling is practically considered the completion point of a quenching process. Therefore, rewet temperature has been considered the most important parameter for the engineering design of cryogenic thermal management systems. As most of the previous correlations for predicting the Leidenfrost temperature and the rewet temperature have been basically developed for water, they are shown to disagree with recent liquid nitrogen pipe chilldown experiments in upward and downward flow directions over a wide range of flow rates, pressures, and degrees of inlet subcooling. In addition to a complete review of the literature, two new correlations are presented in this work, one based on bubble growth and another based on the theoretical maximum limit of superheat. Each correlation performs well over the entire data set.

liquid nitrogen

Synchronized Eruptions on Io: Possible Evidence of Interconnected Subsurface Magma Reservoirs

On 27 December 2024, Juno's JIRAM (Jovian InfraRed Auroral Mapper) instrument observed an unprecedented volcanic event in Io's southern hemisphere, covering a vast region of ∼65,000 km2, near 73°S, 140°E. Within the imaged region, only one hot spot was previously known (Pfd454). This feature was earlier estimated to cover an area of 300 km2 with a total power output of 34 GW. JIRAM results show that the region produces a power output of 140–260 TW, over 1,000 times higher than earlier estimates and likely exceeding the brightest eruption ever recorded on Io, that of Surt in 2001 (∼80 TW). Three adjacent hot spots also exhibited dramatic power increases: P139, PV18, and an unnamed feature south of the main one that surged to ∼1 TW, placing all of them among the top 10 most powerful hot spots observed on Io. A temperature analysis of the features supports the simultaneous onset of these brightenings and suggests a single eruptive event propagating beneath the surface across hundreds of kilometers; this is the first time this has been observed on Io. This in turn would imply a connection among the hotspots' magma reservoirs, while other nearby hotspots that have been known to be active in the recent past, such as Kurdalagon Patera, appear unaffected. The simultaneity supports models of massive, interconnected magma reservoirs. The topology of these regional magma systems may resemble that of a large-scale sponge, in which the massive reservoirs are the pores, interconnected through a largely solid outer shell.

A Mura