Search NASA⌕ Search

SEARCH · Search NASA

Results for “Form factors”

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

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

At least 433 records · Page 24

Trapping of gas mixtures by amorphous water ice

Our studies on gas trapping in amorphous water ice at 24-100 K were extended, by using mixtures of CH4, CO, N2, and Ar, rather than single gases. In 1:1 gas:(water vapor) mixtures, the competition among these gases on the available sites in the ice showed that the trapping capacity for the various gases is determined not only by the structure and dynamics of the ice, but is also influenced by the gas itself. Whereas at 24-35 K all four gases are trapped in the ice indiscriminantly, at 50-75 K there is a clear enhancement, in the order of CH4 > CO > N2 > or approximately Ar. This order is influenced by the gas-water interaction energy, the size of the trapped gas atom or molecule, the type of clathrate-hydrate formed (I or II) and, possibly, other factors. It seems that the gas can be trapped in the amorphous ice in several different locations, each being affected in a different way by the deposition temperature and gas composition. Once a gas atom or molecule is trapped in a specific location, it is predestined to emerge in one of eight different temperature ranges, which are associated with changes in the ice. The experimentally observed enhancements, together with the findings on the gas composition of comet Halley, might enable an estimation of the gas composition in the region of comet formation.

NASA Discipline Number 52-10↗

On-line technique for measuring stable oxygen and hydrogen isotopes from microliter quantities of water

Detailed here is a method for extracting and analyzing oxygen and hydrogen isotopes from 10 microL-sized water samples. Based on the traditional CO2-H2O equilibration technique, the oxygen isotope exchange reaction is done exclusively in sealed 6-mm (o.d.) Pyrex tubes at 25 degrees C, with full isotope exchange completed in at least 28 h. Using the same water sample employed in the 18O equilibration, D/H extractions are done in separate sealed 6-mm (o.d.) Pyrex tubes by reaction with Zn at 450 degrees C to form H2(g). Provided that a correction factor is applied to 18O analyses, accuracy and precision for both 18O and D/H are comparable to standard techniques using much larger samples.

NASA Center JSC↗

Evolution of large shield volcanoes on Venus

We studied the geologic history, topographic expression, and gravity signature of 29 large Venusian shield volcanoes with similar morphologies in Magellan synthetic aperture radar imagery. While they appear similar in imagery, 16 have a domical topographic expression and 13 have a central depression. Typical dimensions for the central depression are 150 km wide and 500 m deep. The central depressions are probably not calderas resulting from collapse of a shallow magma chamber but instead are the result of a corona-like sagging of a previously domical volcano. The depressions all have some later volcanic filling. All but one of the central depression volcanoes have been post-dated by geologic features unrelated to the volcano, while most of the domical volcanoes are at the top of the stratigraphic column. Analysis of the gravity signatures in the spatial and spectral domains shows a strong correlation between the absence of post-dating features and the presence of dynamic support by an underlying plume. We infer that the formation of the central depressions occurred as a result of cessation of dynamic support. However, there are some domical volcanoes whose geologic histories and gravity signatures also indicate that they are extinct, so sagging of the central region apparently does not always occur when dynamic support is removed. We suggest that the thickness of the elastic lithosphere may be a factor in determining whether a central depression forms when dynamic support is removed, but the gravity data are of insufficient resolution to test this hypothesis with admittance methods.

Herrick, Robert R.↗

The Atmospheric Dynamics of Alpha Tau (K5 III) -- Clues to Understanding the Magnetic Dynamo

Using HST/GHRS, HST/STIS and FUSE archival data for (alpha) Tau and the CHIANTI spectroscopic code, we have derived line shifts, volumetric emission measures, and plasma density estimates, and calculated filling factors for a number of UV lines forming between 10,000 K and 300,000 K in the outer atmosphere of this red giant star. The data suggest the presence of low-temperature extended regions and high-temperature compact regions, associated with magnetically open and closed structures in the stellar atmosphere, respectively. The signatures of UV lines from Alpha Tau can be consistently understood via a model of upward-traveling Alfven waves in a gravitationally stratified atmosphere. These wakes cause non-thermal broadening in UV lines due to unresolved wave motions and downward plasma motions in compact magnetic loops heated by resonant .4lf\en wave heating. We discuss implications of this interpretation for understanding the nature of magnetic dynamos operating in late-type giants.

Carpenter Kenneth G.↗

The Atmospheric Dynamics of alpha Tau (K5 III) - Clues to Understanding the Magnetic Dynamo in Late-Type Giant Stars

