SEARCH · Engineering Papers

Results for “Fluid Properties”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Variable Property Flow in Rectangular Ducts with Repeated Rectangular Rib Roughness

PHOENICS was used to simulate turbulent flow of supercritical hydrogen through asymmetrically heated rectangular ducts with one side rib-roughened, in order to provide design data for active cooling systems of hypersonic aircraft where the hydrogen fuel is used as coolant. Variable property and conjugate heat transfer effects were considered. The computations using PHOENICS required provision of property sub-routines for supercritical hydrogen, a low-Reynolds number form of the k-epsilon turbulence model, a rough-wall wall function and a procedure for calculating proper diffusion coefficients at the fluid-solid interface. The effect of each design variable such as width and height of duct and wall thickness, on pressure drop and maximum wall temperature was investigated. PHOENICS version 1.4 was used on an IBM 3090 computer.

Youn, B.↗

NASA NTRS1982

The flux ejection dynamo with small diffusivity. I - Basic properties. II - Illustrative examples

In the first part of this paper, the flux ejection dynamo in an infinitely broad convective layer of finite depth is treated in the limit of electrical conductivity, where the convective layer upper surface is open to empty space while the lower surface is closed to the passage of fluid and field. The idealized convective overturning used allows an exact description of the symmetry and distribution of the magnetic field carried with the fluid. The total flux of the mean horizontal field grows linearly with the number of convective cycles. A discussion is given of the boundary conditions to be applied to an astronomical body whose flux ejection dynamo operates at its surface. In the second part, examples are presented which illustrate the consequences of reverse flux ejection from the surface of a convective layer of conducting fluid. Reverse flux ejection generally has the opposite effect of magnetic buoyancy, burying the fields rather than bringing them through the surface. Reverse flux injection at the surface of an alpha-omega dynamo profoundly alters the character of the solutions of dynamo equations.

Parker, E. N.↗

NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

NASA NTRSDate not supplied

Cavitation Effects on the Structural Dynamics of Turbomachinery Components

The structural integrity of inducer and impeller blades in rocket engine turbomachinery must be evaluated in the face of complex excitation mechanisms including fluctuating pressures due to cavitation. Cavitation occurs when the local fluid pressure drops below the vapor pressure, causing the formation of vapor-filled bubbles. Cavitation can exist to various extents within the typical operating range of rocket engine turbopumps. Despite recent progress towards reducing uncertainties in structural dynamic models of turbomachinery components, the extent to which pump cavitation affects the structural dynamic properties (i.e., natural frequencies, damping, and mode shapes) of inducer and impeller blades remains largely unknown. To study the structural dynamic effects of cavitating flows, experiments are conducted in a high-speed water tunnel. The test articles mimic a single inducer blade that is unwrapped from the hub, thus resembling a low aspect ratio cantilevered plate. To induce cavitation, the plates are oriented at a non-zero angle of attack and have rough strips near their leading edges. To control the extent of cavitation, the static pressure of the water tunnel is controlled between roughly 0.5 and 1.5 atm. A high-speed camera is used to quantify the coverage of the cavitations heet. The sheet coverage is then correlated to changes in fluid-added mass and hydroelstic damping. In addition, a corresponding two-way coupled CFD/FE model is under development for test/analysis correlation.

Ben Davis↗

NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

NASA NTRS1967

Dielectric properties of cryogenic liquids.

Electric breakdown, dielectric constant and dissipation factors of cryogenic liquid nitrogen, hydrogen and helium

LIQUID HYDROGEN↗

NASA NTRS1991

Physical properties of cryovolcanic brines: Applications to the evolution of Ganymede

Carbonaceous chondrites contain abundant veins of water soluble salts, including carbonates and hydrated sulfates of Mg, Ca, Na, Ni, and Fe. These constitute over 1/4 of the mass of the meteorite Orgueil. Magnesium sulfate is the most abundant salt, constituting nearly half the mass of all salt components combined (anhydrous), and 73 pct. of the highly water soluble salts. The assumption that icy satellites and asteroids contain rock compositionally similar to carbonaceous chondrites suggests that salts may be important in the cryoigneous evolution of icy satellites and asteroids. Ordinary chondrites, an alternative rock component of icy satellites, lack abundant salts, although their anhydrous silicate assemblages are unstable with respect to water and would react to produce salts upon initial melting of ice. Some basic physiochemical properties are reviewed of likely cryovolcanic brines and how the existence of soluble salts in Ganymede might affect its structure and evolution is considered. Observations indicate late stage (post heavy bombardment and post tectonic) volcanism on Ganymede. The highly fluid character of Ganymedian volcanism is consistent with extrusions of either water or salt water brines.

