Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multi-phase Flow”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Experimental and Computational Investigation of Flow Boiling in Microgravity

This study explores use of Computational Fluid Dynamics (CFD) to predict near-saturated flow boiling of FC-72 in microgravity. The computational method employs transient analysis to predict detailed interfacial behavior and heat transfer characteristics along a rectangular channel heated along two opposite walls. Predicted results are validated against experimental temperature measurements and high-speed video images captured during a series of parabolic aircraft maneuvers for three sets of operating conditions which include variations of both mass velocity and wall heat flux. The computational method is based on the multi-phase volume of fluid (VOF) model, which is combined with appropriate phase change and turbulence models, and accounts for both shear-lift force on bubbles and conjugate heat transfer along the heating walls. A key advantage of the CFD method is ability to capture details that are very difficult to measure experimentally, including detailed spatial variations of bubble shape, void fraction, mixture fluid temperature, liquid velocity, and vapor velocity, results for which are presented for each of the three test cases. Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with experiment. Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations. Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method. Predicted wall temperature is fairly uniform in the middle of the heated length but increases in the entrance region, due to sensible heat transfer in the subcooled liquid, and decreases toward the exit, mostly because of flow acceleration resulting from increased void fraction.

two-phase cooling; flow boiling; CFD; microgravity↗
DOE OSTI2025

Machine learning modeling and model predictive control of a closed-circuit reverse osmosis system

Closed-circuit reverse osmosis (CCRO) offers a flexible and energy-efficient alternative to conventional reverse osmosis by operating in a semi-batch mode that recycles brine, enabling higher recovery rates and reduced specific energy consumption (SEC). However, developing accurate, system-level dynamic models for CCRO remains challenging due to its nonlinear, multi-phase operation and sensitivity to variable feed water conditions. Traditional modeling approaches, such as NARMAX (nonlinear autoregressive moving average with exogenous inputs), often struggle to generalize across varying inlet feed concentrations, necessitating frequent parameter re-estimation and limiting their utility for real-time control applications. To address these limitations, we developed a long short-term memory (LSTM) neural network model trained on an extensive experimental data set from a CCRO pilot plant. The model accepts three inputs, feed flow rate, recirculation flow rate, and initial feed conductivity, and predicts three key outputs: reject conductivity, feed pump power draw, and recirculation pump power draw. We validated the LSTM model against experimental data, demonstrating its ability to distinguish between different feed conductivities and adapt to variable flow rates. Subsequently, we incorporated the LSTM model within a nonlinear model predictive control (MPC) scheme and conducted closed-loop simulations to optimize the integrated SEC (iSEC). In conclusion, the results project up to a 6% reduction in iSEC by using MPC to optimize performance over the entire experiment duration, without requiring any random excitation for data collection or parameter re-estimation.

Desalination↗
NASA NTRSDate not supplied

Overview of Plume-Surface Interaction Data from Subscale Inert Gas Testing at NASA MSFC Test Stand 300 Vacuum Facilities

Maturation and validation of predictive modeling capabilities for plume-surface interaction, and the effects caused by the impingement of rocket engine exhaust on surface regolith, require relevant, highly characterized ground test data. The multi-phase nature of this problem poses challenges to ground testing across physical scaling, test environments and techniques, and diagnostics. NASA has developed a Physics Focused Ground Test (PFGT) of plume-surface interaction to address this need. The PFGT is a subscale, inert-gas experiment at Mars- and Lunar-relevant ambient pressure conditions that will provide critical data for model validation of cratering, ejecta, and plume physics. In this test, the ambient pressure in a vacuum chamber will be varied, along with the height of a supersonic nozzle above a surface, and the mass flow rate through the nozzle. Testing will be performed on a range of regolith simulants of varying composition and will also include instrumented impingement plate tests in place of regolith simulant to characterize gas-only plume physics. Novel high-speed visual diagnostics will capture temporal data on crater formation and ejecta dynamics, and planar laser-induced fluorescence diagnostics will capture plume structure and behavior. This paper will present an overview of the data collected during this test – for both the cratering and ejecta testing and the gas-only plume physics testing. Emphasis will be placed on trends, effects, and interactions observed in the data across the test parameter space.

