Engineering PapersSearch

Engineering topics

Palmer, Grant

Publications and source records attributed to Palmer, Grant.

At least 19 records

Modeling Heatshield Erosion due to Dust Particle Impacts for Martian Entries

Heatshield design for spacecraft entering the atmosphere of Mars may be affected by the presence of atmospheric dust. Particle impacts with sufficient kinetic energy can cause spallation damage to the heatshield that must be estimated. The dust environment in terms of particle size distribution and number density can be inferred from ground-based or atmospheric observations at Mars. Using a Lagrangian approach, the particle trajectories through the shock layer can be computed using a set of coupled ordinary differential equations. The dust particles are small enough that non-continuum effects must be accounted for when computing the drag coefficient and heat transfer to the particle surface. Surface damage correlations for impact crater diameter and penetration depth are presented for fused-silica, AVCOAT, Shuttle tiles, cork, and Norcoat Liège. The cork and Norcoat Liège correlations are new and were developed in this study. The modeling equations presented in this paper are applied to compute the heatshield erosion due to dust particle impacts on the ExoMars Schiaparelli entry capsule during dust storm conditions.

Palmer, Grant

Modeling Heatshield Erosion Due to Dust Particle Impacts for Martian Entries

Heatshield design for spacecraft entering the atmosphere of Mars may be affected by the presence of atmospheric dust. Particle impacts with sufficient kinetic energy can cause spallation damage to the heatshield that must be estimated. The dust environment in terms of particle size distribution and number density can be inferred from ground-based or atmospheric observations at Mars. Using a Lagrangian approach, the particle trajectories through the shock layer can be computed using a set of coupled ordinary differential equations. The dust particles are small enough that non-continuum effects must be accounted for when computing the drag coefficient and heat transfer to the particle surface. Surface damage correlations for impact crater diameter and penetration depth are presented for fused-silica, AVCOAT, Shuttle tiles, cork, and Norcoat Liège. The cork and Norcoat Liège correlations are new and were developed in this study. The modeling equations presented in this paper are applied to compute the heatshield erosion due to dust particle impacts on the ExoMars Schiaparelli entry capsule during dust storm conditions.

Palmer, Grant

Heatshield for Extreme Entry Environment Technology (HEEET) Thermal Protection System (TPS)

Starting in 2013 and completing in 2019, the Heatshield for Extreme Entry Environment Technology (HEEET) project has been working to mature a 3-D Woven Thermal Protection System (TPS) to Technical Readiness Level (TRL) 6 to support future NASA missions to destinations with extreme entry environments such as Venus, Saturn, Uranus, Neptune and high-speed sample return missions to Earth. A key aspect of the project has been the building and testing of a 1-meter base diameter Engineering Test Unit (ETU) representative of what could be used for a Saturn probe. This paper provides a high-level overview of the HEEET project including manufacturing and testing of the ETU for structural model verification, establish system capability and verify manufacturing workmanship.

Extreme Entry Environment

Heatshield Erosion Due to Dust Particle Impacts on the Schiaparelli Capsule During Martian Entry

One of the unique features of spacecraft entering the atmosphere of Mars is the possibility of a major dust storm occurring during the entry. Design of the thermal protection system (TPS) for Mars missions have to account for the possibility of dust erosion when estimating the thickness of the TPS. Because weight is always a critical factor in designing entry vehicles, accurate assessment of dust erosion is necessary to avoid over-design of the TPS.This study will present computational results of heatshield erosion due to dust particle impacts on the ExoMars Schiaparelli capsule if it had encountered a dust storm during its October 2016 entry. A one-way coupling approach is used where particle trajectories are computed based on an underlying CFD flow solution. Based on a distribution of particle sizes ranging from 1 to 18 microns in diameter, the Icarus material response solver predicted approximately 1 mm of TPS heatshield erosion due to particle impacts, which was 40% of the value due to material charring.

thermal protection system

Heatshield Erosion Due to Dust Particle Impacts on the Schiaparelli Capsule During Martian Entry

On October 2016, a capsule known as Schiaparelli, part of the European Space Agency (ESA) ExoMars mission, entered the Martian atmosphere. Measurements taken during the Schiaparelli descent will be used to validate computational models used to design the thermal protection system (TPS) of future Mars missions. One of the unique features of Schiaparelli entry was the possibility of a major dust storm occurring during the entry. Major dust storms are unpredictable but more likely during the Northern Autumn timeframe. In 2001, for example, regional dust storms merged into a global dust storm that blanketed much of the planet. Even though Schiaparelli did not enter during a major dust storm, future Mars missions will have to account for the possibility of dust erosion (depending on the time of year) when estimating the thickness of the TPS. Because weight is always a critical factor in designing entry vehicles, accurate assessment of dust erosion is necessary to avoid over-design of the TPS. This study will present computational results of heatshield erosion due to dust particle impacts on the Schiaparelli capsule if it had encountered a dust storm during entry.

