Engineering PapersSearch

Engineering topics

Chen, Yih-Kanq

Publications and source records attributed to Chen, Yih-Kanq.

At least 19 records

Arcjet Ablation of Stony and Iron Meteorites

A test campaign was conducted placing meteorites in the 60 MW plasma Arcjet Interaction Heating Facility at NASA Ames Research Center, with the aim to achieve flight-relevant conditions for asteroid impacts in Earth's atmosphere and to provide insight into how meteoritic materials respond to extreme entry heating environments. The test conditions at heat flux of 4000 W/m2 and 140 kPa stagnation pressure are comparable to those experienced by a 30-meter diameter asteroid moving at 20 km/s velocity at 65 km altitude in the Earth's atmosphere. Test objects were a stony type H5 ordinary chondrite (Tamdakht) and an iron type IAB-MG meteorite (Campo Del Cielo), and included the terrestrial analogs Dense Flood Basalt and Fused Silica. All samples were exposed for only a few seconds in the plasma stream. Significant melt flow and vaporization was observed for both the stony and iron meteorites during exposure. Mass loss from spallation of fragments was also observed. Vapor emitted atomic lines from alkali metals and iron, but did not emit the expected MgO molecular band emissions. The meteoritic melts flowed more rapidly, indicating lower viscosity, than those of Fused Silica. The surface recession was mapped. The effective heat of ablation derived from this showed that ablation under these conditions occurred in the melt-dominated regime. Ablation parameters have an effect on ground damage estimates. A bias in ablation parameters towards the melt-dominated regime would imply that impacting asteroids survive to lower altitude, and therefore could possibly have airbursts with a larger ground damage footprint.

Arcjet

Arcjet Ablation of Stony and Iron Meteorites

A test campaign was conducted placing meteorites in the 60 MW plasma Arcjet Interaction Heating Facility at NASA Ames Research Center, with the aim to achieve flight-relevant conditions for asteroid impacts in Earth's atmosphere and to provide insight into how meteoritic materials respond to extreme entry heating environments. The test conditions at heat flux of 4000 W/m2 and 140 kPa stagnation pressure are comparable to those experienced by a 30-meter diameter asteroid moving at 20 km/s velocity at 65 km altitude in the Earth's atmosphere. Test objects were a stony type H5 ordinary chondrite (Tamdakht) and an iron type IAB-MG meteorite (Campo Del Cielo), and included the terrestrial analogs Dense Flood Basalt and Fused Silica. All samples were exposed for only a few seconds in the plasma stream. Significant melt flow and vaporization was observed for both the stony and iron meteorites during exposure. Mass loss from spallation of fragments was also observed. Vapor emitted atomic lines from alkali metals and iron, but did not emit the expected MgO molecular band emissions. The meteoritic melts flowed more rapidly, indicating lower viscosity, than those of Fused Silica. The surface recession was mapped. The effective heat of ablation derived from this showed that ablation under these conditions occurred in the melt-dominated regime. Ablation parameters have an effect on ground damage estimates. A bias in ablation parameters towards the melt-dominated regime would imply that impacting asteroids survive to lower altitude, and therefore could possibly have airbursts with a larger ground damage footprint.

Asteroid

Ablation and Heating During Atmospheric Entry and Its Effect on Airburst Risk