Chad J. Eberhart↗
NASA NTRS2021

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗
NASA NTRS2021

A Study into Validating A Coupled Method of Characteristics And Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

DAC↗
NASA NTRS2008

Aerodynamic Design of a Dual-Flow Mach 7 Hypersonic Inlet System for a Turbine-Based Combined-Cycle Hypersonic Propulsion System

A new hypersonic inlet for a turbine-based combined-cycle (TBCC) engine has been designed. This split-flow inlet is designed to provide flow to an over-under propulsion system with turbofan and dual-mode scramjet engines for flight from takeoff to Mach 7. It utilizes a variable-geometry ramp, high-speed cowl lip rotation, and a rotating low-speed cowl that serves as a splitter to divide the flow between the low-speed turbofan and the high-speed scramjet and to isolate the turbofan at high Mach numbers. The low-speed inlet was designed for Mach 4, the maximum mode transition Mach number. Integration of the Mach 4 inlet into the Mach 7 inlet imposed significant constraints on the low-speed inlet design, including a large amount of internal compression. The inlet design was used to develop mechanical designs for two inlet mode transition test models: small-scale (IMX) and large-scale (LIMX) research models. The large-scale model is designed to facilitate multi-phase testing including inlet mode transition and inlet performance assessment, controls development, and integrated systems testing with turbofan and scramjet engines.

Sanders, Bobby W.↗
NASA NTRS1999

Combustion of Metals in Reduced-Gravity and Extra Terrestrial Environments

The combustion of metals is a field with important practical applications in rocket propellants, high-temperature flames, and material synthesis. Also, the safe operation of metal containers in high-pressure oxygen systems and with cryogenic fuels and oxidizers remains an important concern in industry. The increasing use of metallic components in spacecraft and space structures has also raised concerns about their flammability properties and fire suppression mechanisms. In addition, recent efforts to embark on unmanned and manned planetary exploration, such as on Mars, have also renewed the interest in metal/carbon-dioxide combustion as an effective in situ resource utilization technology. In spite of these practical applications, the understanding of the combustion properties of metals remains far behind that of the most commonly used fuels such as hydrocarbons. The lack of understanding is due to the many problems unique to metal- oxidizer reactions such as: low-temperature surface oxidation prior to ignition, heterogeneous reactions, very high combustion temperatures, product condensation, high emissivity of products, and multi-phase interactions. Very few analytical models (all neglecting the influence of gravity) have been developed to predict the burning characteristics and the flame structure details. Several experimental studies attempting to validate these models have used small metal particles to recreate gravity-free conditions. The high emissivity of the flames, rapid reaction, and intermittent explosions experienced by these particles have made the gathering of any useful information on burning rates and flame structure very difficult. The use of a reduced gravity environment is needed to clarify some of the complex interactions among the phenomena described above. First, the elimination of the intrusive buoyant flows that plague all combustion phenomena is of paramount importance in metal reactions due to the much higher temperatures reached during combustion. Second, a low-gravity environment is absolutely essential to remove the destructive effect of gravity on the shape of a molten metal droplet in order to study a spherically symmetric condition with large bulk samples. The larger size of the spherical metal droplet and the longer burning times available in reduced gravity extend the spatial and temporal dimensions to permit careful probing of the flame structure and dynamics. Third, the influence of the radiative heat transfer from the solid oxides can be studied more carefully by generating a stagnant spherical shell of condensed products undisturbed by buoyancy.

Branch, M.C.↗
DOE OSTI2026

Development and Application of High-Fidelity Models for Heterogeneous CO2 Frost Formation