Kargel, Jeffrey S.↗

NASA NTRS1994

A random distribution reacting mixing layer model

A methodology for simulation of molecular mixing and the resulting velocity and temperature fields has been developed. The ideas are applied to the flow conditions present in the NASA Lewis Planar Reacting Shear Layer (PRSL) facility, and results compared to experimental data. A gaussian transverse turbulent velocity distribution is used in conjunction with a linearly increasing time scale to describe the mixing of different regions of the flow. Equilibrium reaction calculations are then performed on the mix to arrive at a new species composition and temperature. Velocities are determined through summation of momentum contributions. The analysis indicates a combustion efficiency of the order of 80 percent for the reacting mixing layer, and a turbulent Schmidt number of 2/3. The success of the model is attributed to the simulation of large-scale transport of fluid. The favorable comparison shows that a relatively quick and simple PC calculation is capable of simulating the basic flow structure in the reacting and non-reacting shear layer present in the facility given basic assumptions about turbulence properties.

Jones, Richard A.↗

NASA NTRS1994

Recycling of volatiles at subduction zones: Noble gas evidence from the Tabar-Lihir-Tanga-Feni arc of papua New Guinea

Convergent margin processes play an important but poorly understood role in the distribution of terrestrial volatile species. For example, subduction processes filter volatiles from the subducting package, thereby restricting their return to the mantle. In addition, once extracted from the downgoing slab, volatiles become an essential component in the petrogenesis of island arc magmas. The noble gases, with their systematic variation in physical properties and diversity of radiogenic isotopes, should carry a uniquely valuable record of these processes. However, thus far studies of noble gases in arc volcanics have achieved only limited success in this regard. Subduction-related lavas and geothermal fluids carry (3)He/(4)He ratios equal to or slightly lower than those found in the depleted upper mantle source of mid-ocean ridge basalts. Apparently slab-derived helium (which should have (3)He/(4)He much less than MORB) is extensively diluted by MORB-like helium from the mantle wedge, making it difficult to use helium as a tracer of convergent margin processes. Interpretation of the heavier noble gases (Ne-Ar-Kr-Xe) in arc lavas has also proven difficult, because the lavas carry low noble gas concentrations and hence are subject to pervasive atmospheric contamination. The low noble gas concentrations may be a consequence of degassing in the high level magma chambers characteristic of arc stratovolcanos. We have recently initiated a project to better constrain the behavior of volatiles in subduction zones through geochemical studies of the tectonically unusual volcanoes of the Tabar-Lihir-Tanga-Feni (TLTF) arc in the Bismarck Archipelago, Papua New Guinea.

Farley, Kenneth↗

NASA NTRS2003

Crystal Growth of ZnSe and Related Ternary Compound Semiconductors by Vapor Transport

The objective of the project is to determine the relative contributions of gravity-driven fluid flows to the compositional distribution, incorporation of impurities and defects, and deviation from stoichiometry observed in the crystals grown by vapor transport as results of buoyancy-driven convection and growth interface fluctuations caused by irregular fluid-flows. ZnSe and related ternary compounds, such as ZnSeS and ZnSeTe, were grown by vapor transport technique with real time in situ non-invasive monitoring techniques. The grown crystals were characterized extensively to correlate the grown crystal properties with the growth conditions. Significant effects of gravity vector orientation on the growth crystal morphology and point defect distribution were observed.

Source record↗

NASA NTRS1992

Effect of oral calcium and calcium + fluoride treatments on mouse bone properties during suspension

The bone effects of oral dosages of calcium chloride with or without supplementary sodium fluoride were assessed in antiorthostatically suspended mice. Two calcium dosages were used to replace half (3.1 mM) or all(6.3 mM) of the dietary calcium lost due to reduced food intake by the suspended mice. Two groups of 6.3 mM CaCl2-treated mice were additionally treated with 0.25 or 2.5 mM NaF. The results indicate that supplementation of the mouse drinking water with calcium salts prevents bone changes induced by short-term suspension, while calcium salts in combination with fluoride are less effective as fluoride dosage increases. However, the calcium supplements change the relationship between the femur mechanical properties and the mineral composition of the bone. Because of this, it appears that oral calcium supplements are effective through a mechanism other than simple dietary supplementation and may indicate a dependence of bone consistency on systemic and local fluid conditions.

Non-NASA Center↗

NASA NTRS2004

Scalable High Performance Computing: Direct and Large-Eddy Turbulent Flow Simulations Using Massively Parallel Computers

