Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LAVA”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 613 records · Page 34

Computational Support of 9x7 Wind Tunnel Test of Sonic Boom Models with Plumes

NASA and its industry partners are performing studies of supersonic aircraft concepts with low sonic boom pressure signatures. The interaction of the nozzle jet flow with the aircrafts' aft components is typically where the greatest uncertainly in the pressure signature is observed with high-fidelity numerical simulations. An extensive wind tunnel test was conducted in February 2016 in the NASA Ames 9- by 7- Foot Supersonic Wind Tunnel to help address the nozzle jet effects on sonic boom. Five test models with a variety of shock generators of differing waveforms and strengths were tested with a convergent-divergent nozzle for a wide range of nozzle pressure ratios. The LAVA unstructured flow solver was used to generate first CFD comparisons with the new experimental database using best practice meshing and analysis techniques for sonic boom vehicle design for all five different configurations. LAVA was also used to redesign the internal flow path of the nozzle and to better understand the flow field in the test section, both of which significantly improved the quality of the test data.

9x7 Wind Tunnel↗

Experimental Insights into the Geochemistry of Mercury

With the recent estimate of Mercury's surface composition from the X-Ray Spectrometer and Gamma-Ray Spectrometer that were onboard NASA's MErcury Surface, Space ENvironment, Geochemistry and Ranging (MESSENGER) spacecraft, we now have our first opportunity to directly investigate the compositions of lavas from the planet Mercury and indirectly investigate the chemical make-up of its interior. Results from MESSENGER showed exotic surface compositions with more than 3 wt% sulfur in some lavas and relatively low amounts of iron (less than 3 wt%) across the surface. These striking features are consistent with magmatism occurring under highly reducing conditions which has an impact on the thermal and chemical evolution of a planetary body. Here we'll explore the geochemical evolution of Mercury through a series of experimental studies and discuss the implications of low oxygen fugacity on elemental behavior and magmatic processes.

Vander Kaaden, Kathleen E.↗

Using GES DISC Data to Study Kilauea Volcano of 2018

Kilauea volcano in Hawaii which erupted in early May 2018 injected massive amount of SO2 and ash into the atmosphere. The lava flow during the eruption destroyed many home and neighborhoods. The SO2 plume during the eruption of Kilauea volcano is analyzed from May to August 2018 using multiple satellite products such as Level 2 TROPspheric Monitoring Instrument (TROPOMI) and Level 3 Ozone Monitoring Instrument (OMI) from the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC). GES DISC hosts multi-disciplinary Earth science data sets that can be used to analyze natural disasters, such as the Kilauea volcano. Additionally, GES DISC's Giovanni tool can be used to visualize these data. We acquired OMI through the subsetting function, which is processed by the GES DISC in-house developed backend software Level3/4 Regrider and Subsetter (L34RS) and TROPOMI using OPeNDAP.Data from the OMI OMSO2e product showed elevated levels of SO2 amounts during the eruption between May to August 2018. Similarly, ground-based stations at Hawaii Volcanoes National Park recorded higher SO2 concentrations during the same time period. This study uses wind direction from Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) to analyze the transport and dispersion of SO2 plume and map lava flows from the volcano using thermal images from Visible Infrared Imaging Radiometer Suite (VIIRS). Furthermore, satellite observations combined with socioeconomic and public health data are used to analyze its impact in public health.

KC, Binita↗

Computational Simulations of a Mach 0.745 Transonic Truss-Braced Wing Design

A joint effort between the NASA Ames and Langley Research Centers was undertaken to analyze the Mach 0.745 variant of the Boeing Transonic Truss-Braced Wing (TTBW) Design. Two different flow solvers, LAVA and USM3D, were used to predict the TTBW flight performance. Sensitivity studies related to mesh resolution and numerical schemes were conducted to define best practices for this type of geometry and flow regime. Validation efforts compared the numerical simulation results of various modeling methods against experimental data taken from the NASA Ames 11-foot Unitary Wind Tunnel experimental data. The fidelity of the computational representation of the wind tunnel experiment, such as utilizing a porous wall boundary condition to model the ventilated test section, was varied to examine how different tunnel effects influence CFD predictions. LAVA and USM3D results both show an approximate 0.5 angle of attack shift from experimental lift curve data. This drove an investigation that revealed that the trailing edge of the experimental model was rounded in comparison to the CAD model, due to manufacturing tolerances, which had not been accounted for in the initial simulations of the experiment. Simulating the TTBW with an approximation of this rounded trailing-edge reduces error by approximately 60%. An accurate representation of the tested TTBW geometry, ideally including any wing twists and deflections experienced during the test under various loading conditions, will be necessary for proper validation of the CFD.

