Balasubramaniam, R. — Page 2

NASA NTRS2009

Carbothermal Processing of Lunar Regolith Using Methane

The processing of lunar regolith for the production of oxygen is a key component of the In-Situ Resource Utilization plans currently being developed by NASA. Among various candidate processes, the modeling of oxygen production by hydrogen reduction, molten salt electrolysis, and carbothermal processing are presently being pursued. In the carbothermal process, a portion of the surface of the regolith in a container is heated by exposure to a heat source such as a laser beam or a concentrated solar heat flux, so that a small zone of molten regolith is established. The molten zone is surrounded by solid regolith particles that are poor conductors of heat. A continuous flow of methane is maintained over the molten regolith zone. Our model is based on a mechanism where methane pyrolyzes when it comes in contact with the surface of the hot molten regolith to form solid carbon and hydrogen gas. Carbon is deposited on the surface of the melt, and hydrogen is released into the gas stream above the melt surface. We assume that the deposited carbon mixes in the molten regolith and reacts with metal oxides in a reduction reaction by which gaseous carbon monoxide is liberated. Carbon monoxide bubbles through the melt and is released into the gas stream. Oxygen is produced subsequently by (catalytically) processing the carbon monoxide downstream. In this paper, we discuss the development of a chemical conversion model of the carbothermal process to predict the rate of production of carbon monoxide.

Balasubramaniam, R.↗

NASA NTRS2009

Development and Validation of a Model for Hydrogen Reduction of JSC-1A

Hydrogen reduction of lunar regolith has been proposed as a viable technology for oxygen production on the moon. Hydrogen reduces FeO present in the lunar regolith to form metallic iron and water. The water may be electrolyzed to recycle the hydrogen and produce oxygen. Depending upon the regolith composition, FeO may be bound to TiO2 as ilmenite or it may be dispersed in glassy substrates. Some testing of hydrogen reduction has been conducted with Apollo-returned lunar regolith samples. However, due to the restricted amount of lunar material available for testing, detailed understanding and modeling of the reduction process in regolith have not yet been developed. As a step in this direction, hydrogen reduction studies have been carried out in more detail with lunar regolith simulants such as JSC-1A by NASA and other organizations. While JSC-1A has some similarities with lunar regolith, it does not duplicate the wide variety of regolith types on the moon, for example, it contains almost no ilmenite. Nonetheless, it is a good starting point for developing an understanding of the hydrogen reduction process with regolith-like material. In this paper, a model utilizing a shrinking core formulation coupled with the reactor flow is described and validated against experimental data on hydrogen reduction of JSC-1A.

Hegde, U.↗

NASA NTRS2009

Development of a Reactor Model for Chemical Conversion of Lunar Regolith

Lunar regolith will be used for a variety of purposes such as oxygen and propellant production and manufacture of various materials. The design and development of chemical conversion reactors for processing lunar regolith will require an understanding of the coupling among the chemical, mass and energy transport processes occurring at the length and time scales of the overall reactor with those occurring at the corresponding scales of the regolith particles. To this end, a coupled transport model is developed using, as an example, the reduction of ilmenite-containing regolith by a continuous flow of hydrogen in a flow-through reactor. The ilmenite conversion occurs on the surface and within the regolith particles. As the ilmenite reduction proceeds, the hydrogen in the reactor is consumed, and this, in turn, affects the conversion rate of the ilmenite in the particles. Several important quantities are identified as a result of the analysis. Reactor scale parameters include the void fraction (i.e., the fraction of the reactor volume not occupied by the regolith particles) and the residence time of hydrogen in the reactor. Particle scale quantities include the time for hydrogen to diffuse into the pores of the regolith particles and the chemical reaction time. The paper investigates the relationships between these quantities and their impact on the regolith conversion. Application of the model to various chemical reactor types, such as fluidized-bed, packed-bed, and rotary-bed configurations, are discussed.

Hegde, U.↗

NASA NTRS2009

Analysis of Thermal and Reaction Times for Hydrogen Reduction of Lunar Regolith