Using HST/GHRS, HST/STIS and FUSE archival data for alpha Tau and the CHIANTI spectroscopic code, we have derived line shifts, volumetric emission measures, and plasma density estimates, and calculated filling factors for a number of UV lines forming between 10,000 K and 300,000 K in the outer atmosphere of this red giant star. The data suggest the presence of low-temperature extended regions and high-temperature compact regions, associated with magnetically open and closed structures in the stellar atmosphere, respectively. The signatures of UV lines from alpha Tau can be consistently understood via a model of upward-traveling Alfven waves in a gravitationally stratified atmosphere. These waves cause non-thermal broadening in UV lines due to unresolved wave motions and downward plasma motions in compact magnetic loops heated by resonant Alfven wave heating.

Carpenter, Kenneth G.↗

The Atmospheric Dynamics of alpha Tau (K5 III) - Clues to Understanding the Magnetic Dynamo in Late-Type Giant Stars

Using HST/GHRS, HST/STIS and FUSE archival data for alpha Tau and the CHIANTI spectroscopic code, we have derived line shifts, volumetric emission measures, and plasma density estimates, and calculated filling factors for a number of UV lines forming between 10,000 K and 300,000 K in the outer atmosphere of this red giant star. The data suggest the presence of low-temperature extended regions and high-temperature compact regions, associated with magnetically open and closed structures in the stellar atmosphere, respectively. The signatures of UV lines from alpha Tau can be consistently understood via a model of upward-traveling Alfven waves in a gravitationally stratified atmosphere. These waves cause non-thermal broadening in UV lines due to unresolved wave motions and downward plasma motions in compact magnetic loops heated by resonant Alfven wave heating.

Carpenter, Kenneth G.↗

Simulated Altitude Performance of Combustors for the Westinghouse 24C Jet Engine I-24C-2 Combustor

A Westinghouse 24C-2 combustor was investigated at conditions simulating operation of the 24C Jet engine at zero ram over ranges of altitude and engine speed. The investigation was conducted to determine the altitude operational limits, that is, the maximum altitude for various engine speeds at which an average combustor-outlet gas temperature sufficient for operation of the jet engine could be obtained. Information was also obtained regarding the character of the flames, the combustion efficiency, the combustor-outlet gas temperature and velocity distributions, the extent of afterburning, the flow characteristics of the fuel manifolds, the combustor inlet-to-outlet total-pressure drop, and the durability of the combustor basket. The results of the investigation indicated that the altitude operational limits for zero ram decreased from 12,000 feet at an engine speed of 4000 rpm to a minimum of 9000 feet at 6000 rpm, and thence increased to 49,000 feet at 12,000 rpm.. At altitudes below the operational limits, flames were essentially steady, but, at altitudes above the operational limits, flames were often cycling and either blew out or caused violent explosions and vibrations. At conditions on the altitude operational limits the type of combustion varied from steady to cycling with increasing fuel-air ratio and the reverse occurred with decreasing fuel-air ratio. In the region of operation investigated, the combustion efficiency ranged from 75 to 95 percent at altitudes below the operational limits and dropped to 55 percent or less at some altitudes above the operational limits. The deviations in the local combustor-outlet gas temperatures were within +20 to -30 percent of the mean combustor temperature rise for inlet-air temperatures at the low end of the range investigated, but became more uneven (up to +/-100 percent) with increasing inlet-air temperatures. The distribution of the combustor-outlet gas velocity followed a similar trend. Practically no afterburning downstream of the combustor outlet occurred. At conditions of high inlet-air temperature several factors indicated that fuel vapor or air formed in the fuel manifolds and adversely affected combustion. The combustor inlet-to-outlet total-pressure drop can be correlated as a function of the ratio of the combustion-air inlet density to outlet density and of the inlet dynamic pressure. The walls of the combustor basket were warped and burned during 29 hours of operation.

Manganiello, Eugene J.↗

The Atmospheric Dynamics of Alpha Tau (K5 III) - Clues to Understanding the Magnetic Dynamo in Late-Type Giant Stars

Using HST/GHRS, HST/STIS and FUSE archival data for a Tau and the CHIANTI spectroscopic code, we have derived line shifts, volumetric emission measures, and plasma density estimates, and calculated filling factors for a number of UV lines forming between 10,000 K and 300,000 K in the outer atmosphere of this red giant star. The data suggest the presence of low-temperature extended regions and high-temperature compact regions, associated with magnetically open and closed structures in the stellar atmosphere, respectively. The signatures of UV lines from a Tau can be consistently understood via a model of upward-traveling Alfv6n waves in a gravitationally stratified atmosphere. These waves cause nonthermal broadening in UV lines due to unresolved wave motions and downward plasma motions in compact magnetic loops heated by resonant Alfven wave heating.