Maldonado, Daniel↗

Computational Simulations of a Mach 0.745 Transonic Truss-Braced Wing Design

A joint effort between the NASA Ames and Langley Research Centers was undertaken to analyze the Mach 0.745 variant of the Boeing Transonic Truss-Braced Wing (TTBW) Design. Two different flow solvers, LAVA and USM3D, were used to predict the TTBW flight performance. Sensitivity studies related to mesh resolution and numerical schemes were conducted to define best practices for this type of geometry and flow regime. Validation efforts compared the numerical simulation results of various modeling methods against experimental data taken from the NASA Ames 11-foot Unitary Wind Tunnel experimental data. The fidelity of the computational representation of the wind tunnel experiment, such as utilizing a porous wall boundary condition to model the ventilated test section, was varied to examine how different tunnel effects influence CFD predictions. LAVA and USM3D results both show an approximate 0.5o angle of attack shift from experimental lift curve data. This drove an investigation that revealed that the trailing edge of the experimental model was rounded in comparison to the CAD model, due to manufacturing tolerances, which had not been accounted for in the initial simulations of the experiment. Simulating the TTBW with an approximation of this rounded trailing-edge reduces error by approximately 60%. An accurate representation of the tested TTBW geometry, ideally including any wing twists and deflections experienced during the test under various loading conditions, will be necessary for more thorough validation of the CFD.

Maldonado, Daniel↗

Improvements in Simulating a Mach 0.80 Transonic Truss-Braced Wing Configuration using the Spalart-Allmaras and k-ω SST Turbulence

Teams from the NASA Ames Research Centers (ARC) and Langley Research Center (LaRC) have been working on validating their computational fluid dynamics (CFD) results for the Boeing Mach 0.80 Transonic Truss-Braced Wing (TTBW) configuration. Experimental data used for the validation were gathered from a test conducted of a 4.5% scale Mach 0.80 TTBW model in the NASA Ames Research Center 11-by 11-Foot Transonic Wind Tunnel. The CFD simulations were initially run with both LAVA and USM3D Mixed Element solvers utilizing the Spalart-Allmaras (SA) turbulence model. A discrepancy was observed between CFD and experimental loads and moments ranging on average from 0.047 to 0.063 for lift coefficient, 16.5 to 27 in drag counts, and -0.02 to -0.032 for pitching moment, varying with the solver used. With introducing the refactored version of LAVA these values dropped to 0.034 to 0.047 for lift coefficient, 7.74 to 16.5 in drag counts, and lastly -0.012 to -0.02 for pitching moment. Based on findings from Boeing, who observed an improved comparison to experimental data when using the k-ω Shear Stress Transport (SST) turbulence model, the NASA teams conducted simulations with SST to investigate turbulence modeling sensitivities. OVERFLOW and USM3D V6 solvers were used for the comparison of SST and SA simulations. CFD results using the SST turbulence model demonstrate an improvement in matching with experimental CL values, reducing the discrepancies seen by 0.021-0.055 (62-93% reduction in discrepancy respective to the OVERFLOW and USM3D V6). The SST model has varying effects on CD, based on the solver / grid paradigm; for USM3D this prediction is only improved at higher angles of attack (above the cruise design point) resulting in a reduced discrepancy of 8.8–25 drag counts However, at the mid-to-low angles of attack, SST increases the CD discrepancy by 4-55 counts, increasing at lower angles of attack. For OVERFLOW this discrepancy only exists at the low angles of attack and maxes out at 10 counts of drag. Above α = 1.5° OVERFLOW SST reduces discrepancy by 9.3 drag counts on average.

AATT↗

Calcium Isotope Constraints on Recycled Carbonates in Subduction-related Magmas