System analysis of oxygen production by hydrogen reduction of lunar regolith has shown the importance of the relative time scales for regolith heating and chemical reaction to overall performance. These values determine the sizing and power requirements of the system and also impact the number and operational phasing of reaction chambers. In this paper, a Nusselt number correlation analysis is performed to determine the heat transfer rates and regolith heat up times in a fluidized bed reactor heated by a central heating element (e.g., a resistively heated rod, or a solar concentrator heat pipe). A coupled chemical and transport model has also been developed for the chemical reduction of regolith by a continuous flow of hydrogen. The regolith conversion occurs on the surfaces of and within the regolith particles. Several important quantities are identified as a result of the above analyses. Reactor scale parameters include the void fraction (i.e., the fraction of the reactor volume not occupied by the regolith particles) and the residence time of hydrogen in the reactor. Particle scale quantities include the particle Reynolds number, the Archimedes number, and the time needed for hydrogen to diffuse into the pores of the regolith particles. The analysis is used to determine the heat up and reaction times and its application to NASA s oxygen production system modeling tool is noted.

Hegde, U.↗

NASA NTRS2008

Vapor-Gas Bubble Evolution and Growth in Extremely Viscous Fluids Under Vacuum

Formation of vapor and gas bubbles and voids is normal and expected in flow processes involving extremely viscous fluids in normal gravity. Practical examples of extremely viscous fluids are epoxy-like filler materials before the epoxy fluids cure to their permanent form to create a mechanical bond between two substrates. When these fluids flow with a free liquid interface exposed to vacuum, rapid bubble expansion process may ensue. Bubble expansion might compromise the mechanical bond strength. The potential sources for the origin of the gases might be incomplete out-gassing process prior to filler application; regasification due to seal leakage in the filler applicator; and/or volatiles evolved from cure reaction products formed in the hardening process. We embarked on a study that involved conducting laboratory experiments with imaging diagnostics in order to deduce the seriousness of bubbling caused by entrained air and volatile fluids under space vacuum and low gravity environment. We used clear fluids with the similar physical properties as the epoxy-like filler material to mimic the dynamics of bubbles. Another aspect of the present study was to determine the likelihood of bubbling resulting from dissolved gases nucleating from solution. These experimental studies of the bubble expansion are compared with predictions using a modified Rayleigh- Plesset equation, which models the bubble expansion.

Kizito, John↗

NASA NTRS2008

Carbothermal Processing of Lunar Regolith Using Methane

The processing of lunar regolith for the production of oxygen is a key component of the In-Situ Resource Utilization plans currently being developed by NASA. Among various candidate processes, the modeling of oxygen production by hydrogen reduction, molten salt electrolysis, and carbothermal processing are presently being pursued. In the carbothermal process, a portion of the surface of the regolith in a container is heated by exposure to a heat source such as a laser beam or a concentrated solar heat flux, so that a small zone of molten regolith is established. The molten zone is surrounded by solid regolith particles that are poor conductors of heat. A continuous flow of methane is maintained over the molten regolith zone. Our model is based on a mechanism where methane pyrolyzes when it comes in contact with the surface of the hot molten regolith to form solid carbon and hydrogen gas. Carbon is deposited on the surface of the melt, and hydrogen is released into the gas stream above the melt surface. We assume that the deposited carbon mixes in the molten regolith and reacts with metal oxides in a reduction reaction by which gaseous carbon monoxide is liberated. Carbon monoxide bubbles through the melt and is released into the gas stream. Oxygen is produced subsequently by (catalytically) processing the carbon monoxide downstream. In this paper, we discuss the development of a chemical conversion model of the carbothermal process to predict the rate of production of carbon monoxide.

Balasubramaniam, R.↗

NASA NTRS2008

Analysis of Thermal and Reaction Times for Hydrogen Reduction of Lunar Regolith

