Search NASA⌕ Search

SEARCH · Search NASA

Results for “multiphase modeling”

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 109 records · Page 6

CFD Jet Mixing Model Validation Against Zero-Boil-Off Tank (ZBOT) Microgravity Experiment

The Zero-Boil-Off Tank (ZBOT) experiments are a series of small scale two-phase fluids physics experiments conducted aboard the International Space Station (ISS) in 2017. A goal of these experiments is to study various fundamental fluid flow, heat and mass transport, and phase change phenomena that control storage tank pressurization and pressure control in microgravity. In this paper, results of a state-of-the-art two-phase CFD (Computational Fluid Dynamics) storage tank model are presented and validated against jet mixing data provided by microgravity ZBOT experiments. Jet-ullage interaction was studied in the experiment for different jet speeds representing different flow regimes, and tank fill levels. Effects of various physical parameters, such as, jet angle, gravitational acceleration and surface tension value were investigated. Computational parameters, such as turbulence and surface tension models and interface capturing scheme, affecting jet-ullage interaction were also studied. Computational results were compared to experimental results for model validation.

Turbulence↗

CFD Jet Mixing Model Validation against Zero-Boil-Off Tank (ZBOT) Microgravity Experiment

The Zero-Boil-Off Tank (ZBOT) experiments are a series of small scale two-phase fluids physics experiments conducted aboard the International Space Station (ISS) in 2017. A goal of these experiments is to study various fundamental fluid flow, heat and mass transport, and phase change phenomena that control storage tank pressurization and pressure control in microgravity. In this paper, results of a state-of-the-art two-phase CFD storage tank model are presented and validated against jet mixing data provided by microgravity ZBOT experiments. Jet-ullage interaction was studied in the experiment for different jet speeds representing different flow regimes, and tank fill levels. Effects of various physical parameters, such as, jet angle, gravitational acceleration and surface tension value were investigated. Computational parameters, such as turbulence and surface tension models and interface capturing scheme, affecting jet-ullage interaction were also studied. Computational results were compared to experimental results for model validation.

Mixing↗

Aerosols in Atmospheric Chemistry and Biogeochemical Cycles of Nutrients

Atmospheric aerosols have complex and variable compositions and properties.While scientific interest is centered on the health and climatic effects of atmospheric aerosols, insufficient attention is given to their involvement in multiphase chemistry that alters their contribution as carriers of nutrients in ecosystems. However, there is experimental proof that the nutrient equilibria of both land and marine ecosystems have been disturbed during the Anthropocene period. This review study first summarizes our current understanding of aerosol chemical processing in the atmosphere as relevant to biogeochemical cycles. Then it binds together results of recent modeling studies based on laboratory and field experiments, focusing on the organic and dust components of aerosols that account for multiphase chemistry, aerosol ageing in the atmosphere, nutrient (N, P, Fe) emissions, atmospheric transport, transformation and deposition. The human-driven contribution to atmospheric deposition of these nutrients, derived by global simulations using past and future anthropogenic emissions of pollutants, is put into perspective with regard to potential changes in nutrient limitations and biodiversity. Atmospheric deposition of nutrients has been suggested to result in human-induced ecosystem limitations with regard to specific nutrients. Such modifications favor the development of certain species against others and affect the overall functioning of ecosystems. Organic forms of nutrients are found to contribute to the atmospheric deposition of the nutrients N, P and Fe by 20%-40%, 35%-45% and 7%-18%, respectively. These have the potential to be key components of the biogeochemical cycles since there is initial proof of their bioavailability to ecosystems. Bioaerosols have been found to make a significant contribution to atmospheric sources of N and P, indicating potentially significant interactions between terrestrial and marine ecosystems. These results deserve further experimental and modeling studies to reduce uncertainties and understand the feedbacks induced by atmospheric deposition of nutrients to ecosystems.

dust↗

Highly integrated digital electronic control: Digital flight control, aircraft model identification, and adaptive engine control

