SEARCH · Engineering Papers

Results for “diffusion length”

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 757 records · Page 42

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Structure of Soot-Containing Laminar Jet Diffusion Flames

The structure and soot properties of nonbuoyant and weakly-buoyant round jet diffusion flames were studied, considering ethylene, propane and acetylene burning in air at pressures of 0.125-2.0 atm. Measurements of flame structure included radiative heat loss fractions, flame shape and temperature distributions in the fuel-lean (overfire) region. These measurements were used to evaluate flame structure predictions based on the conserved-scalar formalism in conjunction with the laminar flamelet concept, finding good agreement betweem predictions and measurements. Soot property measurements included laminar smoke points, soot volume function distributions using laser extinction, and soot structure using thermophoretic sampling and analysis by transmission electron microscopy. Nonbuoyant flames were found to exhibit laminar smoke points like buoyant flames but their properties are very different; in particular, nonbuoyant flames have laminar smoke point flame lengths and residence times that are shorter and longer, respectively, than buoyant flames.

Mortazavi, S.↗

NASA NTRS2003

Buoyancy Effects in Fully-Modulated, Turbulent Diffusion Flames

Pulsed combustion appears to have the potential to provide for rapid fuel/air mixing, compact and economical combustors, and reduced exhaust emissions. The objective of this experiment (PuFF, for Pulsed-Fully Flames) is to increase the fundamental understanding of the fuel/air mixing and combustion behavior of pulsed, turbulent diffusion flames by conducting experiments in microgravity. In this research the fuel jet is fully-modulated (i.e., completely shut off between pulses) by an externally controlled valve system. This gives rise to drastic modification of the combustion and flow characteristics of flames, leading to enhanced fuel/air mixing compared to acoustically excited or partially-modulated jets. Normal-gravity experiments suggest that the fully-modulated technique also has the potential for producing turbulent jet flames significantly more compact than steady flames with no increase in exhaust emissions. The technique also simplifies the combustion process by avoiding the acoustic forcing generally present in pulsed combustors. Fundamental issues addressed in this experiment include the impact of buoyancy on the structure and flame length, temperatures, radiation, and emissions of fully-modulated flames.

Hermanson, J. C.↗

NASA NTRS2002

Surface Processes of Faceted Growth

Atomic force microscopy and high precision optical interferometry allow to analyze the processes that are in the core of our present understanding of faceted crystal growth. Some of these processes will be reviewed. Recent experiments suggest that the Gibbs-Thomson Law (GTL) may not be applicable to the weakly fluctuating strongly polygonized steps, with low kink density. Propagation rate of a straight short step segment changes with the segment length much steeper than predicted by GTL. If the step fluctuations are still well developed, the typical transition length may be determined just by the radius of rounded comers of the polygonized step. If fluctuations do not have enough time to develop, the kinetic effects may be essential. Indeed, 'communications' between comers of a short straight step segment may occur and, thus, the GTL may be implemented only via kink exchange. E.g., dissolutions of this short segment at macroscopic equilibrium occurs as follows: The kinks are split from one comer as a result of fluctuations, and annihilate with kinks of the opposite sign generated by another comer. Similarly, an exchange via chain of kinks on a rough step is the mechanism to implement the GT Law. If there is a supersaturation with respect to a large crystal, the kinks not only diffuse along the step, but move back the generating comers. This 'wind' prevents communication between the comers and places upper limit supersaturation only below which GTL is applicable. Steps replace kinks in the 3D case, though, of course, another dimensionality brings about another physics. Another specific problem to think about is if fluctuations are always fast enough to provide enough kinks for linear dependence of step rate on supersaturation. Increasing azimuthal polygonization of growth hillocks is a sign of insufficient fluctuation rate. Morphological stability of vicinal faces is usually considered on the basis of mutually parallel steps. Azimuthal anisotropy and interlacing of step bunches inducing the well known bunch splitting is another challenge. Among other challenging problems are: existence of surface diffusion on the crystal/so1ution interface, step interaction in solutions, achievement of step flow modes from liquids, impurity effects. New phenomena essential in biomacromolecular crystallization is another interesting area.

Chernov, A.↗

NASA NTRS2001

Characterization of a Laminate Flat Plate Diffusion Flame in Microgravity using PIV, Visible and CH Emissions

