Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Contour location”

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 271 records · Page 15

CFD Predictions of Boiling Regime Transitions during Line Chilldown validated against a 1G LN2 Experiment

Introduction Before filling a propellant tank on the ground or in Space, the transfer line between the donor and receiver tanks must be cooled down preferably by sacrificing a minimum amount of the cryogenic fluid. The cryogenic line chill-down process involves a transition between different flow boiling regimes, namely, film boiling, transition film boiling, and nucleate boiling which are complex and may be quite gravity-dependent. Capturing these boiling phenomena and predicting the transition between them in a CFD framework is new and challenging both for 1g and microgravity applications. Materials & Methods The present work addresses this challenge by employing a two-phase Eulerian approach in the context of a homogeneous fluid mixture together with the Lee phase change model to capture the film boiling regime of the chill-down process using ANSYS Fluent®. The nucleate boiling regime is predicted by incorporating an in-house developed sub-grid model that accounts for bubble nucleation, bubble growth, bubble departure diameter, and their shedding frequency. The sub-grid model is encoded and implemented into Fluent via a user-defined function for the wall-fluid heat flux calculations. The mathematical formulation and numerical implementation of the CFD model are described in detail. The coupled CFD-Subgrid model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Results Numerical simulation results show good agreements between the CFD predictions of the wall temperature evolution, rewetting temperature, and transition between film and nucleate boiling, with the experimental measurements published by Darr et al [2] for several different LN2 flowrates in the vertical pipe orientation. The CFD predictions for the wall temperature distribution indicate a rapid quenching of the wall at two upstream and downstream temperature sensing locations as compared to the experimental measurement. The only tuning parameter in the CFD model is the Lee mass transfer coefficient. The CFD Model predicts the Liedenfrost rewetting temperature in close agreement with the experiment. This marks a transition between stable and transitionary flow boiling regimes. The CFD-predicted boiling curve for the downstream sensor location is also compared against its experimental counterpart and indicates that the model is able to predict all the key temperature and heat flux parameters during the transitions from stable to transitionary film boiling to nucleate boiling in close agreement with the experiment. A sequence of predicted volume fraction, and temperature contours depicting these transitions will be presented.

Evaporation Condensation↗

Determination and representation of electric charge distributions associated with adverse weather conditions

Algorithms are presented for determining the size and location of electric charges which model storm systems and lightning strikes. The analysis utilizes readings from a grid of ground level field mills and geometric constraints on parameters to arrive at a representative set of charges. This set is used to generate three dimensional graphical depictions of the set as well as contour maps of the ground level electrical environment over the grid. The composite, analytic and graphic package is demonstrated and evaluated using controlled input data and archived data from a storm system. The results demonstrate the packages utility as: an operational tool in appraising adverse weather conditions; a research tool in studies of topics such as storm structure, storm dynamics, and lightning; and a tool in designing and evaluating grid systems.

Rompala, John T.↗

Transient liquid-crystal technique used to produce high-resolution convective heat-transfer-coefficient maps

In this transient technique the preheated isothermal model wall simulates the classic one-dimensional, semi-infinite wall heat transfer conduction problem. By knowing the temperature of the air flowing through the model, the initial temperature of the model wall, and the surface cooling rate measured at any location with time (using the fast-response liquid-crystal patterns recorded on video tape), the heat transfer coefficient can be calculated for the color isothermal pattern produced. Although the test was run transiently, the heat transfer coefficients are for the steady-state case. The upstream thermal boundary condition was considered to be isothermal. This transient liquid-crystal heat-transfer technique was used in a transient air tunnel in which a square-inlet, 3-to-1 exit transition duct was placed. The duct was preheated prior to allowing room temperature air to be suddenly drawn through it. The resulting isothermal contours on the duct surfaces were revealed using a surface coating of thermochromic liquid crystals that display distinctive colors at particular temperatures. A video record was made of the temperature and time data for all points on the duct surfaces during each test. The duct surfaces were uniformly heated using two heating systems: the first was an automatic temperature-controlled heater blanket completely surrounding the test duct like an oven, and the second was an internal hot-air loop through the inside of the test duct. The hot-air loop path was confined inside the test duct by insulated heat dams located at the inlet and exit ends of the test duct. A recirculating fan moved hot air into the duct inlet, through the duct, out of the duct exit, through the oven, and back to the duct inlet. The temperature nonuniformity of the test duct model wall was held very small. Test results are reported for two inlet Reynolds numbers of 200,000 and 1,150,000 (based on the square-inlet hydraulic diameter) and two free-stream turbulence intensities of about 1 percent, which is typical of wind tunnels, and up to 20 percent (using a grid), which is typical of real engine conditions.