System analysis of oxygen production by hydrogen reduction of lunar regolith has shown the importance of the relative time scales for regolith heating and chemical reaction to overall performance. These values determine the sizing and power requirements of the system and also impact the number and operational phasing of reaction chambers. In this paper, a Nusselt number correlation analysis is performed to determine the heat transfer rates and regolith heat up times in a fluidized bed reactor heated by a central heating element (e.g., a resistively heated rod, or a solar concentrator heat pipe). A coupled chemical and transport model has also been developed for the chemical reduction of regolith by a continuous flow of hydrogen. The regolith conversion occurs on the surfaces of and within the regolith particles. Several important quantities are identified as a result of the above analyses. Reactor scale parameters include the void fraction (i.e., the fraction of the reactor volume not occupied by the regolith particles) and the residence time of hydrogen in the reactor. Particle scale quantities include the particle Reynolds number, the Archimedes number, and the time needed for hydrogen to diffuse into the pores of the regolith particles. The analysis is used to determine the heat up and reaction times and its application to NASA s oxygen production system modeling tool is noted.

Hegde, U.↗

NASA NTRS2007

Development of a Reactor Model for Chemical Conversion of Lunar Regolith

Lunar regolith will be used for a variety of purposes such as oxygen and propellant production and manufacture of various materials. The design and development of chemical conversion reactors for processing lunar regolith will require an understanding of the coupling among the chemical, mass and energy transport processes occurring at the length and time scales of the overall reactor with those occurring at the corresponding scales of the regolith particles. To this end, a coupled transport model is developed using, as an example, the reduction of ilmenite-containing regolith by a continuous flow of hydrogen in a flow-through reactor. The ilmenite conversion occurs on the surface and within the regolith particles. As the ilmenite reduction proceeds, the hydrogen in the reactor is consumed, and this, in turn, affects the conversion rate of the ilmenite in the particles. Several important quantities are identified as a result of the analysis. Reactor scale parameters include the void fraction (i.e., the fraction of the reactor volume not occupied by the regolith particles) and the residence time of hydrogen in the reactor. Particle scale quantities include the time for hydrogen to diffuse into the pores of the regolith particles and the chemical reaction time. The paper investigates the relationships between these quantities and their impact on the regolith conversion. Application of the model to various chemical reactor types, such as fluidized-bed, packed-bed, and rotary-bed configurations, are discussed.

Hedge, uday↗

NASA NTRS2006

Two Phase Flow Modeling: Summary of Flow Regimes and Pressure Drop Correlations in Reduced and Partial Gravity

The purpose of this report is to provide a summary of state-of-the-art predictions for two-phase flows relevant to Advanced Life Support. We strive to pick out the most used and accepted models for pressure drop and flow regime predictions. The main focus is to identify gaps in predictive capabilities in partial gravity for Lunar and Martian applications. Following a summary of flow regimes and pressure drop correlations for terrestrial and zero gravity, we analyze the fully developed annular gas-liquid flow in a straight cylindrical tube. This flow is amenable to analytical closed form solutions for the flow field and heat transfer. These solutions, valid for partial gravity as well, may be used as baselines and guides to compare experimental measurements. The flow regimes likely to be encountered in the water recovery equipment currently under consideration for space applications are provided in an appendix.

Balasubramaniam, R.↗

NASA NTRS2004

The formation of spikes in the displacement of miscible fluids

We report on experiments in which a more viscous fluid displaces a less viscous one in a vertical cylindrical tube. These experiments were performed using silicone oils in a vertical pipette of small diameter. The more viscous fluid also had a slightly larger density than the less viscous fluid. In the initial configuration, the fluids were at rest, and the interface was nominally flat. A dye was added to the more viscous fluid for ease of observation of the interface between the fluids. The flow was initiated by pumping the more viscous fluid into the less viscous one. The displacement velocity was such that the Reynolds number was smaller than unity and the Peclet number for mass transfer between the fluids was large compared to unity. For upward displacement of the more viscous fluid from an initially stable configuration, an axisymmetric finger was observed under all conditions. However, a needle-shaped spike was seen to propagate from the main finger in many cases, similar to that observed by Petitjeans and Maxworthy for the displacement of a more viscous fluid by a less viscous one.

Rheology/methods↗

NASA NTRS2004

Thermocapillary migration of a drop: an exact solution with Newtonian interfacial rheology and stretching/shrinkage of interfacial area elements for small Marangoni numbers