Motivated by fire safety concerns and the advent of long-term micro-gravity facilities, a cooperative program has been developed to study the mechanisms and material properties that control flow assisted (co-current) flame spread. This program has used as a common fire scenario a reacting steady-state boundary layer. Preliminary studies explored the aerodynamics of a reacting boundary layer by simulating a condensed fuel by means of a gas burner. Stability curves for ethane air flames were obtained and different burning regimes were identified. An important feature of this study was the independent identification of the different mechanisms leading to the instability of the flow. It was observed that fuel injection velocity and thermal expansion independently contributed to the separation of the flow at the leading edge of the burner. The occurrence of separation resulted in complex three-dimensional flow patterns that have a dominant effect on critical fire safety parameters such as the stand-off distance and flame length. This work was extended to a solid fuel (PMMA) leading to a Sounding Rocket experiment (Mini-Texus-6). The solid phase showed similar flow patterns, mostly present at low flow velocities (<100 mm/s) but the results clearly demonstrated that the thermal balance at the pyrolyzing fuel surface is the dominant mechanism that controls both stand-off distance and flame length. This thermal balance could be described in a global manner by means of a total mass transfer or "B" number. This "B" number incorporates surface re-radiation, radiative feedback and in-depth heat conduction as first prescribed by Emmons. The mass transfer number becomes the single parameter that determines the evolution of these fire safety variables (flame length, stand-off distance) and therefore can be used as a ranking criterion to assess the flammability of materials. The particular configuration is representative of the NASA upward flame spread test (Test 1) therefore this approach can be used in the interpretation of the results obtained from this test. Nevertheless, complete validation of this approach has not been fully achieved due, mainly because all the measurements necessary to compare with the theoretical predictions have not been obtained. Following these studies two different directions have been taken. The first attempts to elucidate the details of the gas phase combustion reaction and the associated flow field by means of quantitative and qualitative measurements. The second approach, a more practical one, is to apply this methodology to the assessment of material flammability. The former is currently being conducted with a gas burner because it allows for easier control and longer experimentation time. The results obtained so far will be presented in more detail. The latter is a new program therefore only a brief summary of the objectives will be presented.

Joulain, P.↗

NASA NTRS1976

One-dimensional numerical study of charged particle trajectories in turbulent electrostatic wave fields

The paper describes a numerical simulation of electron trajectories in weak random electric fields under conditions that are approximately true for Langmuir waves whose wavelength is much longer than the Debye length. Two types of trajectory calculations were made: (1) the initial particle velocity was made equal to the mean phase velocity of the waves, or (2) it was equal to 0.7419 times the mean velocity of the waves, so that the initial velocity differed substantially from all phase velocities of the wave spectrum. When the autocorrelation time is much greater than the trapping time, the particle motion can change virtually instantaneously from one of three states - high-velocity, low-velocity, or trapped state - to another. The probability of instantaneous transition from a high- or low-velocity state becomes small when the difference between the particle velocity and the mean phase velocity of the waves becomes high in comparison to the trapping velocity. Diffusive motion becomes negligible under these conditions also.

Graham, K. N.↗

NASA NTRS1993

Absence of magnetic trapping on closed field lines at Neptune

The large offset of Neptune's magnetic field in an offset tilted dipole (OTD) approximation allows some fraction of magnetospheric field lines to have monotonically increasing field magnitude over their entire length in the magnetosphere. The usual magnetic trapping is impossible on such field lines, meaning that all charged particles on these field lines will precipitate after mirroring at most once. These field lines can be thought of as the part of the global magnetic anomaly where no particles on the corresponding drift shells remain trapped for an entire rotation of the planet. In this paper we examine both the morphology and size of the affected surface area in the OTD2 and O8 magnetic field models. For regions where there are indications that part of the aurora could be due to monotonic field line precipitation we analyze whether these field lines could be sufficiently populated to produce auroral signatures by considering photoelectron production rates and radial diffusion. We conclude that monotonic field line precipitation may contribute to the observed Neptune aurora but does not explain the global phenomenon.

Paranicas, C.↗

NASA NTRS2017

Vortex Generators in a Streamline-Traced, External-Compression Supersonic Inlet