thermal protection system

Material Response Analysis of a Titan Entry Heatshield

Accurate calculation of thermal protection material response is critical to the vehicle design for missions to the Saturn moon Titan. In this study, Icarus, a three-dimensional, unstructured, finite-volume material response solver under active development at NASA Ames Research Center, is used to compute the in-depth material response of the Huygens spacecraft along its November 11 entry trajectory. The heatshield analyzed in this study consists of a five-layer stack-up of Phenolic Impregnated Carbon Ablator (PICA), aluminum honeycomb, adhesive, and face sheetmaterials. During planetary entry, the PICA outer layer is expected to undergo pyrolysis. A surface energy balance boundary condition that captures both time- and spatial-variance of surface properties during entry is used in the simulation.

material response

Overview of Heatshield for Extreme Entry Environment Technology (HEEET) Project

The objective of the Heatshield for Extreme Entry Environment Technology (HEEET) projects is to mature a 3-D Woven Thermal Protection System (TPS) to Technical Readiness Level (TRL) 6 to support future NASA missions to destinations such as Venus and Saturn. Destinations that have extreme entry environments with heat fluxes up to 5000 watts per square centimeter and pressures up to 5 atmospheres, entry environments that NASA has not flown since Pioneer-Venus and Galileo. The scope of the project is broad and can be split into roughly four areas, Manufacturing/Integration, Structural Testing and Analysis, Thermal Testing and Analysis and Documentation. Manufactruing/Integration covers from raw materials, piece part fabrication to final integration on a 1-meter base diameter 45-degree sphere cone Engineering Test Unit (ETU). A key aspect of the project was to transfer as much of the manufacturing technology to industry in preparation to support future mission infusion. The forming, infusion and machining approaches were transferred to Fiber Materials Inc. and FMI then fabricated the piece parts from which the ETU was manufactured. The base 3D-woven material consists of a dual layer weave with a high density outer layer to manage recession in the system and a lower density, lower thermal conductivity inner layer to manage the heat load. At the start of the project it was understood that due to weaving limitations the heat shield was going to be manufactured from a series of tiles. And it was recognized that the development of a seam solution that met the structural and thermal requirements of the system was going to be the most challenging aspect of the project. It was also recognized that the seam design would drive the final integration approach and therefore the integration of the ETU was kept in-house within NASA. A final seam concept has been successfully developed and implemented on the ETU and will be discussed. The structural testing and analysis covers from characterization of the different layers of the infused material as functions of weave direction and temperature, to sub-component level testing such as 4-pt bend testing at sub-ambient and elevated temperature. ETU test results are used to validate the structural models developed using the element and sub-component level tests. Given the seam has to perform both structurally and aerothermally during entry a novel 4-pt bend test fixture was developed allowing articles to be tested while the front surface is heated with a laser. These tests are intended to establish the system's structural capability during entry. A broad range of aerothermal tests (arcjet tests) are being performed to develop material response models for predicting the required TPS thickness to meet a mission's needs and to evaluate failure modes. These tests establish the capability of the system and assure robustness of the system during entry. The final aspect of the project is to develop a comprehensive Design and Data Book such that a future mission will have the information necessary to adopt the technology. This presentation will provide an overview and status of the project and describe the status of the tehnology maturation level for the inner and outer planet as well as earth entry sample return missions.

Ellerby, Don

Verification of the Icarus Material Response Tool

Due to the complex physics encountered during reentry, material response solvers are used for two main purposes: improve the understanding of the physical phenomena; and design and size thermal protection systems (TPS). Icarus, is a three dimensional, unstructured material response tool that is intended to be used for design while maintaining the flexibility to easily implement physical models as needed. Because TPS selection and sizing is critical, it is of the utmost importance that the design tools be extensively verified and validated before their use. Verification tests aim at insuring that the numerical schemes and equations are implemented correctly by comparison to analytical solutions and grid convergence tests.

ablation

Development of an Unstructured, Three-Dimensional Material Response Design Tool