Large meteoroids and asteroids entering the atmosphere endure tremendous heating from the shock heated air, and thereby lose a significant fraction of their mass during atmospheric entry a process known as ablation. The predicted evolution of the asteroids mass as it passes through the atmosphere can affect both the predicted energy deposition profile relevant to an airburst event, or the residual mass that strikes the ground in the case of an impact event. This presentation is divided roughly into two parts. In the first part, an overview of traditional models for heat transfer and ablation that are historically used in the meteor physics community is presented, and the validity in the asteroid entry regime discussed. Sensitivity analyses performed using the recently developed Fragment-Cloud Model (FCM) will be presented which show illustrate the range of sizes and entry parameters for which the predicted asteroid threat is most sensitive to the models for ablation and heat transfer. The second part of the presentation shall focus on recent work done under NASAs Asteroid Threat Assessment Project (ATAP) to develop new models for heat transfer and ablation using high-fidelity numerical simulation in concert with state-of-the-art experiments. Coupled computational fluid dynamics (CFD)radiation transport simulations preformed using the state-of-the-art entry modeling tools at NASA show that, for large meteoroids and asteroids, there can significant attenuation of the heat transfer to the surface (95 in some cases) by the products of ablation. In addition to the heat transfer, new models for the material response and ablation of asteroidal material have been developed [cite]. In the current work, we present finding from recent novel experiments performed in the arc jet facility at NASA Ames, which allows us to, in part, simulate the extreme environment experienced by the asteroid during entry. Briefly, the experimental set-up was comprised of a 1.5 conical article of machined H5 chondrite, which was exposed to a high-enthalpy flow resulting in approximately 4 kWcm2 of heating to the surface. A still frame capture from high-speed video taken during this experiment can be seen in Figure 1. In this figure, we can observe some of the major mechanisms for meteoroid ablation, such as melt flow, spallation (mechanical removal of material), and vaporization. Major findings from this, and other experiments will be discussed, as well progress on utilizing the data from the experiments to inform and develop improved models for ablation.

Asteroids

Thermal Ablation Modeling for Silicate Materials

A general thermal ablation model for silicates is proposed. The model includes the mass losses through the balance between evaporation and condensation, and through the moving molten layer driven by surface shear force and pressure gradient. This model can be applied in the ablation simulation of the meteoroid and the glassy ablator for spacecraft Thermal Protection Systems. Time-dependent axisymmetric computations are performed by coupling the fluid dynamics code, Data-Parallel Line Relaxation program, with the material response code, Two-dimensional Implicit Thermal Ablation simulation program, to predict the mass lost rates and shape change. The predicted mass loss rates will be compared with available data for model validation, and parametric studies will also be performed for meteoroid earth entry conditions.

Ablation

Thermal Ablation Modeling for Silicate Materials

A thermal ablation model for silicates is proposed. The model includes the mass losses through the balance between evaporation and condensation, and through the moving molten layer driven by surface shear force and pressure gradient. This model can be applied in ablation simulations of the meteoroid or glassy Thermal Protection Systems for spacecraft. Time-dependent axi-symmetric computations are performed by coupling the fluid dynamics code, Data-Parallel Line Relaxation program, with the material response code, Two-dimensional Implicit Thermal Ablation simulation program, to predict the mass lost rates and shape change. For model validation, the surface recession of fused amorphous quartz rod is computed, and the recession predictions reasonably agree with available data. The present parametric studies for two groups of meteoroid earth entry conditions indicate that the mass loss through moving molten layer is negligibly small for heat-flux conditions at around 1 MW/cm(exp. 2).

Ablators

Radiant Heating Tests on Meteoritic Material

A new effort geared toward modeling the physics of meteor entry and break-up is underway at NASA Ames Research Center. This is part of a broader interdisciplinary effort on providing physics-based risk assessment models for potentially hazardous objects. As part of the entry modeling task we are seeking to improve our understanding of, among other things, the ablation of meteoric material during high speed entry into earths atmosphere. Meteoroid entry differs greatly in some key respects from spacecraft entry modeling. First, the aerothermal environment at these high velocities (18 km/s) is dominated by radiation. Second, meteoroids less than, say, 50m in size, will likely lose a significant portion of their mass during the high-speed atmospheric entry process due to vaporization as well as melting and spallation. The mass of the object, in turn, directly affects the amount of energy deposited in the atmosphere, and therefore the amount of damage done. Thus it is important for us to understand and be able to model this process in greater detail in order to assess the hazard posed by these objects.In this presentation, we first give an overview of the simple ablation models that are typically used in meteor entry calculations. Additionally, we will present a new model which utilizes a similar approach to what is typically done for spacecraft TPS response modeling. This uses an equilibrium assumption near the surface to compute the ablation rate, as is done in heritage material response codes. Next we describe the radiant heating experiment which uses a high-powered laser to emulate the radiation dominated heating environment experienced by meteoroids during atmospheric entry. The facility the Laser Hardened Material Evaluation Laboratory (LHMEL) permits us to expose samples of meteoritic material to heating rates in excess of 100kW/sq.cm. An overview of the experimental set-up and test plan for the initial exploratory campaign at this facility will be given. Then we present both qualitative and quantitative results from this initial test series. Comparisons between predicted and measured ablation rates suggest that there may be significant blockage of the incident beam by the ablation plume. Furthermore, melt is shown to be a significant ablation mechanics, even at high heating rates (16 kW/sq.cm). Finally, comparisons between the phenomenology of the ablation of terrestrial rocks namely, basalt -- to that of meteorites show very different behavior. This is shown to likely be due in part to the effect of composition on the melt viscosity.