Hippensteele, Steven A.↗

Snow Dunes: A Controlling Factor of Melt Pond Distribution on Arctic Sea Ice

The location of snow dunes over the course of the ice-growth season 2007/08 was mapped on level landfast first-year sea ice near Barrow, Alaska. Landfast ice formed in mid-December and exhibited essentially homogeneous snow depths of 4-6 cm in mid-January; by early February distinct snow dunes were observed. Despite additional snowfall and wind redistribution throughout the season, the location of the dunes was fixed by March, and these locations were highly correlated with the distribution of meltwater ponds at the beginning of June. Our observations, including ground-based light detection and ranging system (lidar) measurements, show that melt ponds initially form in the interstices between snow dunes, and that the outline of the melt ponds is controlled by snow depth contours. The resulting preferential surface ablation of ponded ice creates the surface topography that later determines the melt pond evolution.

arctic↗

Jet-induced ground effects on a parametric flat-plate model in hover

The jet-induced forces generated on short takeoff and vertical landing (STOVL) aircraft when in close proximity to the ground can have a significant effect on aircraft performance. Therefore, accurate predictions of these aerodynamic characteristics are highly desirable. Empirical procedures for estimating jet-induced forces during the vertical/short takeoff and landing (V/STOL) portions of the flight envelope are currently limited in accuracy. The jet-induced force data presented significantly add to the current STOVL configurations data base. Further development of empirical prediction methods for jet-induced forces, to provide more configuration diversity and improved overall accuracy, depends on the viability of this STOVL data base. The data base may also be used to validate computational fluid dynamics (CFD) analysis codes. The hover data obtained at the NASA Ames Jet Calibration and Hover Test (JCAHT) facility for a parametric flat-plate model is presented. The model tested was designed to allow variations in the planform aspect ratio, number of jets, nozzle shape, and jet location. There were 31 different planform/nozzle configurations tested. Each configuration had numerous pressure taps installed to measure the pressures on the undersurface of the model. All pressure data along with the balance jet-induced lift and pitching-moment increments are tabulated. For selected runs, pressure data are presented in the form of contour plots that show lines of constant pressure coefficient on the model undersurface. Nozzle-thrust calibrations and jet flow-pressure survey information are also provided.

Wardwell, Douglas A.↗

Mean winds of the upper middle atmosphere (60-110 km): A global distribution from radar systems (M.F., METEOR, VHF)

Since the development of the last CIRA in 1972, the number of radars providing winds in the upper middle atmosphere has increased significantly. These systems fill the data gap between 60 km and 110 km. The radars include medium frequency (MF) radars or partial reflection systems giving data from 60/70 to 100/110 km; meteor radars, 80 to 110 km, and M.S.T. radars operating as meteor radars. Data from 12 locations are shown, which represent a good Northern Hemispheric (NH) North American chain, an Oceanian chain which is mainly in the Southern Hemisphere (SH), and some Western Europe data. Generally tidal oscillations have been removed from days or groups of days, and the remaining mean winds and longer period oscillations plotted as height-time contours. Composite cross sections from the years 1978 to 1982 were formed where possible so that only the major temporal features remain.