Carpenter, Kenneth G.↗

Life Beyond the Planet of Origin and Implications for the Search for Life on Mars

Outer space is vast, cold, devoid of matter, radiation filled with essentially no gravity. These factors present an environmental challenge for any form of life. Earth's biosphere has evolved for more than 3 billion years shielded from the hostile environment of outer space by the protective blanket of the atmosphere and magnetosphere. Space is a nutritional wasteland with no liquid water and readily available organic carbon. Moving beyond a life's planet of origin requires a means for transport, the ability to withstand transport, and the ability to colonize, thrive and ultimately evolve in the new environment. Can life survive beyond its home planet? The key to answering this question is to identify organisms that first have the ability to withstand space radiation, space vacuum desiccation and time in transit, and second the ability to grow in an alien environment. Within the last 60 years space technology allowed us to transport life beyond Earth's protective shield so we may study, in situ, their responses to selected conditions of space. To date a variety of microbes ranging from viruses, to Bacteria, to Archaea, to Eukarya have been tested in the space environment. Most died instantly, but not all. These studies revealed that ultraviolet radiation is the near-term lethal agent, while hard radiation is the long-term lethal agent when the organism is shielded from ultraviolet radiation. In fact, bacterial spores, halophilic cyanobacteria and Archaea as well as some lichens survive very well if protected from ultraviolet radiation [1]. Some microbes, then, may be able to survive the trip in outer space to Mars on a spacecraft or in a meteorite. Once on Mars can a terrestrial microbe survive? Although the conditions on Mars are not as harsh as those in space, they are not hospitable for a terrestrial microbe. Studies, however, have shown that certain microbes that can survive in space for several years may also be able to survive on Mars if protected from ultraviolet radiation [1]. Laboratory simulation experiments using a mock-up of the Phoenix lander have shown that microbes transported to the surface of Mars on a spacecraft come off the spacecraft and mix into the Martian regolith [2]. Additionally, studies simulating Martian dust storms demonstrate that microbes can survive in the Martian wind blown dust and be scattered across the Martian surface away from the spacecraft. Would these microbes that may survive on Mars metabolize and propagate? Growth requires liquid water, a carbon source and an energy source. Survival on Mars also requires protection from ultraviolet radiation. In the cold, dry environment of Mars the probability of microbial metabolism and growth at or just beneath the surface is extremely low. Although the probability is low, Mars may be contaminated with potentially live terrestrial organisms. In light of that statistic we must be extremely diligent and cautious in our search for Martian life. If we are not cautious we may find life on Mars and it may be a contaminant from Earth.

Survival in Space↗

Towards a General Equation for the Survival of Microbes Transferred between Solar System Bodies

It should be possible to construct a general equation describing the survival of microbes transferred between Solar System bodies. Such an equation will be useful for constraining the likelihood of transfer of viable organisms between bodies throughout the lifetime of the Solar System, and for refining Planetary Protection constraints placed on future missions. We will discuss the construction of such an equation, present a plan for definition of pertinent factors, and will describe what research will be necessary to quantify those factors. Description: We will examine the case of microbes transferred between Solar System bodies as residents in meteorite material ejected from one body (the "intial body") and deposited on another (the "target body"). Any microbes transferred in this fashion will experience four distinct phases between their initial state on the initial body, up to the point where they colonize the target body. Each of these phases features phenomena capable of reducing or exterminating the initial microbial population. They are: 1) Ejection: Material is ejected from the initial body, imparting shock followed by rapid desiccation and cooling. 2) Transport: Material travels through interplanetary space to the target body, exposing a hypothetical microbial population to extended desiccation, irradiation, and temperature extremes. 3) Infall: Material is deposited on the target body, diminishing the microbial population through shock, mass loss, and heating. 4) Adaptation: Any microbes which survive the previous three phases must then adapt to new chemophysical conditions of the target body. Differences in habitability between the initial and target bodies dominate this phase. A suitable general-form equation can be assembled from the above factors by defining the initial number of microbes in an ejected mass and applying multiplicitive factors based on the physical phenomena inherent to each phase. It should be possible to present the resulting equation in terms of initial ejection mass, ejection shock magnitude, transfer time, initial microbial load and/or other terms and generate graphs defining the number of surviving microbes. The general form of the equation is: x(sub f) = x(sub i) f(sub1) f(sub 2) f(sub 3) f(sub 4) Where x(sub f) is the final number of microbes to survive transfer, x(sub i) is the initial population prior to ejection, and f(sub 1-4) are mortality factors for the four phases described above. Among other considerations, f(sub 1) will vary with respect to impact shock magnitude and f(sub 2) will be time-dependent. Considerable research has been performed to date to quantify the survival rates of various microbes in response to portions of these four phases, both as vegetative cells and/or spores. Results indicate that many species tend to respond differently to the pertinent mortality factors, especially in the case of extremophiles. Therefore, a complete equation will include species-specific responses to the mortality factors.

