SEARCH · Engineering Papers

Results for “thermal diffusion”

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,549 records · Page 86

All sources NASA NTRS DOE OSTI

NASA NTRS2018

NASA Plum Brook Station In-Space Propulsion Facility Test Stand Characterization Hot Fire Test

A test facility modification to enable small scale altitude propulsion testing at the NASA Glenn Research Center's In-Space Propulsion (ISP) Facility was verified with a hot fire test campaign. As the facility's primary steam supply system undergoes refurbishment, the alternate facility configuration, known as the "vacuum accumulator" mode, would enable rocket engine testing up to 10,000 lbf thrust. The NASA Johnson Space Center developed the vehicle for the verification test campaign: the Integrated Cryogenic Propulsion Test Article (ICPTA). Constructed primarily from assets of the former Morpheus Project, the ICPTA provided an integrated liquid oxygen (LOX) / liquid methane (LCH4) propulsion system including a 2,800 lbf thrust main engine. The ISP Facility's vacuum accumulator configuration leveraged the large test volume of the facility and a diffuser insert to maintain altitude conditions. During hot fire, the ICPTA main engine "started" the diffuser insert constructed for the test campaign. As a result, the test chamber upstream of the diffuser insert remained at altitude conditions throughout the hot fire. Upon engine shut down, a backflow deflector mitigated blow back into the test chamber by restricting the mass flow and redirecting it away from the test article. The test campaign successfully characterized the performance of the vacuum accumulator configuration. In addition, it provided an opportunity to collect data for an integrated LOX / LCH4 propulsion system in an altitude and thermal vacuum environment.

Rocket Engine Testing↗

NASA NTRS1980

Advanced methods for preparation and characterization of infrared detector materials

Mercury cadmium telluride crystals were prepared by the Bridgman method with a wide range of crystal growth rates and temperature gradients adequate to prevent constitutional supercooling under diffusion-limited, steady state, growth conditions. The longitudinal compositional gradients for different growth conditions and alloy compositions were calculated and compared with experimental data to develop a quantitative model of the crystal growth kinetics for the Hg(i-x)CdxTe alloys, and measurements were performed to ascertain the effect of growth conditions on radial compositional gradients. The pseudobinary HgTe-CdTe constitutional phase diagram was determined by precision differential thermal analysis measurements and used to calculate the segregation coefficient of Cd as a function of x and interface temperature. Computer algorithms specific to Hg(1-x)CdxTe were developed for calculations of the charge carrier concentrations, charge carrier mobilities, Hall coefficient, optical absorptance, and Fermi energy as functions of x, temperature, ionized donor and acceptor concentrations, and neutral defect concentrations.

Lehoczky, S. L.↗

NASA NTRS1981

Electronic properties of deep-level defects in proton irradiated AlGaAs-GaAs solar cells

Deep level transient spectroscopy and capacitance voltage techniques as well as analysis of the forward current voltage (I-V) characteristics and SEM-EIC data were carried out for proton irradiated GaAs solar cells over a wide range of proton energies and proton fluences. Defect and recombination parameters such as defect energy levels and density, carrier capture cross sections and lifetimes as well as diffusion lengths in the undoped n-GaAs LPE layers were determined. Good correlation between these defect parameters and solar cell performance parameters was obtained for GaAs solar cells irradiated by 200 and 290 KeV protons. It was found that 200 to 290 KeV protons will produce the most defects and damages to the GaAs solar cell structure used. The influence of the low temperature (200 to 400 C) periodic thermal annealing on the deep level defects and the performance of the 200 KeV proton irradiated cells is discussed.

Li, S. S.↗

NASA NTRS1994

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields

