Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational pyrolysis”

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.

85 records · Page 5

Effects of Non-Equilibrium Chemistry and Darcy-Forchheimer Flow of Pyrolysis Gas for a Charring Ablator

The Fully Implicit Ablation and Thermal Response code, FIAT, simulates pyrolysis and ablation of thermal protection materials and systems. The governing equations, which include energy conservation, a three-component decomposition model, and a surface energy balance, are solved with a moving grid. This work describes new modeling capabilities that are added to a special version of FIAT. These capabilities include a time-dependent pyrolysis gas flow momentum equation with Darcy-Forchheimer terms and pyrolysis gas species conservation equations with finite-rate homogeneous chemical reactions. The total energy conservation equation is also enhanced for consistency with these new additions. Parametric studies are performed using this enhanced version of FIAT. Two groups of analyses of Phenolic Impregnated Carbon Ablator (PICA) are presented. In the first group, an Orion flight environment for a proposed Lunar-return trajectory is considered. In the second group, various test conditions for arcjet models are examined. The central focus of these parametric studies is to understand the effect of pyrolysis gas momentum transfer on PICA material in-depth thermal responses with finite-rate, equilibrium, or frozen homogeneous gas chemistry. Results are presented, discussed, and compared with those predicted by the baseline PICA/FIAT ablation and thermal response model developed by the Orion Thermal Protection System Advanced Development Project.

Ablation↗

Simplified jet-A kinetic mechanism for combustor application

Successful modeling of combustion and emissions in gas turbine engine combustors requires an adequate description of the reaction mechanism. For hydrocarbon oxidation, detailed mechanisms are only available for the simplest types of hydrocarbons such as methane, ethane, acetylene, and propane. These detailed mechanisms contain a large number of chemical species participating simultaneously in many elementary kinetic steps. Current computational fluid dynamic (CFD) models must include fuel vaporization, fuel-air mixing, chemical reactions, and complicated boundary geometries. To simulate these conditions a very sophisticated computer model is required, which requires large computer memory capacity and long run times. Therefore, gas turbine combustion modeling has frequently been simplified by using global reaction mechanisms, which can predict only the quantities of interest: heat release rates, flame temperature, and emissions. Jet fuels are wide-boiling-range hydrocarbons with ranges extending through those of gasoline and kerosene. These fuels are chemically complex, often containing more than 300 components. Jet fuel typically can be characterized as containing 70 vol pct paraffin compounds and 25 vol pct aromatic compounds. A five-step Jet-A fuel mechanism which involves pyrolysis and subsequent oxidation of paraffin and aromatic compounds is presented here. This mechanism is verified by comparing with Jet-A fuel ignition delay time experimental data, and species concentrations obtained from flametube experiments. This five-step mechanism appears to be better than the current one- and two-step mechanisms.

Lee, Chi-Ming↗

TPSAS-NF1676L-21180-DND

Ablative materials are used as the Thermal Protection System (TPS) for space vehicles that enter a planetary atmosphere. During the entry phase, the vehicle experiences high rates of aerodynamic heating. These ablative materials, acting as a heat shield, protect the vehicle from the harsh environment and ensure a safe entry. Ablative materials reduce the incoming heat flux through various phenomena such as near-surface oxidation, pyrolysis chemical reaction, vaporization, and other erosive processes, collectively known as ablation. Spallation is one such ablative process, in which solid particles are ejected from the material into the flow field. The many conditions leading to spallation, however, are not well understood and cannot be described accurately by simple phenomenological models. However, the spallation phenomenon might impact the aerodynamic heating rates of re-entry vehicles using ablative heat shields. For example, the solid particles ejected from the surface can modify the near surface chemical processes and flow field. To investigate spallation effects, a code was developed to compute the dynamics of spalled particles. The code uses a finite-rate chemistry model to study the chemical interactions of the particles with the flow field. The spallation code is dynamically coupled to a CFD solver that models the time accurate hypersonic flow field around an ablative sample. An experimental investigation was also performed at the NASA Langley HYMETS facility. HYMETS is a 400 kW arcjet that can simulate hypersonic flight and Earth or Martian entry conditions. Using high speed cameras, an IR camera, NIR and UV-VIS spectrometers, and other diagnostic instruments, test data on various graphite, LI-900, FiberForm, and PICA test samples were obtained, using target heat flux conditions of 100, 200, and 400 W/cm 2 . From the test campaign, measurements of particle velocity, acceleration and trajectory are obtained.