This final report contains reports of research related to the tasks "Scalable High Performance Computing: Direct and Lark-Eddy Turbulent FLow Simulations Using Massively Parallel Computers" and "Devleop High-Performance Time-Domain Computational Electromagnetics Capability for RCS Prediction, Wave Propagation in Dispersive Media, and Dual-Use Applications. The discussion of Scalable High Performance Computing reports on three objectives: validate, access scalability, and apply two parallel flow solvers for three-dimensional Navier-Stokes flows; develop and validate a high-order parallel solver for Direct Numerical Simulations (DNS) and Large Eddy Simulation (LES) problems; and Investigate and develop a high-order Reynolds averaged Navier-Stokes turbulence model. The discussion of High-Performance Time-Domain Computational Electromagnetics reports on five objectives: enhancement of an electromagnetics code (CHARGE) to be able to effectively model antenna problems; utilize lessons learned in high-order/spectral solution of swirling 3D jets to apply to solving electromagnetics project; transition a high-order fluids code, FDL3DI, to be able to solve Maxwell's Equations using compact-differencing; develop and demonstrate improved radiation absorbing boundary conditions for high-order CEM; and extend high-order CEM solver to address variable material properties. The report also contains a review of work done by the systems engineer.

Morgan, Philip E.↗

NASA NTRS2026

Spacecraft Radiator Protection from Ionizing Radiation, Dust, and Excessive Heat Loss

Under a Phase II SBIR project funded by NASA Johnson Space Center (Contract No. 80NSSC25C0088), Analytical Scientific Products LLC (ASP) has been developing an actively controlled louver to protect spacecraft radiators from degradation due to exposure to various types of environmental conditions. Of particular interest are ionizing radiation during spacecraft transit through the Van Allen belts, dust during spacecraft landing and surface operations on the moon, and excessive heat loss during the long lunar night especially near the poles where the local ambient temperatures can drop below -200°C. Exposure to ionizing radiation and dust can degrade the optical properties of the radiator coating that in turn reduces its ability to reject excess heat from the spacecraft into the environment. Exposure to the extremely low temperature conditions during the lunar night near the poles can freeze the radiator fluids that can compromise the integrity of the radiator. Passive louvers constructed using thick aluminum vanes are used currently to protect spacecraft radiators from some of the above effects, but they have several drawbacks: (i) their high aerial density makes it impractical to scale them to protect the much larger spacecraft radiators needed for future manned missions to the moon and beyond, (ii) the bimetallic actuators used to open and close the vanes in passive louvers rely on external temperature alone and so the louver cannot offer protection from dust and ionizing radiation when the ambient temperatures are high, and (iii) the bimetallic actuators need time scales of the order of hours to open and close. We have addressed the above limitations of passive louvers through the development of a low aerial density and rapidly actuating actively controlled louver. It is constructed from a custom alloy that offers much higher levels of protection against the ionizing radiation, dust and excessive heat loss at a fraction of weight compared to passive louvers while providing opening and closing time scales of the order of a second. Our modular design makes it easy to scale the system up or down depending on the application. We have recently constructed a 31 inch × 31 inch module of this louver and tested its functionality and effectiveness. These tests have shown that the louver can be opened and closed in less than 2 s. It reduces the transmission of ionizing radiation by 78% to 100% (depending on the radiation source), dust transmission by more than 93%, and heat loss by more than 97%. We are currently preparing to evaluate this device at the Johnson Space Center’s cryogenic vacuum chamber where it can be subjected to simulated lunar surface conditions. This paper provides a detailed discussion of the test designs as well as the data from tests that demonstrate the effectiveness in reducing the transmission of ionizing radiation, dust, and heat under laboratory conditions.

Radiator↗

NASA NTRSDate not supplied

Spacecraft Radiator Protection from Ionizing Radiation, Dust, and Excessive Heat Loss

Radiator↗

NASA NTRS1992

Theoretical study of thermodynamic properties and reaction rates of importance in the high-speed research program

One of the primary goals of NASA's high-speed research program is to determine the feasibility of designing an environmentally safe commercial supersonic transport airplane. The largest environmental concern is focused on the amount of ozone destroying nitrogen oxides (NO(x)) that would be injected into the lower stratosphere during the cruise portion of the flight. The limitations placed on NO(x) emission require more than an order of magnitude reduction over current engine designs. To develop strategies to meet this goal requires first gaining a fundamental understanding of the combustion chemistry. To accurately model the combustor requires a computational fluid dynamics approach that includes both turbulence and chemistry. Since many of the important chemical processes in this regime involve highly reactive radicals, an experimental determination of the required thermodynamic data and rate constants is often very difficult. Unlike experimental approaches, theoretical methods are as applicable to highly reactive species as stable ones. Also our approximation of treating the dynamics classically becomes more accurate with increasing temperature. In this article we review recent progress in generating thermodynamic properties and rate constants that are required to understand NO(x) formation in the combustion process. We also describe our one-dimensional modeling efforts to validate an NH3 combustion reaction mechanism. We have been working in collaboration with researchers at LeRC, to ensure that our theoretical work is focused on the most important thermodynamic quantities and rate constants required in the chemical data base.