In this paper we analyze the effects of the following phenomena associated with the thermocapillary migration of a drop. The first is the influence of Newtonian surface rheology of the interface and the second is that of the energy changes associated with stretching and shrinkage of the interfacial area elements, when the drop is in motion. The former occurs because of dissipative processes in the interfacial region, such as when surfactant molecules are adsorbed at the interface in sufficient concentration. The interface is typically modeled in this instance by ascribing to it a surface viscosity. This is a different effect from that of interfacial tension gradients arising from surfactant concentration gradients. The stretching and shrinkage of interfacial area elements leads to changes in the internal energy of these elements that affects the transport of energy in the fluids adjoining the interface. When an element on the interface is stretched, its internal energy increases because of the increase in its area. This energy is supplied by the neighboring fluids that are cooled as a consequence. Conversely, when an element on the interface shrinks, the adjoining fluids are warmed. In the case of a moving drop, elements of interfacial area are stretched in the forward half of the drop, and are shrunk in the rear half. Consequently, the temperature variation on the surface of the drop and its migration speed are modified. The analysis of the motion of a drop including these effects was first performed by LeVan in 1981, in the limit when convective transport of momentum and energy are negligible. We extend the analysis of LeVan to include the convective transport of momentum by demonstrating that an exact solution of the momentum equation is obtained for an arbitrary value of the Reynolds number. This solution is then used to calculate the slightly deformed shape of the drop from a sphere.

Rheology/methods↗

NASA NTRS2003

Thermocapillary Convection Due to a Stationary Bubble - A Paradox

We analyze the velocity and temperature fields at steady state due to thermocapillary convection around a gas bubble that is stationary in a liquid. A linear temperature field is imposed in the undisturbed liquid. Our interest is in investigating the effect of convective transport of momentum and energy on the velocity and temperature fields. We assume the pertinent physical properties to be constant, and that buoyant convection is negligible. Suitably defined Reynolds and Marangoni numbers are assumed to be small compared with unity. When both the Reynolds and Marangoni numbers are set equal to zero, a solution can be found. In this solution, far from the bubble, the velocity field decays as the inverse of the distance from the bubble, and the disturbance temperature field decays as the inverse of the square of this distance. We now attempt to obtain a solution when the Reynolds number is zero, but the Marangoni number is small, but non-zero, by a perturbation expansion in the Marangoni number. When the temperature field is expanded in a regular perturbation series in the Marangoni number, we show that the problem for the first correction field is ill-posed. The governing equation for this perturbation field contains an inhomogeneity, and the corresponding particular solution neither decays far from the bubble, nor can be canceled by a homogeneous solution. Additional information is included in the original extended abstract.

Balasubramaniam, R.↗

NASA NTRS2002

Marangoni Effects of a Drop in an Extensional Flow: The Role of Surfactant Physical Chemistry

While the changes in stresses caused by surfactant adsorption on non-deforming interfaces have been fairly well established, prior to this work, there were few studies addressing how surfactants alter stresses on strongly deforming interfaces. We chose the model problem of a drop in a uniaxial extensional flow to study these stress conditions To model surfactant effects at fluid interfaces, a proper description of the dependence of the surface tension on surface concentration, the surface equation of state, is required. We have adopted a surface equation of state that accounts for the maximum coverage limit; that is, because surfactants have a finite cross sectional area, there is an upper bound to the amount of surfactant that can adsorb in a monolayer. The surface tension reduces strongly only when this maximum coverage is approached. Since the Marangoni stresses go as the derivative of the surface equation of state times the surface concentration gradient, the non-linear equation of state determines both the effect of surfactants in the normal stress jump, (which is balanced by the product of the mean curvature of the interface times the surface tension), and the tangential stress jump, which is balanced by Marangoni stresses. First, the effects of surface coverage and intermolecular interactions among surfactants which drive aggregation of surfactants in the interface were studied. (see Pawar and Stebe, Physics of Fluids).

Stebe, Kathleen J.↗

NASA NTRS2001

The Dynamics of Miscible Fluids: A Space Flight Experiment (MIDAS)