Stern, Eric

Evaluation of Finite-Rate Gas/Surface Interaction Models for a Carbon Based Ablator

Two sets of finite-rate gas-surface interaction model between air and the carbon surface are studied. The first set is an engineering model with one-way chemical reactions, and the second set is a more detailed model with two-way chemical reactions. These two proposed models intend to cover the carbon surface ablation conditions including the low temperature rate-controlled oxidation, the mid-temperature diffusion-controlled oxidation, and the high temperature sublimation. The prediction of carbon surface recession is achieved by coupling a material thermal response code and a Navier-Stokes flow code. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and Ablation Program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting full Navier-Stokes equations using Data Parallel Line Relaxation method. Recession analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities with heat fluxes ranging from 45 to 1100 wcm2 are performed and compared with data for model validation. The ablating material used in these arc-jet tests is Phenolic Impregnated Carbon Ablator. Additionally, computational predictions of surface recession and shape change are in good agreement with measurement for arc-jet conditions of Small Probe Reentry Investigation for Thermal Protection System Engineering.

Ablation

Physics-Based Modeling of Meteor Entry and Breakup

A new research effort at NASA Ames Research Center has been initiated in Planetary Defense, which integrates the disciplines of planetary science, atmospheric entry physics, and physics-based risk assessment. This paper describes work within the new program and is focused on meteor entry and breakup. Over the last six decades significant effort was expended in the US and in Europe to understand meteor entry including ablation, fragmentation and airburst (if any) for various types of meteors ranging from stony to iron spectral types. These efforts have produced primarily empirical mathematical models based on observations. Weaknesses of these models, apart from their empiricism, are reliance on idealized shapes (spheres, cylinders, etc.) and simplified models for thermal response of meteoritic materials to aerodynamic and radiative heating. Furthermore, the fragmentation and energy release of meteors (airburst) is poorly understood. On the other hand, flight of human-made atmospheric entry capsules is well understood. The capsules and their requisite heatshields are designed and margined to survive entry. However, the highest speed Earth entry for capsules is less than 13 km/s (Stardust). Furthermore, Earth entry capsules have never exceeded diameters of 5 m, nor have their peak aerothermal environments exceeded 0.3 atm and 1 kW/cm2. The aims of the current work are: (i) to define the aerothermal environments for objects with entry velocities from 13 to greater than 20 km/s; (ii) to explore various hypotheses of fragmentation and airburst of stony meteors in the near term; (iii) to explore the possibility of performing relevant ground-based tests to verify candidate hypotheses; and (iv) to quantify the energy released in airbursts. The results of the new simulations will be used to anchor said risk assessment analyses. With these aims in mind, state-of-the-art entry capsule design tools are being extended for meteor entries. We describe: (i) applications of current simulation tools to spherical geometries of diameters ranging from 1 to 100 m for an entry velocity of 20 km/s and stagnation pressures ranging from 1 to 100 atm; (ii) the influence of shape and departure of heating environment predictions from those for a simple spherical geometry; (iii) assessment of thermal response models for silica subject to intense radiation; and (iv) results for porosity-driven gross fragmentation of meteors, idealized as a collection of smaller objects. Lessons learned from these simulations will be used to help understand the Chelyabinsk meteor entry up to its first point of fragmentation.

Aerothermodynamics

Physics-Based Modeling of Meteor Entry and Breakup