The highly integrated digital electronic control (HIDEC) program at NASA Ames Research Center, Dryden Flight Research Facility is a multiphase flight research program to quantify the benefits of promising integrated control systems. McDonnell Aircraft Company is the prime contractor, with United Technologies Pratt and Whitney Aircraft, and Lear Siegler Incorporated as major subcontractors. The NASA F-15A testbed aircraft was modified by the HIDEC program by installing a digital electronic flight control system (DEFCS) and replacing the standard F100 (Arab 3) engines with F100 engine model derivative (EMD) engines equipped with digital electronic engine controls (DEEC), and integrating the DEEC's and DEFCS. The modified aircraft provides the capability for testing many integrated control modes involving the flight controls, engine controls, and inlet controls. This paper focuses on the first two phases of the HIDEC program, which are the digital flight control system/aircraft model identification (DEFCS/AMI) phase and the adaptive engine control system (ADECS) phase.

Baer-Riedhart, Jennifer L.↗

A Design Tool for Liquid Rocket Engine Injectors

A practical design tool for the analysis of flowfields near the injector face has been developed and used to analyze the Fastrac engine. The objective was to produce a computational design tool which was detailed enough to predict the interactive effects of injector element impingement angles and points and the momenta of the individual orifice flows. To obtain a model which could be used to simulate a significant number of individual orifices, a homogeneous computational fluid dynamics model was developed. To describe liquid and vapor sub- and super-critical flows, the model included thermal and caloric equations of state which were valid over a wide range of pressures and temperatures. A homogeneous model was constructed such that the local state of the flow was determined directly, i.e. the quality of the flow was calculated. Such a model does not identify drops or their distribution, but it does allow the flow along the injector face and into the acoustic cavity to be predicted. It also allows the film coolant flow to be accurately described. The initial evaluation of the injector code was made by simulating cold flow from an unlike injector element and from a like-on-like overlapping fan (LOL) injector element. The predicted mass flux distributions of these injector elements compared well to cold flow test results. These are the same cold flow tests which serve as the data base for the JANNAF performance prediction codes. The flux distributions 1 inch downstream of the injector face are very similar; the differences were somewhat larger at further distances from the faceplate. Since the cold flow testing did not achieve good mass balances when integrations across the entire fan were made, the CFD simulation appears to be reasonable alternative to future cold flow testing. To simulate the Fastrac, an RP-1/LOX combustion model must be chosen. This submodel must be relatively simple to accomplish three-dimensional, multiphase flow simulations. Single RP-1 pyrolysis and partial oxidation steps were chosen and the combustion was completed with the wet CO mechanism. Soot was also formed with a single global reaction. To validate the combustion submodel, global data from gas generator tests and from subscale motor test were used to predict qualitatively correct mean molecular weights, temperature, and soot levels. Because such tests do not provide general kinetics rates, the methodology is not necessarily appropriate for other than rocket type flows conditions. Soot predictions were made so that radiation heating to the motor walls can be made. These initial studies of the Fastrac were for a small region close to the injector face and chamber wall which included a segment of the acoustic cavity. The region analyzed includes 11 individual orifice holes to represent the LOL elements and the H2 film coolant holes. Typical results of this simulation are shown in Figure 1. At this point the only available test data to verify the predictions are temperatures measured in the acoustic cavity. These temperatures are in reasonable agreement at about 2000R (1111 K). Future work is expected to include improving the computational efficiency or the CFD model and/or using more computer capacity than the single Pentium PC with which these simulations were made.

Farmer, Richard C.↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

Modeling Nitrogen Oxides in the Lower Stratosphere