A Martin↗

Formation, Structure and Properties of Amorphous Carbon Char from Polymer Materials in Extreme Atmospheric Reentry Environments

Amorphous carbonaceous char produced from the pyrolysis of polymer solids has many desirable properties for ablative heat shields for space vehicles. Molecular dynamics simulations are presented to study the transformation of the local atomic structure from virgin polymer to a dense, disordered char [1]. Release of polymer hydrogen is found to be critical to allow the system to collapse into a highly coordinated char structure. Mechanisms of the char formation process and the morphology of the resulting structures are elucidated. Thermal conductivity and mechanical response of the resulting char are evaluated [2]. During reenty, the optical response and oxidative reactivity of char are also important properties. Results of ab initio computations of char optical functions [3] and char reactivity [4] are also presented.

Lawson, John W.↗

Thermal Response Modeling System for a Mars Sample Return Vehicle

A multi-dimensional, coupled thermal response modeling system for analysis of hypersonic entry vehicles is presented. The system consists of a high fidelity Navier-Stokes equation solver (GIANTS), a two-dimensional implicit thermal response, pyrolysis and ablation program (TITAN), and a commercial finite-element thermal and mechanical analysis code (MARC). The simulations performed by this integrated system include hypersonic flowfield, fluid and solid interaction, ablation, shape change, pyrolysis gas eneration and flow, and thermal response of heatshield and structure. The thermal response of the heatshield is simulated using TITAN, and that of the underlying structural is simulated using MARC. The ablating heatshield is treated as an outer boundary condition of the structure, and continuity conditions of temperature and heat flux are imposed at the interface between TITAN and MARC. Aerothermal environments with fluid and solid interaction are predicted by coupling TITAN and GIANTS through surface energy balance equations. With this integrated system, the aerothermal environments for an entry vehicle and the thermal response of the entire vehicle can be obtained simultaneously. Representative computations for a flat-faced arc-jet test model and a proposed Mars sample return capsule are presented and discussed.

Chen, Y.-K.↗

Thermal Response Modeling System for a Mars Sample Return Vehicle

A multi-dimensional, coupled thermal response modeling system for analysis of hypersonic entry vehicles is presented. The system consists of a high fidelity Navier-Stokes equation solver (GIANTS), a two-dimensional implicit thermal response, pyrolysis and ablation program (TITAN), and a commercial finite element thermal and mechanical analysis code (MARC). The simulations performed by this integrated system include hypersonic flowfield, fluid and solid interaction, ablation, shape change, pyrolysis gas generation and flow, and thermal response of heatshield and structure. The thermal response of the heatshield is simulated using TITAN, and that of the underlying structural is simulated using MARC. The ablating heatshield is treated as an outer boundary condition of the structure, and continuity conditions of temperature and heat flux are imposed at the interface between TITAN and MARC. Aerothermal environments with fluid and solid interaction are predicted by coupling TITAN and GIANTS through surface energy balance equations. With this integrated system, the aerothermal environments for an entry vehicle and the thermal response of the entire vehicle can be obtained simultaneously. Representative computations for a flat-faced arc-jet test model and a proposed Mars sample return capsule are presented and discussed.

Chen, Y.-K.↗

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.↗

An Approximate Ablative Thermal Protection System Sizing Tool for Entry System Design

A computer tool to perform entry vehicle ablative thermal protection systems sizing has been developed. Two options for calculating the thermal response are incorporated into the tool. One, an industry-standard, high-fidelity ablation and thermal response program was integrated into the tool, making use of simulated trajectory data to calculate its boundary conditions at the ablating surface. Second, an approximate method that uses heat of ablation data to estimate heat shield recession during entry has been coupled to a one-dimensional finite-difference calculation that calculates the in-depth thermal response. The in-depth solution accounts for material decomposition, but does not account for pyrolysis gas energy absorption through the material. Engineering correlations are used to estimate stagnation point convective and radiative heating as a function of time. The sizing tool calculates recovery enthalpy, wall enthalpy, surface pressure, and heat transfer coefficient. Verification of this tool is performed by comparison to past thermal protection system sizings for the Mars Pathfinder and Stardust entry systems and calculations are performed for an Apollo capsule entering the atmosphere at lunar and Mars return speeds.