Vortex generators within a streamline-traced, external-compression supersonic inlet for Mach 1.66 were investigated to determine their ability to increase total pressure recovery and reduce total pressure distortion. The vortex generators studied were rectangular vanes arranged in counter-rotating and co-rotating arrays. The vane geometric factors of interest included height, length, spacing, angle-of-incidence, and positions upstream and downstream of the inlet terminal shock. The flow through the inlet was simulated numerically through the solution of the steady-state, Reynolds-averaged Navier-Stokes equations on multi-block, structured grids using the Wind-US flow solver. The vanes were simulated using a vortex generator model. The inlet performance was characterized by the inlet total pressure recovery and the radial and circumferential total pressure distortion indices at the engine face. Design of experiments and statistical analysis methods were applied to quantify the effect of the geometric factors of the vanes and search for optimal vane arrays. Co-rotating vane arrays with negative angles-of-incidence positioned on the supersonic diffuser were effective in sweeping low-momentum flow from the top toward the sides of the subsonic diffuser. This distributed the low-momentum flow more evenly about the circumference of the subsonic diffuser and reduced distortion. Co-rotating vane arrays with negative angles-of-incidence or counter-rotating vane arrays positioned downstream of the terminal shock were effective in mixing higher-momentum flow with lower-momentum flow to increase recovery and decrease distortion. A strategy of combining a co-rotating vane array on the supersonic diffuser with a counter-rotating vane array on the subsonic diffuser was effective in increasing recovery and reducing distortion.

computational fluid dynamics↗

NASA NTRS2021

An Encounter with the Ion and Electron Diffusion Regions at a Flapping and Twisted Tail Current Sheet

We analyze data returned by the Magnetospheric Multiscale mission (MMS) constellation during a rapid (~1.5 s) traversal of a flapping and reconnecting current sheet (CS) in the near-Earth magnetotail (X ~ -20 RE). The CS was highly tilted, with its normal pointing strongly duskward. Its extreme thinness was confirmed by a curvature analysis of the magnetic field lines. The event was associated with a guide field of 8% of the reconnecting components. From the pitch angle distributions of low-energy electrons we infer that the crossing occurred earthward of the X-line. Traveling practically normal to the CS, MMS encountered an ion diffusion region (IDR) in which was embedded an electron diffusion region (EDR). IDR signatures included breaking of the ion frozen-in condition in the presence of Hall B and E fields. EDR signatures included a strong out of-plane current associated with a super-Alfvenic electron jet, evidence of positive energy transfer, and a temperature anisotropy (TeII > Te1) which disappeared at the field reversal. Derived scale sizes normal to the CS are: ~6.9 de (EDR) and ~0.4 di (IDR) (40 and 100 km). We estimate the average dimensionless reconnection rate as 0.077 + 0.050. The observations and inferences are supported by PIC numerical simulations. We find very good agreement in the reconnection rates. We also discuss the effects of asymmetries in the density, temperature and magnetic field strength on the Hall fields and length of the outflow jets. The event is associated with a substorm onset which began 7 min after the MMS observations.

C J Farrugia↗

NASA NTRS1986

Computer simulation of a wind tunnel test section with discrete finite-length wall slots

A computer simulation of a slotted wind tunnel test section which includes a discrete, finite-length wall slot representation with plenum chamber constraints and accounts for the nonlinear effects of the dynamic pressure of the slot outflow jet and of the low energy of slot inflow air was developed. The simulation features were selected to be those appropriate for the intended subsequent use of the simulation in a wall interference assessment procedure using sparsely located wall pressure measurements. It is demonstrated that accounting for slot discreteness is important in interpreting wall pressure measured between slots, and that accounting for nonlinear slot flow effects produces significant changes in tunnel-induced velocity distributions and, in particular, produces a longitudinal component of tunnel-induced velocity due to model lift. A characteristic mode of tunnel flow interaction with constraints imposed by the plenum chamber and diffuser entrance is apparent in simulation results and is derived analytically through a simplified analysis.

Kemp, W. B., Jr.↗

NASA NTRS1995

Convergence properties of finite-difference hydrodynamics schemes in the presence of shocks