This talk will focus on the status of current understanding (not a historical review) as regards modeling nitrogen oxides (NOy) in the lower stratosphere (LS). The presentation will be organized around three major areas of process understanding: 1) NOy sources, sinks, and transport to the LS, 2) NOy species partitioning, and 3) polar multiphase processes. In each area, process topics will be identified with an estimate of the degree of confidence associated with their representation in numerical models. Several exotic and/or speculative processes will also be discussed. Those topics associated with low confidence or knowledge gaps, weighted by their prospective importance in stratospheric chemical modeling, will be collected into recommendations for further study. Suggested approaches to further study will be presented for discussion.

Kawa, S. Randy↗

CFD Simulation of Ignition Overpressure and Plume-Water Interaction in the SLS Scale Model Acoustics Test

The solid rocket boosters (SRBs) on the Space Launch System (SLS) generate powerful ignition overpressure (IOP) waves that can potentially damage the vehicle and surrounding ground structures. Based on historical Shuttle experience, the IOP will be mitigated by utilizing an IOP/Sound Suppression (SS) water system which sprays water in and around the SRB and RS-25 engine plumes. However, the complex multiphase physics of the resulting plume-water interaction can lead to undesirable effects, such as high velocity water spray onto the engine nozzles and SRB nozzle side loads. Prior to launch, insight into this environment can be achieved through scale model testing and predictive computational fluid dynamics (CFD) modeling. The Scale Model Acoustic Test (SMAT) was conducted at Marshall Space Flight Center to provide the data necessary to develop the SLS liftoff acoustic environments by using a 5% scale SLS vehicle, fully integrated with a scaled IOP/SS water system, mobile launcher (ML), and tower. NASA engineers have put forth significant effort in validating the Loci/CHEM-Real Fluids multiphase CFD solver on a wet hold down test from the SMAT series. In this paper, the IOP acoustics and plume-water interactions predicted by the CFD model are compared and validated against the experiment, enabling relatively inexpensive, rapid turnaround for multiphase launch environment analyses to support future flights and design evolutions of SLS and other launch vehicles.

Computational Fluid Dynamics↗

Research on turbojet combustor pollutant formation

The research for developing a modular computer program for describing gas-turbine combustor flow and chemical fields is reported. The status of investigations is presented for three key elements of the unified model: turbulent mixing with swirl, turbulent mixing with swirl including finite-rate kinetics, and multiphase recirculation zone analysis

Source record↗

Radiation/convection coupling in rocket motor and plume analysis

A method for describing radiation/convection coupling to a flow field analysis was developed for rocket motors and plumes. The three commonly used propellant systems (H2/O2, RP-1/O2, and solid propellants) radiate primarily as: molecular emitters, non-scattering small particles (soot), and scattering larger particles (Al2O3), respectively. For the required solution, the divergence of the radiation heat flux was included in the energy equation, and the local, volume averaged intensity was determined by a solution to the radiative transfer equation. A rigorous solution to this problem is intractable, therefore, solution methods which use the ordinary and improved differential approximation were developed. This radiation model was being incorporated into the FDNS code, a Navier-Stokes flowfield solver for multiphase, turbulent combusting flows.

Saladino, A. J.↗

Contribution of the Recent AUSM Schemes to the Overflow Code: Implementation and Validation

We shall present results of a recent collaborative effort between the authors attempting to implement the numerical flux scheme, AUSM+ and its new developments, into a widely used NASA code, OVERFLOW. This paper is intended to give a thorough and systematic documentation about the solutions of default test cases using the AUSNI+ scheme. Hence we will address various aspects of numerical solutions, such as accuracy, convergence rate, and effects of turbulence models, over a variety of geometries, speed regimes. We will briefly describe the numerical schemes employed in the calculations, including the capability of solving for low-speed flows and multiphase flows by employing the concept of numerical speed of sound. As a bonus, this low Mach number formulations also enhances convergence to steady solutions for flows even at transonic speed. Calculations for complex 3D turbulent flows were performed with several turbulence models and the results display excellent agreements with measured data.

Liou, Meng-Sing↗

Shock waves and shock tubes; Proceedings of the Fifteenth International Symposium, Berkeley, CA, July 28-August 2, 1985