A preliminary verification and validation of a new material response model is presented. This model, Icarus, is intended to serve as a design tool for the thermal protection systems of re-entry vehicles. Currently, the capability of the model is limited to simulating the pyrolysis of a material as a result of the radiative and convective surface heating imposed on the material from the surrounding high enthalpy gas. Since the major focus behind the development of Icarus has been model extensibility, the hope is that additional physics can be quickly added. The extensibility is critical since thermal protection systems are becoming increasing complex, e.g. woven carbon polymers. Additionally, as a three-dimensional, unstructured, finite-volume model, Icarus is capable of modeling complex geometries as well as multi-dimensional physics, which have been shown to be important in some scenarios and are not captured by one-dimensional models. In this paper, the mathematical and numerical formulation is presented followed by a discussion of the software architecture and some preliminary verification and validation studies.

Material Response Model

Heatshield for Extreme Entry Environment Technology (HEEET) Development Status

The Heatshield for Extreme Entry Environment Technology (HEEET) Project is a NASA STMD (Space Technology Mission Directorate) and SMD (Science Mission Directorate) co-funded effort. The goal is to develop and mission infuse a new ablative Thermal Protection System that can withstand extreme entry. It is targeted to support NASA's high priority missions, as defined in the latest decadal survey, to destinations such as Venus and Saturn in-situ robotic science missions. Entry into these planetary atmospheres results in extreme heating. The entry peak heat-flux and associated pressure are estimated to be between one and two orders of magnitude higher than those experienced by Mars Science Laboratory or Lunar return missions. In the recent New Frontiers community announcement NASA has indicated that it is considering providing an increase to the PI (Principal Investigator) managed mission cost (PIMMC) for investigations utilizing the Heatshield for Extreme Entry Environment Technology (HEEET) and in addition, NASA is considering limiting the risk assessment to only their accommodation on the spacecraft and the mission environment.

Thermal Protection System

Arc Jet Facility Test Condition Predictions Using the ADSI Code

The Aerothermal Design Space Interpolation (ADSI) tool is used to interpolate databases of previously computed computational fluid dynamic solutions for test articles in a NASA Ames arc jet facility. The arc jet databases are generated using an Navier-Stokes flow solver using previously determined best practices. The arc jet mass flow rates and arc currents used to discretize the database are chosen to span the operating conditions possible in the arc jet, and are based on previous arc jet experimental conditions where possible. The ADSI code is a database interpolation, manipulation, and examination tool that can be used to estimate the stagnation point pressure and heating rate for user-specified values of arc jet mass flow rate and arc current. The interpolation is performed in the other direction (predicting mass flow and current to achieve a desired stagnation point pressure and heating rate). ADSI is also used to generate 2-D response surfaces of stagnation point pressure and heating rate as a function of mass flow rate and arc current (or vice versa). Arc jet test data is used to assess the predictive capability of the ADSI code.

fluid dynamics

Uncertainty Determination for Aeroheating in Uranus and Saturn Probe Entries by the Monte Carlo Method

The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.

Palmer, Grant

Coupled Fluids-Radiation Analysis of a High-Mass Mars Entry Vehicle

The NEQAIR line-by-line radiation code has been incorporated into the DPLR Navier-Stokes flow solver such that the NEQAIR subroutines are now callable functions of DPLR. The coupled DPLR-NEQAIR code was applied to compute the convective and radiative heating rates over high-mass Mars entry vehicles. Two vehicle geometries were considered - a 15 m diameter 70-degree sphere cone configuration and a slender, mid-L/D vehicle with a diameter of 5 m called an Ellipsled. The entry masses ranged from 100 to 165 metric tons. Solutions were generated for entry velocities ranging from 6.5 to 9.1 km/s. The coupled fluids-radiation solutions were performed at the peak heating location along trajectories generated by the Traj trajectory analysis code. The impact of fluids-radiation coupling is a function of the level of radiative heating and the freestream density and velocity. For the high-mass Mars vehicles examined in this study, coupling effects were greatest for entry velocities above 8.5 km/s where the surface radiative heating was reduced by up 17%. Generally speaking, the Ellipsled geometry experiences a lower peak radiative heating rate but a higher peak turbulent convective heating rate than the MSL-based vehicle.

Palmer, Grant

Numerical Prediction of Radiation Measurements Taken in the X2 Facility for Mars and Titan Gas Mixtures

