- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- 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
Journal of Chemistry
Volume 2013 (2013), Article ID 836187, 11 pages
Numerical Modeling Chemical Vapor Infiltration of SiC Composites
Institute of Solid Mechanics, Karlsruhe Institute of Technology (KIT), Kaiserstrasse 12, Room 310, Building 10.91, 76131 Karlsruhe, Germany
Received 29 June 2012; Accepted 22 October 2012
Academic Editor: Ekaterina Tsipis
Copyright © 2013 Yaochan Zhu and Eckart Schnack. 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.
The multiphase field model for chemical vapor infiltration (CVI) of SiC/SiC composites is developed in this study, thereby to reproduce the microstructure evolution during CVI process and to achieve better understanding of the effect of process parameters (e.g., temperature, pressure, etc.) on the final product. In order to incorporate the thermodynamics of methyltrichlorosilane (MTS) pyrolysis into phase field model framework, the reduced chemical reaction mechanism is adopted. The model consists of a set of nonlinear partial differential equations by coupling Ginzburg-Landau type phase field equations with mass balance equations (e.g., convection-diffusion equation) and the modified Navier-Stokes equations which accounts for the fluid motion. The microstructure of preferential codeposition of Si, SiC under high ratio of H2 to MTS is simulated and the potential risk of blockage of the premature pores during isothermal CVI process is predicted. The competitive growth mechanism of SiC grains is discussed and the formation process of potential premature pore blockage is reproduced.
The lightweight SiC fiber-reinforced SiC matrix composite (SiC/SiC) is of great importance in industry because of its excellent thermophysical and mechanical properties, for example, high strength, light weight, and excellent chemical-thermal resistance [1, 2]. It is expected to be served as structural and functional components, especially in extreme environments, such as for gas turbine, spacecrafts, and nuclear fusion reactors. This composite is now industrially widely produced by chemical vapor infiltration (CVI) process, a unique method for producing continuous fiber ceramic composite which belongs to chemical vapor deposition (CVD) on the internal surface of a porous preform. Among the different precursors for the deposition of SiC, methyltrichlorosilane (MTS) is the most frequently used precursor due to the equivalent ratio of Si to C, wide deposition temperature range, and low toxicity [3, 4]. If MTS is the only precursor for manufacturing of the SiC matrix by infiltration, SiC with excess Si, SiC with excess C, and stoichiometric SiC can be deposited by pyrolysis of MTS as shown in Figure 1.
Extensive experimental works have been devoted to investigate the influence of the infiltration conditions on SiC deposition, for example, the temperature, pressure, and the kind of precursor. Many studies also concentrated on the morphology of the deposited matrix. Most of them are focused on the composition of the deposited matrix in terms of the C/Si atomic ratio caused by macroscopic infiltration conditions [5–7]. Considering that the experimental investigation is time consuming and costly, more and more studies are now focused on numerically modeling chemical vapor infiltration process with the development of the computational technology and advanced numerical algorithms. Based on the conception of porosity evolution, both the chemical reaction mechanism and mass transport are taken into account in the model [8–10].
Besides the macroscale porosity, the morphology of the deposited matrix and microscale components play a critical role in determining the final properties of the composites. There are mainly two kinds of methodologies to describe the microstructure evolution on mesoscale: sharp interface method and diffusive interface method. The sharp interface method describes the interface between different phases or species explicitly by the mathematical analytical model, which always results in the complicated refinement of the mesh near free moving boundary during numerical discretization and is nearly impractical 3D cases. As an example of diffusive interface method, the level-set method uses the convection equation to implicitly capture the evolution of microstructures, which may be numerically instable. Recently, the phase field method has become an important and extremely versatile technique for modeling microstructure evolution on mesoscale . The phase field method inherits the advantages of the level-set method and the diffusive region is introduced to replace the sharp interface in the sharp interface method. The phase field order parameter is introduced to denote the thermodynamics state there. On the other hand, the convection diffusion equation is applied to describe the evolution of the interface to avoid the numerical problems in the level-set method. In the past two decades, the phase-field method has exhibited its potential in different applications.
Until now, few studies investigating the morphological phenomena during the CVI process have been published. Palmer [12, 13] has studied the morphological instabilities during CVD process. Jin et al. [14, 15] tried to apply the level-set method to predict the evolution of the microstructures during CVI process. There is still a challenge in the morphological study of the CVI process. In the our previous research [16, 17], the phase field model has been successfully developed and applied in predicting the microstructure evolution during CVI process of C/C and SiC/SiC composites in two dimensions. Due to the codeposition of condensed solid phases -SiC(s) with Si(s) and/or C(graphite), a multiphase field model is formulated in this paper to model the CVI process of SiC fiber-reinforced SiC matrix composite (SiC/SiC).
2. Multiphase Field Model of the Chemical Vapor Infiltration
2.1. Phase Field Equations
In phase field model, the spatial diffusive distribution of phase field order parameter is used to denote the gas-solid interface implicitly instead of the explicit sharp interface (shown in Figure 2). Therefore the spatiotemporal behavior of phase field order parameter reproduces the deposited SiC microstructure evolution. Due to the fact that the codeposition of condensed solid phases occurs under specific processing conditions during the SiC CVI process, the multiphase field order parameters are adopted to describe not only the gas-solid but also the solid-solid interfaces. The existence of different bulk phase is represented by phase field order parameters ( is the number of phases or species) to be 1, otherwise 0 in the rest phases, and in the interfaces between and the other phases. The sum of these phase field order parameters should be conversed as the following equation: here denotes the multi-component gaseous phase and , , denote silicon, silicon carbide, and carbon, respectively.
Benefited from the equivalent ratio of silicon to carbon, methyltrichlorosilane (, MTS) is adopted as precursor in the isothermal chemical vapour infiltration (ICVI) process of SiC/SiC composites. During the process, MTS as precursor and hydrogen () as catalytic gas and as diluting gas are pumped into the reactor. The deposition production (SiC with Si and/or C) from MTS depends on the temperature, pressure and the ratio of . In the previous work [16, 18], the lumped reaction is adopted to predict the dissociation of MTS using the Arrhenius’ equation for the dependence of the reaction rate on temperature. For this lumped reaction, it still seems not being well understood which is indicated by the following facts: the reaction order with respect to pyrolysis of MTS has been reported varying from −3 to 2.5 and activation energies from 40 kJmol−1 to 400 kJmol−1  depending on the infiltration condition. Now experimental research increasingly focuses on the reduced mechanism of thermal decomposition of MTS. Instead of simplified morphological system free energy description under lumped reaction in previous research, the reduced mechanism permits a detailed description which is necessary for constructing phase field model. On the other side, only several species are considered in the reduced mechanism instead of hundreds of species in equilibrium gas-phase reaction thermodynamics calculation of CVD/CVI system.
In the following, based on the work by Zhang and Huttinger , Allendof and Melius , and Fitzer and Kehr  a model of both homogeneous gas phase chemistry and heterogeneous surface chemistry is presented. Here and species are assumed to be the carbon sources, while and are assumed to be the silicon sources for the deposition of C, Si, and SiC.
Dissociation reaction of MTS is as follows: C-bearing species is as follows: Si-bearing species is as follows: Formation of the depositions ( indicates the adsorption of the intermediate species on free surface site) is as follows: The above chemical reactions are governed by thermodynamics and kinetics. The thermodynamics defines the driving force which indicates the potential direction that the reaction will proceed. The kinetics defines the transport process and subsequently determines the rate-control step during chemical vapor deposition. Both of these two aspects of information should be introduced into the multiphase field model. The free energy changing for a reaction is defined as where is the stoichiometric coefficient of species in reaction (negative for reactants, positive for products). () is the ideal gas constant. is the reaction quotient which is related to the activity and partial pressure of species ; is the standard free energy of formation of species at temperature and 1 atm as where with coefficient being obtained for each species from the thermodynamic database in the NASA equilibrium code. The standard thermodynamic data, that is, the enthalpies and the entropies at 298 K can be found in JANAF thermochemical tables . The thermodynamics properties of individual substances can be refered to .
The Ginzburg-Landau type of system free energy can be expressed as here is the gradient energy coefficient and is the height of the energy barrier. is dependent on the interface tensor and the interface thickness . Anisotropy of the interface energy is incorporated in the model by the dependence of on the orientation  (, ). The Lagrange multiplier is introduced to account for that the sum of phase field in the system is conserved. The first item with the gradient of phase fields accounts for energy penalty for the formation of interfaces/boundaries. The second item is the energy barrier to prevent phase transition occurring randomly. is the chemical free energy density of local bulk phase, which can be interpolated by free energy density of the coexisting phases using the weight function : here the weight function takes the form to ensure the high numeric efficiency. Thereby any point in the system is assumed to be a mixture of phases with fraction of for phase. In previous work [16, 18], the is phenomenologically expressed by for gas phase and for solid phase, where , , and are phenomenological constants. These two expressions also indicate that the driving force of chemical reaction is change of the Gibb’s free energy which is accordant with thermodynamic principle; however, they fail to indicate that that Gibb’s free energy is highly dependent on temperature and pressure, not only on the concentration of species.
The evolution of phase order parameter fields follows the time-dependent Ginzburg-Landau equations (TGDL) assuming minimization of total system free energy as follow: where is the identification function to locate interface between phase and phase. By interface field theory , where the step function if and otherwise. The interface field theory assumes that the junction is composed of the multipairwise interfaces and thereby the Lagrange multiplier limit is canceled. is the interface mobility coefficient which includes kinetics information of the reaction should come from experiments. is the number of phases really coexisting in a given point and can be expressed as
2.2. Mass Transport
In the following, the mass transport equation in the multiphase field model will be deduced. The modified mass transport of the species in the process of ICVI can be stated in the compact form of a time-dependent convection-diffusion problem: where () is the molar concentration of gaseous species . indicates that the diffusion of each species in solid phases is disregarded due to its low diffusivity. () is the velocity vector. () is the diffusivity which is a function of effective binary diffusion coefficient in gas. () is the effective chemical reaction rate. In this paper the Milke equation  is used to calculate the effective binary diffusion coefficient : where is the molar fraction of component . is the binary diffusion coefficient of species and , which can be calculated by Chapman-Enskog equation: here () is the molecular weight of the species , () is the Avogadro number, (Pa) is the pressure of the gas mixture, is the collision diameter of species and , and is the dimensionless collision integral for species , . Using mixture rule, then the diffusivity can be expressed as: where an artificial small diffusivity is assumed for solid phase in order to prevent from being zero in the diffusive interface near solid phase side.
Two kinds of reactions occur during the ICVI process of MTS: homogeneous gaseous reaction and heterogeneous surface reaction. The effective chemical reaction rate can be regarded as a combination of the homogeneous gaseous reaction and heterogeneous surface reaction with the reaction intensities in gas phase and in the diffuse interface, respectively. It is assumed that the heterogeneous process ends in the full-deposited solid region (shown in Figure 2). To account for this property of the process, a step function is introduced in present work: Based on the above function, we can describe the effective chemical reaction rate in the diffuse domain of species as with as the homogeneous reaction rate (gaseous reaction rate) of species and as the heterogeneous reaction rate (surface reaction rate) of species .
2.3. Fluid Motion
In the previous work [16, 18], the solid phase is assumed to be Newtonian fluid with a large viscosity. Thereby the system consists of two Newtonian fluids: gas fluid and solid fluid, while, in this work, the solid-fluid system is assumed to be rigid and stationary. Assigning a zero-velocity to the solid, the calculation for solid phase is avoided instead of still solving the equations in previous work. Here a modified Navier-Stokes equations is employed to describe the evolution of the gaseous fluid motion: here () is the velocity. () is the density (the same density assumed for different gaseous species for simplicity), and () is the viscosity. is the unit metric and is the inversion sign. is the parameter depending on the gradient energy coefficient . The last item in (27) is the force item, which is chosen such that the no slip condition can be accurately reproduced regardless of the diffusive interface thickness .
2.4. Multiphase Field Model of CVI
2.5. Finite Element Implementation for Multiphase Field Model
The numerical discretization for the phase field formulation consists of two parts: temporal discretization and spatial discretization. In this work, we use the implicit Euler finite difference approximation for temporal discretization and Galerkin weighted residual method for spatial discretization. The implicit Euler format discretization is of the advantage of stability and loose limit on time step:
In the following, we will propose the weak formulations for (28). Based on the Green formula, the solutions of phase field model should satisfy the following equations for the arbitrary test functions , , , and : with as the outward unit normal.
Here the spatial discretization of (32) is proposed as an example. Defining the approximation , where is vector composed by shape functions, , is the node number in element. Choosing shape functions as test functions and assuming the boundary conditions , then we have
Notice that the coefficient matrix is large, sparse, and nonsymmetric. Here the parallel generalized minimal residual method (GMRES) is employed as an iterative linear solver.
3. Result and Discussion
3.1. Initial and Boundary Conditions
In this paper, the evolution of the microstructure in unidirectional fibre-reinforced composites is studied. The geometric model is shown in Figure 3. For three-dimensional case, the one-fiber or the multifiber model can be used. Due to the symmetry along the radial section of the fiber, this three-dimensional model often can be reduced to a 2D model as illustrated in Figure 3. To accomplish the multiphase field model equations system, appropriate initial and boundary conditions should be supplemented. For the one-fiber or multifiber model in 3D case, the boundaries are composed of (a) the inlet face and outlet face which are normal to the fiber direction; (b) the artificial wall faces with same physical properties which are parallel to the fiber direction; (c) the gas/solid boundary which is especially required by fluid motion for no slip boundary condition is defined by .
For the multiphase field order parameters , there is no flux at the inlet, outlet, and wall, which can be described by the Neumann boundaries:
For the species variable , they are limited by mass conservation which leads to either a specification of the concentration boundary conditions with Dirichlet type or a prescription of the mass fluxes with generalized Neumann type. The precursor MTS imported at the inlet is assumed to be a constant (Dirichlet type) and there are no mass fluxes for the other species (Neumann type). Compared with convection, the influence of the diffusion at the outlet can be neglected (Neumann type). On artificial wall, there is no mass flux across the boundaries due to the close system, the symmetric boundary conditions are implemented (Neumann type). Then the boundary conditions for species take the following forms:
Here is the molar concentration of MTS in the imported gas. It can be calculated by the ideal gas law ( is the partial pressure of MTS):
For the variables () in Navier-Stokes equations, the velocity along fiber direction on the inlet is constant while those normal to the fiber direction are zero. The no slip boundary is applied on the gas/solid boundary. On the outlet, the total pressure is fixed as a constant. It is assumed there are no viscous effects hence no boundary layer develops on the wall. Then we have where and is a tangential vector to the boundary.
The initial conditions for variables are straightforward: if th phase exists; for not MTS species; and .
The similar procedure can be applied to the 2D case, except that the bottom wall in 3D case now is occupied by fiber as shown in the section of Figure 3. Therefore the corresponding boundary conditions are and . Notice that it is not necessary to set up boundary condition for the variables in Navier-Stokes equations (), because the actual boundary for these variables come to be the interface with , namely, gas/solid boundary . The parameters related to model are listed in Table 1.
3.2. 2D and 3D Simulation
The numerical simulation is performed based on the software “COMSOL Multiphysics” with the relative tolerance and absolute tolerance . The distributed parallelization calculation is implemented on high performance computing (HPC) clusters “InstitutsCluster” in Karlsruhe Institute of Technology. Each node contains two Quad-core Intel Xeon processors which run at a clock speed of 2.667 GHz and each node has 16 GB of main memory, 4 local disks (250 GB each).
Figure 4 shows the microstructure evolution and corresponding MTS concentration field during SiC deposition. It can be found that the growth of solid phase nearby precursor inlet side is much faster than that nearby outlet side, which implies that the risk of blockage of the premature pores increases. Also we can find that the C (graphite) phase tends to be eliminated by SiC phase because of its overgrowth. This is accordant with experiments [30, 31] that only codeposition of Si(s) and SiC(s) exists while C (graphite) phase is absent in under the high ratios of to MTS (>104) as in this study. As time passes by, concentration of MTS nearby the inlet side is much higher than that of the rest, and this distribution favours the growth of the solid phase nearby inlet side. The concentration of MTS at section A and B at different moment is shown in Figure 5. The step increasing in the curves along section B indicates the location of gas/solid interface, with left part in solid phase and the right part in gas phase. It shows that as time goes on, the along both these two sections decreases; the gradient decreasing of along section A is obvious which is mainly caused by chemical reaction as being explained later; there barely exists change of along section B because the effect of diffusion is weak comparing to that caused by chemical reaction.
By changing the reaction diffusion equation, the effect of fluid flow and chemical reaction on distribution is investigated. Three cases are set up: (a) only diffusion is considered; (b) constant fluid flow is considered in mass transport equation based on case a; (c) chemical reaction in mass transport equation is considered based on case a. Simulation shows that the fluid flow in case b changes MTS concentration field as in case a and tends to reduce the concentration gradient as shown in Figure 6. This homogenization of MTS is in favour of the even deposition to reduce the potential pore blockages which cause the porosity in matrix. After chemical reaction is introduced into the diffusion equation as in case c the concentration nearby inlet side is quickly decreased to small value nearby outlet side of the reactor, which indicates that MTS is highly reactive and is prone to decomposing. This inhomogeneous MTS distribution impedes the even deposition to achieve highly densified production. Because the decomposition of MTS is highly dependent on temperature; the thermogradient CVI (TGCVI) is frequently employed in practice to improve the SiC deposition environment. The phase field modeling TGCVI process will be a part of our further research as scheduled in renewal proposal. From this simulation, it can be found that the flow is in favour of the homogeneous distribution of MTS while this concentration distribution is mainly dominated by chemical reaction.
Under appreciated infiltration conditions, the single phase SiC is deposited. Figure 7(a) shows the growth of multi-SiC grains with different initial seed size. It is found that as deposition proceeds, the SiC grains with small size tend to be eliminated by those surrounding grains by the competitive growth. The coarse SiC grains in final microstructure are accordant with experimental results which show that the tiny grain on the deposition surface is gradually coarsen. The 3D simulation of SiC grains growth with uniform seeds is shown in Figure 7(b). The SiC deposition on SiC fiber surface under the constant fluid is simulated as shown in Figure 8. Herein the fibers are unidirectionally preformed as bundle. The gas flow is along -axis direction, from left to right. The slices of fiber at different moment show clearly how the premature pore blockage comes into being with narrow entrance and wide end.
4. Summary and Conclusion
The multiphase field model is formulated in this paper to predict the microstructures evolution during the CVI process of SiC. By this model the codeposition behavior of the multisolid phases is reproduced. The model consists of phase field equations, mass conservation equations, and the modified Navier-stokes equations. Thereby the complex coupling is accounted among the interprocesses of chemical reaction, deposition, diffusion, and fluid flow. A new form of the system free energy density is constructed based on chemical reaction free energy changing, instead of phenomenological description of the dependence of free energy on the concentration of the gaseous reactants in previous research. Both homogeneous process and heterogeneous process are included in the model through process intensities. The modified Navier-Stokes equations for fluid motion is proposed by using phase field order parameter as coupling variable, which successfully avoids the assumption for solid phase as Newtonian fluids; thereby potential computation load is reduced. Based on the analysis of linear system which arises from the finite element discretization of the model, the parallelized implement is applied. The microstructure simulation shows the potential risk of blockage of the premature pores during CVI process. It can be found that the diffusion effect tends to cause the homogeneous distribution of MTS and meanwhile the flow is in favour of it; however, this concentration distribution of MTS is mainly dominated by chemical reaction. Because the pyrolysis of MTS is highly sensitive to temperature, the thermogradient CVI (TGCVI) is expected to improve the SiC deposition environment aiming at the full densification of SiC preform. The quick decreasing of MTS from inlet side to the outlet side is accordant with the experiment investigation. The competitive growth mechanism is reproduced during SiC grain deposition and result shows that the SiC grains with small size tend to be eliminated by those surrounding grains. Its 3D simulation illustrates the formation process of the premature pore blockage characterized by narrow entrance and wide end.
The financial support by the German Research Foundation (DFG Schn 245/31-1) is gratefully acknowledged.
- H. G. Wulz, U. Trabandt, and T. Schmid, “High temperature ceramic matrix composites iii,” Key Engineering Materials, vol. 164, pp. 1445–1449, 1999.
- J. I. Kim, W. J. Kim, D. J. Choi, and J. Y. Park, “Deposition of compositionally graded SiC/C layers on C–C composites by low pressure chemical vapor deposition,” Journal of Nuclear Materials, vol. 307–311, no. 2, pp. 1084–1087, 2002.
- K. Minato and K. Fukuda, “Structure of chemically vapour deposited silicon carbide for coated fuel particles,” Journal of Materials Science, vol. 23, no. 2, pp. 699–706, 1988.
- J. H. Oh, B. J. Oh, D. J. Choi, G. H. Kim, and H. S. Song, “The effect of input gas ratio on the growth behavior of chemical vapor deposited SiC films,” Journal of Materials Science, vol. 36, no. 7, pp. 1695–1700, 2001.
- D. A. Streitwieser, Kinetic investigation of the Chemical Vapor Infiltration and Reaction (CVI-R) process for the production of SiC and TiC biomorphic ceramics from paper preforms [Ph.D. thesis], University of Erlangen-Nürnberg, 2004.
- F. Kobayashi, K. Ikawa, and K. Iwamoto, “Formation of carbon-excess SiC from pyrolysis of CH3SiCl3,” Journal of Crystal Growth, vol. 28, no. 3, pp. 395–396, 1975.
- A. Josiek and F. Langlais, “Kinetics of CVD of stoichiometric and Si-excess SiC in the system MTS/H2 at medium decomposition of MTS,” Chemical Vapor Deposition, vol. 2, no. 4, pp. 141–146, 1996.
- T. Tago, M. Kawase, Y. Ikuta, and K. Hashimoto, “Numerical simulation of the thermal-gradient chemical vapor infiltration process for production of fiber-reinforced ceramic composite,” Chemical Engineering Science, vol. 56, no. 6, pp. 2161–2170, 2001.
- X. Wei, L. Cheng, L. Zhang, Y. Xu, and Q. Zeng, “Numerical simulation for fabrication of C/SiC composites in isothermal CVI reactor,” Computational Materials Science, vol. 38, no. 2, pp. 245–255, 2006.
- W. G. Zhang and K. J. Huettinger, “Cvd of sic from methyltrichlorosilane. part I: deposition rates,” Chemical Vapor Deposition, vol. 1677, 172 pages, 2001.
- N. Moelans, B. Blanpain, and P. Wollants, “An introduction to phase-field modeling of microstructure evolution,” Computer Coupling of Phase Diagrams and Thermochemistry, vol. 32, no. 2, pp. 268–294, 2008.
- B. J. Palmer and R. G. Gordon, “Local equilibrium model of morphological instabilities in chemical vapor deposition,” Thin Solid Films, vol. 158, no. 2, pp. 313–341, 1988.
- B. J. Palmer and R. G. Gordon, “Kinetic model of morphological instabilities in chemical vapor deposition,” Thin Solid Films, vol. 177, no. 1-2, pp. 141–159, 1989.
- S. Jin, X. Wang, and T. L. Starr, “Model for front evolution with a nonlocal growth rate,” Journal of Materials Research, vol. 14, no. 10, pp. 3829–3832, 1999.
- S. Jin and X. Wang, “Robust numerical simulation of porosity evolution in chemical vapor infiltration. II. Two-dimensional anisotropic fronts,” Journal of Computational Physics, vol. 179, no. 2, pp. 557–577, 2002.
- E. Schnack, F. W. Wang, and A. J. Li, “Phase-field model for the chemical vapor infiltration of silicon carbide,” Journal of the Electrochemical Society, vol. 157, no. 7, pp. D377–D386, 2010.
- A. Ekhlakov, S. Dimitrov, T. A. Langhoff, and E. Schnack, “Phase-field model for deposition of pyrolytic carbon,” Communications in Numerical Methods in Engineering, vol. 24, no. 12, pp. 2139–2154, 2008.
- E. Schnack, F. W. Wang, T. A. Langhoff, and A. J. Li, “Modeling and simulation of composites in the design process,” in Proceedings of the 7th International Symposium on Tools and Methods in Competitive Engineering, 2008.
- W. G. Zhang and G. K. Huttinger, “Chemical vapor deposition of sic from methyltrichlorosilane part ii: Composition of the gas phase and the deposit,” Advanced Materials, vol. 7, pp. 173–181, 2001.
- M. D. Allendorf and C. F. Melius, “Theoretical study of thermochemistry of molecules in the silicon-carbon-chlorine-hydrogn system,” Journal of Physical Chemistry, vol. 97, no. 3, pp. 720–728, 1993.
- E. Fitzer and D. Kehr, “Carbon, carbide and silicide coatings,” Thin Solid Films, vol. 39, no. C, pp. 55–67, 1976.
- M. W. Chase, J. L. Curnutt, R. A. McDonald, and A. N. Syverud, “Nist-janaf thermochemical talbles,” Journal of Physical and Chemical Reference Data, vol. 9, 1998.
- L. V. Gurvich, I. V. Veyts, and C. B. Alcock, Thermodynamic Properties of Individual Substances, Hemisphere Publishing, 1989.
- J. J. Eggleston, G. B. McFadden, and P. W. Voorhees, “A phase-field model for highly anisotropic interfacial energy,” Physica D, vol. 150, no. 1-2, pp. 91–103, 2001.
- I. Steinbach, F. Pezzolla, B. Nestler et al., “A phase field concept for multiphase systems,” Physica D, vol. 94, no. 3, pp. 135–147, 1996.
- R. H. Perry, D. W. Green, and J. O. Maloney, Perry's Chemical Engineers' Handbook, McGraw-Hill, 7th edition, 1997.
- C. Beckermann, H. J. Diepers, I. Steinbach, A. Karma, and X. Tong, “Modeling melt convection in phase-field simulations of solidification,” Journal of Computational Physics, vol. 154, no. 2, pp. 468–496, 1999.
- V. Arnautu and C. Morosanu, “Numerical approximation for the phase field transition system,” International Journal of Computer Mathematics, vol. 62, pp. 209–221, 1995.
- M. R. Dorr, J. L. Fattebert, M. E. Wickett, J. F. Belak, and P. E. A. Turchi, “A numerical algorithm for the solution of a phase-field model of polycrystalline materials,” Journal of Computational Physics, vol. 229, no. 3, pp. 626–641, 2010.
- F. Loumagne, F. Langlais, and R. Naslain, “Experimental kinetic study of the chemical vapour deposition of SiC-based ceramics from CH3SiCl3H2 gas precursor,” Journal of Crystal Growth, vol. 155, no. 3-4, pp. 198–204, 1995.
- G. D. Papasouliotis and S. V. Sotirchos, “Experimental study of atmospheric pressure chemical vapor deposition of silicon carbide from methyltrichlorosilane,” Journal of Materials Research, vol. 14, no. 8, pp. 3397–3409, 1999.