Manson, A. H.↗

Kalman Filtering and RTS Smoothing for Arc-Jet Sample Edge Tracking

Accurate and temporally consistent measurements of test-article recession are required to characterize material response during arc-jet ablation experiments. However, image-based boundary measurements are often affected by segmentation noise, brightness variations, and frame-to-frame variability. This work extends arcjetCV [1] with filtering methods for tracking the evolving boundaries of hemispherical and ISO-Q test articles. Two boundary-tracking approaches based on Kalman filtering [2] were implemented. The first applies a point-wise Kalman filter followed by a Rauch–Tung–Striebel smoother [3] to individual boundary-point locations. The second applies the same filtering and smoothing framework to a reduced set of geometry-dependent shape parameters. The point-wise method reduces local frame-to-frame fluctuations while preserving spatial details along the detected boundary. The shape-parameter method provides a compact and geometrically constrained estimate of the sample contour. Figure 1 illustrates the two approaches for a hemispherical test article. Both methods improve temporal consistency and support more robust estimation of surface recession from arc-jet imagery. The two filtering approaches provide complementary representations of boundary evolution and improve the reliability of image-based recession measurements during arc-jet experiments.

recession measurement↗

Stripping of gas and dust from the elliptical galaxy M86

Past observations of the x ray morphology of M86 have revealed that the galaxy is experiencing ram-pressure stripping due to its large velocity (1500 km s(-1)) relative to the intracluster medium of Virgo (Forman et al. 1979, Fabian, Schwartz, and Forman 1980). Observations indicate that the x ray emitting gas in the plume of M86 is still being produced from the continual heating of gas and dust stripped from nearer the galaxy's center. Researchers obtained two-dimensional Infrared Astronomy Satellite (IRAS) images of M86 which have revealed that there are two spatially separated regions of emission, one at 60 microns and the other at 100 microns of the IRAS wavebands. The 100 microns emission, presumably from cool dust (at approximately 20 K), appears to be located near the center of the galaxy together with HI (detected by Bregman, Roberts and Giovanelli 1988), while the 60 microns emission appears to lie more than 3 arcminutes away from the optical center in a direction slightly south of the center of the plume. Optical images produced by scanning U.K. Schmidt plates, reveal asymmetric isophotal contours along the major axis of the galaxy (first reported by Nulsen and Carter in 1987, which they propose as excess emission due to star formation). This excess optical emission is co-incident with the direction of the 60 micron infra-red emission.

White, D. A.↗

A high-frequency, secondary instability of crossflow vortices that leads to transition

Three-dimensional boundary-layer transition experiments are currently being conducted on a 45 deg swept wing in the Arizona State University Unsteady Wind Tunnel. Crossflow-dominated transition is produced via a model with contoured end liners to simulate infinite swept-wing flow. Fixed-wavelength, stationary and travelling crossflow vortices are observed. The frequencies of the most amplified travelling waves are in agreement with linear-stability theory; however travelling waves at frequencies an order of magnitude higher than predicted are also observed near transition. Near the transition location, the distorted boundary layer is due to the stationary crossflow vortex and is subject to a Rayleigh-type instability in the stream direction. As a result, a high-frequency secondary instability is detected in the transition region and spatial relations of the process are well documented by coupled use of flow visualization and hot-wire measurements.

Kohama, Yasuaki↗

The local distribution of Na I interstellar gas