The primary accomplishments of the project were as follows: (1) From an overall standpoint, the primary accomplishment of this research was the development of a complete gasdynamic-radiatively coupled nonequilibrium viscous shock layer solution method for axisymmetric blunt bodies. This method can be used for rapid engineering modeling of nonequilibrium re-entry flowfields over a wide range of conditions. (2) Another significant accomplishment was the development of an air radiation model that included local thermodynamic nonequilibrium (LTNE) phenomena. (3) As part of this research, three electron-electronic energy models were developed. The first was a quasi-equilibrium electron (QEE) model which determined an effective free electron temperature and assumed that the electronic states were in equilibrium with the free electrons. The second was a quasi-equilibrium electron-electronic (QEEE) model which computed an effective electron-electronic temperature. The third model was a full electron-electronic (FEE) differential equation model which included convective, collisional, viscous, conductive, vibrational coupling, and chemical effects on electron-electronic energy. (4) Since vibration-dissociation coupling phenomena as well as vibrational thermal nonequilibrium phenomena are important in the nonequilibrium zone behind a shock front, a vibrational energy and vibration-dissociation coupling model was developed and included in the flowfield model. This model was a modified coupled vibrational dissociation vibrational (MCVDV) model and also included electron-vibrational coupling. (5) Another accomplishment of the project was the usage of the developed models to investigate radiative heating. (6) A multi-component diffusion model which properly models the multi-component nature of diffusion in complex gas mixtures such as air, was developed and incorporated into the blunt body model. (7) A model was developed to predict the magnitude and characteristics of the shock wave precursor ahead of vehicles entering the Earth's atmosphere. (8) Since considerable data exists for radiating nonequilibrium flow behind normal shock waves, a normal shock wave version of the blunt body code was developed. (9) By comparing predictions from the models and codes with available normal shock data and the flight data of Fire II, it is believed that the developed flowfield and nonequilibrium radiation models have been essentially validated for engineering applications.

Carlson, Leland A.↗

NASA NTRS2003

A 2-D Piston Effect Solution for the Relaxation in the MISTE Cell

When a large density stratification is no longer a problem in a microgravity environment, one would like to increase the sample size in order to increase the signal-to-noise ratio for a specific heat measurement. To reduce the equilibration time associated with the large sample size, we designed a cylindrical cell containing a stack of plates that separate the bulk fluid into 60 equally thin layers. To understand the thermal behavior of the whole cell, we analyzed the thermal behavior of a 2-D composite system of a cylindrical near-critical fluid layer in contact with a cylindrical copper plate. In this 2-D analysis, the circumference boundary of the two cylindrical layers is subjected to a step temperature change. The solution of this 2-D composite system includes the piston effect that speeds up the equilibration in the near-critical fluid layer and the pure diffusion in the copper plate. The results of this analysis indicate that the characteristic length for the equilibration of the stacked cell is determined by an effective thickness of a single fluid layer instead of the total height of the cylindrical cell.

Weilert, Mark↗

NASA NTRS2011

MODIS On-orbit Calibration Uncertainty Assessment

MODIS has 20 reflective solar bands (RSB) and 16 thermal emissive bands (TEB). Compared to its heritage sensors, MODIS was developed with very stringent calibration uncertainty requirements. As a result, MODIS was designed and built with a set of on-board calibrators (OBC), which allow key sensor performance parameters and on-orbit calibration coefficients to be monitored and updated. In terms of its calibration traceability, MODIS RSB calibration is reflectance based using an on-board solar diffuser (SD) and the TEB calibration is radiance based using an on-board blackbody (BB). In addition to on-orbit calibration coefficients derived from its OBC, calibration parameters determined from sensor pre-launch calibration and characterization are used in both the RSB and TEB calibration and retrieval algorithms. This paper provides a brief description of MODIS calibration methodologies and an in-depth analysis of its on-orbit calibration uncertainties. Also discussed in this paper are uncertainty contributions from individual components and differences due to Terra and Aqua MODIS instrument characteristics and on-orbit performance.

Chiang, Vincent↗

NASA NTRS1976

Similarity theory of the buoyantly interactive planetary boundary layer with entrainment

A similarity model is developed for the vertical profiles of turbulent flow variables in an entraining turbulent boundary layer of arbitrary buoyant stability. In the general formulation the vertical profiles, internal rotation of the velocity vector, discontinuities or jumps at a capping inversion and bulk aerodynamic coefficients of the boundary layer are given by solutions to a system of ordinary differential equations in the similarity variable. To close the system, a formulation for buoyantly interactive eddy diffusivity in the boundary layer is introduced which recovers Monin-Obukhov similarity near the surface and incorporates a hypothesis accounting for the observed variation of mixing length throughout the boundary layer. The model is tested in simplified versions which depend only on roughness, surface buoyancy, and Coriolis effects by comparison with planetary-boundary-layer wind- and temperature-profile observations, measurements of flat-plate boundary layers in a thermally stratified wind tunnel and observations of profiles of terms in the turbulent kinetic-energy budget of convective planetary boundary layers. On balance, the simplified model reproduced the trend of these various observations and experiments reasonably well, suggesting that the full similarity formulation be pursued further.