Calcium isotope ratios are readily mass-fractionated by low-temperature depositional processes that can be recorded in sedimentary rocks and, therefore, have the potential to track distinct geochemical signatures of recycled materials in mantle-derived igneous rocks. In this study, I report calcium isotopic compositions for well-characterized lavas from the Central American volcanic arc that exhibit a range of trace element and radiogenic isotope signatures that have been used to indicate variable amounts of subducted marine carbonate and hemipelagic sedimentary rocks, fresh carbonatite lavas from Oldoinyo Lengai Volcano, Tanzania, and an intrusive carbonatite clast erupted 12.8 ka from the Laacher See Volcano, Germany. I observed no calcium isotope evidence for recycled sedimentary rocks in the Central American arc magmas or in the Oldoinyo Lengai carbonatite volcanic rocks. They all exhibit isotopic compositions similar to rocks dominated by a primitive, mantle-like, bulk silicate Earth (BSE) composition. The exception in this work is the calcium isotope composition measured in the intrusive Laacher See carbonatite that is resolvable from BSE (δ^44 Ca = 0.0). Although this calcium isotopic signature might be related to ancient carbonate recycling, the magnitude and relatively light calcium isotopic composition (δ^44 Ca = –0.4) of this carbonatite could also reflect an origin that involved partial melting of altered lithospheric mantle, from which isotopically heavy Ca-bearing fluids have been lost. The decoupled signatures between trace elements and their radiogenic isotopes and those recorded by calcium isotope data are important because they likely reflect different sources and processes, and demonstrate that sediment subduction is not a bulk mixing process.

Ca isotopes↗

Basaltic Fissure Types on Earth: Suitable Analogs to Evaluate the Origins of Volcanic Terrains on the Moon and Mars?