Langhoff, Stephen↗

NASA NTRS2007

Mesh Deformation Based on Fully Stressed Design: The Method and Two-Dimensional Examples

Mesh deformation in response to redefined boundary geometry is a frequently encountered task in shape optimization and analysis of fluid-structure interaction. We propose a simple and concise method for deforming meshes defined with three-node triangular or four-node tetrahedral elements. The mesh deformation method is suitable for large boundary movement. The approach requires two consecutive linear elastic finite-element analyses of an isotropic continuum using a prescribed displacement at the mesh boundaries. The first analysis is performed with homogeneous elastic property and the second with inhomogeneous elastic property. The fully stressed design is employed with a vanishing Poisson s ratio and a proposed form of equivalent strain (modified Tresca equivalent strain) to calculate, from the strain result of the first analysis, the element-specific Young s modulus for the second analysis. The theoretical aspect of the proposed method, its convenient numerical implementation using a typical linear elastic finite-element code in conjunction with very minor extra coding for data processing, and results for examples of large deformation of two-dimensional meshes are presented in this paper. KEY WORDS: Mesh deformation, shape optimization, fluid-structure interaction, fully stressed design, finite-element analysis, linear elasticity, strain failure, equivalent strain, Tresca failure criterion

Hsu, Su-Yuen↗

NASA NTRS1970

The dynamics of oceanic fronts. Part 1: The Gulf Stream

The establishment and maintenance of the mean hydrographic properties of large scale density fronts in the upper ocean is considered. The dynamics is studied by posing an initial value problem starting with a near surface discharge of buoyant water with a prescribed density deficit into an ambient stationary fluid of uniform density. The full time dependent diffusion and Navier-Stokes equations for a constant Coriolis parameter are used in this study. Scaling analysis reveals three independent length scales of the problem, namely a radius of deformation or inertial length scale, Lo, a buoyance length scale, ho, and a diffusive length scale, hv. Two basic dimensionless parameters are then formed from these length scales, the thermal (or more precisely, the densimetric) Rossby number, Ro = Lo/ho and the Ekman number, E = hv/ho. The governing equations are then suitably scaled and the resulting normalized equations are shown to depend on E alone for problems of oceanic interest. Under this scaling, the solutions are similar for all Ro. It is also shown that 1/Ro is a measure of the frontal slope. The governing equations are solved numerically and the scaling analysis is confirmed. The solution indicates that an equilibrium state is established. The front can then be rendered stationary by a barotropic current from a larger scale along-front pressure gradient. In that quasisteady state, and for small values of E, the main thermocline and the inclined isopycnics forming the front have evolved, together with the along-front jet. Conservation of potential vorticity is also obtained in the light water pool. The surface jet exhibits anticyclonic shear in the light water pool and cyclonic shear across the front.

Kao, T. W.↗

NASA NTRS1991

Planetary origin, evolution, and structure

Three areas of recent and ongoing research are presented. The first area is giant planet heatflows. Conventional wisdom attributes the heatflow of the giant planets to the gradual loss of primordial heat, except in the case of Saturn where helium separation is evidently occurring. There are two problems with this picture: (1) the observed helium abundance of Saturn's atmosphere is so low that Jupiter must also be differentiating helium since its internal entropy cannot be much higher than Saturn; and (2) the heatflow of Neptune (not to mention Uranus) is too high to be consistent with adiabatic cooling from an initial hot state. A self-consistent solution to these two problems is presented. The second area covered is that of the despinning protogiant planets. Modeling of the possible despinning of these protoplanets by hydromagnetic torques was performed and the model results are discussed. The third area covered is how Titan hides its ocean. Until recently, the favored picture of Titan's surface was a roughly kilometer-thick ethane/methane ocean, presumably global in extent with at most a few outcroppings of dry land. The depth of the ocean is well constrained by observed atmospheric properties, and the constraints on subaerial topography are obtained indirectly from tidal considerations. A different picture of Titan's surface was pursued which was motivated by the perspective that methane on Titan should more properly be considered as a magmatic fluid. In this picture, methane is stored subsurface in magma chambers fed from deep-seated sources of methane, most probably due to the high pressure breakdown of methane clathrate. Other aspects of this model of Titan are presented.

Stevenson, David J.↗