Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “AEROTHERMODYNAMICS”

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 739 records · Page 41

Progress Towards the CFD Vision 2030

Outline: NASA Motivation; CFD Vision 2030 Study Summary - The Vision - Roadmap - Validation Experiments - CFD 2030 Activities Since Report Release; T3 Technical Challenges Supporting CFD Vision 2030 - Revolutionary Computational Aerosciences (RCA) [Aerodynamics - Aerothermodynamics - Aerostructures - Aeroacoustics - Propulsion Integration] - Multidisciplinary Design, Analysis, and Optimization (MDAO) - Combustion Modeling - Innovative Measurements (Proposed); Implementation Plan, Going Forward; Summary.

Vision 2030↗

Radiative Heating of Large Meteoroids During Atmospheric Entry

A high-fidelity approach for simulating the aerothermodynamic environments of meteor entries was developed, which allows the commonly assumed heat transfer coefficient of 0.1 to be assessed. This model uses chemically reacting computational fluid dynamics (CFD), coupled with radiation transport and surface ablation. Coupled radiation accounts for the impact of radiation on the flowfield energy equations, while coupled ablation explicitly models the injection of ablation products within the flowfield and radiation simulations. For a meteoroid with a velocity of 20 km/s, coupled radiation is shown to reduce the stagnation point radiative heating by over 60%. The impact of coupled ablation (with coupled radiation) is shown to provide at least a 70% reduction in the radiative heating relative to cases with only coupled radiation. This large reduction is partially the result of the low ionization energies of meteoric ablation products relative to air species. The low ionization energies of ablation products, such as Mg and Ca, provide strong photoionization and atomic line absorption in regions of the spectrum that air species do not. MgO and CaO are also shown to provide significant absorption. Turbulence is shown to impact the distribution of ablation products through the shock-layer, which results in up to a 100% increase in the radiative heating downstream of the stagnation point. To create a database of heat transfer coefficients, the developed model was applied to a range of cases. This database considered velocities ranging from 14 to 20 km/s, altitudes ranging from 20 to 50 km, and nose radii ranging from 1 to 100 m. The heat transfer coefficients from these simulations are below 0.045 for the range of cases, for both laminar and turbulent, which is significantly lower than the canonical value of 0:1. When the new heat transfer model is applied to a Tunguska-like 15 Mt entry, the effect of the new model is to lower the height of burst by up to 2 km, depending on assumed entry angle. This, in turn, results in a significantly larger ground damage footprint than when the canonical heating assumption is used.

Christopher O Johnston↗

Turbulent Aeroheating Measurements on a 7-deg Half-Angle Sphere-Cone in a High-Enthalpy CO2 Expansion Tunnel

A database of heating and pressure measurements on a 7-deg half-angle cone in a highenthalpy expansion tunnel in CO2 has been generated to support development and validation of computational models to be employed in the design of future Mars missions. Laminar, transitional, and turbulent simulations were performed at the test conditions for comparisons with the data. Close agreement was obtained for both fully-laminar and fully turbulent conditions. For the remaining transitional/turbulent conditions, agreement to within, or slightly more than, the estimated experimental uncertainty was demonstrated. The influence of transition intermittency and transition length models on predicted heating levels was demonstrated, as were differences in turbulent heating predictions generated using various algebraic, one-equation, and two-equation turbulence models. These comparisons provide some measure of confidence in turbulent simulation capabilities; however, because the data were not obtained on a relevant entry vehicle geometry, it is not possible to fully quantify computational uncertainties for the definition of Mars mission aerothermodynamic environments at this time

Brian R Hollis↗

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Comparisons are made between potential energy surfaces (PES) for N2 + N and N2 + N2 collisions and between rate coefficients for N2 dissociation that were computed using the quasiclassical trajectory method (QCT) on these PESs. For N2 + N we compare the Laganas empirical LEPS surface with one from NASA Ames Research Center based on ab initio quantum chemistry calculations. For N2 + N2 we compare two ab initio PESs (from NASA Ames and from the University of Minnesota). These use different methods for computing the ground state electronic energy for N4, but give similar results. Thermal N2 dissociation rate coefficients, for the 10,000K-30,000K temperature range, have been computed using each PES and the results are in excellent agreement. Quasi-stationary state (QSS) rate coefficients using both PESs have been computed at these temperatures using the Direct Molecular Simulation of Schwartentruber and coworkers. The QSS rate coefficients are up to a factor of 5 lower than the thermal ones and the thermal and QSS values bracket the results of shock-tube experiments. We conclude that the combination of ab initio quantum chemistry PESs and QCT calculations provides an attractive approach for the determination of accurate high-temperature rate coefficients for use in aerothermodynamics modeling.