We investigate the asymptotic convergence of finite-difference schemes for the Euler equations when the limiting solution contains shocks. The Lax-Wendroff theorem guarantees that certain conservative schemes converge to correct, physically valid solutions. We focus on two one-dimensional operator-split schemes with explicit artificial-viscosity terms. One, an internal-energy scheme, does not satisfy the assumptions of Lax-Wendroff; the other, a conservative total-energy scheme, does. With viscous lengths chosen proportional to the grid size, we find that both schemes converge to their zero-grid-size limits at the theoretically expected rate, but only the conversative scheme converges toward correct solutions of the inviscid fluid equations. We show that the difference in their behaviors results directly from the presence of shocks in the limiting solution. Empirically, we find that when the viscous lenghts tend toward zero more slowly than the grid size, however the nonconservative scheme also converges toward correct solutions. We characterize the asymptotic behavior of the total-energy scheme in a particular problem in which a shock forms. As the grid is refined, a Cauchy error approaches the expected rate of change slowly. We show that the changes in the artificial viscosity alter the diffusion of small-amplitude waves. The differences associated with such waves make the dominant contribution to the Cauchy error. We formulate an analytic model to relate the rate of approach to the effect of varying diffusion in waves and find quantitative agreement with our numerical results.

Kimoto, Paul A.↗

NASA NTRS2006

Carbon-Nanotube-Carpet Heat-Transfer Pads

Microscopic thermal-contact pads that include carpet-like arrays of carbon nanotubes have been invented for dissipating heat generated in integrated circuits and similarly sized single electronic components. The need for these or other innovative thermal-contact pads arises because the requisite high thermal conductances cannot be realized by scaling conventional macroscopic thermal-contact pads down to microscopic sizes. Overcoming limitations of conventional thermal-contact materials and components, the carbon-nanotube thermal-contact pads offer the high thermal conductivities needed to accommodate the high local thermal power densities of modern electronic circuits, without need for large clamping pressures, extreme smoothness of surfaces in contact, or gap-filling materials (e.g., thermally conductive greases) to ensure adequate thermal contact. Moreover, unlike some conventional thermal-contact components, these pads are reusable. The figure depicts a typical pad according to the invention, in contact with a rough surface on an electronic component that is to be cooled. Through reversible bending and buckling of carbon nanotubes at asperities on the rough surface, the pad yields sufficiently, under relatively low contact pressure, that thermal contact is distributed to many locations on the surface to be cooled, including valleys where contact would not ordinarily occur in conventional clamping of rigid surfaces. Hence, the effective thermal-contact area is greater than that achievable through scaling down of a macroscopic thermal-contact pad. The extremely high longitudinal thermal conductivities of the carbon nanotubes are utilized to conduct heat away from potential hot spots on the surface to be cooled. The fibers protrude from a layer of a filler material (Cu, Ag, Au, or metal-particle- filled gels), which provides both mechanical support to maintain the carbon nanotubes in alignment and thermal conductivity to enhance the diffusion of concentrated heat from the nanotubes into the larger adjacent volume of a heat sink. The array of carbon nanotubes, the filler material, and the heat sink are parts of a unitary composite structure that is fabricated as follows: 1. Using techniques that have been reported previously, the array of substantially perpendicularly oriented carbon nanotubes is grown on a metal, silicon, or other suitable thermally conductive substrate that is intended to become the heat sink. 2. By means of chemical vapor deposition, physical vapor deposition, plasma deposition, ion sputtering, electrochemical deposition, or casting from a liquid phase, some or all of the interstitial volume between carbon nanotubes is filled with the aforementioned layer of mechanically supporting, thermally conductive material. 3. To cause the carbon nanotubes to protrude the desired length from the filler material, an outer layer of filler is removed by mechanical polishing, chemical mechanical polishing, wet chemical etching, electrochemical etching, or dry plasma etching.

Li, Jun↗

NASA NTRS1985

Generation of Z mode radiation by diffuse auroral electron precipitation

The generation of Z mode waves by diffuse auroral electron precipitation is investigated assuming that a loss cone exists in the upgoing portion of the distribution due to electron interactions with the atmosphere. The waves are generated at frequencies above, but very near, the local electron cyclotron frequency omega(e) and at wave normal angles larger than 90 deg. In agreement with Hewitt et al. (1983), the group velocity is directed downward in regions where the ratio of the upper hybrid frequency omega(pe) to Omega(e) is less than 0.5, so that Z mode waves excited above a satellite propagate toward it and away from the upper hybrid resonance. Z mode waves are excited in a frequency band between Omega(e) and about 1.02 Omega(e), and with maximum growth rates of about 0.001 Omega(e). The amplification length is about 100 km, which allows Z mode waves to grow to the intensities observed by high-altitude satellites.