Dec, John A.↗

An Approximate Ablative Thermal Protection System Sizing Tool for Entry System Design

A computer tool to perform entry vehicle ablative thermal protection systems sizing has been developed. Two options for calculating the thermal response are incorporated into the tool. One, an industry-standard, high-fidelity ablation and thermal response program was integrated into the tool, making use of simulated trajectory data to calculate its boundary conditions at the ablating surface. Second, an approximate method that uses heat of ablation data to estimate heat shield recession during entry has been coupled to a one-dimensional finite-difference calculation that calculates the in-depth thermal response. The in-depth solution accounts for material decomposition, but does not account for pyrolysis gas energy absorption through the material. Engineering correlations are used to estimate stagnation point convective and radiative heating as a function of time. The sizing tool calculates recovery enthalpy, wall enthalpy, surface pressure, and heat transfer coefficient. Verification of this tool is performed by comparison to past thermal protection system sizings for the Mars Pathfinder and Stardust entry systems and calculations are performed for an Apollo capsule entering the atmosphere at lunar and Mars return speeds.

Dec, John A.↗

Integrated Thermal Response Tool for Earth Entry Vehicles

A system is presented for multi-dimensional, fully-coupled thermal response modeling of hypersonic entry vehicles. The system consists of a two-dimensional implicit thermal response, pyrolysis and ablation program (TITAN), a commercial finite-element thermal and mechanical analysis code (MARC), and a high fidelity Navier-Stokes equation solver (GIANTS). The simulations performed by this integrated system include hypersonic flow-field, fluid and solid interaction, ablation, shape change, pyrolysis gas generation and flow, and thermal response of heatshield and structure. The thermal response of the ablating and charring heatshield material is simulated using TITAN, and that of the underlying structural is simulated using MARC. The ablating heatshield is treated as an outer boundary condition of the structure, and continuity conditions of temperature and heat flux are imposed at the interface between TITAN and MARC. Aerothermal environments with fluid and solid interaction are predicted by coupling TITAN and GIANTS through surface energy balance equations. With this integrated system, the aerothermal environments for an entry vehicle and the thermal response of both the heatshield and the structure can be obtained simultaneously. Representative computations for a proposed blunt body earth entry vehicle are presented and discussed in detail.

Chen, Y.-K.↗

Optimization of thermal protection systems for the space shuttle vehicle. Volume 1: Final report

A study performed to continue development of computational techniques for the Space Shuttle Thermal Protection System is reported. The resulting computer code was used to perform some additional optimization studies on several TPS configurations. The program was developed in Fortran 4 for the CDC 6400, and it was converted to Fortran 5 to be used for the Univac 1108. The computational methodology is developed in modular fashion to facilitate changes and updating of the techniques and to allow overlaying the computer code to fit into approximately 131,000 octal words of core storage. The program logic involves subroutines which handle input and output of information between computer and user, thermodynamic stress, dynamic, and weight/estimate analyses of a variety of panel configurations. These include metallic, ablative, RSI (with and without an underlying phase change material), and a thermodynamic analysis only of carbon-carbon systems applied to the leading edge and flat cover panels. Two different thermodynamic analyses are used. The first is a two-dimensional, explicit precedure with variable time steps which is used to describe the behavior of metallic and carbon-carbon leading edges. The second is a one-dimensional implicity technique used to predict temperature in the charring ablator and the noncharring RSI. The latter analysis is performed simply by suppressing the chemical reactions and pyrolysis of the TPS material.

Source record↗

Thermodynamic Analysis of Boundary Layer Species During Ablation