Fries, M.↗

The optimal fiber volume fraction and fiber-matrix property compatibility in fiber reinforced composites

Although the question of minimum or critical fiber volume fraction beyond which a composite can then be strengthened due to addition of fibers has been dealt with by several investigators for both continuous and short fiber composites, a study of maximum or optimal fiber volume fraction at which the composite reaches its highest strength has not been reported yet. The present analysis has investigated this issue for short fiber case based on the well-known shear lag (the elastic stress transfer) theory as the first step. Using the relationships obtained, the minimum spacing between fibers is determined upon which the maximum fiber volume fraction can be calculated, depending on the fiber packing forms within the composites. The effects on the value of this maximum fiber volume fraction due to such factors as fiber and matrix properties, fiber aspect ratio and fiber packing forms are discussed. Furthermore, combined with the previous analysis on the minimum fiber volume fraction, this maximum fiber volume fraction can be used to examine the property compatibility of fiber and matrix in forming a composite. This is deemed to be useful for composite design. Finally some examples are provided to illustrate the results.

Pan, Ning↗

Full Life-Cycle Defect Management Assessment: Initial Inspection Data Collection Results and Research Questions for Further Study

It is often the case in software projects that when schedule and budget resources are limited, the Verification and Validation (V&V) activities suffer. Fewer V&V activities can be afforded and moreover, short-term challenges can result in V&V activities being scaled back or dropped altogether. As a result, too often the default solution is to save activities for improving software quality until too late in the life-cycle, relying on late-term code inspections followed by thorough testing activities to reduce defect counts to acceptable levels. As many project managers realize, however, this is a resource-intensive way of achieving the required quality for software. The Full Life-cycle Defect Management Assessment Initiative, funded by NASA s Office of Safety and Mission Assurance under the Software Assurance Research Program, aims to address these problems by: Improving the effectiveness of early life-cycle V&V activities to make their benefits more attractive to team leads. Specifically, we focus on software inspection, a proven method that can be applied to any software work product, long before executable code has been developed; Better communicating this effectiveness to software development teams, along with suggestions for parameters to improve in the future to increase effectiveness; Analyzing the impact of early life-cycle V&V on the effectiveness and cost required for late life-cycle V&V activities, such as testing, in order to make the tradeoffs more apparent. This white paper reports on an initial milestone in this work, the development of a preliminary model of inspection effectiveness across multiple NASA Centers. This model contributes toward reaching our project goals by: Allowing an examination of inspection parameters, across different types of projects and different work products, for an analysis of factors that impact defect detection effectiveness. Allowing a comparison of this NASA-specific model to existing recommendations in the literature regarding how to plan effective inspections. Forming a baseline model which can be extended to incorporate factors describing: the numbers and types of defects that are missed by inspections; how such defects flow downstream through software development phases; how effectively they can be caught by testing activities in the late stages of development. The model has been implemented in a prototype web-enabled decision-support tool which allows developers to enter their inspection data and receive feedback based on a comparison against the model. The tool also allows users to access reusable materials (such as checklists) from projects included in the baseline. Both the tool itself and the model underlying it will continue to be extended throughout the remainder of this initiative. As results of analyzing inspection effectiveness for defect containment are determined, they can be shared via the tool and also via updates to existing training courses on metrics and software inspections. Moreover, the tool will help satisfy key CMMI requirements for the NASA Centers, as it will enable NASA to take a global view across peer review results for various types of projects to identify systemic problems. This analysis can result in continuous improvements to the approach to verification.

Shull, Forrest↗

A method for measuring precipitation parameters and raindrop size distributions using radar reflectivity and optical extinction

A method of determining precipitation parameters from two remotely measurable quantities, the radar reflectivity factor and the optical extinction, is described. The raindrop size spectrum is approximated by a two-parameter exponential form; when these parameters are evaluated in terms of the radar reflectivity factor and the optical extinction, an exponential spectrum is obtained that is generally in very good agreement with the observed size spectrum. Other calculated precipitation parameters, such as rainfall rate and liquid water content, which are derived from the exponential approximation, also agree with experimental data. It is indicated that other combinations of two remote measurables can also be used to obtain more accurate estimates of precipitation parameters than can be obtained by the use of an empirical relationship.