We present high-resolution absorption measurements (lambda/Delta lambda approximately 75,000) of the interstellar Na I D lines at 5890 A toward 80 southern hemisphere early-type stars located in the local interstellar medium (LISM). Combining these results with other sodium measurements taken from the literature, we produce galactic maps of the distribution of neutral sodium column density for a total of 293 stars generally lying within approximately 250 pc of the Sun. These maps reveal the approximate shape of the mid-plane contours of the rarefied region of interstellar space termed the Local Bubble. Its shape is seen as highly asymmetric, with a radius ranging from 30 to 300 pc, and with an average radius of 60 pc. Similar plots of the Galactic mid-plane distribution of sources emitting extreme ultraviolet radiation show that they also trace out similar contours of the Local Bubble derived from Na I absorption measurements. We conclude that the Local Bubble absorption interface can be represented by a hydrogen column density, Nu(sub ETA) = 2 x 10(exp 19) cm(exp -2), which explains both the local distribution of Na I absorption and the observed galactic distribution of extreme ultraviolet sources. The derived mid-plane contours of the Bubble generally reproduce the large-scale features carved out in the interstellar medium by several nearby galactic shell structures.

Welsh, B. Y.↗

Low-Cost Preform and Molding Processes

The entry cost for prototyping a composite component for manufacture using automated, high rate processes is prohibitively expensive in many cases, especially for small business, where tooling costs may be several $100k. Discussions with industry also indicate that many small companies, tier 1 and 2 suppliers, have an interest to mold composite parts but do not want to deal with the capital cost, material handling issues, and labor associated with dry fiber preforming operations. While the molders may locate near the end user for logistics reasons, it may be more cost effective for the performer to remain regional and invest in capital equipment to support preform automation, thus keeping costs to a minimum. This project was designed to explore and demonstrate several options to meet these industry needs. Dry fiber preforming approaches were evaluated which allow for low pressure resin infusion, single sided tooling options such a vacuum assisted resin transfer molding (VARTM) or low pressure resin transfer molding (RTM-light). Unlike sheet molding compound, SMC compression molding where typical molding pressures of 1000 psi are required to push material into the desired location; positioning of a dry fiber preform into the desired location on the tool allows for low molding pressures of 10-50 psi. Lower molding pressures allow for use of low cost, additive fabrication of polymeric tooling. Polymeric tooling is suitable for rapid part prototyping and limited production. Dry fiber preforming approaches evaluated included use of commercial chopped strand mat, robotic chopper gun deposition, and continuous fiber preform augmentation using tailored fiber placement (TFP). Use of chopped strand mat does not require a robotic deposition method, however a cutting table is generally required and there is typically 20-30% scrap generation. While various fiber areal weights are available, the preform is not readily optimized for minimal fiber use or weight savings. In contrast, a robotic chopper gun approach allows for localized deposition where fiber is required to meet structural requirements. The robotic method is highly automated and minimizes fiber scrap, however the capital cost of the equipment and engineering labor for programming can result in higher preform cost compared to chopped strand mat in certain cases depending on preform complexity. Dry fiber preforming using the robotic chopper gun method allows for creation of three dimensional forms. This approach may be ideal for molding in-house, or if the preforms stack together densely to allow for efficient shipping. Applications evaluated for this program considered trade-off between fabrication of a fully 3D preform versus production of a flat preform which is designed to readily drape into the final desired shape. Such a preform design greatly simplifies robotic programming and requires no specialized tooling. The flat preforms are easily stacked and shipped to the final molding location. Flat preforms are much easier to augment with TFP continuous fiber to provide local reinforcement. The demonstration and evaluation of these preforming and tooling methods were completed on three component applications. The first application was a battery box cover for an electric vehicle which was highly three dimensional. The second demonstrator article was comprised of complex contours and was used to demonstrate the use of TFP and RTM-light molding process. The third demonstration article was the roof of an operator’s cab for large construction equipment. The roof is relatively flat however it is comprised of complex changes in thickness which clearly demonstrate the advantage of robotic chopper gun approach as compared to using numerous preform layers of chopped strand mat. The cost trades for the various preforming methods are summarized to help guide the reader as to preforming method considerations. Finally, these demonstrations all used glass fiber roving. A fourth, exploratory task was added to evaluate the ability to make preforms using Zoltek’s carbon fiber split tow roving. We were able to adapt the chopper gun to make flat preforms for laminate testing, but further development effort would be required to make suitable preforms.