Carbon America has developed a cryogenic carbon capture technology ("FrostCC") that separates CO2 from point source emissions by solidifying it at cold temperatures through preferential desublimation. Cooling is achieved through a series of interlinked compression, heat exchange, and expansion operations. In the current system, frosting of CO2 happens in heat exchangers, followed by CO2 recovery in a separate extraction step. In this work, multiphysics computational fluid dynamics (CFD) models are developed and validated for compressible and low Mach flows to simulate the formation of solid CO2 in flue gas flowing in a heat exchanger geometry. The models track the mass transfer rate of CO2 from gas phase to solid phase, heat released from desublimation, and the evolution of the solid CO2 layer. Simulations are used to answer scientific questions related to the angle of heat exchanger pipes, where buoyancy effects from flow velocity and pipe orientation influence CO2 frosting. Results show that upwardly angled pipes produce notably different flow structures compared to horizontal or vertical configurations, and that carbon capture efficiency correlates with buoyancy effects for pipe angles within plus or minus 23 degrees of horizontal.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Determining the extent of potential fugitive fluid migration from geologic carbon storage in hydrocarbon-bearing reservoirs: Insights from one-dimensional numerical modeling

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted R f . The area within R f comprises regions of substantially elevated pressure and free-phase CO 2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H 2 O, CO 2 , CH 4 , C 4 H 10 , and C 10 H 22 . The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO 2 . The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO 2 -EOR fields being considered for dedicated CO 2 storage.

CO₂ plume migration↗
NASA NTRS1997

[Tail Plane Icing]

The Aviation Safety Program initiated by NASA in 1997 has put greater emphasis in safety related research activities. Ice-contaminated-tailplane stall (ICTS) has been identified by the NASA Lewis Icing Technology Branch as an important activity for aircraft safety related research. The ICTS phenomenon is characterized as a sudden, often uncontrollable aircraft nose- down pitching moment, which occurs due to increased angle-of-attack of the horizontal tailplane resulting in tailplane stall. Typically, this phenomenon occurs when lowering the flaps during final approach while operating in or recently departing from icing conditions. Ice formation on the tailplane leading edge can reduce tailplane angle-of-attack range and cause flow separation resulting in a significant reduction or complete loss of aircraft pitch control. In 1993, the Federal Aviation Authority (FAA) and NASA embarked upon a four-year research program to address the problem of tailplane stall and to quantify the effect of tailplane ice accretion on aircraft performance and handling characteristics. The goals of this program, which was completed in March 1998, were to collect aerodynamic data for an aircraft tail with and without ice contamination and to develop analytical methods for predicting the effects of tailplane ice contamination. Extensive dry air and icing tunnel tests which resulted in a database of the aerodynamic effects associated with tailplane ice contamination. Although the FAA/NASA tailplane icing program generated some answers regarding ice-contaminated-tailplane stall (ICTS) phenomena, NASA researchers have found many open questions that warrant further investigation into ICTS. In addition, several aircraft manufacturers have expressed interest in a second research program to expand the database to other tail configurations and to develop experimental and computational methodologies for evaluating the ICTS phenomenon. In 1998, the icing branch at NASA Lewis initiated a second multi-phase research program for tailplane icing (TIP II) to develop test methodologies and tailplane performance and handling qualities evaluation tools. The main objectives of this new NASA/Industry/Academia collaborative research programs were: (1) define and evaluate a sub-scale wind tunnel test methodology for determining tailplane performance degradation due to icing. (2) develop an experimental database of tailplane aerodynamic performance with and without ice contamination for a range of tailplane configurations. Wind tunnel tests were planned with representative general aviation aircraft, i.e., the Learjet 45, and a twin engine low speed aircraft. This report summarizes the research performed during the first year of the study, and outlines the work tasks for the second year.

Source record↗
NASA NTRS2013

Numerical Simulations of Instabilities in Single-Hole Office Elements

An orifice element is commonly used in liquid rocket engine test facilities either as a flow metering device, a damper for acoustic resonance or to provide a large reduction in pressure over a very small distance in the piping system. While the orifice as a device is largely effective in stepping down pressure, it is also susceptible to a wake-vortex type instability that generates pressure fluctuations that propagate downstream and interact with other elements of the test facility resulting in structural vibrations. Furthermore in piping systems an unstable feedback loop can exist between the vortex shedding and acoustic perturbations from upstream components resulting in an amplification of the modes convecting downstream. Such was the case in several tests conducted at NASA as well as in the Ariane 5 strap-on P230 engine in a static firing test where pressure oscillations of 0.5% resulted in 5% thrust oscillations. Exacerbating the situation in cryogenic test facilities, is the possibility of the formation of vapor clouds when the pressure in the wake falls below the vapor pressure leading to a cavitation instability that has a lower frequency than the primary wake-vortex instability. The cavitation instability has the potential for high amplitude fluctuations that can cause catastrophic damage in the facility. In this paper high-fidelity multi-phase numerical simulations of an orifice element are used to characterize the different instabilities, understand the dominant instability mechanisms and identify the tonal content of the instabilities.