Ulbrich, C. W.↗

An implicit factored scheme for the compressible Navier-Stokes equations

An implicit finite difference scheme is developed for the numerical solution of the compressible Navier-Stokes equations in conservation-law form. The algorithm is second-order-time accurate, noniterative, and spatially factored. In order to obtain an efficient factored algorithm, the spatial cross-derivatives are evaluated explicitly. However, the algorithm is unconditionally stable and, although a three-time-level scheme, requires only two-time-levels of data storage. The algorithm is constructed in a 'delta' form (i.e., increments of the conserved variables and fluxes) that provides a direct derivation of the scheme and leads to an efficient computational algorithm. In addition, the delta form has the advantageous property of a steady-state (if one exists) independent of the size of the time step. Numerical results are presented for a two-dimensional shock boundary-layer interaction problem.

Beam, R. M.↗

Ultrahigh Temperature Ceramics for Thermal Protection of Next Generation Space Vehicles

Materials with improved properties are needed for thermal protection of next generation space vehicles. Sharp leading edges on these vehicles will have to withstand exposure to high temperatures (> 2200 C or 4000 F) and severe thermal cycling in both neutral and oxidizing environments. These extreme conditions will require materials that possess superior oxidation resistance, low creep, and excellent thermal shock properties. This presentation will first discuss the system requirements for thermal protection of advanced space vehicles and then show how they are driving development of new materials systems. The presentation will focus on ultrahigh temperature ceramics (UHTCs) that are promising candidates for such applications. ZrB2 and HfB2 and composites of those ceramics with SiC are two particular families of UHTCs that are currently under development for sharp leading edges. These ceramics are appealing because their melting temperatures are 3245 C (5873 F) for ZrB2 and 3380 C (6116 F) for HfB2 and because they may form protective, oxidation resistant coatings in use. The mechanical properties of the UHTCs are strongly dependent on phase purity and the processing route used to make them, both of which factors are being actively investigated. For example, oxide impurities could form glassy grain boundary phases that soften at high temperatures and make the ceramic susceptible to creep deformation. Results from scanning and transmission electron microscopic studies of the UHTCs have shown how their processing can be improved to give better properties. This presentation will discuss the UHTC characterization results in some detail, focusing particularly on the structure and composition of the ceramic grain boundaries. The presentation will conclude with some remarks on how the properties of these promising UHTCs can be further improved and how they might be made more economically.

Loehman, R. E.↗

A Summary of Test and Analysis Results from a Second Lift+Cruise Full-Scale Drop Test

The realization of advanced air mobility markets is enabling new forms of transportation to take shape in the United States and around the world. Though currently in development, as these markets mature, new types of vertical take-off and landing (VTOL) vehicles have been undergoing development for use. There are many factors which must be addressed prior to these types of vehicles becoming viable alternative forms of transportation in these markets. These factors include incorporation into the existing airspaces, the logistics of operating in urban environments, along with numerous factors associated with safety and reliability. To address some of the safety aspects associated with the development of these new types of vehicles, NASA has been conducting research into the performance of an example electric VTOL (eVTOL) aircraft as a part of the Revolutionary Vertical Lift Technology (RVLT) project. Over the course of this research, many aspects including the development of energy absorbing components, the evaluation of seating systems, the development of advanced finite element material model systems and the acquisition of full-scale vehicle impact data were investigated. The report will discuss aspects related to the acquisition of full-scale vehicle data which occurred in the form of a full-scale impact test conducted in the Summer of 2025. This test was on a NASA designed Lift+Cruise composite cabin test article and represented a partial capstone in the entirety of previous eVTOL research conducted for the project. In this test, a variety of experiments were included in order to investigate the effect of a full-scale environment on the experiment results. In parallel, the development of a computational impact model to simulate the full-scale test will be discussed in this report. A model of the Lift+Cruise test article was developed utilizing data collected from previous sub- and full-scale test data and then simulated in the current test environment. The model development, its use in pre-test predictions, and its use in post-test correlation will all be presented. This report will present the test data acquired from the Lift+Cruise test and document several of the results obtained. One intended result is to determine the effect of a complex full-scale crash impact on the identification of occupant injury risk within seat and vehicle designs. A second intended result is to determine whether high-fidelity models can be used with some confidence in the prediction of test events and can allow for additional test cases to be simulated without the need of having to conduct additional tests. The overall goal of the test is to provide the community with data that can be used for design, development or certification efforts, along with providing data on what an example eVTOL crash incident could entail.

energy storage systems↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