A new research effort at NASA Ames Research Center has been initiated in Planetary Defense, which integrates the disciplines of planetary science, atmospheric entry physics, and physics-based risk assessment. This paper describes work within the new program and is focused on meteor entry and breakup. Over the last six decades significant effort was expended in the US and in Europe to understand meteor entry including ablation, fragmentation and airburst (if any) for various types of meteors ranging from stony to iron spectral types. These efforts have produced primarily empirical mathematical models based on observations. Weaknesses of these models, apart from their empiricism, are reliance on idealized shapes (spheres, cylinders, etc.) and simplified models for thermal response of meteoritic materials to aerodynamic and radiative heating. Furthermore, the fragmentation and energy release of meteors (airburst) is poorly understood. On the other hand, flight of human-made atmospheric entry capsules is well understood. The capsules and their requisite heat shields are designed and margined to survive entry. However, the highest speed Earth entry for capsules is 13 kms (Stardust). Furthermore, Earth entry capsules have never exceeded diameters of 5 m, nor have their peak aerothermal environments exceeded 0.3 atm and 1 kWcm2. The aims of the current work are: (i) to define the aerothermal environments for objects with entry velocities from 13 to 20 kms; (ii) to explore various hypotheses of fragmentation and airburst of stony meteors in the near term; (iii) to explore the possibility of performing relevant ground-based tests to verify candidate hypotheses; and (iv) to quantify the energy released in airbursts. The results of the new simulations will be used to anchor said risk assessment analyses.With these aims in mind, state-of-the-art entry capsule design tools are being extended for meteor entries. We describe: (i) applications of current simulation tools to spherical geometries of diameters ranging from 1 to 100 m for an entry velocity of 20 kms and stagnation pressures ranging from 1 to 100 atm; (ii) the influence of shape and departure of heating environment predictions from those for a simple spherical geometry; (iii) assessment of thermal response models for silica subject to intense radiation; and (iv) results for porosity-driven gross fragmentation of meteors, idealized as a collection of smaller objects. Lessons learned from these simulations will be used to help understand the Chelyabinsk meteor entry up to its first point of fragmentation.

Atmospheric Entry

Implicit Coupling Approach for Simulation of Charring Carbon Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver with nonequilibrium gas/surface interaction for simulation of charring carbon ablators can be performed using an implicit approach. The material thermal response code used in this study is the three-dimensional version of Fully Implicit Ablation and Thermal response program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation method. Coupling between the material response and flow codes is performed by solving the surface mass balance in flow solver and the surface energy balance in material response code. Thus, the material surface recession is predicted in flow code, and the surface temperature and pyrolysis gas injection rate are computed in material response code. It is demonstrated that the time-lagged explicit approach is sufficient for simulations at low surface heating conditions, in which the surface ablation rate is not a strong function of the surface temperature. At elevated surface heating conditions, the implicit approach has to be taken, because the carbon ablation rate becomes a stiff function of the surface temperature, and thus the explicit approach appears to be inappropriate resulting in severe numerical oscillations of predicted surface temperature. Implicit coupling for simulation of arc-jet models is performed, and the predictions are compared with measured data. Implicit coupling for trajectory based simulation of Stardust fore-body heat shield is also conducted. The predicted stagnation point total recession is compared with that predicted using the chemical equilibrium surface assumption

Ablation

Effects of Nonequilibrium Chemistry and Darcy-Forchheimer Pyrolysis Flow for Charring Ablator

The fully implicit ablation and thermal response code 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 code. 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. Two groups of parametric studies of the phenolic impregnated carbon ablator are performed. 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 material in-depth thermal responses with finite-rate, equilibrium, or frozen homogeneous gas chemistry. Results indicate that the presence of chemical nonequilibrium pyrolysis gas flow does not significantly alter the in-depth thermal response performance predicted using the chemical equilibrium gas model.

Chen, Yih-Kanq

Effect of Non-Equilibrium Surface Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using non-equilibrium gas/surface interaction model provides time-accurate solutions for the multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and AblatioN Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas mass conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between the material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of nonequilibrium gas/surface interaction chemistry between air and the carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kanq

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

Validation of a Three-Dimensional Ablation and Thermal Response Simulation Code