Hoffert, M. I.↗

NASA NTRS1987

Temperature characteristics of silicon space solar cells and underlying parameters

Silicon space cells, 2 cm x 2 cm, with 10 ohm-cm p-base resistivity, 8-mil base thickness, and no back-surface fields have been investigated over the temperature range from 301 to 223 K by measurements of dark forward and reverse current-voltage characteristics and current-voltage relations under illumination. From dark forward bias data, the first and second diode saturation currents, I01 and I02, are determined and hence the base diffusion length and lifetime of minority carriers as functions of temperature. Lifetime increases exponentially with temperature and is explained by a Shockley-Read-Hall model with deep recombination levels 0.245 eV above the valence band. The I02 variation with temperature follows the Sah-Noyce-Shockley-Choo model except at low temperature where extra transitions raise the value above the predicted level. Reverse bias current at low voltage is a thermally assisted tunneling process via deep levels which are observed in base recombination at higher temperatures. The tunneling effects tend to become independent of temperature in the low-temperature region. These results demonstrate the ability to deduce basic parameters such as lifetime from simple measurements and show that back-surface fields offer no advantage at temperatures below 230 K. The analysis also explains the fall in lifetimes observed as the base conductivity increases, attributing it to native defects (perhaps carbon-oxygen-vacancy complexes) rather than the concentration of base dopant.

Anspaugh, B. E.↗

NASA NTRS1995

ROSAT PSPC and HRI observations of the composite starburst/Seyfert 2 galaxy NGC 1672

The nearby barred spiral galaxy NGC 1672 has been observed with the Position Sensitive Proportional Counter (PSPC) and High Resolution Imager (HRI) instruments on board the ROSAT X-ray satellite. NGC 1672 is thought to have an obscured Seyfert nucleus, and it has strong starburst activity as well. Three bright X-ray sources with luminosities 1-2 x 10(exp 40) erg/s are clearly identified with NGC 1672. The strongest lies at the nucleus, and the other two lie at the ends of NGC 1672's prominent bar, locations that are also bright in H alpha and near-infrared images. The nuclear source is resolved by the HRI on about the scale of the recently identified nuclear ring, and one of the sources at the ends of the bar is also probably resolved. The X-ray spectrum of the nuclear source is quite soft, having a Raymond-Smith plasma temperature of approximately equals 0.7 keV and little evidence for intrinsic absorption. The ROSAT band X-ray flux of the nuclear source appears to be dominated not by X-ray binary emission but rather by diffuse gas emission. The absorption and emission properties of the sources, as well as their spatial extents, lead us to models of superbubbles driven by supernovae. However, the large density and emission measure of the nuclear X-ray source stretch the limits that can be comfortably accommodated by these models. We do not detect direct emission from the putative Seyfert nucleus, although an alternative model for the nuclear source is thermal emission from gas that is photoionized by a hidden Seyfert nucleus. The spectra of the other two X-ray sources are harder than that of the nuclear source, and have similar difficulties with regard to superbubble models.

Brandt, W. N.↗

NASA NTRS2001

Observations of O VI Emission from the Diffuse Interstellar Medium

We report the first Far Ultraviolet Spectroscopic Explorer (FUSE) measurements of diffuse O(VI) (lambda lambda 1032,1038) emission from the general diffuse interstellar medium outside of supernova remnants or superbubbles. We observed a 30 arcsec x 30 arcsec region of the sky centered at l = 315.0 deg and b = -41.3 deg. From the observed intensities (2930 +/- 290 (random) +/- 410 (systematic) and 1790 +/- 260 (random) +/- 250 (systematic) photons/sq cm/s/sr in the 1032 and 1038 angstrom emission lines, respectively), derived equations, and assumptions about the source location, we calculate the intrinsic intensity, electron density, thermal pressure, and emitting depth. The intensities are too large for the emission to originate solely in the Local Bubble. Thus, we conclude that the Galactic thick disk and lower halo also contribute. High velocity clouds are ruled out because there are none near the pointing direction. The calculated emitting depth is small, indicating that the O(VI)-bearing gas fills a small volume. The observations can also be used to estimate the cooling rate of the hot interstellar medium and constrain models. The data also yield the first intensity measurement of the C(II) 3s 2S(1/2) to 2p 2P(3/2) emission line at 1037 angstrom and place upper limits on the intensities of ultraviolet line emission from C(I), C(III), Si(II), S(III), S(IV), S(VI), and Fe(III).