We propose a space flight experiment to study the dynamics of miscible interfaces. A less viscous fluid displaces one of higher viscosity within a tube. The two fluids are miscible in all proportions. An intruding "finger" forms that occupies a fraction of the tube. As time progresses diffusion at the interface combined with flow induced straining between the two fluids modifies the concentration and velocity distributions within the whole tube. Also, under such circumstances it has been proposed that the interfacial stresses could depend on the local concentration gradients (Korteweg stresses) and that the divergence of the velocity need not be zero, even though the flow is incompressible. We have obtained reasonable agreement for the tip velocity between numerical simulations (that ignored the Korteweg stress and divergence effects) and physical experiments only at high Peelet Numbers. However at moderate to low Pe agreement was poor. As one possibility we attributed this lack of agreement to the disregard of these effects. We propose a space experiment to measure the finger shape, tip velocity, and the velocity and concentration fields. From intercomparisons between the experiment and the calculations we can then extract values for the coefficients of the Korteweg stress terms and confirm or deny the importance of these stresses.

Maxworthy, T.↗

NASA NTRS2001

The Dynamics of Miscible Fluids: A Space Flight Experiment (MIDAS)

We propose a space flight experiment to study the dynamics of miscible interfaces. A less viscous fluid displaces one of higher viscosity within a tube. The two fluids are miscible in all proportions. An intruding "finger" forms that occupies a fraction of the tube. As time progresses diffusion at the interface combined with flow induced straining between the two fluids modifies the concentration and velocity distributions within the whole tube. Also, under such circumstances it has been proposed that the interfacial stresses could depend on the local concentration gradients (Korteweg stresses) and that the divergence of the velocity need not be zero, even though the flow is incompressible. We have obtained reasonable agreement for the tip velocity between numerical simulations (that ignored the Korteweg stress and divergence effects) and physical experiments only at high Peclet Numbers. However at moderate to low Pe agreement was poor. As one possibility we attributed this lack of agreement to the disregard of these effects. We propose a space experiment to measure the finger shape, tip velocity, and the velocity and concentration fields. From intercomparisons between the experiment and the calculations we can then extract values for the coefficients of the Korteweg stress terms and confirm or deny the importance of these stresses.

Maxworthy, T.↗

NASA NTRS2000

Quantum Corrections to the 'Atomistic' MOSFET Simulations

We have introduced in a simple and efficient manner quantum mechanical corrections in our 3D 'atomistic' MOSFET simulator using the density gradient formalism. We have studied in comparison with classical simulations the effect of the quantum mechanical corrections on the simulation of random dopant induced threshold voltage fluctuations, the effect of the single charge trapping on interface states and the effect of the oxide thickness fluctuations in decanano MOSFETs with ultrathin gate oxides. The introduction of quantum corrections enhances the threshold voltage fluctuations but does not affect significantly the amplitude of the random telegraph noise associated with single carrier trapping. The importance of the quantum corrections for proper simulation of oxide thickness fluctuation effects has also been demonstrated.

Asenov, Asen↗

NASA NTRS2000

Effect of Single-Electron Interface Trapping in Decanano MOSFETs: A 3D Atomistic Simulation Study

We study the effect of trapping/detrapping of a single-electron in interface states in the channel of n-type MOSFETs with decanano dimensions using 3D atomistic simulation techniques. In order to highlight the basic dependencies, the simulations are carried out initially assuming continuous doping charge, and discrete localized charge only for the trapped electron. The dependence of the random telegraph signal (RTS) amplitudes on the device dimensions and on the position of the trapped charge in the channel are studied in detail. Later, in full-scale, atomistic simulations assuming discrete charge for both randomly placed dopants and the trapped electron, we highlight the importance of current percolation and of traps with strategic position where the trapped electron blocks a dominant current path.

Asenov, Asen↗

NASA NTRS2000

Random Telegraph Signal Amplitudes in Sub 100 nm (Decanano) MOSFETs: A 3D 'Atomistic' Simulation Study

In this paper we use 3D simulations to study the amplitudes of random telegraph signals (RTS) associated with the trapping of a single carrier in interface states in the channel of sub 100 nm (decanano) MOSFETs. Both simulations using continuous doping charge and random discrete dopants in the active region of the MOSFETs are presented. We have studied the dependence of the RTS amplitudes on the position of the trapped charge in the channel and on the device design parameters. We have observed a significant increase in the maximum RTS amplitude when discrete random dopants are employed in the simulations.

Asenov, Asen↗