36 MATERIALS SCIENCE↗

Pressure distribution on a 1- by 3-meter semispan wing at sweep angles from 0 deg to 40 deg in subsonic flow

A 1- by 3-meter semispan wing of taper ratio 1.0 with NACA 0012 airfoil section contours was tested in the Langley V/STOL tunnel to measure the pressure distribution at five sweep angles, 0 deg, 10 deg, 20 deg, 30 deg, and 40 deg, through an angle-of-attack range from -6 deg to 20 deg. The pressure data are presented as plots of pressure coefficients at each static-pressure tap location on the wing. Flow visualization wing-tuft photographs are also presented for a wing of 40 deg sweep. A comparison between theory and experiment using two inviscid theories and a viscous theory shows good agreement for pressure distributions, normal forces, and pitching moments for the wing at 0 deg sweep.

Yip, L. P.↗

Space Station Reboost with Electrodynamic Tethers

This paper presents the results of a study of an electrodynamic tether system to reboost the International Space Station (ISS). One recommendation is to use a partially bare tether for electron collection. Locations are suggested as to where the tether system is to be attached at the space station. The effects of the tether system on the microgravity environment may actually be beneficial, because the system can neutralize aerodrag during quiescent periods and, if deployed from a movable boom, can permit optimization of laboratory positioning with respect to acceleration contours. Alternative approaches to tether deployment and retrieval are discussed. It is shown that a relatively short tether system, 7 km long, operating at a power level of 5 kW could provide cumulative savings or over a billion dollars during a 10-year period ending in 2012. This savings is the direct result of a reduction in the number or nights that would otherwise be required to deliver propellant for reboost, with larger cost savings for higher tether usage. In addition to economic considerations, an electrodynamic tether promises a practical backup system that could ensure ISS survival in the event of an (otherwise) catastrophic delay in propellant delivery.

Vas, Irwin E.↗

Coupled elastic response of open section laminated composite beams subject to generalized beam loading

A generalized beam theory for arbitrary open section laminated composite beams subjected to generalized beam loading is developed using a strength of materials approach. At present, the theory is limited to statistically determinate beams. Solution of the general system of equations yields strain and curvature distributions as functions of the axial coordinate x and the contour coordinate s. Average stresses on a given cross section or ply stresses in a local coordinate system are determined using these distributions. Data obtained for thin-walled open section composite beams reveal that elastic coupling occurs among curvatures, shearing strains, and axial strain, depending on the laminate stiffness parameters. Laminate ply layup does not affect the location of the shear center in the thin-wall approximation, but may have a significant impact for thicker walled open sections. The average axial stress distribution is found to be distinctly different from the isotropic distribution, but the average shear stress distribution is the same as that in the isotropic case.

Zvarick, Albert G.↗

Nonlinear Görtler Vortices and Their Secondary Instability in a Hypersonic Boundary Layer

Nonlinear development of the Görtler instability over a concave surface gives rise to a highly distorted inflectional flow field in the boundary layer that leads to both wall-normal and spanwise gradients in the flow. Such nonlinear structures are susceptible to strong, high-frequency secondary instabilities that may lead to the onset of laminar-turbulent transition. The present numerical study uses direct numerical simulations and linear secondary instability theory to investigate finite amplitude Görtler vortices and their secondary instability characteristics, respectively, in the hypersonic flow over an axisymmetric cone with a concave aft body. To complement previous studies in the literature wherein the Görtler instability was usually studied for a flat plate and initiated at some upstream location by imposing an eigenfunction as the inflow condition or by blowing and suction at the wall, the present investigation is focused on fully realizable Görtler instability that is excited by an azimuthally periodic array of surface protuberances. Furthermore, while the previous work had mostly focused on the secondary instability of Görtler vortices with cross-plane velocity contours that resembled bell-shaped structures, the present results confirm that fully developed mushroom structures also exist in the hypersonic regime when the Görtler vortex amplitude is sufficiently large. Computations further reveal that the dominant modes of secondary instability in these mushroom-shaped structures correspond to an antisymmetic (i.e., sinuous) “stem” mode that concentrates within the strong, nearly wall-normal internal shear layers surrounding the stem regions underneath the caps of the mushroom structures. Additionally, there exist a multitude of other significantly unstable secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes originating from the Mack mode instability play an important role during the nonlinear breakdown process.