A detailed survey is presented of shock tube experiments, theoretical developments, and applications being carried out worldwide. The discussions explore shock tube physics and the related chemical, physical and biological science and technology. Extensive attention is devoted to shock wave phenomena in dusty gases and other multiphase and heterogeneous systems, including chemically reactive mixtures. Consideration is given to techniques for measuring, visualizing and theoretically modeling flowfield, shock wave and rarefaction wave characteristics. Numerical modeling is explored in terms of the application of computational fluid dynamics techniques to describing flowfields in shock tubes. Shock interactions and propagation, in both solids, fluids, gases and mixed media are investigated, along with the behavior of shocks in condensed matter. Finally, chemical reactions that are initiated as the result of passage of a shock wave are discussed, together with methods of controlling the evolution of laminar separated flows at concave corners on advanced reentry vehicles.

Bershader, D.↗

Development and testing of a unique carousel wind tunnel to experimentally determine the effect of gravity and the interparticle force on the physics of wind-blown particles

In the study of planetary aeolian processes the effect of gravity is not readily modeled. Gravity appears in the equations of particle motion along with the interparticle forces but the two are not separable. A wind tunnel that perimits multiphase flow experiments with wind blown particles at variable gravity was built and experiments were conducted at reduced gravity. The equations of particle motion initiation (saltation threshold) with variable gravity were experimentally verified and the interparticle force was separated. A uniquely design Carousel Wind Tunnel (CWT) allows for the long flow distance in a small sized tunnel since the test section if a continuous loop and develops the required turbulent boundary layer. A prototype model of the tunnel where only the inner drum rotates was built and tested in the KC-135 Weightless Wonder 4 zero-g aircraft. Future work includes further experiments with walnut shell in the KC-135 which sharply graded particles of widely varying median sizes including very small particles to see how interparticle force varies with particle size, and also experiments with other aeolian material.

Leach, R. N.↗

Pressure-Fed LOX/LCH4 Reaction Control System for Spacecraft: Transient Modeling and Thermal Vacuum Hotfire Test Results