Ahuja, Vineet↗
NASA NTRS1996

Modeling of Thermal Performance of Multiphase Nuclear Fuel Cell Under Variable Gravity Conditions

A unique numerical method has been developed to model the dynamic processes of bulk evaporation and condensation processes, associated with internal heat generation and natural convection under different gravity levels. The internal energy formulation, for the bulk liquid-vapor phase change problems in an encapsulated container, was employed. The equations, governing the conservation of mass, momentum and energy for both phases involved in phase change, were solved. The thermal performance of a multiphase uranium tetra-fluoride fuel element under zero gravity, micro-gravity and normal gravity conditions has been investigated. The modeling yielded results including the evolution of the bulk liquid-vapor phase change process, the evolution of the liquid-vapor interface, the formation and development of the liquid film covering the side wall surface, the temperature distribution and the convection flow field in the fuel element. The strong dependence of the thermal performance of such multiphase nuclear fuel cell on the gravity condition has been revealed. Under all three gravity conditions, 0-g, 10(exp -3)-g, and 1-g, the liquid film is formed and covers the entire side wall. The liquid film covering the side wall is more isothermalized at the wall surface, which can prevent the side wall from being over-heated. As the gravity increases, the liquid film is thinner, the temperature gradient is larger across the liquid film and smaller across the vapor phase. This investigation provides valuable information about the thermal performance of multi-phase nuclear fuel element for the potential space and ground applications.

Ding, Z.↗
NASA NTRS2004

Simulations of Cavitating Cryogenic Inducers

Simulations of cavitating turbopump inducers at their design flow rate are presented. Results over a broad range of Nss, numbers extending from single-phase flow conditions through the critical head break down point are discussed. The flow characteristics and performance of a subscale geometry designed for water testing are compared with the fullscale configuration that employs LOX. In particular, thermal depression effects arising from cavitation in cryogenic fluids are identified and their impact on the suction performance of the inducer quantified. The simulations have been performed using the CRUNCH CFD[R] code that has a generalized multi-element unstructured framework suitable for turbomachinery applications. An advanced multi-phase formulation for cryogenic fluids that models temperature depression and real fluid property variations is employed. The formulation has been extensively validated for both liquid nitrogen and liquid hydrogen by simulating the experiments of Hord on hydrofoils; excellent estimates of the leading edge temperature and pressure depression were obtained while the comparisons in the cavity closure region were reasonable.

Dorney, Dan↗
NASA NTRS2005

Simulations of inducers at Low-Flow Off-Design Conditions

The ability to accurately model details of inlet back flow for inducers operating at low-flow, off-design conditions is evaluated. A sub-scale version of a three- bladed liquid hydrogen inducer tested in water with detailed velocity and pressure measurements is used as a numerical test bed. Under low-flow, off-design conditions the length of the separation zone as well as the swirl velocity magnitude was under predicted with a standard k-E model. When the turbulent viscosity coefficient was reduced good comparison was obtained at all the flow conditions examined with both the magnitude and shape of the profile matching well with the experimental data taken half a diameter upstream of the leading edge. The velocity profiles and incidence angles at the leading edge itself were less sensitive to the back flow length predictions indicating that single-phase performance predictions may be well predicted even if the details of flow separation modeled are incorrect. However, for cavitating flow situations the prediction of the correct swirl in the back flow and the pressure depression in the core becomes critical since it leads to vapor formation. The simulations have been performed using the CRUNCH CFD@ code that has a generalized multi-element unstructured framework and an advanced multi-phase formulation for cryogenic fluids. The framework has been validated rigorously for predictions of temperature and pressure depression in cryogenic fluid cavities and has also been shown to predict the cavitation breakdown point for inducers at design conditions.