Shelton, R. L.↗

NASA NTRS2010

Ablation, Thermal Response, and Chemistry Program for Analysis of Thermal Protection Systems

In previous work, the authors documented the Multicomponent Ablation Thermochemistry (MAT) and Fully Implicit Ablation and Thermal response (FIAT) programs. In this work, key features from MAT and FIAT were combined to create the new Fully Implicit Ablation, Thermal response, and Chemistry (FIATC) program. FIATC is fully compatible with FIAT (version 2.5) but has expanded capabilities to compute the multispecies surface chemistry and ablation rate as part of the surface energy balance. This new methodology eliminates B' tables, provides blown species fractions as a function of time, and enables calculations that would otherwise be impractical (e.g. 4+ dimensional tables) such as pyrolysis and ablation with kinetic rates or unequal diffusion coefficients. Equations and solution procedures are presented, then representative calculations of equilibrium and finite-rate ablation in flight and ground-test environments are discussed.

Milos, Frank S.↗

NASA NTRS2009

Multiscale Computer Simulation of Tensile and Compressive Strain in Polymer- Coated Silica Aerogels

While the low thermal conductivities of silica aerogels have made them of interest to the aerospace community as lightweight thermal insulation, the application of conformal polymer coatings to these gels increases their strength significantly, making them potentially useful as structural materials as well. In this work we perform multiscale computer simulations to investigate the tensile and compressive strain behavior of silica and polymer-coated silica aerogels. Aerogels are made up of clusters of interconnected particles of amorphous silica of less than bulk density. We simulate gel nanostructure using a Diffusion Limited Cluster Aggregation (DLCA) procedure, which produces aggregates that exhibit fractal dimensions similar to those observed in real aerogels. We have previously found that model gels obtained via DLCA exhibited stress-strain curves characteristic of the experimentally observed brittle failure. However, the strain energetics near the expected point of failure were not consistent with such failure. This shortcoming may be due to the fact that the DLCA process produces model gels that are lacking in closed-loop substructures, compared with real gels. Our model gels therefore contain an excess of dangling strands, which tend to unravel under tensile strain, producing non-brittle failure. To address this problem, we have incorporated a modification to the DLCA algorithm that specifically produces closed loops in the model gels. We obtain the strain energetics of interparticle connections via atomistic molecular statics, and abstract the collective energy of the atomic bonds into a Morse potential scaled to describe gel particle interactions. Polymer coatings are similarly described. We apply repeated small uniaxial strains to DLCA clusters, and allow relaxation of the center eighty percent of the cluster between strains. The simulations produce energetics and stress-strain curves for looped and nonlooped clusters, for a variety of densities and interaction parameters.

Good, Brian↗

NASA NTRS1975

Further numerical calculations of the circulation of the atmosphere of Venus

Results of two-dimensional simulations of the deep circulation of Venus are presented. They prove that the high surface temperature can only be explained by the greenhouse effect, and that Goody and Robinson's (1966) dynamical model is not valid. Analytical radiative equilibrium solutions for a semi-grey atmosphere, both with and without an internal heat source, are presented. It is shown that the greenhouse effect is sufficient to produce the high surface temperature if the thermal optical depth is much greater than 100 and the ratio of solar to thermal optical depths does not exceed 0.005. A two-dimensional version of a three-dimensional model is used to test the validity of the new mechanism proposed by Gierasch (1975) to explain the 4-day circulation. Numerical experiments with large horizontal viscosities failed to show strong zonal velocities even for the case of large Prandtl numbers. It is observed that the dissipation of angular momentum introduced by the strong horizontal diffusion more than compensates for the upward transport of angular momentum due to the Hadley cell.

De Rivas, E. K.↗

NASA NTRS1992

