Recent Advancements in Fully Implicit Numerical Methods for Hypersonic Reacting Flows: With Application to Reentry and Arcjet Modeling
No abstract available
SEARCH · Engineering Papers
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.
No abstract available
This paper presents the numerical simulations of confined three-dimensional coaxial water jets. The objectives are to validate the newly proposed nonlinear turbulence models of momentum and scalar transport, and to evaluate the newly introduced scalar APDF and DWFDF equation along with its Eulerian implementation in the National Combustion Code (NCC). Simulations conducted include the steady RANS, the unsteady RANS (URANS), and the time-filtered Navier-Stokes (TFNS); both without and with invoking the APDF or DWFDF equation. When the APDF (ensemble averaged probability density function) or DWFDF (density weighted filtered density function) equation is invoked, the simulations are of a hybrid nature, i.e., the transport equations of energy and species are replaced by the APDF or DWFDF equation. Results of simulations are compared with the available experimental data. Some positive impacts of the nonlinear turbulence models and the Eulerian scalar APDF and DWFDF approach are observed.
The high altitude of Jupiter's Great Red Spot (GRS) may enhance the upward flux of gaseous ammonia (NH3 ) into the high troposphere, where NH3 molecules can be photodissociated and initiate a chain of chemical reactions with downwelling acetylene molecules (C2H2 ). These reactions, experimentally studied earlier by (Ferris and Ishikawa [1987] Nature 326, 777-778) and (Ferris and Ishikawa [1988] J. Amer. Chem. Soc. 110, 4306-4312), produce chromophores that absorb in the visible and ultraviolet regions. In this work we photolyzed mixtures of NH3 and C2H2 using ultraviolet radiation with a wavelength of 214 nm and measured the spectral transmission of the deposited films in the visible region (400-740 nm). From these transmission data we estimated the imaginary indices of refraction. Assuming that ammonia grains at the top of the GRS clouds are coated with this material, we performed layered sphere and radiative transfer calculations to predict GRS reflection spectra. Comparison of those results with observed and previously unreported Cassini visible spectra and with true-color images of the GRS show that the unknown GRS chromophore is spectrally consistent with the coupled NH3-C2H2 photochemical products produced in our laboratory experiments. Using high-resolution mass spectrometry and infrared spectroscopy we infer that the chromophore-containing residue is composed of aliphatic azine, azo, and diazo compounds.
Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure. As one of the nation's largest rocket testing facilities, the performance and design limitations of diffusers are of great interest to NASA's Stennis Space Center. This paper describes a series of tests conducted on four diffuser configurations to better understand the effects of inlet geometry and throat area on starting behavior and boundary layer separation. The diffusers were tested for a duration of five seconds with a 1455-pound thrust, LO2/GH2 thruster to ensure they each reached aerodynamic steady state. The effects of a water spray ring at the diffuser exits and a water-cooled deflector plate were also evaluated. Static pressure and temperature measurements were taken at multiple axial locations along the diffusers, and Computational Fluid Dynamics (CFD) simulations were used as a tool to aid in the interpretation of data. The hot combustion products were confirmed to enable the diffuser start condition with tighter second throats than predicted by historical cold-flow data or the theoretical normal shock method. Both aerodynamic performance and heat transfer were found to increase with smaller diffuser throats. Spray ring and deflector cooling water had negligible impacts on diffuser boundary layer separation. CFD was found to accurately capture diffuser shock structures and full-flowing diffuser wall pressures, and the qualitative behavior of heat transfer. However, the ability to predict boundary layer separated flows was not consistent.
Aerothermodynamics and hypersonic flows involve complex multi-disciplinary physics, including finite-rate gas-phase kinetics, finite-rate internal energy relaxation, gas-surface interactions with finite-rate oxidation and sublimation, transition to turbulence, large-scale unsteadiness, shock-boundary layer interactions, fluid-structure interactions, and thermal protection system ablation and thermal response. Many of the flows have a large range of length and time scales, requiring large computational grids, implicit time integration, and large solution run times. The University of Minnesota NASA US3D code was designed for the simulation of these complex, highly-coupled flows. It has many of the features of the well-established DPLR code, but uses unstructured grids and has many advanced numerical capabilities and physical models for multi-physics problems. The main capabilities of the code are described, the physical modeling approaches are discussed, the different types of numerical flux functions and time integration approaches are outlined, and the parallelization strategy is overviewed. Comparisons between US3D and the NASA DPLR code are presented, and several advanced simulations are presented to illustrate some of novel features of the code.
Hypersonic boundary layer transition is critical to the design of all hypersonic vehicles due to its effect on the heat transfer into the vehicle surface and potential drag enhancement or reduction during reentry. Boundary layer transition and boundary layer stability analysis under hypersonic conditions has been studied for decades, yet there is ample room for improved accuracy and further investigations into the relevant phenomena. In this work, we present a recent implementation of chemical equilibrium, finite-rate chemistry, and thermochemical nonequilibrium capabilities into LASTRAC, an existing well-established boundary-layer stability analysis code. Verification against existing numerical results in the literature are presented. LASTRAC was previously able to address calorically perfect flows. By using solutions of the Parabolized Stability Equations (PSE) with chemical and thermal nonequilibrium, we are able to investigate the effects of chemical and thermal nonequilibrium on a variety of phenomena including stationary crossflow instability on a swept wing and 2nd mode instabilities over a wedge.
The Thermal Protection System (TPS) is a significant portion of the mass of reentry vehicles, planetary probes, and Martian entry vehicles. Reducing this mass has benefits in terms of decreased fuel requirements and increased payload; however, due to the high risk and uncertainty, the TPS or heat shields are designed conservatively by assuming fully turbulent flow. Laminar flow results in reduced heat flux, and improved transition prediction has the potential to reduce TPS mass and uncertainty in aerothermodynamic predictions. Limited previous research exists examining the problem of transition prediction for Martian atmospheric entry, with studies available on transition on the Mars Science Laboratory TPS. Although transition was demonstrated in wind-tunnel tests, uncertainty in the transition location resulted in a TPS designed for fully turbulent flow, and therefore greater mass than required for a partially laminar condition. In this work, we extend the boundary layer stability code LASTRAC, recently modified to include chemical and thermal nonequilibrium capabilities, to a model of the Martian atmosphere. LASTRAC provides Parabolized Stability Equations (PSE) as well as Linear Stability Theory (LST) to predict the stability of a boundary layer and transition with semi-empirical e(sup N) methods. Results included in this work compare disturbance growth characteristics between air and Martian atmosphere at similar non-dimensional freestream conditions on a simple flat plate geometry. Both chemical nonequilibrium and thermochemical nonequilibrium, as well as both PSE and LST, are used.
The Thermal Protection System (TPS) is a significant portion of the mass of reentry vehicles, planetary probes, and Martian entry vehicles. Reducing this mass has benefits in terms of decreased fuel requirements and increased payload; however, due to the high risk and un-certainty, the TPS or heat shields are designed conservatively by assuming fully turbulent flow. Laminar flow results in reduced heat flux, and improved transition prediction has the potential to reduce TPS mass and uncertainty in aerothermodynamic predictions. Limited previous research exists examining the problem of transition prediction for Martian atmospheric entry, with studies available on transition on the Mars Science Laboratory TPS. Although transition was demonstrated in wind-tunnel tests, uncertainty in the transition location resulted in a TPS designed for fully turbulent flow, and therefore greater mass than required for a partially-laminar condition. In this work, we extend the boundary layer stability code LASTRAC, recently modified to include chemical and thermal nonequilibrium capabilities, to a model of the Martian atmosphere. LASTRAC provides Parabolized Stability Equations (PSE) as well as Linear Stability Theory (LST) to predict the stability of a boundary layer and transition with semi-empirical eN methods. Results included in this work compare disturbance growth characteristics between air and Martian atmosphere at similar nondimensional freestream conditions on a simple flat plate geometry. Both chemical nonequilibrium and thermochemical nonequilibrium, as well as both PSE and LST, are used.
Explore the source record for details and available documents.
Project Overview This project is focused on developing technologies that enable construction of a Lunar Safe Haven type structure on the moon. A full scale architectural and structural design will be completed based on lunar conditions and Artemis mission needs. A scaled structure will be 3D printed in simulated lunar environments. Technical Approach Polymer bound regolith composite materials will be developed and characterized for lunar additive construction applications. A full scale architectural and structural design will be completed based on the latest scientific data on lunar environments including radiation, meteoroids, thermal conditions, reduced gravity, and moonquakes. The project will culminate with a demonstration of additive construction in simulated lunar vacuum, thermal, regolith, and UV conditions. Results/Summary Preliminary material specifications have been developed and an initial batch has been produced. A first iteration of the architectural and structural designs have been completed. Hardware for 3D printing in simulated lunar conditions is under development. Contributing Partners AI Space Factory – Winner of the NASA 3D Printed Habitat Centennial Challenge Infusion and Transition Plan On-surface construction of infrastructure is a critical capability of the sustainable phase of the Artemis Program and eventual sustainable human presence on Mars. It provides the ability to create protective shelters on-demand using local resources. The 3D-printing technology will achieve TRL 6 in 1g/vacuum tests with thermal/UV exposure in this project. CLPS missions will be pursued for infusion to demonstrate vertical construction with lunar regolith and conditions. Modeling and simulation will be used to define the pilot-scale shelter construction mission for long-term exposure on the lunar surface to validate the system for Artemis operational deployment.
While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].
Advancing the TRL of technologies developed during the NASA 3D Printed Habitat Centennial Challenge. Focused on the design and construction of a lunar unpressurized radiation protection shelter called LINA.
Upper stage rocket engines are optimized for operation at the edge of space. Ground testing often requires simulation of high-altitude environments to prevent damage from off-design operation. Passive supersonic diffusers are an effective means of establishing the partial vacuum needed around the nozzle exit. The low cost of diffusers relative to active vacuum sources (ejectors, pumps, etc.) provides a perpetual incentive to improve their performance and expand the envelope of passively testable engines. Second-throat diffusers have been the de facto standard for rocket testing since the 1950s, and numerous studies have focused on parametric optimization of the topology. Yet, even the best second-throat designs may require more driving pressure than an upper stage engine can provide. To overcome conventional limits, engineers at NASA’s Stennis Space Center (SSC) adapted an obscure cold-flow diffuser topology to the extremes of immersion in supersonic combustion products. Aerodynamic experimentation was conducted at SSC’s E-3 test stand to quantify the performance of prototypical hot-fire spike diffuser hardware and enable direct comparison to prior states of the art. A LO X /GH 2 thruster was contoured to approximate the regeneratively-cooled nozzle of an upper-stage engine and fired into 28 downstream diffuser configurations: 16 second-throat, 2 centerbody, and 10 spike. Temperature and pressure were recorded along the test article walls for a range of chamber pressures spanning 3.5 - 5.0 MPa. The second-throat configurations with the lowest start and unstart pressure ratios were shown to outperform equivalent systems in literature and taken as baselines for evaluation of advanced diffuser performance. Compared to these exemplars, centerbody diffusers decreased start pressure ratio by 11% but increased unstart pressure ratio by 15%. Spike diffusers reduced start and unstart pressure ratios by 24% and 9%, respectively.
Goal: To demonstrate construction of a civil engineered, unpressurized protective shelter in simulated Lunar conditions using regolith polymer composites in a Fused Filament Fabrication (FFF) additive process
The ability to recover oxygen from metabolic carbon dioxide is crucial for crewed long duration missions in which the atmosphere must be recycled. Two physiochemical methods are primarily employed for this purpose: the Bosch process and the Sabatier reaction. This paper focuses on the Bosch process. It utilizes hydrogen in the presence of a catalyst to form water and solid carbon out of metabolic carbon dioxide. The solid carbon product deposits onto the catalyst and eventually causes reactor over pressurization issues. Multiple methods are currently being investigated to regenerate the catalyst’s surface and break up the pressure causing carbon clogs. This paper explores whether water and mechanical abrasion would work to remove Bosch carbon by testing multiple washing conditions. The different washing conditions involved agitation intensity, heating, duration, and volume of water. It was found that it is possible to break up the clumps of carbon formed in the reactor through agitation with room temperature water in one hour and shows promise for completely cleaning the carbon off the beads’ surface.
Static tests were made on the compatibility of liquid fluorine with several nonmetallic materials at -3200 F and at pressures of 0 and 1500 pounds per square inch gage. The results are compared with those from previous work with gaseous fluorine at the same pressures, but at atmospheric temperature. In general, although environmental effects were not always consistent, reactivity was least with the low-temperature, low-pressure liquid fluorine. Reactivity was greatest with the warm, high-pressure gaseous fluorine. None of the liquids and greases tested was found to be entirely suitable for use in fluorine systems. Polytrifluorochloroethylene and N-43, the formula for which is (C4F9)3N, did not react with liquid fluorine at atmospheric pressure or 1500 pounds per square inch gage under static conditions, but they did react when injected into liquid fluorine at 1500 pounds per square inch gage; they also reacted with gaseous fluorine at 1500 pounds per square inch gage. While water did not react with liquid fluorine at 1500 pounds per square inch gage, it is known to react violently with fluorine under other conditions. The pipe-thread lubricant Q-Seal did not react with liquid fluorine, but did react with gaseous fluorine at 1500 pounds per square inch gage. Of the solids, ruby (Al2O3) and Teflon did not react under the test conditions. The results show that the compatibility of fluorine with nonmetals depends on the state of the fluorine and the system design.
A tool that would be useable in both conventional and self-reacting friction stir welding (FSW) has been proposed. The tool would embody both a prior tooling concept for self-reacting FSW and an auto-adjustable pin-tool (APT) capability developed previously as an augmentation for conventional FSW. Some definitions of terms are prerequisite to a meaningful description of the proposed tool. In conventional FSW, depicted in Figure 1, one uses a tool that includes (1) a rotating shoulder on top (or front) of the workpiece and (2) a rotating pin that protrudes from the shoulder into the depth of the workpiece. The main axial force exerted by the tool on the workpiece is reacted through a ridged backing anvil under (behind) the workpiece. When conventional FSW is augmented with an APT capability, the depth of penetration of the pin into the workpiece is varied in real time by a position- or force-control system that extends or retracts the pin as needed to obtain the desired effect. In self-reacting (also known as self-reacted) friction stir welding (SR-FSW), there are two rotating shoulders: one on top (or front) and one on the bottom (or back) of the workpiece. In this case, a threaded shaft protrudes from the tip of the pin to beyond the back surface of the workpiece. The back shoulder is held axially in place against tension by a nut on the threaded shaft. The main axial force exerted on the workpiece by the tool and front shoulder is reacted through the back shoulder and the threaded shaft, back into the FSW machine head, so that a backing anvil is no longer needed. A key transmits torque between the bottom shoulder and the threaded shaft, so that the bottom shoulder rotates with the shaft. A tool for SRFSW embodying this concept was reported in "Mechanism for Self-Reacted Friction Stir Welding" (MFS-31914), NASA Tech Briefs, Vol. 28, No. 10 (October 2004), page 53. In its outward appearance, the proposed tool (see Figure 2) would fit the above description of an SR-FSW tool. In this case, the FSW machine would have an APT capability and the pin would be modified to accept a bottom shoulder. The APT capability could be used to vary the distance between the front and back shoulders in real time to accommodate process and workpiece-thickness variations. The tool could readily be converted to a conventional FSW tool, with or without APT capability, by simply replacing the modified pin with a conventional FSW pin.