quasiclassical trajectory method↗

Boundary Layer Stability and Transition in a Chemically Reacting Martian Atmosphere Using LASTRAC

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.

Kline, H. L.↗

Estimating Drag and Heating Coefficients for Hollow Reentry Objects in Transitional Flow Using DSMC

In NASA’s Object Reentry Survival Analysis Tool (ORSAT), aerodynamic drag and aerothermal heating coefficients are computed for each of the free-molecular, continuum, and transitional flow regimes using analytical and semi-analytical methods. These methods are typically limited to convex, blunt objects (such as spheres) and are applied to other objects such as boxes and cylinders using multiplicative “shape factors” to account for the different behavior. Previous literature has analyzed the aerodynamic and aerothermodynamic properties of flow around sharp-edged objects like boxes and cylinders in transitional flow, though only those objects with solid external boundaries. However, many reentry objects we have encountered in real spacecraft have been hollow (i.e., with the potential to allow flow through them). We present here preliminary results from analyses performed using the NASA Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) on hollow cylinders and boxes (with varying wall thickness-diameter ratio).

Marichalar, Jeremiah J.↗

Effect of Out-Gassing on the Onset of Hypersonic Boundary Layer Transition

Prediction and control of the onset of transition and the associated variation in aerothermodynamic parameters in high-speed flows is key to optimize the performance and design of Thermal Protection Systems (TPS) of next-generation aerospace vehicles [1]. Boundary Layer Transition (BLT) characteristics can influence the surface heating budget determining the TPS thickness and consequently its weight penalty. Ablative heatshields are designed to alleviate the high heat flux at the surface through pyrolysis of their polymeric matrix and subsequent fiber ablation [2]. Pyrolysis leads to out-gassing and non-uniform ablation lead to surface roughness, both of which are known to influence the transition process. An ablator impacts BLT through three main routes: gas injecting into the boundary layer from the wall, changing the surface heat transfer due to wall-flow chemical reactions, and modifying surface roughness [3]. In preparation to Mars 2020 mission post-flight analysis, the predictive transition capability has been initiated toward hard-coupling porous material response analysis and aerothermal environment calculation.

Karimi, Mona↗

Effect of Distributed Patch of Smooth Roughness Elements on Transition in a High-Speed Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of hypersonic vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, smooth-shaped roughness elements on the laminar boundary layer over a 7-degree half-angle, circular cone for flow conditions corresponding to a selected trajectory point from the ascent phase of the HIFiRE-1 flight experiment. For peak-to-valley roughness heights of 50 percent or less in comparison with the thickness of the unperturbed boundary layer, the computations converge to a stationary flow, suggesting that the flow is globally stable. Analysis of convective instabilities in the wake of the roughness patch indicates two dominant families of unstable disturbances, namely, a high frequency mode that corresponds to Mack mode waves modified by the wake and a lower frequency mode that corresponds to shear layer instabilities associated with the streaks in the roughness wake. Even though the peak growth rate of the later mode is more than 35 percent greater than the peak growth rates of the Mack modes, the latter modes achieve higher amplification ratios, and hence, are likely to dominate the onset of transition, which is estimated to occur slightly later than that in the unperturbed, i.e., smooth surface boundary layer. Additional computations are performed to investigate the effects of various roughness patch configurations on a Mach 3.5 flat plate boundary layer, to help guide an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at NASA Langley Research Center. In this case, the cumulative reinforcement of basic state distortion over the length of the roughness patch is predicted to yield a significantly earlier transition than that over a smooth plate or a plate with a shorter length roughness patch.

Meelan Choudhari↗

Effects of Wall Cooling on Supersonic Modes in High-Enthalpy Hypersonic Boundary Layers over a Cone