Results of a zonally truncated three-dimensional model of the Venus middle atmosphere

Although the equatorial rotational speed of the solid surface of Venus is only 4 m s(exp-1), the atmospheric rotational speed reaches a maximum of approximately 100 m s(exp-1) near the equatorial cloud top level (65 to 70 km). This phenomenon, known as superrotation, is the central dynamical problem of the Venus atmosphere. We report here the results of numerical simulations aimed at clarifying the mechanism for maintaining the equatorial cloud top rotation. Maintenance of an equatorial rotational speed maximum above the surface requires waves or eddies that systematically transport angular momentum against its zonal mean gradient. The zonally symmetric Hadley circulation is driven thermally and acts to reduce the rotational speed at the equatorial cloud top level; thus wave or eddy transport must counter this tendency as well as friction. Planetary waves arising from horizontal shear instability of the zonal flow (barotropic instability) could maintain the equatorial rotation by transporting angular momentum horizontally from midlatitudes toward the equator. Alternatively, vertically propagating waves could provide the required momentum source. The relative motion between the rotating atmosphere and the pattern of solar heating, which as a maximum where solar radiation is absorbed near the cloud tops, drives diurnal and semidiurnal thermal tides that propagate vertically away from the cloud top level. The effect of this wave propagation is to transport momentum toward the cloud top level at low latitudes and accelerate the mean zonal flow there. We employ a semispectral primitive equation model with a zonal mean flow and zonal wavenumbers 1 and 2. These waves correspond to the diurnal and semidiurnal tides, but they can also be excited by barotropic or baroclinic instability. Waves of higher wavenumbers and interactions between the waves are neglected. Symmetry about the equator is assumed, so the model applies to one hemisphere and covers the altitude range 30 to 110 km. Horizontal resolution is 1.5 deg latitude, and vertical resolution is 1.5 km. Solar and thermal infrared heating, based on Venus observations and calculations drive the model flow. Dissipation is accomplished mainly by Rayleigh friction, chosen to produce strong dissipation above 85 km in order to absorb upward propagating waves and limit extreme flow velocities there, yet to give very weak Rayleigh friction below 70 km; results in the cloud layer do not appear to be sensitive to the Rayleigh friction. The model also has weak vertical diffusion, and very weak horizontal diffusion, which has a smoothing effect on the flow only at the two grid points nearest the pole.

Newman, M.↗

NASA NTRS1993

Centrifugally driven winds from protostellar disks. I - Wind model and thermal structure

The thermal structure of a wind that is centrifugally driven from the surface of a protostellar disk is studied. A generalized version of the Blandford and Payne self-similar wind model is introduced, and the temperature and ionization distributions in the outflow are investigated. For the evolution of atomic winds, the heat equation and the rate equations that describe the ionization and excitation state of hydrogen are solved self-consistently. Ambipolar diffusion is found to be a robust mechanism for heating the gas. In the more powerful outflows, molecular hydrogen is collisionally dissociated close to the disk surface, and hydrogen is mainly atomic within a few astronomical units from the central source. It is also demonstrated that these outflows have enough momentum to lift dust grains from the disk surface.

Safier, Pedro N.↗

NASA NTRS1998

Relationship of Topside Ionospheric Ion Outflows to Auroral Forms and Precipitation, Plasma Waves, and Convection Observed by Polar

The POLAR satellite often observes upflowing ionospheric ions (UFIs) in and near the aurora] oval on southern perigee (approx. 5000 km altitude) passes. We present the UFI features observed by the thermal ion dynamics experiment (TIDE) and the toroidal imaging mass angle spectrograph (TIMAS) in the dusk-dawn sector under two different geomagnetic activity conditions in order to elicit their relationships with auroral forms, wave emissions, and convection pattern from additional POLAR instruments. During the active interval, the ultraviolet imager (UVI) observed a bright discrete aurora on the duskside after the substorm onset and then observed a small isolated aurora form and diffuse auroras on the dawnside during the recovery phase. The UFIs showed clear conic distributions when the plasma wave instrument (PWI) detected strong broadband wave emissions below approx. 10 kHz, while no significant auroral activities were observed by UVI. At higher latitudes, the low-energy UFI conics gradually changed to the polar wind component with decreasing intensity of the broadband emissions. V-shaped auroral kilometric radiation (AKR) signatures observed above -200 kHz by PWI coincided with the region where the discrete aurora and the UFI beams were detected. The latitude of these features was lower than that of the UFI conics. During the observations of the UFI beams and conics, the lower-frequency fluctuations observed by the electric field instrument were also enhanced, and the convection directions exhibited large fluctuations. It is evident that large electrostatic potential drops produced the precipitating electrons and discrete auroras, the UFI beams, and the AKR, which is also supported by the energetic plasma data from HYDRA. Since the intense broadband emissions were also observed with the UFIs, the ionospheric ions could be energized transversely before or during the parallel acceleration due to the potential drops.