An integrated cryogenic liquid oxygen, liquid methane (LOX/LCH4) reaction control system (RCS) was tested at NASA Glenn Research Center's Plum Brook Station in the Spacecraft Propulsion Research Facility (B-2) under vacuum and thermal vacuum conditions. The RCS is a subsystem of the Integrated Cryogenic Propulsion Test Article (ICPTA), a pressure-fed LOX/LCH4 propulsion system composed of a single 2,800 lbf main engine, two 28 lbf RCS engines, and two 7 lbf RCS engines. Propellants are stored in four 48 inch diameter 5083 aluminum tanks that feed both the main engine and RCS engines in parallel. Helium stored cryogenically in a composite overwrapped pressure vessel (COPV) flows through a heat exchanger on the main engine before being used to pressurize the propellant tanks to a design operating pressure of 325 psi. The ICPTA is capable of simultaneous main engine and RCS operation. The RCS engines utilize a coil-on-plug (COP) ignition system designed for operation in a vacuum environment, eliminating corona discharge issues associated with a high voltage lead. There are two RCS pods on the ICPTA, with two engines on each pod. One of these two engines is a heritage flight engine from Project Morpheus. Its sea level nozzle was removed and replaced by an 85:1 nozzle machined using Inconel 718, resulting in a maximum thrust of 28 lbf under altitude conditions. The other engine is a scaled down version of the 28 lbf engine, designed to match the core and overall mixture ratios as well as other injector characteristics. This engine can produce a maximum thrust of 7 lbf with an 85:1 nozzle that was additively manufactured using Inconel 718. Both engines are film-cooled and capable of limited duration gas-gas and gas-liquid operation, as well as steady-state liquid-liquid operation. Each pod contains one of each version, such that two engines of the same thrust level can be fired as a couple on opposite pods. The RCS feed system is composed of symmetrical 3/8 inch lines that tap off of the main propellant manifold to send LOX and LCH4 outboard to the RCS pods. A Thermodynamic Vent System (TVS) is used to condition propellants at each pod by venting through an orifice and then routing the cold expansion products back through tubing that is welded along a large portion of the main RCS feed lines. Prior to final installation on the ICPTA, the RCS engines were tested in a small vacuum chamber at the Johnson Space Center (JSC) Energy Systems Test Area (ESTA) to verify functionality of the new COP ignition system and check out operation of the vacuum nozzles. After engine-level testing, the RCS engines were installed on the vehicle and a series of integrated hot-fire tests were performed at JSC consisting of various pulsing and steady-state firings as well as integrated main engine/RCS operation. The ICPTA was then integrated into the Plum Brook B-2 facility for vacuum and thermal/vacuum testing. Testing in the B-2 facility was composed of multiple thermal and pressure environments. The first set of tests were performed under ambient temperature and altitude pressure conditions. These tests consisted of a range of minimum impulse bit (MIB) pulsing sequences with low duty cycle, analogous to a coast phase in which the RCS is primarily used for station keeping. The primary goal of this sequence is to understand how propellant conditions were effected without an active TVS. In this scenario, consistent gas-gas operation is desirable since it results in a smaller MIB and more efficient propellant consumption. Multiple skin thermocouples are mounted on the feedlines, in addition to a submerged thermocouple on each commodity, in order to gather thermal data on the system. Higher duty cycle pulsing tests were then performed, analogous to an ascent or landing mission phase. The primary goal of this sequence was to examine how well the engines self-conditioned without active TVS when starting from a quiescent state. The TVS was then activated during some tests to demonstrate the capability to quickly condition the engines for higher pulsing demand scenarios. A thermocouple at the TVS outlet allows for the calculation of energy absorbed by the vented propellant. Lastly, tests with longer pulses and multiple engines firing either in sequence or simultaneously were run in order to gather transient system response data on waterhammer. Six total high-speed pressure transducers are installed on the RCS system, one sensor at the end of each propellant manifold line on the pods, and one at the tap-off location for each commodity. This will allow for the accurate characterization of waterhammer in the system under various propellant conditions and firing sequences. Other instrumentation for this test series includes nozzle throat thermocouples, chamber pressure measurement, heat soakback measurement, and tank wall plume impingement temperature measurement. The next set of tests were performed to demonstrate simultaneous main engine and RCS operation. Data from this test will be used to examine if there is any change to nominal operation of the RCS as a result of feed system interaction or other phenomenon. Some of these tests began under high vacuum conditions (target ambient pressure less than 1x10(exp -3) torr) and others began at altitude conditions. The last set of tests were performed with the B-2 cold wall active. Under these tests, many of the same low duty cycle MIB tests were repeated in order to characterize how propellant conditions changed with the lower heat leak. In this scenario the RCS manifold experiences much less heat leak, resulting in a change to how well the engines self-condition. As a result, an increase in maximum waterhammer pressures and a change in natural frequency of the system was expected due to higher density propellants. The lower heat leak should also result in a change to the MIB pulse profile, and data will be examined to understand how MIB repeatability is affected in the different operating environments. Parallel to the test efforts, a set of transient model development efforts were made to predict RCS performance. The primary effort was aimed at producing a SINDA/FLUINT model to predict propellant conditioning up to the engine inlet as a function of different environmental and operating parameters, with the goal of predicting chamber pressure, TVS performance, and propellant consumption over time. Preliminary results for this effort will be presented in comparison with test data. Additional modeling efforts were made using SINDA/FLUINT to predict waterhammer in the system since the software is capable of handling multiphase transient fluid dynamics. These results will be compared with the high-speed pressure transducer test data for validation purposes.

Atwell, Matthew J.↗

Radar scene generation for tactical decision aids