The 3dFIAT code simulates pyrolysis, ablation, and shape change of thermal protection materials and systems in three dimensions. The governing equations, which include energy conservation, a three-component decomposition model, and a surface energy balance, are solved with a moving grid system to simulate the shape change due to surface recession. This work is the first part of a code validation study for new capabilities that were added to 3dFIAT. These expanded capabilities include a multi-block moving grid system and an orthotropic thermal conductivity model. This paper focuses on conditions with minimal shape change in which the fluid/solid coupling is not necessary. Two groups of test cases of 3dFIAT analyses of Phenolic Impregnated Carbon Ablator in an arc-jet are presented. In the first group, axisymmetric iso-q shaped models are studied to check the accuracy of three-dimensional multi-block grid system. In the second group, similar models with various through-the-thickness conductivity directions are examined. In this group, the material thermal response is three-dimensional, because of the carbon fiber orientation. Predictions from 3dFIAT are presented and compared with arcjet test data. The 3dFIAT predictions agree very well with thermocouple data for both groups of test cases.

Chen, Yih-Kanq

Computational Analysis of Arc-Jet Wedge Tests Including Ablation and Shape Change

Coupled fluid-material response analyses of arc-jet wedge ablation tests conducted in a NASA Ames arc-jet facility are considered. These tests were conducted using blunt wedge models placed in a free jet downstream of the 6-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. The fluid analysis includes computational Navier-Stokes simulations of the nonequilibrium flowfield in the facility nozzle and test box as well as the flowfield over the models. The material response analysis includes simulation of two-dimensional surface ablation and internal heat conduction, thermal decomposition, and pyrolysis gas flow. For ablating test articles undergoing shape change, the material response and fluid analyses are coupled in order to calculate the time dependent surface heating and pressure distributions that result from shape change. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator. Effects of the test article shape change on fluid and material response simulations are demonstrated, and computational predictions of surface recession, shape change, and in-depth temperatures are compared with the experimental measurements.

Goekcen, Tahir

Ablation, Thermal Response, and Chemistry Program for Analysis of Thermal Protection Systems

In previous work, the authors documented the Multicomponent Ablation Thermochemistry (MAT) and Fully Implicit Ablation and Thermal response (FIAT) programs. In this work, key features from MAT and FIAT were combined to create the new Fully Implicit Ablation, Thermal response, and Chemistry (FIATC) program. FIATC is fully compatible with FIAT (version 2.5) but has expanded capabilities to compute the multispecies surface chemistry and ablation rate as part of the surface energy balance. This new methodology eliminates B' tables, provides blown species fractions as a function of time, and enables calculations that would otherwise be impractical (e.g. 4+ dimensional tables) such as pyrolysis and ablation with kinetic rates or unequal diffusion coefficients. Equations and solution procedures are presented, then representative calculations of equilibrium and finite-rate ablation in flight and ground-test environments are discussed.

Milos, Frank S.

Validation of PICA Ablation and Thermal-Response Model at Low Heat Flux

Phenolic Impregnated Carbon Ablator (PICA) was the forebody heatshield material on the Stardust sample-return capsule and is also a primary candidate material for the Mars Science Lander (MSL), the Orion Crew Module, and the SpaceX Dragon vehicle. As part of the heatshield qualification for Orion, physical and thermal properties of virgin and charred PICA were measured, and an ablation and thermal response model was developed. We validated the model by comparing it with recession and temperature data from stagnation arcjet tests conducted over a wide range of stagnation heat flux of 107 to 1102 W/sq cm. The effect of orthotropic thermal conductivity was evident in the thermal response of the arcjet models. In general, model predictions compared well with the data; however, the uncertainty of the recession prediction was greatest for heat fluxes below 200 W/sq cm. More recent MSL testing focused on the low heat flux regime of 45 to 250 W/sq cm. The new results confirm the recession uncertainty, especially for pressures below 6 kPa. In this work we focus on improving the model predictions for MSL and Orion tests below 250 W/sq cm.

Milos, Frank S.

Arc Jet Testing of Carbon Phenolic for Mars Sample Return and Future NASA Missions

The objective of the Mars Sample Return (MSR) Mission is to return a sample of MArtian soil to Earth. The Earth Entry Vehicle (EEV) brings te samples through the atmosphere to the ground.The program aims to: Model aerothermal environment during EEV flight; On the basis of results, select potential TPS materials for EEV forebody; Fabricate TPS materials; Test the materials in the arc jet environment representative of predicted flight environment;Evaluate material performance; Compare results of modeling predictions with test results.

Laub, Bernard