Hirahara, M.↗

NASA NTRS1999

Modelling of Spherical Gas Bubble Oscillations and Sonoluminescence

The discovery of single-bubble sonoluminescence has led to a renewed interest in the forced radial oscillations of gas bubbles. Many of the more recent studies devoted to this topic have used several simplifications in the modelling, and in particular in accounting for liquid compressibility and thermal processes in the bubble. In this paper the significance of these simplifications is explored by contrasting the results of Lohse and co-workers with those of a more detailed model. It is found that, even though there may be little apparent difference between the radius-versus time behaviour of the bubble as predicted by the two models, quantities such as the spherical stability boundary and the threshold for rectified diffusion are affected in a quantitatively significant way. These effects are a manifestation of the subtle dependence upon dissipative processes of the phase of radial motion with respect to the driving sound field. The parameter space region, where according to the theory of Lohse and co-workers, sonoluminescence should be observable, is recalculated with the new model and is found to be enlarged with respect to the earlier estimate. The dependence of this parameter region on sound frequency is also illustrated.

Prosperetti, A.↗

NASA NTRS2003

CO2 Acquisition Membrane (CAM)

The objective of CAM is to develop, test, and analyze thin film membrane materials for separation and purification of carbon dioxide (CO2) from mixtures of gases, such as those found in the Martian atmosphere. The membranes are targeted toward In Situ Resource Utilization (ISRU) applications that will operate in extraterrestrial environments and support future unmanned and human space missions. A primary application is the Sabatier Electrolysis process that uses Mars atmosphere CO2 as raw material for producing water, oxygen, and methane for rocket fuel and habitat support. Other applications include use as an inlet filter to collect and concentrate Mars atmospheric argon and nitrogen gases for habitat pressurization, and to remove CO2 from breathing gases in Closed Environment Life Support Systems (CELSS). CAM membrane materials include crystalline faujasite (FAU) zeolite and rubbery polymers such as silicone rubber (PDMS) that have been shown in the literature and via molecular simulation to favor adsorption and permeation of CO2 over nitrogen and argon. Pure gas permeation tests using commercial PDMS membranes have shown that both CO2 permeance and the separation factor relative to other gases increase as the temperature decreases, and low (Delta)P(Sub CO2) favors higher separation factors. The ideal CO2/N2 separation factor increases from 7.5 to 17.5 as temperature decreases from 22 C to -30 C. For gas mixtures containing CO2, N2, and Ar, plasticization decreased the separation factors from 4.5 to 6 over the same temperature range. We currently synthesize and test our own Na(+) FAU zeolite membranes using standard formulations and secondary growth methods on porous alumina. Preliminary tests with a Na(+) FAU membrane at 22 C show a He/SF6 ideal separation factor of 62, exceeding the Knudsen diffusion selectivity by an order of magnitude. This shows that the membrane is relatively free from large defects and associated non-selective (viscous flow) transport mechanisms. The Membrane Test Facility (MTF) has been developed to measure membrane permeance over a wide range of temperature and pressure. The facility uses two volume compartments separated by the membrane that are instrumented to measure temperature, delta pressure across the membrane, and gas composition. A thermal shroud supports and encloses the membrane, and provides temperature control. Methods were developed to determine membrane permeance using the first order decay of the pressure difference between the sealed compartments, using the total pressure for pure gases, and partial pressure of each species in gas mixtures. The technique provides an end-to-end measurement of gas permeance that includes concentration polarization effects. Experiments have shown that in addition to membrane permeance properties, the geometry and design of associated structures play an important role in how membrane systems will function on Mars.

Mason, Larry W.↗