High-enthalpy hypersonic boundary-layer transition plays an important role in many entry/descent vehicles. Aerothermodynamic performance of these vehicles strongly depends on the transition location on the surface. However, detailed transition flow physics in these chemically reacting boundary layers are poorly understood and transition estimates during the design phase rely heavily on empirically-derived transition criteria such as 𝑹𝒆(𝜽)/𝑴(𝒆). One of the most intriguing characteristics in hypersonic boundary layers is the presence of unstable supersonic modes, first identified in the 1990s. Due to a recent surge in hypersonic applications, there has been renewed interest in studying the flow physics of supersonic modes using either theory or direct numerical simulations. This paper investigates the rise of supersonic modes in a high-enthalpy hypersonic flow over a 𝟓° half-angle cone at various wall temperatures using both quasiparallel linear stability theory and linear PSE. It was found that supersonic modes exist in all wall temperature conditions including the adiabatic wall case. Cooler wall temperature causes the second Mack mode to become an unstable supersonic mode naturally downstream of the upper-branch neutral location when the wall is sufficiently cooled. In terms of the integrated growth, the second mode is still the dominant mode. Nonetheless, supersonic modes can cause additional series of relatively weaker growth than that of the second mode beyond the peak amplitude location. According to the present linear PSE results, contrary to what was speculated in the literature, supersonic-mode pressure disturbance structure radiated into the freestream is nonacoustic in nature and the formation of unstable supersonic modes is mainly associated with the synchronization of phase speed between the instability and acoustic waves in a nonp ayer, not due to nonlinear modal interaction as suggested in the literature.

Chang, Chau-Lyan↗

Boundary Layer Stability and Transition in a Chemically Reacting Martian Atmosphere using LASTRAC

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.

Heather L Kline↗

TPSAS-NF1676L-13013-DND

Description, capabilities and utilization of the NASA Langley Aerothermodynamics Laboratory (LAL) are presented. The LAL consists of four hypersonic, blow-down-to-vacuum wind tunnels that collectively provide a range of Mach number from 6 to 10 (18 simulation), unit Reynolds number from 0.05 to 8 million per foot and, most importantly for blunt configurations, normal shock density ratio from 4 to 12. These wide ranges of hypersonic simulation parameters are due, in part, to the use of two different test gases (air and tetrafluoromethane, CF4), thereby making the facilities unique. The complex represents half of the conventional (as opposed to impulse) -type hypersonic wind tunnels operational in this country. LAL facilities are used to assess and optimize the hypersonic aerodynamic performance and aeroheating characteristics of aerospace vehicle concepts and to provide benchmark aerodynamic/aeroheating data for generating the flight aerodynamic databook and final design of the thermal protection system (TPS). Modifications and enhancements of LAL hardware components and instrumentation have been pursued to increase capability, reliability, and productivity in support of programmatic goals. The facilities are utilized to generate essentially all of the experimental hypersonic aerodynamic and aeroheating information for high-priority, fast-paced Agency programs and were a major contributor to the Shuttle Columbia Accident Investigation, space shuttle Return-to-Flight program and Hypersonic Thermodynamic Infrared Measurements (HYTHIRM) project. The LAL has also contributed to the development of the Orion Crew Exploration Vehicle, Ares launch system, Mars Science Laboratory, Genesis Sample Return Capsule and supported DOD programs including HIFiRE, X-37, Falcon HTV-2 and X-51A.

John R Micol↗

Analytical and Experimental Demonstration of an Alternate Mixing Performance Metric for High-Speed Fuel Mixing Studies

To experimentally assess the fuel/air mixing performance of high-speed fuel injectors, one-dimensional metrics that quantify the degree of mixing completeness downstream of the fuel injection location are required. The most accurate assessment of mixing performance is achieved with the mixing efficiency parameter. In order to experimentally determine the mixing efficiency parameter, the spatial distributions of both mass flux and fuel mass fraction must be measured. In-stream gas sampling techniques are commonly used to measure the fuel mass fraction distribution; however, the mass flux distribution is not easily determined because it requires the measurement of three independent aerothermodynamic variables in addition to the gas composition. Therefore, to experimentally determine the mixing efficiency parameter, the spatial distributions of four independent properties must be measured, with each property generally requiring its own unique probe. Because of this difficulty, it is commonly assumed that alternate metrics, which rely solely on the fuel distribution, are good indicators of mixing performance. However, since these alternate metrics do not provide a mass flux-weighted measure of mixing completeness, they can lead to incorrect conclusions being drawn about the mixing performance of the studied fuel injector configuration. Recognizing this shortcoming, this work proposes two new alternate mixing performance metrics that are easier to obtain than the mixing efficiency parameter. The analytical development of the new metrics, as well as their application to relevant CFD and experimental data of high-speed fuel injector configurations, is presented in this work. For two different experiments, the new metrics are shown to provide an excellent representation of the true mass flux-weighted mixing performance, unlike the traditionally-used alternatives. The results presented herein suggest that the new metrics can serve as accurate surrogates for mixing efficiency in future high-speed fuel/air mixing studies.

