- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
International Journal of Spectroscopy
Volume 2013 (2013), Article ID 764321, 9 pages
Plasma Wind Tunnel Investigation of European Ablators in Nitrogen/Methane Using Emission Spectroscopy
Institute of Space Systems (IRS), University of Stuttgart, Pfaffenwaldring 29, 70569 Stuttgart, Germany
Received 19 March 2013; Revised 16 May 2013; Accepted 21 May 2013
Academic Editor: Jozef Kaiser
Copyright © 2013 Ricarda Wernitz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
For atmospheric reentries at high enthalpies ablative heat shield materials are used, such as those for probes entering the atmosphere of Saturn’s moon Titan, such as Cassini-Huygens in December, 2004. The characterization of such materials in a nitrogen/methane atmosphere is of interest. A European ablative material, AQ60, has been investigated in plasma wind tunnel tests at the IRS plasma wind tunnel PWK1 using the magnetoplasma dynamic generator RD5 as plasma source in a nitrogen/methane atmosphere. The dimensions of the samples are 45 mm in length with a diameter of 39 mm. The actual ablator has a thickness of 40 mm. The ablator is mounted on an aluminium substructure. The experiments were conducted at two different heat flux regimes, 1.4 MW/m2 and 0.3 MW/m2. In this paper, results of emission spectroscopy at these plasma conditions in terms of plasma species’ temperatures will be presented, including the investigation of the free-stream species, N2 and , and the major erosion product C2, at a wavelength range around 500 nm–600 nm.
For atmospheric entries concerning celestial bodies such as Saturn’s moon Titan, a suitable shielding of the atmospheric entry vehicle from the high heat fluxes generated during such entries is required. For such high speed entries, ablative materials are usually selected as heat shield materials [1, 2]. Thus, on the one hand, the characteristics of the ablators chosen need to be determined. On the other hand, the processes occurring during ablation are currently not well known and understood. Here, basic investigations into the processes in the boundary layer are required in order to facilitate a better understanding thereof and to enable the choice of and improving a suitable heat shield material. Optical plasma investigation such as emission spectroscopy is well suited to this task, as it enables the identification of the predominant plasma species over the spatial region of interest, providing a great amount of spectral information within the restricted measurement interval of less than 50 s, as well as their temporal progression.
In this paper, therefore, emission spectroscopic results of four representative measurements at two heat flux regimes and two measurement distances to the probe surface each will be presented. These measurements are both resolved in time and space in the vertical direction perpendicular to the plasma jet axis. Results are presented in terms of plasma species vibration and rotation temperatures obtained fitting a numerically obtained spectrum to the measured one. A temporal as well as a spatial profile of these temperatures is presented for each test case.
An investigation into surface recession or in-depth specimen temperatures is not included in this paper.
2. Test Facility, Probe Geometry, and Plasma Conditions
Four experiments in all were conducted in the IRS plasma wind tunnel facility PWK1 using the magneto-plasma dynamic generator RD5 [3, 4]. The vacuum tank used for PWK 1 is a 6 m long steel tank with a diameter of 2 m and a double-wall cooling. It is closed with a hemispherical part which is connected to the vacuum system and protected against heat by water-cooled copper shields. On the other side the vacuum chamber can be opened by moving a plane cover plate on a guide rail. The plasma source is not located in the vacuum tank but flanged on a conical part of the plate. The cone-shaped element of the end plate enables the plasma source to be fixed at that point. The whole plasma jet range is accessible by optical methods. The facility is equipped with a 4-axis positioning system on which the different probes and the specimen support system can be mounted. This allows the simulation of parts of the reentry trajectories by automatically moving the specimen in different plasma flow regimes. A certain point on a reentry trajectory is defined by the ambient pressure, the axial distance to the generator, the arc current, and the mass flow. Together, these values result in adequate enthalpy, total pressure, and temperature or heat flux. Optical glass windows allow pyrometric temperature measurements on the front side of the specimen at distances from the plasma source between about 50 mm and 1 m. Moreover, optical measurements perpendicular to the plasma jet axis are possible through three movable flanges of three optical glasses each, which are located on both sides of the tank opposite to each other and on the top. The electric power is supplied by a current-regulated thyristor rectifier consisting of six identical units supplying 1 MW each. These may be connected in series or parallel, thus varying the desired output level of current, voltage, and power. The current ripple is less than 0.5%. The maximum current is 48 kA supplied at 125 V, and the maximum voltage is 6 kV at a current of 1 kA. Self-field MPGs of different sizes have been developed over the past decades at IRS; they are employed successfully in the IRS plasma wind tunnels PWK 1 and PWK 2. The nozzle-type MPG plasma generators consist of two coaxial electrodes, separated by water-cooled neutrodes. The nozzle exit, which is also a water-cooled segment, forms the anode. The cathode, made of 2% thoriated tungsten, is mounted in the center of the plenum chamber. The arc is ignited by a Paschen breakdown. The current passes through the expansion nozzle from the cathode tip to the end of the nozzle. The test gas is dissociated and partly ionized. To avoid oxidation of the cathode, only the nitrogen component of the test gas is fed in along the cathode into the plenum chamber, heated up by the arc, and accelerated partly by thermal expansion and partly by electromagnetic forces due to the self-induced magnetic field in the nozzle. The magnitude of the magnetic acceleration force strongly depends on the current level of operation. With this MPG, the oxygen needed for the duplication of high enthalpy air flows is fed in radially at a high velocity at the supersonic part of the nozzle but still within the arc region. The various gas injection points enable the operation of the MPG with different gas mixtures. So the investigation of entry maneuvers into the atmospheres of other celestial bodies such as Titan and Mars, containing CH4 and CO2, is possible .
Special efforts have been made to minimize the erosion of the plasma generator. To avoid a spotty arc attachment on the anode, which would cause a contamination of the plasma flow, a small amount of argon is injected tangentially along the anode contour. This method has been shown experimentally to eliminate anode erosion. Within the whole region of operation the only contaminating part of the MPG is the cathode. The 2% thoriated tungsten cathode reaches more than 3000 K during steady-state operation. The high temperature and the low work function of the cathode result in a diffuse arc attachment and, consequently, a very low cathode erosion rate at the order of sublimation. The very low cathode erosion results in operation periods of the MPG of hundreds of hours.
For the MPG generator RD5 with a nozzle exit diameter of 125 mm, the mass flows for air are between 0.3 g/s and 50 g/s at current levels between 200 A to 4 kA and power levels of 40 kW to 1 MW, while the average specific enthalpy at the nozzle exit varies between 2 MJ/kg and 150 MJ/kg.
For testing large TPS structures an MPG generator with an enlarged nozzle was built. The nozzle diameter of RD7 is 320 mm, which enables the investigation of TPS materials and substructures up to 400 mm in diameter .
For the experimental entry simulations for Titan, a gas mixture of 98.5% N2 and 1.5% CH4 was chosen. Here, two test conditions have been chosen, at heat fluxes of 1.4 MW/m², at a testing time of 25 s, and of 0.3 MW/m², at a testing time of 50 s; see Table 1. For each test condition two tests have been performed using the common European ablative material AQ60. This material is a standard ablation-based heat shield material for atmospheric entry capsules entering atmospheres under modest heat loads. It was, for example, used as heat shield material for the European Huygens Titan entry probe and is an isotropic material made of 10% (volume) glass fiber and 6% phenolic resin. The material has a very low density (0.3 g/cm3) and very good insulation properties; see also . The dimensions of the samples are approximately mm in length with a diameter of the ablative material of mm. The actual ablator has a thickness of mm. The ablator is mounted on an aluminium substructure and installed into a water-cooled copper probe with a diameter of 80 mm; see Figure 1.
Figure 1, however, depicts the constellation of the sample holder with the steady state calorimetric heat flux sensor inserted. This sensor has the same diameter as the ablation samples and allows the determination of copper-based heat fluxes including the simultaneous measurement of the Pitot pressure (see the bore in the center of the front view of the plasma probe system). Heat fluxes and total pressures were correspondingly measured in previously performed calibration tests using the combined Pitot pressure and copper heat flux probe setup with a diameter of 80 mm and the insert for heat flux determination of 40 mm diameter.
Table 1 shows the facility conditions and, in addition, the individual parameters of conditions C1 and C2.
Surface temperatures given in Table 1 (bottom) were obtained using the linear pyrometer LP3 at an operating wavelength of 726.5 nm, which covers a temperature range up to 3000°C and even higher. This pyrometer is manufactured by KE, Germany.
Mach numbers have been estimated according to Herdrich and Petkow  using adiabatic coefficients of : This formula is valid for the subsonic case, only. However, as shown in Table 1, the plasma condition is clearly subsonic and the boundary layer is artificially broadened by the flat face probe design, enabling measurements at the axial distances of 3 mm and 5 mm investigated here.
Local specific enthalpies are generally estimated using Marvin and Pope’s  method, employing an effective probe nose radius as determined by Boisen and Curtiss  of 2.7 times the probe radius. However, for a methane/nitrogen atmosphere, the factor employed by Pope is not known; hence no enthalpy calculation for the plasma conditions investigated here could be performed.
3. Experimental Setup for Emission Spectroscopy
The experimental setup for the emission spectroscopic investigation of the plasma jet and ablative measurements is displayed in Figure 2. A monochromator/spectrograph in Czerny-Turner configuration, ACTON SpectraPro 2750 (by Princeton Instruments), with a focal length of , an aperture ratio of , and both 300 L/mm and 600 L/mm gratings are used. A micrometer controlled side entrance slit is included in the spectrometer. For data acquisition, a CCD camera (ANDOR Newton DU920N-OE) is used. Measurements are conducted in imaging mode. Thus, spectra are composed of 255 separate spectra recorded for each row of the CCD and provide additional information on the vertical spatial distribution of the plasma (radial intensity profile of the plasma jet) and the ablative process.
Light emitted from an area mm or mm in front of the probe surface is focused onto the spectrometer entrance slit via two mirrors, a plane mirror reflecting the emission of the plasma jet onto a concave mirror with a focal length of mm and a ratio of . For plasma characterization at wavelengths above 450 nm, a suitable filter is included in the setup to avoid influence of second-order spectra, directly in front of the spectrometer entrance slit. Additionally, an aperture has been installed in front of the optical setup, reducing the intensity and the influence of aberrations.
Test parameters including spectrometer grids, 300 L/mm and 600 L/mm, exposure times (1 s at a measurement frequency of 1.1874 Hz), and wavelength ranges are given in Table 2. A wavelength range between approximately 525–575 nm has been chosen, as in this range, C2, a well-known product in ablation processes using carbon containing ablators, is a strong radiator.
Spectra were taken mm and mm in front of the probe surface. Here, by placing a light source with a 0.2 mm aperture in the measurement region and moving it in 0.2 mm increments in the axial direction, one obtains an axial intensity distribution. Fitting a Gaussian to these intensities yields an area of 1.5 mm (axial) width at half maximum, which the radiation incident on the spectrometer originates from.
In the vertical direction, 51 mm of the plasma jet core has been observed and is distributed over 255 pixel lines, yielding a maximum possible vertical resolution of 0.2 mm/pixel. Effectively, vertical resolution is 0.6 mm at half maximum. This vertical resolution was obtained by placing a light source with a 0.2 mm aperture in front of the optical setup and determining the width at half maximum of the resulting peak on the camera’s CCD.
Summarizing, all spectra are an integration of an area of mm2 at half maximum, taken mm and mm in front of the sample.
Spectra were obtained using 300 L/mm and 600 L/mm grids at a maximum theoretical spectral resolution of 0.12 nm/pixel and 0.06 nm/pixel, respectively. The wavelength calibration has been conducted comparing lines of a conventional mercury lamp to pixel numbers and using a 3rd degree polynomial fit. Intensity calibration was conducted using an integrating sphere (Ulbricht sphere). It consists of a hollow spherical cavity with its interior covered with a diffuse white reflective coating, with small holes for entrance and exit port(s). Therefore, it provides a uniform diffusing effect. Light rays incident on any point on the inner surface are distributed equally to all other points. Hence, effects of the original direction of the light are minimized. An integrating sphere may be thought of as a radiation diffuser which preserves power but destroys spatial information. It is typically used with some light source and a detector for optical power measurement. Correspondingly, measurements of power flux can be performed without inaccuracy caused by the directional characteristics of the source. Therefore, the sphere creates a reference radiation source.
4. Emission Spectroscopic Free-Stream Characterization
Before analyzing ablation species and processes, the free-stream needs to be characterized for both plasma conditions at 1.4 MW/m² and 0.3 MW/m².
For the nitrogen/methane gas composition used, CN, N2, and have been identified as the predominant plasma species, CN, resulting from the methane fraction of the gas mixture.
In Figures 3 and 4 two observed spectra are shown. In addition, emission spectra have been simulated using the PARADE  (plasma radiation database) program. PARADE is the plasma radiation database that was developed by Fluid Gravity Engineering and the Institute of Space Systems (IRS) for ESA. The currently available version can be obtained upon request from ESTEC (within member state entities). The code can be either used independently (as in this work) or in a coupled mode that is coupled to an adequate CFD code. A Monte Carlo-based radiation transport code is used for coupling (HERTA). Examples are simulated spectra for entry capsules such as STARDUST  or the simulation of test cases within ESA’s radiation working group .
Rotational and vibrational temperatures at the plasma jet center have been extracted, computing the optimal fit of a simulated spectrum to a measured one. Since, for the investigation of the free-stream temperatures at 3 mm and 5 mm distance to the probe surface the difference in the plasma composition is very small, results for those two distances for the free-stream are not presented separately.
All spectra have been normalized with their respective maxima in order to facilitate comparability between the molecular bands of the simulated and measured spectra.
In Figure 3, for heat fluxes of 1.4 MW/m², N2, , and CN have been identified as the predominant emitters. At the plasma jet center (plasma jet axis), CN shows vibrational temperatures of approximately K, and rotational temperatures of K. Temperatures for N2 range from K and K and for from K and K.
In Figure 4 for heat fluxes of 0.3 MW/m² and CN are the predominant species; see Figure 4. For plasma condition C2 no N2 emission was detected, which may however be due to the emission of N2 being very weak in comparison to and CN, which are strong radiators at the wavelength range (290 nm–410 nm) investigated for the free-stream. At the plasma jet center, CN shows vibrational temperatures of K and rotational temperatures of K. Temperatures for range from K and K.
5. Emission Spectroscopic Results for 1.4 MW/m2
For the study of ablation, experiments concentrated on the emission of a common ablative plasma species, C2. Here, the wavelength range chosen for all detected spectra facilitates the observation of a Swan band of C2, 1, 2 and 0, 1.
Swan bands were originally a characteristic of the spectra of carbon stars, comets, and of burning hydrocarbon fuels. They are named after the Scottish physicist William Swan who first studied the spectral analysis of radical carbon C2 in 1856.
The wavelength range allows for the observation of C2 without disturbing influence of other plasma species. Particularly CN was to be avoided, as it is a strong radiator and ablation product, but also, here, present in the free-stream.
For plasma condition C1, emission spectra of C2 were taken using a 600 L/mm grid at a wavelength range of 525 nm–575 nm at a distance of mm to the probe surface. A second test has been performed at a distance of mm to the probe surface using a 300 L/mm grid at a wavelength range of approximately 500 nm–600 nm. Plasma species temperatures were obtained by fitting simulated spectra to measured ones, while changing rotational and vibrational temperatures of the predominant species (CN and ) and adjusting the density quotient of N2//C2. The temperature of N2 does not show any significant influence on the spectra; thus they were set to the equilibrium temperatures K. Table 3 shows an overview outcome for the C2 temperatures for 1.4 MW/m2 heat flux at the distanc 3 mm and 5 mm.
5.1. mm Axial Distance to Probe Surface
At a heat flux of 1.4 MW/m², the investigation of during ablation yielded constant temperatures of K and K perpendicular to the plasma jet axis along the probe surface.
Although considerably weaker compared to the emission band of CN, parts of the C2 Swan band could be well observed during the first second. Figure 5 exemplary shows one observed spectrum and its respective optimum PARADE fit.
Leaving the density quotient of N2//C2 = 3000/70/1 constant, temperatures of C2 increase with perpendicular distance to the probe center from K and K at the probe center to K and K at a distance of 32.4 mm to the probe center. However, it is interesting to note that simulated and measured spectra are also in acceptable agreement if the density of C2 is halved at that distance.
Additional to the spatial distribution of C2 emission along the probe surface, C2 temperatures and concentrations (in terms of emission intensity) show a temporal progression. C2 emission intensity decreases with time. After approximately 5 s of exposure of the probe to the plasma, C2 intensity in the molecular band observed in this study is almost nonexistent. Reliable simulations could therefore only be obtained until 3 s into ablation. Here, C2 emission at 32.4 mm vertical distance to the probe symmetry axis was too weak to yield reasonable simulation results.
Summarizing, at plasma condition C1 at 5 mm axial distance to the probe surface, showed constant temperatures of K and K along the vertical axis, as well as over time.
C2 temperatures at the probe center ranged between K and K at the onset of ablation, 1 s after first exposing the probe to the plasma, and K and K after 3 s of exposure to the plasma. C2 temperatures increased slightly along the vertical (probe) axis to K and K at a distance of 32.4 mm to the probe center. This increase of temperatures along the probe surface was also observed after 2 s and 3 s, respectively, of exposure to the plasma; see Figure 5.
During 3 s the density quotient of the simulated molecules was adapted from approximately N2//C2 (1 s) = 3000/70/1 to N2//C2 (3 s) = 3000/70/0.4.
5.2. 3 mm Axial Distance to Probe Surface
At an axial distance of 3 mm to the probe surface, at the average heat flux of 1.4 MW/m², the investigation of during ablation also yielded constant temperatures of K and K perpendicular to the plasma jet axis along the probe surface. These temperatures were constant for all radial distances to the probe center investigated, as well as constant in time.
At a density quotient of N2//C2 = 5500/70/1, rotational and vibrational temperatures of C2 increased along the vertical distance to the probe center from K and K at the probe center to K and K at 32.4 mm distance to the probe center; see Figure 6.
At 3 mm axial distance to the probe surface C2 also showed a transient temporal behavior, temperatures increasing slightly with time towards K and K (at 32.4 mm to the probe center and N2//C2 = 5500/70/1), while particle density quotients decreased to N2//C2 = 8500/90/0.9.
6. Emission Spectroscopic Results for 0.3 MW/m2
For plasma condition C2, emission spectra of C2 were taken using a 600 L/mm grid at a wavelength range of 525 nm–575 nm at an axial distance of mm to the probe surface. A second test has been performed using the same spectrometer grid at the same wavelength range at a distance of mm to the probe surface.
6.1. mm Axial Distance to Probe Surface
At an average heat flux of 0.3 MW/m² at 5 mm axial distance to the probe surface the absolute emission intensity of C2 was very weak and only spectra measured during the first 2 s of the experiment could reasonably be compared to simulated spectra.
Here, for both density quotients used (N2//C2 = 3000/90/0.5, 4000/90/0.9, and 5000/90/0.5), at all radial distances to the probe center investigated (0 mm and 10.8 mm), both and C2 yielded constant temperatures for both time frames (after 1 s and 2 s of exposure to the plasma). For a rotational and vibration temperatures of K and K and for C2 of K and K were obtained.
6.2. mm Axial Distance to Probe Surface
At 3 mm distance to the probe surface at an average heat flux of 0.3 MW/m² to the probe, C2 emission has been observed for more than 9 s, the longest time span of all four experiments.
For the numerical simulations three different densities were used; see Table 4, where the temperature of N2 was again kept at a constant K, and yielded temperatures of K and K for all vertical distances and times investigated.
For all densities, times (1 s to 8.3 s) and vertical distances (0 mm and 10.8 mm to the probe center) were investigated; C2 yielded a vibrational temperature of 6000 K, except for the last spectrum investigated 9.5 s into ablation, where a vibration temperature of 5500 K has been obtained; see Table 5. The respective rotational temperatures are 5000 K and 4500–7000 K. Rotational temperatures at 10.8 mm vertical distance are generally slightly higher than at the probe center at all density quotients investigated, see also Figures 7 and 8.
To ascertain the species identified Table 6 highlights the predominant molecular species, their respective transitions, and band head positions.
In the present study, C2 emission exhibits a transient behavior with time. Several explanations are feasible for this phenomenon: Pyrolysis, surface/gas phase chemistry, and spallation. A numerical rebuild by Fluid Gravity Engineering  indicates that the temporal, transient behaviour of C2 emission may be due to spallation along the samples’ surface.
Pyrolysis is ruled out, as in the rebuild; the pyrolysis rate increases with time, and C2 is not a common species as pyrolysis gas.
Considering gas phase/surface chemistry, surface temperatures are too low for C2 to be produced by sublimation and C2 is not a common decomposition product of pyrolysis gases.
Thus, as a hypothesis, spallation is likely to be the cause for C2 transient behaviour. Here, small carbon particles are assumedly ejected into the flow due to the pressure built up of pyrolysis gases in the initial stages of testing. These particles are then transported and heated in the flow, which has been characterized by emission spectroscopy and is therefore known to have temperatures higher than approximately 5000 K. Due to these high temperatures, C2 may therefore be produced by sublimation from the spallated particles. As the test proceeds, fissures are formed allowing pyrolysis gases to escape and hence decreasing the rate at which particles are ejected into the flow.
For a nitrogen/methane gas mixture all in all four ablative probes, consisting of the European ablative material AQ60, have been investigated at two different heat flux regimes (1.4 MW/m² and 0.3 MW/m²) at two different axial distances to the probe surface (5 mm and 3 mm) each. Emission spectra of C2 were taken during these tests, yielding temporal and spatial information on rotational and vibration temperature distributions.
After first exposure of the probe to the plasma, rotational temperatures of C2 showed an increase for approximately 5 s, see Figure 9, and are then stable at the vertical distance of 10.8 mm, whereas at the probe center, rotational temperatures decreased slightly after 5 s.
For both investigated heat fluxes, the temperatures of C2 showed a spatial distribution in the vertical direction along the probe surface. Mostly, temperatures increased slightly with distance to the center point, which may be due to the particles being heated by the flow.
Furthermore, C2 emission in the observed band showed a transient behavior, intensity decreasing with time, with no visible emission observed after a maximum of 9 s (0.3 MW/m², 3 mm vertical distance to probe center). In a very short bullet type summary the following can be derived.(i)Ablative heat shield materials are characterized in N2/CH4 plasma as missions to Titan are planned. (ii)Plasma wind tunnel tests were performed at different heat flux levels, 1.4 MW/m² and 0.3 MW/m². (iii)Emission spectroscopy analysis for free-stream and boundary layer (N2, and C2) is performed.(iv)C2 emission behaves transient and numerical rebuild implies spallated particles’ sublimation. (v)As time proceeds, formed fissures promote pyrolysis gas escape. Particles’ ejection is decreased.
Conflict of Interests
There is no conflict of interests.
The authors would like to thank Nelli Eswein, Andreas Steinbeck, and Hannes Fulge for their support during plasma wind tunnel tests. The authors would also very much like to thank ESA for their support within the Aerothermochemistry project, especially for the procurement of the test specimen. Furthermore, The authors want to thank Fluid Gravity Engineering Ltd. for providing numerical simulations within the Aerothermochemistry project, giving a possible explanation for C2 temporal behaviour. This work was supported by the German Research Foundation (DFG) within the funding programme Open Access Publishing.
- W. Röck and M. Auweter-Kurtz, “Experimental investigation of huygens entry into the titan atmosphere within a plasma wind tunnel,” in Proceedings of the 30th AIAA Thermophysics Conference, San Diego, Calif, USA, June 1995.
- W. Röck, M. Auweter-Kurtz, P. Dabalà, H. Früholz, H. Habiger, and S. Laure, “Experimental simulation of the entry of huygens into the titan atmosphere for the thermal protection qualification,” in Proceedings of the 44th Congress of the International Astronautical Federation, Graz, Austria, October 1993.
- G. Herdrich, M. Fertig, and S. Löhle, “Experimental simulation of high enthalpy planetary entries,” The Open Plasma Physics Journal, vol. 2, no. 15, pp. 150–164, 2009.
- G. Herdrich, S. Löhle, M. Auweter-Kurtz, P. Endlich, M. Fertig, and S. Pidan, “IRS ground-testing facilities: thermal protection system development, code validation and flight experiment development,” in Proceedings of the 24th AIAA Aerodynamic Measurement Technology and ground Testing Conference, Paper AIAA 2004-2596, pp. 2004–2596, Portland, Ore, USA, June–July 2004.
- G. Herdrich and D. Petkow, “High-enthalpy, water-cooled and thin-walled ICP sources characterization and MHD optimization,” Journal of Plasma Physics, vol. 74, no. 3, pp. 391–429, 2008.
- J. G. Marvin and R. B. Pope, “Laminar convective heating and ablation in the mars atmosphere,” AIAA Journal, vol. 5, no. 2, pp. 240–248, 1967.
- J. C. Boisen and H. A. Curtiss, “An experimental investigation of blunt body stagnation point velocity gradient,” ARS Journal, vol. 29, pp. 130–135, 1959.
- A. J. Smith, A. Wood, J. Dubois, M. Fertig, and B. Pfeiffer, “Plasma Radiation Database PARADE V2. 2,” Tech. Rep. TR28/96 3, Fluid Gravity Engineering, 2006.
- H. Liebhart, R. Wernitz, G. Herdrich et al., “Advances for radiation modeling for earth re-entry in parade: application to the STARDUST atmospheric entry,” in Proceedings of the 43rd AIAA Thermophysics Conference (AIAA' 12), pp. 25–28, New Orleans, La, USA, June 2012.
- H. Liebhart, U. Bauder, G. Herdrich, S. Fasoulas, and H. P. Röser, “Numerical modeling of radiative and convective heat flux for entry flights in CO2 containing atmospheres,” in Proceedings of the 43rd AIAA Thermophysics Conference, AIAA-2012-3195, pp. 25–28, New Orleans, La, USA, June.
- J. Beck, “Validation of aerothermal chemistry models for reentry applications: executive summary,” ESA Report and Final Presentation TTS052-11, 2011.