Basaltic eruptive fissures of the Great Rift and surroundings on the eastern Snake River Plain of Idaho, USA, and selected volcanic features in Hawai`i, Iceland and northern Africa were surveyed for their relevancy as planetary analogs. Evaluated during field investigations and in satellite imagery for structures, physiography, and geologic setting, fissures were categorized into four broad types: (1) simple, monogenetic fissures with obvious volcanic constructs or deposits, (2) monogenetic fissures now obscured by low shields or relatively large cones, (3) polygenetic volcanic rift zones with multiple vents and deposits, and (4) compound regional fissure systems or dike swarms that comprise major rift zones or large volcanic terrains. Using this classification as an initial base, we surveyed imagery of volcanic features for likely fissure vents in two major geologic settings on the Moon: floor-fractured craters (FFCs) and mare and cryptomare provinces. Two major regions on Mars, the volcanic plains around Alba Mons and the greater Tharsis region, were also surveyed for fissure types and volcanic associations of fissure-like features. The planetary surveys suggest that the proposed classification provides a suitable analog starting point to interpret structures associated with fissure systems on the Moon and Mars. With few exceptions, our survey indicates that each of the studied terrains exhibits a dominant fissure type. Type 1 fissures, most with pyroclastic deposits, prevail in lunar FFCs and mare-like regions; whereas type 2 fissures are ubiquitous in the Tharsis region of Mars and a few exist on the Moon as low shields. Type 3 volcanic rift zones are not common on either the Moon or Mars, although they might become evident in future work on chemically evolved terrains. Type 4 fissures are inferred in mare terrains, often represented as the extensions of major linear rille networks or rimae, with possibly complex dike swarms that were buried beneath voluminous mare basalt lava flows. Likewise, numerous flood lavas on Mars are possibly associated with now-obscured or difficult to define type 4 fissure systems.

Scott S. Hughes↗

Wall Modeled Large Eddy Simulations for NASA’s Jet Noise Consensus Database of Single-Flow, Round, Convergent Jets

A campaign of wall-modeled large-eddy simulations (WMLES) using structured curvilinear overlapping grids has been performed with the Launch Ascent and Vehicle Aerodynamics(LAVA) computational fluid dynamics (CFD) software to predict jet noise for single-stream axisymmetric round jets. The simulations address the new Prediction Uncertainty Reduction(PUR) technical challenge within the context of NASA’s Commercial Supersonic Technology(CST) project. The goal of PUR is to quantify and reduce uncertainties from scale-resolving simulations to assess noise characteristics of next generation quiet supersonic commercial jets during takeoff and landing conditions where the noise from the exhaust jet dominates. The focus of this effort is to generate a simulation database for single-stream axisymmetric round nozzles at several conditions both for static (no ambient co-flow), which is the focus of this article, and in-flight (flight stream co-flow) conditions, which are beyond the current scope. Nine different flow conditions ranging in jet exit Mach number from 0.38 to 1.0 with nozzle temperature ratios (NTR) from 0.84 to 2.7 have been conducted. Details of the structured overset grids, numerical discretization and wall-model are provided. Near-field comparisons to PIV data show great agreement for both velocity and normal stresses, however a systematic TKE overshoot at the nozzle exit is seen in the lip line shear-layer. A permeable Ffowcs Williams Hawkings (FWH) surface, enclosing the jet, is used to predict far-field noise from the simulated flow-field. Comparison of CFD predictions to microphone array measurements demonstrate excellent agreement within the resolved frequency range. A systematic under-prediction of far-aft observer angles larger than 150 degrees has been observed across all simulations. We achieved a cost reduction of an order of magnitude for these WMLES compared to an earlier study of this configuration due to algorithmic and software improvements. The accuracy of the results and short turnaround time demonstrate that WMLES within the LAVA framework is a cost-effective approach for jet noise predictions that could soon be incorporated into the design cycle of jet noise reduction technologies.

CST↗

Chronology of Steno Crater Basalt

The Apollo Next Generation Sample Analysis (ANGSA) program was designed to investigate specially curated samples returned by the Apollo 15 and 17 missions [1]. As part of this effort, we are examining a suite of Apollo 17 basalts collected from the rim of Steno Crater, to investigate the magmatic and post-magmatic volatile histories of these samples. Here we report initial 207Pb-206Pb chronology of accessory phases for two steno crater basalts, which will ultimately be paired with (40)Ar*/(39)Ar chronology and (38)Ar exposure ages to constrain crystallization and surface histories. The chronology of these samples will help determine if the samples derive from the one lava sequence or represent multiple lava flows, and will allow us to place petrologic and volatile data into the wider context of the magmatic suite [1,2]. Note: Extended abstract on document.

S. J. Pomeroy↗

Predicting SLS Launch Environment using a Novel Multiphase Formulation

Powerful acoustic waves generated during ignition of launch vehicles may be dangerous to the vehicle, its payload, or the surrounding structures. The water-based Ignition Overpressure and Sound Suppression (IOP/SS) system at Kennedy Space Center’s (KSC) Launch Complex 39B (LC-39B) will be used to protect the Space Launch System (SLS) from the acoustic vibrations generated during launch. The IOP/SS system uses enormous amounts of water to dampen and attenuate these sound waves. To better understand the launch environment risks and to study the effectiveness of the IOP/SS system it is desirable to have time-accurate unsteady simulations of the vehicle ignition with water-based sound suppression. This paper presents results obtained with a novel, high-order accurate, and robust numerical method designed for simulating compressible multiphase flows. A positivity-preserving finite difference scheme is utilized which is formally high-order accurate and also provably robust. Robustness is critical due to the extreme nature of the flow which exhibits highly nonlinear shock and rarefaction waves interacting with liquid-gas interfaces with density ratios of the order of 1000:1. Furthermore, the high-order accuracy (and the high resolution property) is desirable for predicting wave phenomena like IOP waves since the signal can be resolved accurately and propagated long distances with fewer grid points. This finite-difference method was developed using NASA’s Launch, Ascent, and Vehicle Aerodynamics (LAVA) Cartesian immersed boundary framework. We present a validation case by applying our solver to the SLS Scale Model Acoustic Test (SMAT). The SLS SMAT is a well-instrumented 5% scale model test meant to represent the SLS at NASA KSC’s LC-39B pad. Scale IOP tests were performed with and without the sound suppression water and included many sensors which recorded the pressure waves produced during ignition. For this validation case we conduct two simulations, likewise with and without sound suppression water, and compare the SLS SMAT pressure sensor signals with our numerical signals at identical locations. Following this validation case we present a study of the SLS launch environment to examine engineering safety concerns about the mobile launch pad. Engineers at KSC redesigned the main flame deflector at LC-39B anticipating the increased loads from the SLS and to repair damage from prior Shuttle missions. This deflector redesign made use of surface pressure and temperature data from LAVA full-scale SLS simulations without the sound suppression system. The engineers were questioning the possibility of increased pressure loads on the underside of the mobile launcher due to the water in the flame trench. Based on the results established in our simulations of the SLS SMAT, we performed updated calculations for SLS at LC-39B with and without water systems active to assess the readiness of the launch pad for Artemis I launch. Our results show that the IOP/SS system is effective at reducing the overpressure signal and overall sound pressure levels felt by the vehicle and additionally that the pressure loads experienced by the mobile launcher (ML) during engine startup is not increased by the presence of water.

EGS↗

Evidence for Early Extension and Pressure Drop Related to Magma Plumbing in Noctis Labyrinthus (Mars)

Noctis Labyrinthus is a complex district located between the western side of Valles Marineris’ large outflow channel and the Tharsis volcanic plateau. This region is characterized by a system of steep-walled canyons and an inter-connected network of scalloped pits that manifest in various orientations and size. The formation history of this area has been a topic of substantial discussion in the past. Previous studies have proposed a karst landscape with associated caves and water-related processes [1,2,3] a network of lava tubes [4,5], or a volcano-tectonic activity as a driving process for the formation of Noctis Labyrinthus [6,7,8,9]. However, none of these models have been confirmed or entirely accepted; although, the erosional nature of Noctis Labyrinthus is supported by the presence of channelized lava flows still visible on its floor [5]. In this work, we use MOLA and HRSC data to produce a new structural map of the region, we analyse fault systems (chronology and orientation), the relationships between faults and pit chains, and we finally propose a deformational model that can explain the sequences of events responsible about the formation of Noctis Labyrinthus.

M. El Yazidi↗

Enabling Supersonic Over-Land Flight Using Computational Modeling

NASA's QueSST mission is interested in exploring the potential of commercial, over-land supersonic flight, which has been banned since 1973 due to the noise intensity of the sonic booms that supersonic vehicles create. NASA has partnered with Lockheed Martin to design and construct an experimental aircraft, the X-59, that will shape the sonic boom to achieve a perceived loudness on the ground that is lower than 75dB. In this talk, we present new developments in the Launch Ascent and Vehicle Aerodynamics (LAVA) CFD solver framework which are tailored to the analysis of supersonic vehicles, like the X-59, and sonic boom propagation. Utilizing a high-order space marching method, adjoint-based grid adaptation, and robust database capabilities have enabled significant cost savings over traditional CFD methods for the same level of accuracy in perceived boom loudness. Example cases demonstrating the toolset are presented, along with comparisons to experimental wind tunnel data. These developments will enable LAVA to quickly and accurately provide simulation data to the project during X-59's first and subsequent flights in 2024 and beyond.

CST↗

Measurement of Surface Deformation Related to the December 2018 Mt. Etna Eruption Using Time-Series Interferometry and Magma Modeling for Hazard Zone Mapping

Mount Etna has erupted several times since it was first formed. Recently, Mount Etna began erupting again over 24–27 December 2018. Because it erupts frequently, Mount Etna should be observed on a frequent basis. From June 2018 to October 2019, 34 and 56 interferometric synthetic aperture radar (InSAR) images were acquired from the ascending and descending tracks of the Sentinel-1 satellite, respectively. We employed the Stanford Method for Persistent Scatterers (StaMPS) and a refined small baseline subset (SBAS) InSAR method to produce a surface deformation time-series map. In the time-series analysis, the phase signal remained unaltered with time. The Okada model was then applied to the result to generate a modeled interferogram, and the Q-LavHA program was run to generate a lava flow prediction model. A direct comparison of the results showed that Persistent Scatterers Interferometry (PSI)-StaMPS and the refined SBAS technique were comparable in terms of the displacement pattern, with slightly different velocity values obtained for individual points. In particular, a velocity range of −25 to 21 cm/yr was obtained from PSI-StaMPS, whereas a range of −30 to 25 cm/yr was obtained from the refined SBAS method. Upon computation of the vertical and east-west displacement components based on ascending and descending track data using both methods, deformation velocities of 51.5 and 52.5 cm/yr in the westerly direction on the western flank of Mount Etna were obtained from PSI-StaMPS and the refined SBAS method, respectively, whereas on the eastern flank, deformation toward the east was estimated to occur at a velocity of 50.1 or 54.2 cm/yr, respectively. PSIStaMPS estimated a vertical deformation velocity of −5.3 to 18.3 cm/yr, whereas the refined SBAS method produced a velocity range of approximately −7 to 19 cm/yr. The interferogram obtained via Okada modeling showed two fault sources in the 2018 Mount Etna eruption and a total volume change of approximately 12.39 × 10 6 m 3 . From the modeling results, a lava flow prediction model was generated using the Q-LavHA program. The approaches described in this study can be used by government officials, authorities, and other decision-makers to monitor and assess the risk of volcanic activity in the region.

Suci Ramayanti↗

Wall-Modeled Large Eddy Simulations of Transonic Buffet Over a Supercritical Airfoil

A series of scale-resolving simulations of flow over the ONERA OAT15A airfoil have been performed at an angle of attack of 3.5◦, just past the onset of buffet. The focus of this study is to document the sensitivity of the wall-modeled large eddy simulation (WMLES) methodology for curvilinear structured overset grids within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework to mesh spacing, mesh distribution, and domain size. A secondary purpose of the study is to compare the results from WMLES to unsteady Reynolds-averaged Navier Stokes (URANS) simulations and hybrid RANS-LES (HRLES) within the same LAVA solver framework. The study provides a unique perspective regarding comparisons between different turbulence modeling approaches, time-integration methods, and computational performance since many of the same numerical routines are used for all three types of simulations. The results are compared with experiments and previous numerical studies of the same geometry and flow conditions.

TTT↗

Test and Recommendation of Flight-forward Resistive Temperature Detector for Resource Prospector Mission

The Resource Prospector Mission (RPM) is an in-situ resource utilization (ISRU) technology demonstration mission planned to launch in 2018. The mission will use the Regolith and Environment Science & Oxygen and Lunar Volatile Extraction (RESOLVE) Payload to prospect for lunar volatiles such as water, oxygen, and carbon dioxide. These compounds will validate ISRU capability. The payload, particularly the Lunar Advanced Volatile Analysis (LAVA) subsystem, requires numerous temperature measurements to accurately control on-board heaters that keep the volatiles in the vapor phase to allow quantification and prevent the clogging of delivery lines. Previous spaceflight missions have proven that Resistive Temperature Detector (RTD) failure impedes mission success. The research resulted in a recommendation for a flight-forward RTD. The recommendation was based on accuracy, consistency, and ease of installation of RTDs procured from IST, QTI, and Honeywell.

LAVA↗

Test and Recommendation of Flight-Forward Resistive Temperature Detector for Resource Prospector Mission

The Resource Prospector Mission (RPM) is an in-situ resource utilization (ISRU) technology demonstration mission planned to launch in 2018. The mission will use the Regolith and Environment Science & Oxygen and Lunar Volatile Extraction (RESOLVE) Payload to prospect for lunar volatiles such as water, oxygen, and carbon dioxide. These compounds will validate ISRU capability. The payload, particularly the Lunar Advanced Volatile Analysis (LAVA) subsystem, requires numerous temperature measurements to accurately control on-board heaters that keep the volatiles in the vapor phase to allow quantification and prevent the clogging of delivery lines. Previous spaceflight missions have proven that Resistive Temperature Detector (RTD) failure impedes mission success. The research resulted in a recommendation for a flight-forward RTD. The recommendation was based on accuracy, consistency, and ease of installation of RTDs procured from IST, QTI, and Honeywell.

LAVA↗

The Role of Cooling in Pahohoe Emplacement on Planetary Surfaces.

Abundant evidence is emerging that many lavas on Mars were emplaced as slow-moving pahoehoe flows. Models for such scenarios contrast sharply with those for steep-sloped applications where gravity is the dominant force. The mode of flow emplacement on low slopes is characterized by toe formation and inflation. In the latter phase of pahoehoe flow emplacement, stagnation, inflation, and toe formation are most closely tied to the final topography, dimensions, and morphologic features. This mode of emplacement is particularly relevant to the low slopes of planetary surfaces such as the plains of Mars, Io and the Moon.

Pahoehoe Lava Flows; Pressure Dependent Simulation↗