Dusenbery, P. B.↗

NASA NTRS1998

The Effects of Ground and Space Processing on the Properties of Organic, Polymeric, and Colloidal Materials

In recent years, a great deal of interest has been directed toward the use of organic materials in the development of high-efficiency optoelectronic and phototonic devices. There is a myriad of possibilities among organic materials which allow flexibility in the design of unique structures with a variety of functional groups. The use of nonlinear optical (NLO) organic materials as thin film wave-guides allows full exploitation of their desirable qualifies by permitting long interaction lengths and large susceptibilities allowing modest power input. There are several methods in use to prepare thin films such as Langmuir-Blodgett (LB) and self-assembly techniques, vapor deposition, growth from sheared solution or melt, and melt growth between glass plates. Organic-based materials have many features that make them desirable for use in optical devices, such as high second-and third-order nonlinearity, flexibility of molecular design, and damage resistance to optical radiation. However, their use in devices has been hindered by processing difficulties for crystals and thin films. We discuss the potential role of microgravity processing of a few organic and polymeric materials. It is of interest to note how materials with second-and third-order NLO behavior may be improved in a diffusion-limited environment and ways in which convection may be detrimental to these materials. We focus our discussion on third-order materials for all-optical switching, and second-order materials for frequency conversion and electrooptics. The goal of minimizing optical loss obviously depends on processing methods. For solution-based processes, such as solution crystal growth and solution photopolymerization, it is well known that thermal and solutal density gradients can initiate buoyancy-driven convection. Resultant fluid flows can affect transport of material to and from growth interfaces and become manifest in the morphology and homogeneity of the growing film or crystal. Likewise, buoyancy-driven convection can hinder production of defect-free, high-quality crystals or films during crystal and film growth by vapor deposition.

Frazier, Donald O.↗

NASA NTRS1994

Changes in intervertebral disc cross-sectional area with bed rest and space flight

STUDY DESIGN. We measured the cross-sectional area of the intervertebral discs of normal volunteers after an overnight rest; before, during, and after 5 or 17 weeks of bed rest; and before and after 8 days of weightlessness. OBJECTIVES. This study sought to determine the degree of expansion of the lumbar discs resulting from bed rest and space flight. SUMMARY OF BACKGROUND DATA. Weightlessness and bed rest, an analog for weightlessness, reduce the mechanical loading on the musculoskeletal system. When unloaded, intervertebral discs will expand, increasing the nutritional diffusion distance and altering the mechanical properties of the spine. METHODS. Magnetic resonance imaging was used to measure the cross-sectional area and transverse relaxation time (T2) of the intervertebral discs. RESULTS. Overnight or longer bed rest causes expansion of the disc area, which reaches an equilibrium value of about 22% (range 10-40%) above baseline within 4 days. Increases in disc area were associated with modest increases in disc T2. During bed rest, disc height increased approximately 1 mm, about one-half of previous estimates based on body height measurements. After 5 weeks of bed rest, disc area returned to baseline within a few days of ambulation, whereas after 17 weeks, disc area remained above baseline 6 weeks after reambulation. After 8 days of weightlessness, T2, disc area, and lumbar length were not significantly different from baseline values 24 hours after landing. CONCLUSIONS. Significant adaptive changes in the intervertebral discs can be expected during weightlessness. These changes, which are rapidly reversible after short-duration flights, may be an important factor during and after long-duration missions.

Non-NASA Center↗

NASA NTRS1999

A Validated All-Pressure Fluid Drop Model and Lewis Number Effects for a Binary Mixture