Hosangadi, Ashvin↗
NASA NTRS2005

Simulations of a Liquid Hydrogen Inducer at Low-Flow Off-Design Flow Conditions

The ability to accurately model details of inlet back flow for inducers operating a t low-flow, off-design conditions is evaluated. A sub-scale version of a three-bladed liquid hydrogen inducer tested in water with detailed velocity and pressure measurements is used as a numerical test bed. Under low-flow, off-design conditions the length of the separation zone as well as the swirl velocity magnitude was under predicted with a standard k-E model. When the turbulent viscosity coefficient was reduced good comparison was obtained a t all the flow conditions examined with both the magnitude and shape of the profile matching well with the experimental data taken half a diameter upstream of the leading edge. The velocity profiles and incidence angles a t the leading edge itself were less sensitive to the back flow length predictions indicating that single-phase performance predictions may be well predicted even if the details of flow separation modeled are incorrect. However, for cavitating flow situations the prediction of the correct swirl in the back flow and the pressure depression in the core becomes critical since it leads to vapor formation. The simulations have been performed using the CRUNCH CFD(Registered Trademark) code that has a generalized multi-element unstructured framework and a n advanced multi-phase formulation for cryogenic fluids. The framework has been validated rigorously for predictions of temperature and pressure depression in cryogenic fluid cavities and has also been shown to predict the cavitation breakdown point for inducers a t design conditions.

Hosangadi, A.↗
NASA NTRS2015

Particle Image Velocimetry Applications Using Fluorescent Dye-Doped Particles

Polystyrene latex sphere particles are widely used to seed flows for velocimetry techniques such as Particle Image Velocimetry (PIV) and Laser Doppler Velocimetry (LDV). These particles may be doped with fluorescent dyes such that signals spectrally shifted from the incident laser wavelength may be detected via Laser Induced Fluorescence (LIF). An attractive application of the LIF signal is achieving velocimetry in the presence of strong interference from laser scatter, opening up new research possibilities very near solid surfaces or at liquid/gas interfaces. Additionally, LIF signals can be used to tag different fluid streams to study mixing. While fluorescence-based PIV has been performed by many researchers for particles dispersed in water flows, the current work is among the first in applying the technique to micron-scale particles dispersed in a gas. A key requirement for such an application is addressing potential health hazards from fluorescent dyes; successful doping of Kiton Red 620 (KR620) has enabled the use of this relatively safe dye for fluorescence PIV for the first time. In this paper, basic applications proving the concept of PIV using the LIF signal from KR620-doped particles are exhibited for a free jet and a twophase flow apparatus. Results indicate that while the fluorescence PIV techniques are roughly 2 orders of magnitude weaker than Mie scattering, they provide a viable method for obtaining data in flow regions previously inaccessible via standard PIV. These techniques have the potential to also complement Mie scattering signals, for example in multi-stream and/or multi-phase experiments.

Petrosky, Brian J.↗
NASA NTRS2009

Findings from the Supersonic Qualification Program of the Mars Science Laboratory Parachute System