Cody R Ground↗

Recent Improvements to the LAURA and HARA Codes

This paper describes recent improvements to the LAURA and HARA codes. LAURA is a CFD code for aerothermodynamics, and HARA evaluates the shock-layer radiation that provides the radiative source term for the flowfield energy equations and radiative heating to a surface. The next release of LAURA and HARA includes a variety of new capabilities. These new capabilities include an automated uncertainty quantification workflow for radiative heat transfer, options for specifying surface roughness and turbulent transition location in the algebraic turbulence models, and improved grid and solution interpolation techniques. Additionally, the computational efficiency of both LAURA and HARA have been improved. Optimization of the MPI communication routines in LAURA are shown to improve the parallel efficiency of the primary flow when running with multiple processes per block, and recent optimization of HARA leverage graphics processing unit (GPU) acceleration in the radiation calculations. Using GPU acceleration of HARA is shown to decrease the cost of the radiation line-of-sight calculation by approximately one order of magnitude for a 10.5 km/s Earth entry simulation.

LAURA HARA CFD 5.6↗

Pterodactyl: Aerodynamic and Aerothermal Modeling for a Symmetric Deployable Earth Entry Vehicle with Flaps

NASA’s Pterodactyl project has investigated a deployable atmospheric entry vehicle integrated with a flap control system that provides precision targeting during reentry. The control system consists of eight flaps mounted at the edge of a heatshield that can deflect in and out of the flow. An aerodynamic and aerothermodynamic analysis process was developed for preliminary vehicle design and aerodynamic database generation using engineering and CFD tools with varying levels of fidelity. The objective of this analysis was to: 1) understand the inherent aerodynamics, 2) provide an aerodynamics database for stability and control analysis, and 3) provide qualitative and quantitative aeroheating analysis for Thermal Protection System modeling of the flaps. A high fidelity Euler code, Cart3D, was used to resolve complex flow features such as secondary shocks and shock impingement. A wide range of supersonic (Mach 2) and hypersonic Mach numbers (up to Mach 40) were tested with an Earth atmosphere model. A process was developed to utilize the adaptive volume mesh generator utility of the perfect gas model of Cart3D to create meshes with high cell efficiency and numerical stability prior to running a 2nd order accurate solution with a real gas model. High-fidelity aerothermal CFD simulations were performed using US3D to further improve the aerothermal analysis on hypersonic flow around a complex entry vehicle shape including viscosity, chemical reactions of air species, vibrational energy, and catalytic surface reactions. In this process, it was found that the flap control system provided multi-axis control that can be utilized for entry precision targeting. Additionally, the Pterodactyl vehicle can achieve up to a trim L/D of 0.2. Finally, the increased fidelity of the aerothermal heating environments revealed that the shear stress contributes to increasing the heating on the flaps.

Deployable Entry Vehicle↗

Pterodactyl: Aerodynamic and Aerothermal Modeling for a Symmetric Deployable Earth Entry Vehicle with Flaps