The Mueller matrix and polarization covariance matrix for polarimetric radar systems was studied. The clutter is modeled by a layer of random permittivity, described by a three-dimensional correlation function, with variance, and horizontal and vertical correlation lengths. A general mixing formula was derived for discrete scatters immersed in a host medium. The results are applicable to general multiphase mixtures. The strong fluctuation theory was used to derive the backscattering cross sections, and was further extended to include higher order co-polarized and cross-polarized moments. A two-layer anisotropic random medium model was developed for the active and passive remote sensing of ice fields. A three-layer random medium model was adopted to study the volume scattering effects for the active and passive microwave remote sensing of snow-covered ice fields. The snow layer was simulated by an isotropic random medium and the ice layer by an anisotropic random medium. The vegetation canopy and snow-covered ice field were studied with a three-layer model, an isotropic random medium layer overlying an anisotropic random medium. The dyadic Green's function of the three-layer medium and the scattered electromagnetic intensities with Born approximation were calculated.

Kong, J. A.↗

CFD Modeling of Cryogenic Chilldown in a Complex Channel under Normal and Low Gravity Conditions

Future NASA architectures have baselined cryogenic propulsion systems as well as cryogenic fluid management to support lunar missions and ultimately to support future missions to Mars. These missions will require chilling hardware down prior to engine restart as well as chilling lines and tanks prior to transferring and refueling these propulsion elements in orbit. In lieu of expensive tests conducted on-orbit, accurate predictive computational models of these chilldown processes can be used to reduce system and propellant mass as well as mission risk. To gain confidence in these computational models, appropriate anchoring and validation to experimental data in a relevant environment needs to be performed. Recent ground and sub-orbital flight experiments conducted by the Japan Aerospace Exploration Agency (JAXA) investigated chilldown of a complex channel resembling a turbopump bearing cavity at low flow rates. This work presents Computational Fluid Dynamics (CFD) model development of the chilldown experiment employing two-phase flow boiling models available in commercial CFD software STAR-CCM+ using Volume of Fluid (VOF) and the traditional Euler-Euler multiphase flow solvers. Comparisons of the numerical and experimental results under normal and low-gravity conditions are presented. An assessment of solid wall temperatures and phase distribution yielded important insights into multiphase solver choice, dependence on gravity environment, and challenges associated with cryogenic flow boiling prediction and validation.

CFD↗

Base Heating Sensitivity Study for a 4-Cluster Rocket Motor Configuration in Supersonic Freestream

In support of launch vehicle base heating and pressure prediction efforts using the Loci-CHEM Navier-Stokes computational fluid dynamics solver, 35 numerical simulations of the NASA TND-1093 wind tunnel test have been modeled and analyzed. This test article is composed of four JP-4/LOX 500 lbf rocket motors exhausting into a Mach 2 - 3.5 wind tunnel at various ambient pressure conditions. These water-cooled motors are attached to a base plate of a standard missile forebody. We explore the base heating profiles for fully coupled finite-rate chemistry simulations, one-way coupled RAMP (Reacting And Multiphase Program using Method of Characteristics)-BLIMPJ (Boundary Layer Integral Matrix Program - Jet Version) derived solutions and variable and constant specific heat ratio frozen flow simulations. Variations in turbulence models, temperature boundary conditions and thermodynamic properties of the plume have been investigated at two ambient pressure conditions: 255 lb/sq ft (simulated low altitude) and 35 lb/sq ft (simulated high altitude). It is observed that the convective base heat flux and base temperature are most sensitive to the nozzle inner wall thermal boundary layer profile which is dependent on the wall temperature, boundary layer s specific energy and chemical reactions. Recovery shock dynamics and afterburning significantly influences convective base heating. Turbulence models and external nozzle wall thermal boundary layer profiles show less sensitivity to base heating characteristics. Base heating rates are validated for the highest fidelity solutions which show an agreement within +/-10% with respect to test data.

Mehta, Manish↗