AP01_45.vp 1 Introduction Thermal plasmas are plasmas in which the thermody- namic state approaches Local Thermodynamic Equilibrium (LTE). Although such plasmas are characterized by a single temperature common to all species and to all their ther- mal movements (translation, rotation, vibration, electronic temperature), the analysis of a thermal plasma system is highly complex. It involves an intricate interplay between fluid dynamics, turbulent transport, thermal radiation, chem- ical reactions and interaction with other phases. This paper deals with the theoretical and experimental analysis of a typical engineering system utilizing thermal plasma – a system for Chemical Vapor Deposition (CVD), in particular deposition of diamond. The system under con- sideration is shown schematically in Figure 1, where the major processes are: (1) heating of the process gas with a DC or an RF plasma torch, (2) expansion of the gas into the reactor and introduction of the diamond growth precursor species, and (3) interaction of the jet with the substrate, and diamond formation. The focus of this paper is on the second stage – the interaction of the plasma flow (plasma jet) with the surround- ing atmosphere and with the substrate located downstream. The DC torch power of the system under consideration is 6.3 kW, the jet power is about 50 % of this, i.e. 3.1 kW. The argon flow rate through the torch is 16 SLM, and the hydro- gen flow rate through the torch is 5 SLM. Methane is used as a diamond growth precursor and is injected through three probes within 3 mm of the torch nozzle exit with a total flow rate of 0.25 SLM. The reactor pressure is 12.5 kPa the substrate stand-off distance is 8 cm. The plasma torch exit temperature is of the order of 5000 K and the torch exit velocity is of the order of 2500 m/s. The system is located and operates at FG Plasma-Oberflächentechnik, Technische Universität Ilmenau, Thüringen, Germany. 2 Laminar flow modeling In this study, it was assumed that the gas exits the plasma torch with known radial profiles of pressure, velocity, enthalpy and chemical composition. These boundary conditions were obtained in two ways: (1) by a simple axial integration of the conservation equations with a 1D model of the torch region, and (2) from the spectroscopic measurements made by Jahn [1] on the plasma system under consideration. Jahn mea- sured the H� lines in the region close to the torch nozzle exit and from the results evaluated the plasma jet temperature. The temperature was evaluated assuming an LTE in the torch exit. In reality, the specification of the inflow boundary condi- tions turns out to be not so much an input but rather one of the outputs of the design analysis, and the selected final profiles are based on the best agreement of the overall simula- tion results with experiments. Assuming the upstream boundary conditions are known, one forms a set of conservation equations that are to be solved for the unknown flow characteristics, temperature/energy, pressure and species concentrations. The set includes [2] the © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 29 Acta Polytechnica Vol. 41 No. 4 – 5/2001 Analysis of a Thermal Plasma Diamond CVD System D. Kolman This paper deals with the analysis of a typical engineering system utilizing thermal plasma – a system for Diamond Chemical Vapor Deposition. It defines the system – a slightly overexpanded plasma jet impinging at a downstream-located substrate, outlines the theoretical description of the system – the Navier-Stokes and species conservation equations, and presents key theoretical results on the major and most troublesome factors influencing diamond deposition – velocity and temperature of the jet. Then, the paper demostrates the necessity to shift from a laminar to a turbulent flow description and compares both results to experiments. An explanation of the remaining discrepancy – insufficient velocity drop in the jet – is attempted. Keywords: modeling, thermal plasma, CVD. Fig. 1: The thermal plasma diamond CVD system continuity equation, two momentum equations, the energy equation, species continuity equations with temperature-de- pendent chemical kinetics, and as a closure the equation of the state. The species set may but does not include ions due to a rather low jet temperature in the region of interest. None of the viscous and/or diffusion effects, and compressibility effects can be excluded. As there is no electric arc passing through the reactor domain there is no need to solve the Maxwell-Boltzmann electromagnetic equations. Thermal dif- fusion terms and radiation are also neglected due to their small values in the situation under consideration (a rather low temperature). The obtained equation set is solved with a finite differ- ence semi-implicit pressure-based algorithm for compressible viscous flows with an arbitrary number of conservation equa- tions [3]. The final form of the program adopts a hybrid convection-diffusion flux approximation at the computa- tional cell interfaces. The 2D computational domain includes the reactor and the substrate and is non-uniformly subdivided into ~ 100 computational cells in both directions with the fin- est spacing in the vicinity of the precursor injection location and at the substrate. At the substrate, no slip boundary conditions are adopted for the momentum equation, the known radial tempera- ture distribution is imposed at the substrate surface, and a source/sink is included in the species conservation equation where the strength of the surface net production rates is given by a system of surface chemical reactions. For this purpose the momentary chemical state of the surface is described via concentrations of appropriately defined surface species (a radical site, a site ended with a hydrogen atom, etc.). Energy release by surface recombination is included in the energy balance. In the present model, the Ar-H2-CH4 kinetics mechanism consists of 34 gas phase species (Ar, H, H2, CH0�4 C2H0�6) and 52 gas phase reactions. The production rates of the individual species are given by the sum of contributions from the indi- vidual gas phase reactions. All temperature and chemical composition-dependent thermodynamic and transport properties are evaluated using the Wilke mixing rule which has been found to be sufficient for the situation under consid- eration. Equilibrium chemical composition is assumed in the torch exit. For a discussion of the surface chemistry see [2]. In all the figures presented, the jet is assumed slightly overex- panded with the nozzle pressure ~ 9000 Pa, axial velocity ~ 2700 m/s and temperature ~ 4150 K. (The supersonic char- acter of the jet, M ~ 2, is supported by the diamond shocks sometimes visible in the actual jet.) After the theoretical model was developed and the results obtained, a set of enthalpy probe measurements [4] was per- formed in the reactor – see Figure 2. These data were used to validate the theoretical results, through a comparison to gain more advanced insight into the deposition process, and to suggest the further course of investigation. In evaluation of the measurements it was assumed that the static pressure of the jet is close to the chamber pressure – a reasonable as- sumption far away from the torch exit (an assumption of this kind is necessary to evaluate the results across the bow shock formed in front of the enthalpy probe in the case of compress- ible high velocity flows). To the disappointment of the author it was found that the experimental data agree rather poorly with the predictions. The main problem is that both the axial velocity and the jet enthalpy are too high and do not fall off as fast as the experi- ment suggests. After a detailed analysis with the experimental researchers the following explanations have been proposed: (1) The spectroscopic data are too high: when computing the temperature from the H� line data a state of LTE was assumed. However, this is not likely to be the case in the torch exit, where electron temperature tends to be ele- vated above the heavy species temperature. Thus the temperature derived from the spectroscopic measure- ments corresponds more to the electron temperature than to the heavy species temperature and should not be taken as a boundary condition for heavy particle jet flow. (2) However, even when starting with a lower nozzle exit tem- perature, 4150 K, the theoretical curves do not decline as fast as the experiments require. This leads to the suspi- cion that the jet is turbulent. The transition Reynolds number for an axisymmetric jet is about 2300. In our case, the torch orifice diameter = 6 mm, plasma velocity 30 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 41 No. 4 – 5/2001 (a) (b) D D A E Fig. 2: (a) Axial jet velocity; (b) jet enthalpy – comparison of simu- lations and experiments ~ 2700 m/s, density ~ 0.1 kg/m3, and viscosity ~ 0.00016 kg/(m�s), and thus the Reynolds number ~ 10000. This value decreases to some extent with increasing torch exit temperature, but always stays above the transition limit. © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 31 Acta Polytechnica Vol. 41 No. 4 – 5/2001 D D D D D D R R R R R R Fig. 3: Contours of the axial velocity and enthalpy for a laminar, Prandtl and k-� turbulent plasma jet 3 Turbulent flow modeling Modeling of the turbulence was approached in two ways: (1) Using Prandtl’s free shear layer model for eddy viscosity [5], � � �t C u� � � � max (1) where C is constant for the particular type of mixing layer, C ~ 0.012 for a round jet, � is local density, � is the shear layer thickness taken as the nozzle radius, and umax is the centerline jet velocity at a particular distance from the torch nozzle. (2) Using the standard low Reynolds number k-� model of turbulence [6, 7]. Figure 2 shows that the character of the solution has changed drastically: One notes that the enthalpy corresponds quite well with the experimental measurements, both qualita- tively and quantitatively. However, the agreement in the axial velocity is still rather poor. In the experimental data, which are unfortunately available only from a distance 55 mm away from the nozzle exit, the axial velocity falls from ~ 2500 m/s at 55 mm to ~ 1200 m/s at 67 mm. In the simulation, the velocity falls sufficiently but much more gradually and over most of the distance to the nozzle. The following remedies and explanations are available at the present time: (1) One can modify the nozzle boundary conditions via an in- crease in the nozzle exit velocity and a decrease in the density, and pull the whole velocity profile up that way (ei- ther by decreasing the torch temperature further below 4000 K, which is unrealistic, or by reducing the jet pres- sure further and thus causing greater overexpansion). However, this fix will not lead to a sudden drop in velocity, as the experiment suggests. (2) The interpretation of the experimental data may be ques- tioned on two counts: the assumption that the jet static presure equals the chamber pressure; and the finding that the pressure ratio in the entalpy probe measurements closer to the torch than 55 mm is so high that it does not allow velocity and enthalpy evaluation for those values. These points will need further clarification with the exper- imental researchers. Finally, an illustrative comparison between individual com- putational cases is given in Figure 3. The two turbulent simula- tions produce much a shorter and much wider jet than the laminar case, as expected (k-� leading to even more extensive entrainment of the ambient atmosphere than the Prandtl model). However, again, no turbulent modeling can repro- duce the sudden velocity drop around 6 cm from the nozzle that the experiments suggest. 4 Summary and conclusions This paper has reported the process of designing and ana- lysing of a typical thermal plasma chemical vapor deposition system. It has shown the development of a theoretical model from a simple laminar model to a more demanding turbulent one, and has pointed out the reasons leading to such change in understanding. Also, close cooperation with the experi- mental researchers has been emphasized. At the current stage, the benefits of the analysis under- taken are fourfold: (1) The simulation has confirmed the su- personic character of the flow; (2) Recirculation of the gas (incl. hydrocarbon precursor) in the region close to the torch exit has been brought to attention; (3) Numerous cross-checks and critical evaluations of the experiments were enforced, specific suggestions of further experiments have been made, some potentially erroneous assumptions in the experimental work have been exposed; and (4) Fruitful discussions have been initiated on various aspects of the deposition process, incl. the precursor flow pattern and its chemical composition, and boundary layer effects on chemistry. Future work should include (1) a critical reevaluation of the enthalpy probe measurements, (2) better suited turbu- lence modeling, and (3) greater attention to the torch exit conditions. Only then trustworthy predictions and further optimization of the deposition process take place. Acknowledgement This project was in part carried out using computa- tional resources provided under Research Grants AV ČR A1057001/120/00 and MSM/216200031. References [1] Jahn, E.: unpublished results, FG Plasma-Oberflächen- technik, Technische Universität Ilmenau, 1997 [2] Kolman, D., Heberlein, J., Pfender, E.: Influence of deposi- tion parameters on diamond thermal plasma chemical vapor deposition with liquid feedstock injection. Diam. Rel. Mat. 7/1998, pp. 794 –801 [3] Patankar, S. V.: Numerical Heat Transfer and Fluid Flow. Hemisphere Publishing, New York, 1980 [4] Schwenk, A., Sember V., Gravelle, D. V., Boulos, M. I., Nutsch, G.: Enthalpy probe measurements in supersonic plasma jets. VIII. Workshop Plasmatechnik, Technische Universität Ilmenau, Germany, Juni 2000 [5] Patankar, S. V.: Computation of Heat Transfer and Fluid Flow. Complementary notes, University of Minnesota [6] Jones, W. P., Launder, B. E.: The prediction of laminariza- tion with a two-equation model of turbulence. Int. J. Heat Mass Transfer 15, pp. 301–314 [7] Lam, C. K. G., Bremhorst, K.: A modified form of the k-� model for predicting wall turbulence. Transactions ASME 103, Sept. 1981, pp. 456 – 459. Ing. David Kolman, PhD. phone: +420 2 2435 7541 e-mail: kolman@marian.fsik.cvut.cz Department of Technical Mathematics Czech Technical University in Prague Faculty of Mechanical Engineering Karlovo nám. 13, 121 35 Praha 2, Czech Republic 32 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 41 No. 4 – 5/2001