Laminar turbulent transition↗

A 1.49 GHz atlas of the IRAS Bright Galaxy Sample

The VLA has been used in its A-, B-, C-, and D-configurations to make 1.49 GHz maps of sources in both the original and revised IRAS Bright Galaxy Samples of strong extragalactic sources selected at a wavelength of 60 microns. Integrated 1.49 GHz flux densities were obtained from the lowest resolution maps, and maps were made with higher resolution so that nearly all of the radio sources have been at least partially resolved. Only NGC 1377 was not detected at 1.49 GHz. An atlas of contour maps, a table of total flux densities plus other radio source parameters, and references to published radio maps are given. Since the infrared and radio continuum brightness distributions of IR-selected galaxies are usually similar, these high-resolution radio maps can be used as substitutes for the unavailable IR maps to indicate the sizes and precise locations of the IR-emitting regions.

Condon, J. J.↗

Comparison of central axis and jet ring coolant supply for turbine disk cooling on a SSME-HPOTP model

The test facility, test methods and results are presented for an experimental study modeling the cooling of turbine disks in the blade attachment regions with multiple impinging jets, in a configuration simulating the disk cooling method employed on the Space Shuttle Main Engine oxygen turbopump. The study's objective was to provide a comparison of detailed local convection heat transfer rates obtained for a single center-supply of disk coolant with those obtained with the present flight configuration where disk coolant is supplied through an array of 19 jets located near the disk outer radius. Specially constructed disk models were used in a program designed to evaluate possible benefits and identify any possible detrimental effects involved in employing an alternate disk cooling scheme. The study involved the design, construction and testing of two full scale rotating model disks, one plane and smooth for baseline testing and the second contoured to the present flight configuration, together with the corresponding plane and contoured stator disks. Local heat transfer rates are determined from the color display of encapsulated liquid crystals coated on the disk in conjunction with use of a computer vision system. The test program was composed of a wide variety of disk speeds, flowrates, and geometrical configurations, including testing for the effects of disk boltheads and gas ingestion from the gas path region radially outboard of the disk-cavity.

Kim, Y. W.↗

Deep Basalt Aquifers in Orcus Patera, Elysium Basin Mars: Perspectives for Exobiology Exploration