NASA’s Pterodactyl project has investigated a deployable atmospheric entry vehicle integrated with a flap control system that provides precision targeting during reentry. The control system consists of eight flaps mounted at the edge of a heatshield that can deflect in and out of the flow. An aerodynamic and aerothermodynamic analysis process was developed for preliminary vehicle design and aerodynamic database generation using engineering and CFD tools with varying levels of fidelity. The objective of this analysis was to: 1) understand the inherent aerodynamics, 2) provide an aerodynamics database for stability and control analysis, and 3) provide qualitative and quantitative aeroheating analysis for Thermal Protection System modeling of the flaps. A high fidelity Euler code, Cart3D, was used to resolve complex flow features such as secondary shocks and shock impingement. A wide range of supersonic (Mach 2) and hypersonic Mach numbers (up to Mach 40) were tested with an Earth atmosphere model. A process was developed to utilize the adaptive volume mesh generator utility of the perfect gas model of Cart3D to create meshes with high cell efficiency and numerical stability prior to running a 2nd order accurate solution with a real gas model. High-fidelity aerothermal CFD simulations were performed using US3D to further improve the aerothermal analysis on hypersonic flow around a complex entry vehicle shape including viscosity, chemical reactions of air species, vibrational energy, and catalytic surface reactions. In this process, it was found that the flap control system provided multi-axis control that can be utilized for entry precision targeting. Additionally, the Pterodactyl vehicle can achieve up to a trim L/D of 0.2. Finally, the increased fidelity of the aerothermal heating environments revealed that the shear stress contributes to increasing the heating on the flaps.

Deployable Entry Vehicle↗

Overview of Dragonfly Entry Aerosciences Measurements (DrEAM)

NASA Ames Research Center (ARC) leads the Dragonfly Entry Aerosciences Measurements (DrEAM) project, which is an aeroshell instrumentation suite on the Dragonfly mission that fulfills the Engineering Science Investigation requirement for the New Frontiers mission. NASA ARC is partnering with NASA Langley Research Center (LaRC) and the German Aerospace Center (DLR) to provide a comprehensive sensor suite, including a DLR-provided Data Acquisition System (DAS). DrEAM will provide key aerothermodynamic data and performance analysis for Dragonfly’s forebody and backshell Thermal Protection System (TPS). Titan’s atmosphere predominantly consists of ni-trogen (~98% by mole) with small amounts of me-thane (~2% by mole) and other trace gases. CN is a strong radiator, and is found in nonequilibrium con-centrations for Titan entry. The accurate modeling of nonequilibrium CN radiation has proven to be a diffi-cult task. Prompted by the Huygens mission, many experimental campaigns and analyses were performed to better understand the aerothermal environments experienced by the probe during Titan entry [1]. How-ever, the Huygens probe carried no heatshield instru-mentation. Therefore, the DrEAM sensor suite will sig-nificantly advance the state-of-the-art not only by documenting the environment and performance of Dragonfly’s entry system but also by making key in situ measurements in Titan’s atmosphere for the first time. Aerothermal environments and TPS response will be measured using sensors whose flight heritage is tak-en from the Mars Entry, Descent, and Landing In-strumentation 2 (MEDLI2) thermocouple plugs and the COMbined Aerothermal and Radiometer Sensor (COMARS) suite [2], with the latter supplied by DLR. The MEDLI2 project used embedded thermocouples to directly measure the in-depth TPS temperature-time history at several locations on the heat shield and backshell of the Mars 2020 entry vehicle. These tem-perature measurements, in turn, can be used to infer surface environments via an inverse analysis proce-dure analogous to that used for MEDLI. For DrEAM, the thermocouple plug subsystem will be known as Dragonfly Sensors for Aero-Thermal Reconstruction (DragSTR). On Schiaparelli, the COMARS suite in-cluded three total surface-mounted heat flux sensors, three pressure sensors, and one radiometer. For DrEAM, the COMARS package will be known as the COmbined Sensor System for Titan Atmosphere (COSSTA). Since the methane concentration in the Titan atmos-phere is directly proportional to the radiative heat flux, the COSSTA measurements will be used to reduce the current uncertainty in the methane volume fraction. Atmospheric density measurements and capsule aero-dynamic data will be obtained through the onboard Inertial Measurement Unit (IMU), supplemented by pressure transducers similar to those used by the MEDLI and MEDLI2 projects. The DrEAM pressure sensors will be known as the Dragonfly Atmospheric Flight Transducers. (DrAFT). The pressure measure-ments, when combined with data from the on-board IMU, will allow for reconstruction of such quantities as vehicle Mach number, freestream density, and atmos-pheric winds. DrAFT measurements will enhance Dragonfly trajectory reconstruction and enable a sepa-ration of the aerodynamics from the atmosphere, as was done for MEDLI [3] and is currently in process for the MEDLI2 flight data set.

EDL↗