In 2012, the Mars Science Laboratory Mission (MSL) will deploy NASA's largest extra-terrestrial parachute, a technology integral to the safe landing of its advanced robotic explorer on the surface. The supersonic parachute system is a mortar deployed 21.5 m disk-gap-band (DGB) parachute, identical in geometric scaling to the Viking era DGB parachutes of the 1970's. The MSL parachute deployment conditions are Mach 2.3 at a dynamic pressure of 750 Pa. The Viking Balloon Launched Decelerator Test (BLDT) successfully demonstrated a maximum of 700 Pa at Mach 2.2 for a 16.1 m DGB parachute in its AV4 flight. All previous Mars deployments have derived their supersonic qualification from the Viking BLDT test series, preventing the need for full scale high altitude supersonic testing. The qualification programs for Mars Pathfinder, Mars Exploration Rover, and Phoenix Scout Missions were all limited to subsonic structural qualification, with supersonic performance and survivability bounded by the BLDT qualification. The MSL parachute, at the edge of the supersonic heritage deployment space and 33% larger than the Viking parachute, accepts a certain degree of risk without addressing the supersonic environment in which it will deploy. In addition, MSL will spend up to 10 seconds above Mach 1.5, an aerodynamic regime that is associated with a known parachute instability characterized by significant canopy projected area fluctuation and dynamic drag variation. This aerodynamic instability, referred to as "area oscillations" by the parachute community has drag performance, inflation stability, and structural implications, introducing risk to mission success if not quantified for the MSL parachute system. To minimize this risk and as an alternative to a prohibitively expensive high altitude test program, a multi-phase qualification program using computation simulation validated by subscale test was developed and implemented for MSL. The first phase consisted of 2% of fullscale supersonic wind tunnel testing of a rigid DGB parachute with entry-vehicle to validate two high fidelity computational fluid dynamics (CFD) tools. The computer codes utilized Large Eddy Simulation and Detached Eddy Simulation numerical approaches to accurately capture the turbulent wake of the entry vehicle and its coupling to the parachute bow-shock. The second phase was the development of fluid structure interaction (FSI) computational tools to predict parachute response to the supersonic flow field. The FSI development included the integration of the CFD from the first phase with a finite element structural model of the parachute membrane and cable elements. In this phase, a 4% of full-scale supersonic flexible parachute test program was conducted to provide validation data to the FSI code and an empirical dataset of the MSL parachute in a flight-like environment. The final phase is FSI simulations of the full-scale MSL parachute in a Mars type deployment. Findings from this program will be presented in terms of code development and validation, empirical findings from the supersonic testing, and drag performance during supersonic operation.

EDL↗
NASA NTRS2022

An Investigation Into Some Important Aspects of Droplet Breakup/Vaporization Behavior Caused By a Gas Flow (Both Nonreacting & Detonative Combustion) in a Shock Tube

As a part of the rotating detonation engine (RDE) technology enablement project at NASA Glenn Research Center (GRC), an effort was undertaken to extend our current computational capabilities of OpenNCC in some important ways with the implementation of a modeling approach to account for the droplet breakup caused by a shock-induced gas motion, & a vaporization model valid over a wide range of pressure conditions encountered in multiphase detonation. With the modified code, a study was undertaken to investigate the individual droplet behavior followed by the passage of a shock front. The study is carried out by tracking a sparse group of droplets to gain some understanding of shock induced droplet behavior under various shock strengths & fuel injector conditions. The study also looks into the effect of randomization involved in determining the droplet breakup outcomes. Over a wide range of sub-critical conditions examined, larger droplets are observed to undergo significant changes in droplet behavior following their breakup. However, smaller droplets (10 µm or less ) remain unaffected by any shock induced breakup. In a follow-on work, we investigated the impact of shock and droplet interaction in a detonation study involving both gaseous as well as gas/liquid (droplet clouds) fuel/air stoichiometric mixtures in a simple 3D shock-tube configuration. The droplet clouds are made up of different initial droplet sizes of either 6, 10, or 30 µm. We also investigated the individual droplet behavior followed by the passage of a detonation front. The results represent conditions that lead to both overdriven and C-J (ChapmanJouguet) detonations. Under both test conditions, most of the droplet vaporization is completed within a short distance (duration) behind the detonation front & well within the region of complete combustion observed in a corresponding equivalent gas-phase fuel/air mixture. The impact of the shock-induced droplet breakup is found to be significant in the calculations involving the 30-µm droplets. Subsequent to the breakup, the drop sizes vary from 1 to 10 µm. Another factor that contributed to the observed rapid vaporization is the result of vaporization taking place under supercritical conditions. The overall detonation properties of various droplet clouds (made up of different initial sizes) are similar to those observed in a corresponding gaseous fuel/air mixture. In the calculation involving a gaseous fuel, the calculated C-J detonation velocity is 1822 m/s involving Jet-A/air and φ = 1. In the overdriven detonation, it is 2044 m/s. In the calculations involving droplet clouds, the corresponding detonation velocities are lower. The impact of increased droplet size is primarily seen in a higher reduction in the detonation velocity.

shock/droplet interaction↗