Thermochemical relaxation behind a normal shock in Mars and Titan gas mixtures is simulated using a CFD solver, DPLR, for a hemisphere of 1 m radius; the thermochemical relaxation along the stagnation streamline is considered equivalent to the flow behind a normal shock. Flow simulations are performed for a Titan gas mixture (98% N2, 2% CH4 by volume) for shock speeds of 5.7 and 7.6 km/s and pressures ranging from 20 to 1000 Pa, and a Mars gas mixture (96% CO2, and 4% N2 by volume) for a shock speed of 8.6 km/s and freestream pressure of 13 Pa. For each case, the temperatures and number densities of chemical species obtained from the CFD flow predictions are used as an input to a line-by-line radiation code, NEQAIR. The NEQAIR code is then used to compute the spatial distribution of volumetric radiance starting from the shock front to the point where thermochemical equilibrium is nominally established. Computations of volumetric spectral radiance assume Boltzmann distributions over radiatively linked electronic states of atoms and molecules. The results of these simulations are compared against experimental data acquired in the X2 facility at the University of Queensland, Australia. The experimental measurements were taken over a spectral range of 310-450 nm where the dominant contributor to radiation is the CN violet band system. In almost all cases, the present approach of computing the spatial variation of post-shock volumetric radiance by applying NEQAIR along a stagnation line computed using a high-fidelity flow solver with good spatial resolution of the relaxation zone is shown to replicate trends in measured relaxation of radiance for both Mars and Titan gas mixtures.

Palmer, Grant

Coupled Fluids-Radiation Analysis of a High-Mass Mars Vehicle

The NEQAIR line-by-line radiation code has been incorporated into the DPLR Navier-Stokes flow solver such that the NEQAIR subroutines are now callable functions of DPLR. The coupled DPLR-NEQAIR code was applied to compute the convective and radiative heating rates over high-mass Mars entry vehicles. Two vehicle geometries were considered - a 15 m diameter 70-degree sphere cone configuration and a slender, mid-L/D vehicle with a diameter of 5 m called an Ellipsled. The entry masses ranged from 100 to 165 metric tons. Solutions were generated for entry velocities ranging from 6.5 to 9.1 km/s. The coupled fluids-radiation solutions were performed at the peak heating location along trajectories generated by the Traj trajectory analysis code. The impact of fluids-radiation coupling is a function of the level of radiative heating and the freestream density and velocity. For the high-mass Mars vehicles examined in this study, coupling effects were greatest for entry velocities above 8.5 km/s where the surface radiative heating was reduced by up 17%. Generally speaking, the Ellipsled geometry experiences a lower peak radiative heating rate but a higher peak turbulent convective heating rate than the MSL-based vehicle.

Palmer, Grant

Creating Web-Based Scientific Applications Using Java Servlets

There are many advantages to developing web-based scientific applications. Any number of people can access the application concurrently. The application can be accessed from a remote location. The application becomes essentially platform-independent because it can be run from any machine that has internet access and can run a web browser. Maintenance and upgrades to the application are simplified since only one copy of the application exists in a centralized location. This paper details the creation of web-based applications using Java servlets. Java is a powerful, versatile programming language that is well suited to developing web-based programs. A Java servlet provides the interface between the central server and the remote client machines. The servlet accepts input data from the client, runs the application on the server, and sends the output back to the client machine. The type of servlet that supports the HTTP protocol will be discussed in depth. Among the topics the paper will discuss are how to write an http servlet, how the servlet can run applications written in Java and other languages, and how to set up a Java web server. The entire process will be demonstrated by building a web-based application to compute stagnation point heat transfer.

Palmer, Grant

Viscosity Coefficient Curve Fits for Ionized Gas Species Grant Palmer

Viscosity coefficient curve fits for neutral gas species are available from many sources. Many do a good job of reproducing experimental and computational chemistry data. The curve fits are usually expressed as a function of temperature only. This is consistent with the governing equations used to derive an expression for the neutral species viscosity coefficient. Ionized species pose a more complicated problem. They are subject to electrostatic as well as intermolecular forces. The electrostatic forces are affected by a shielding phenomenon where electrons shield the electrostatic forces of positively charged ions beyond a certain distance. The viscosity coefficient for an ionized gas species is a function of both temperature and local electron number density. Currently available curve fits for ionized gas species, such as those presented by Gupta/Yos, are a function of temperature only. What they did was to assume an electron number density. The problem is that the electron number density they assumed was unrealistically high. The purpose of this paper is two-fold. First, the proper expression for determining the viscosity coefficient of an ionized species as a function of both temperature and electron number density will be presented. Then curve fit coefficients will be developed using the more realistic assumption of an equilibrium electron number density. The results will be compared against previous curve fits and against highly accurate computational chemistry data.

Palmer, Grant