The focus of this work is on carbon clusters, hydrocarbons and other carbonaceous species which are constituents of pyrolysis gases injected into the boundary layer. Some of these molecules are observed to have absorption spectra in the VUV and UV region that match the emission spectra of atomic nitrogen and oxygen. Hence, they can potentially absorb the radiation impinging on the heat shield. Accurate thermodynamic data is not available for many of these species over the relevant temperature range for reentry. The objective of this work is to determine the thermochemical properties of potential radiation absorbing molecules using accurate ab initio quantum chemistry calculations. We have compiled a database of formation enthalpies, anharmonic vibrational frequencies and rotational constants for the ground and low-lying excited electronic states of these carbonaceous species. Using these thermochemical properties, we compute equilibrium mole fractions and observe that the mole fractions of some species of interest are significant, thus they could affect the radiative heat flux. Furthermore, comparisons of the formation enthalpies and mole fractions are made with data from literature.

heatshields↗

Real-Time X-ray Radiography Diagnostics of Components in Solid Rocket Motors

Solid rocket motors (SRMs) typically use nozzle materials which are required to maintain their shape as well as insulate the underlying support structure during the motor operation. In addition, SRMs need internal insulation materials to protect the motor case from the harsh environment resulting from the combustion of solid propellant. In the nozzle, typical materials consist of high density graphite, carbon-carbon composites and carbon phenolic composites. Internal insulation of the motor cases is typically a composite material with carbon, asbestos, Kevlar, or silica fibers in an ablative matrix such as EPDM or NBR. For both nozzle and internal insulation materials, the charring process occurs when the hot combustion products heat the material intensely. The pyrolysis of the matrix material takes away a portion of the thermal energy near the wall surface and leaves behind a char layer. The fiber reinforcement retains the porous char layer which provides continued thermal protection from the hot combustion products. It is of great interest to characterize both the total erosion rates of the material and the char layer thickness. By better understanding of the erosion process for a particular ablative material in a specific flow environment, the required insulation material thickness can be properly selected. The recession rates of internal insulation and nozzle materials of SRMs are typically determined by testing in some sort of simulated environment; either arc-jet testing, flame torch testing, or subscale SRMs of different size. Material recession rates are deduced by comparison of pre- and post-test measurements and then averaging over the duration of the test. However, these averaging techniques cannot be used to determine the instantaneous recession rates of the material. Knowledge of the variation in recession rates in response to the instantaneous flow conditions during the motor operation is of great importance. For example, in many SRM configurations the recession of the solid propellant grain can drastically alter the flow-field and effect the recession of internal insulation and nozzle materials. Simultaneous measurement of the overall erosion rate, the development of the char layer, and the recession of the char-virgin interface during the motor operation can be rather difficult. While invasive techniques have been used with limited success, they have serious drawbacks. Break wires or make wire sensors can be installed into a sufficient number of locations in the charring material from which a time history of the charring surface can be deduced. These sensors fundamentally alter the local structure of the material in which they are imbedded. Also, the location of these sensors within the material is not known precisely without the use of an X-ray. To determine instantaneous recession rates, real-time X-ray radiography (X-ray RTR) has been utilized in several SRM experiments at PSU. The X-ray RTR system discussed in this paper consists of an X-ray source, X-ray image intensifier, and CCD camera connected to a capture computer. The system has been used to examine the ablation process of internal insulation as well as nozzle material erosion in a subscale SRM. The X-ray source is rated to 320 kV at 10 mA and has both a large (5.5 mm) and small (3.0 mm) focal spot. The lead-lined cesium iodide X-ray image intensifier produces an image which is captured by a CCD camera with a 1,000 x 1,000 pixel resolution. To produce accurate imagery of the object of interest, the alignment of the X-ray source to the X-ray image intensifier is crucial. The image sequences captured during the operation of an SRM are then processed to enhance the quality of the images. This procedure allows for computer software to extract data on the total erosion rate and the char layer thickness. Figure 1 Error! Reference source not found.shows a sequence of images captured during the operation the subscale SRM with the X-ray RTR system. The X-rayTR system, alignment procedure, uncertainty determination, and image analysis process will be discussed in detail in the full manuscript.

Cortopassi, A. C.↗