Direct indicators of shorelines, spillways, and terraces allowed to determine the extent of the Elysium Paleolake between the contour-lines 1000 and 500 m below the Martian datum. The Elysium Paleolake is bordered north by Orcus Patera (14N/181W), which lies west of the Tartarus Montes and Tartarus Colles. The Orcus Patera displays an ellipse-shaped collapsed caldera of 360-km long and 100-km wide. Viking topographic data show that the bottom of the caldera is located at 2500 below the Martian datum, and surrounded by a steep-walled ram art which crest is located at about 0 m elevation. Considering the localization of Orcus Patera in the Elysium paleolake, its altimetry, and the magmatic origin of this caldera, we propose the existence of a paleolake in Orcus Patera generated (a) by juvenile water from magma during the Noachian period, and (b) by intermittent influx of the Elysium Basin from Hesperian to Amazonian. Results are encouraging to consider this site as a potential high-energy source environment for microbial communities. are circumscribed by a 50-km wide lava field mapped as Noachian material. The structure of Orcus Patera represents the record of material erupted from a magmatic reservoir. The caldera is enclosed by steep inner walls (25% measured from topographic data), values which could be in agreement with the presence of a deep magmatic reservoir, as suggested by the typology of Crumpler et.al. The depth of the caldera might be due to the collapse of the magma reservoir, and the release of gases accompanying the magma thermal evolution. Origins of water for the paleolake(s): The water that generated a paleolake in Orcus Patera may have come from two origins: (1) Juvenile water: Plescia and Crips estimated a magma H20 content by weight between 0.5% and 1.5% using for the first value a comparison with terrestrial basalt, and for the second values from a Martian meteorite. The amount of H20 can be estimated by the volume of erupted lava, and the lava content of the caldera. In this study, we adopt a water content of 1%. The total volume of magma that has been contained in the caldera, and the volume of lava contained in the observed lava field is about 110 x 10(exp 6) cubic km, that gives a total volume of 1.10 x 10(exp 6) cubic km of water. The juvenile water expelled by the overpressure within the magma chamber charged with desolved water-vapor may have moved into the crust. The decrease in overburden pressure led to bubble formation. The ascent of these bubbles generated a pressurization of the magma, which was sufficient to fracture the overlaying magma layer, (2) Water from Elysium paleolake. During the Amazonian, the rise of the Elysium paleolake level generated an overspilling that supplied the caldera with water. The southern portion of the crest shows a deep gap 12-km wide at -1500 m elevation, locating the gap between 500 to 1000 in below the assumed water of Elysium paleolake, thus facilitating the influx of Elysium paleolake water into Orcus Patera. Bathymetric calculations give a floor area of 25,500 sq km at -2000 m elevation, and a water volume of 42,000 cubic km, with a lake-level at -1500 m. A substantial amount of water may have percolated through the fractured lava, and part of the volume may have overspilled the northern crest of Orcus Patera to debouch in the Tartarus Montes region. We envision the formation of a subsurface aqueous environment in basaltic rocks at the contact of the two water-source origins, possibly the percolating surface lake water, and more likely the juvenile water. Similarly to terrestrial calderas, Orcus Patera might be surrounded by ring-fractures caused by the collapse of the magma chamber that followed the release of gases. These ring-fractures may have been covered later by sedimentation in the caldera (lacustrine, aeolian, and volcanic), and by mass wasting. The detumescence of the magma in the caldera, and the vesiculation of the juvenile water may have operated simultaneously. Comparatively to terrestrial melts, Martian iron-rich melts are denser. This greater density implies greater effusion rates (eight-times terrestrial values), and larger fissuration widths (two-times terrestrial ones). With increasing vesiculation of magma, the bubbles interact with one-another because there are of similar pressure. They make a magma froth at the contact with the caldera surface, and on the walls of the fractures. In the saturated magma, froth, where the volume ratio of gases-to-liquid is about 4:1, the bubbles form a huge surface area of interconnected spaces. Bubbles near the caldera surface disrupt the magma, and fragmentation takes place, which moves downward through the magma column. On Earth, the bubbles are likely to grow between 1 and 50 mm in diameter due to the difference between the magma surface tension, and the bubble supersaturation pressure. The Martian low-pressure at surface level is likely to accelerate the expansion of the bubbles, and increase their final diameter and number, creating more voids in the magma. The strong magma froth with enclosed juvenile water bubbles interconnected with exsolved gas bubbles constitute a potential geothermal environment for geochemical energy production from basalt and water that does not require excessive temperatures. This process can start at +20C. Similar types of environments have been shown on Earth as potential energy sources for microbial metabolism, and could have provided deep aqueous basaltic niches for possible Martian microorganisms, even geologically recently. During the Amazonian, combination of volcanism and water activity still existed on Mars. Moreover, this type of potential niches open ways for investigation of possible oases of extinct or extant life, not only on paleolakes, and surface hydrothermalism spring areas, but also all large systems of fossae, which combine hydrologic and volcanic activities, and which provide an energy source, and an underground shelter to prevent surface UV bombardment. Additional information contained in the original.

Grin, E. A.↗