The differences between subcritical liquid drop and supercritical fluid drop behavior are discussed. Under subcritical, evaporative high emission rate conditions, a film layer is present in the inner part of the drop surface which contributes to the unique determination of the boundary conditions; it is this film layer which contributes to the solution's convective-diffusive character. In contrast, under supercritical condition as the boundary conditions contain a degree of arbitrariness due to the absence of a surface, and the solution has then a purely diffusive character. Results from simulations of a free fluid drop under no-gravity conditions are compared to microgravity experimental data from suspended, large drop experiments at high, low and intermediary temperatures and in a range of pressures encompassing the sub-and supercritical regime. Despite the difference between the conditions of the simulations and experiments (suspension vs. free floating), the time rate of variation of the drop diameter square is remarkably well predicted in the linear curve regime. The drop diameter is determined in the simulations from the location of the maximum density gradient, and agrees well with the data. It is also shown that the classical calculation of the Lewis number gives qualitatively erroneous results at supercritical conditions, but that an effective Lewis number previously defined gives qualitatively correct estimates of the length scales for heat and mass transfer at all pressures.

Harstad, K.↗

NASA NTRS2021

Microstructural Observations on Advanced Multilayer Protective Coatings for Nuclear Thermal Propulsion

The present research discusses the microstructures of a new multilayered coating concept proposed to solve the “midrange corrosion” problem observed in NbC and ZrC-coated graphite (Gr) nuclear fuels during the NERVA/Rover programs. The concept envisions designing a compliant multilayered coatings architecture to accommodate the thermal expansion mismatch strains between the ZrC outer coating and the Gr/(U,Zr)C fuel matrix along with using a diffusion barrier to minimize carbon and hydrogen diffusion, respectively. Multilayered coatings were deposited on Gr substrates by chemical vapor deposition (CVD). Graphite disks, as well as the inner channels of a 19-hole hexagonal rod 152.4 mm long, were coated by CVD in proof-of concept process demonstration studies. Coated disk specimens were thermally-cycled between ambient and 1900 K for as long as 20 cycles with a 1h hold time at the temperature using both stepped and non-interrupted thermal cycling methods. Detailed cross-sectional microstructural analyses were conducted on the thermally-cycled coated disks. These results showed that the coatings were intact after 20 thermal cycles at 1900 K although the bonding of the overlay Mo coating appeared to be influenced by the oxygen content in the layer. Seven transverse sections were cut along the length from the coated hexagonal rod for microstructural studies. All the nineteen channels were observed to have been coated in all the sections. This paper reports the results of detailed microstructural analyses, compositional line scans and x-ray dot maps conducted on three coated channels from a mid-section of the rod. The measured values of carbon (C) were observed to be exceptionally high in all the specimens, which were attributed to hydrocarbon contamination of the specimen surfaces. It was observed that the oxygen (O) levels were high in all the coating layers even after approximately correcting for the observed high values of C. The source of O contamination of the CVD reactor system was attributed to either the presence of oxychloride impurities in the precursor salts or a pin hole leak in the coating system. It is concluded that high values of O result in poor bonding of the Mo overlay to the ZrC layer. These trials have established that it is possible to deposit multilayered coatings by CVD in 19-hole hexagonal rods 152.4 mm long similar to the NERVA fuel rod design although the quality of the coatings and the uniformity of their thicknesses require further improvement.

Nuclear thermal propulsion, multilayer coatings, m↗

NASA NTRS1997

Microgravity Processing and Photonic Applications of Organic and Polymeric Materials

In recent years, a great deal of interest has been directed toward the use of organic materials in the development of high-efficiency optoelectronic and photonic devices. There is a myriad of possibilities among organics which allow flexibility in the design of unique structures with a variety of functional groups. The use of nonlinear optical (NLO) organic materials such as thin-film waveguides allows full exploitation of their desirable qualities by permitting long interaction lengths and large susceptibilities allowing modest power input. There are several methods in use to prepare thin films, such as Langmuir-Blodgett (LB) and self-assembly techniques, vapor deposition, growth from sheared solution or melt, and melt growth between glass plates. Organics have many features that make them desirable for use in optical devices such as high second- and third-order nonlinearities, flexibility of molecular design, and damage resistance to optical radiation. However, their use in devices has been hindered by processing difficulties for crystals and thin films. In this chapter, we discuss photonic and optoelectronic applications of a few organic materials and the potential role of microgravity on processing these materials. It is of interest to note how materials with second- and third-order nonlinear optical behavior may be improved in a diffusion-limited environment and ways in which convection may be detrimental to these materials.

Frazier, Donald O.↗

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↗