1. Introduction
Composite materials are an important and oftentimes used type of engineering materials with tuneable properties which depend on the materials combined. They are usually made of two or more materials. The first phase, commonly referred to as “reinforcement” is constituted by a solid material, which can be employed in several morphologies such as fibers (short or continuous), particles, or as a co-continuous domain of a solid phase which occupies a fraction of the available volume (foams) [1-7]. The second phase, normally called “matrix,” can be gaseous, liquid, or solid, and is in contact with the reinforcement through a proper interface. Composite materials are very common in nature, where a combination of fibers and matrix is often found. A good example is wood, which is composed of flexible fibers of cellulose embedded in a rigid matrix of lignin. Artificially produced composites are usually classified, according to the nature of their matrix, into metal, ceramic, or polymer matrix composites. Each type of composite presents its advantages and disadvantages, and together they cover a great scope of applications [1, 8-12].
Composite materials manufacturing can be accomplished through several techniques. Each technique is mainly the result of the trade-off between the target properties of the final product, the starting materials properties, and the final cost. In most of the cases, the matrix is in a liquid state during processing. A common action in composite materials manufacturing is “to infiltrate” the reinforcement with the matrix. When a liquid is used, the impregnation of the remaining free volume is highly dependent on how the first phase is placed in the space and on the physical properties of the liquid itself. In the case of polymeric matrix composites (PMC) techniques involving a liquid matrix are the only way to process them, as polymers are in a liquid state before their cross-linking in thermosets, or can be melted if they are thermoplastics. Also, ceramic matrix composites are processed through liquid routes (e.g., sol-gel techniques for oxide composites, or polymer infiltration and pyrolysis (PIP) [13-14] for ceramic matrix composites (CMC) produced by pre-ceramic precursors. In the latter case, where several impregnations are performed, the permeability of the porous materials decreases with the increasing of the final density of the composite under production). A widely utilized technique which stands among ceramic and metallic matrix composites (MMC) manufacturing is the so-called reactive silicon infiltration [15]. Its peculiarity is that, during processing, silicon is molten and by reacting with extra carbon present in the primary phase it yields silicon carbide. This method was not originally conceived to produce composite materials, but for the synthesis of bulk SiC, still, it proved to be a suitable methodology for ceramic matrix composites manufacture as well.
Some industrial products obtained with this infiltration technology are pantographs, connecting rods, loudspeakers, piston rings, heat sinks, electro-brushes, engine cylinders, brake discs, etc. In most of these examples, primary phase (reinforcement) content is above 50% in volume. For reinforcement contents of such magnitude, the most commonly used technique for manufacturing composites is infiltration [1, 5-7, 16-18].
Actually, infiltration processes are not exclusive of materials manufacture, but can be observed in many everyday life activities, such as making coffee, where hot water is infiltrated into a porous bed of coffee particles. Perhaps one of the most important cases where infiltration phenomenon is observed is in the penetration of rainwater into the ground, which ends up filling aquifers that guarantee water supply to population in dry regions where the surface runoff is scarce. Because of the great importance of infiltration in geological sciences, it was in this field where the deepest and most detailed studies of this phenomenon were performed. Many advances in this area of knowledge have been expanded and applied to other scientific fields such as materials science. The clearest example is Darcy’s law which is widely used in the design of composite materials. This empirical law was obtained by Henry Darcy whose original work published in 1856 is entitled: “Les Fontaines Publiques de la Ville de Dijon” (The Public Fountains of the City of Dijon). However, Darcy's law is not the only example that can be found in the literature. The relationship between water saturation and organic compounds in soils of different origins has resulted in the semiempirical equations of Brooks and Corey [19] and Van Genuchten [20], both widely used in composite materials.
The aim of the present chapter is to introduce the reader to the basic concepts of wetting and infiltration, and to combine these concepts to solve the equations governing the infiltration process in several particular cases that often arise in the composite materials manufacturing industry.
2. Experimental
Although it will be discussed in detail in the upcoming sections of the chapter, it is important to clarify that infiltration, in most composite materials manufacturing processes, does not occur spontaneously (in the case of PMCs it does proceed spontaneously, but it is so slow that it is necessary to force it), so the liquid must be forced to infiltrate by applying an external pressure. From an engineering point of view, mainly two technologies have been developed for this purpose: the so-called gas-assisted pressure infiltration (GPI), where a gas is used to push the liquid; and another called squeeze-casting (SC) in which pressure is applied by means of a mechanical device. Each technique has its upsides and downsides, which will be discussed throughout this chapter. The main advantage of GPI over SC is that it allows a better control of physical parameters (pressure and temperature), and admits higher complexity in the manufactured parts. On the other hand, the main advantage of SC over GPI is that higher pressures can be achieved and the automatization of the process becomes easier.
First the GPI technique will be analyzed, and for this purpose, the example of a laboratory equipment developed at the University of Alicante will be used [10, 11, 21-29]. A scheme of the equipment is shown in Figure 1. The equipment consists of a pressure chamber (maximum admissible pressure at 300 °C of 5 MPa) equipped with a resistance furnace, which is thermally isolated from the chamber by the use of conventional refractory material (temperature of chamber wall is always below 100°C). The furnace has a power of 1500 W and can reach temperatures of 1000°C. For measurement and control of furnace temperature, K type thermocouples, connected to a temperature controller, are used. The pressure control is entirely performed by a system of pneumatic and electronic valves connected to a pressure controller. The pressure vessel is made entirely of stainless steel and consists of a cylindrical body with an internal diameter of 30 cm, an outer diameter of 34.6 cm, a height of 45 cm, and a wall thickness of 3 cm. The closing of the camera is made with 12 M16 mounting studs. Both, the electrical feed (to thermocouple input and power cables) and the gas inlet and outlet for pressurizing the chamber are positioned in the body. Additionally, the camera is connected to a primary vacuum pump, which allows the process to be performed in the absence of oxygen (low partial pressure).
The manufacturing process of a composite material using GPI can be described as follows: a preform, which consists of a porous bed of compacted particles, is placed into a quartz (or BN- coated graphite) crucible and a metal ingot is placed on top of it (see Figure 1b). The crucible is then introduced inside the pressure chamber and vacuum is applied during at least 30 min., until a pressure of about 100 Pa is reached. The chamber is heated under vacuum at a very low rate (3°C/min) up to 250°C, to remove humidity and gas adsorbed on particles. Subsequently, further heating is applied up to a temperature 50°C higher than the melting point of the metal at a rate of approximately 5°C/min. The maximum temperature is then maintained for 30 min followed by argon injection to reach the desired infiltration pressure. Once the preform has been infiltrated (approximately 2 min. under pressure) the crucible is pushed towards the bottom of the chamber, which is cold, and directional solidification occurs, ensuring a high-quality composite (null porosity). Figure 2 shows a composite material obtained at the University of Alicante with this process, where it is possible to observe the absence of porosity and an intimate contact between reinforcement and matrix.
At the University of Alicante, another equipment, in this case for squeeze-casting technique [30], was also developed. Figure 3 shows a schematic representation of the main body of this equipment and how it operates. The main feature of this apparatus is that it can work up to 450°C and at a maximum pressure of 150 MPa. To resist these operating conditions, special (W720 Bölher) steel was used to build the body of the equipment. Embedded resistances in the central body heat the system, and temperature is measured and controlled by K type thermocouples connected to a temperature controller.
The manufacturing process of a composite material using SC can be described as follows: The system is first preheated to the operating temperature (between 350 and 450 °C depending on the alloy) to avoid early solidification of the metal matrix during infiltration. The preform is preheated (600 to 700°C) and the alloy is melted (150°C above the melting point) in two auxiliary furnaces. Thereafter, the preform is transferred to the main chamber, and the liquid metal is added. The desired pressure is applied as indicated in Figure 3. After infiltration, the system is cooled down and the manufactured material is removed as indicated in Figure 3.
The main problem with this technique is that the processing is complex because the times and temperatures must be very accurate to keep the metal from solidifying before penetrating into the preform. This disadvantage is also an advantage, since the high velocity of the process generates a finer microstructure and thus improves the mechanical properties of the final material.
Most importantly, it is possible to find in the market, both SC and GPI fully automated devices, but the basic operating principles are the same for the ones developed at University of Alicante with a scale factor between 2 and 10.
3. Wetting basic laws
One of the most relevant characteristics for the production of MMC through liquid processing routes is the wetting of the reinforcement by the liquid metals, which is quantified by measuring the contact angle between the reinforcement and matrix [31]. The contact angle between a flat substrate and a drop of liquid metal can be routinely measured by the method of the sessile drop. However, this measurement is not possible when the solid is in particulate form, in which case the wetting can be studied by means of liquid metal infiltration into compacted particulate samples [10, 11, 21-29].
Considering a solid, flat, non-deformable and chemically homogeneous surface in contact with a nonreactive liquid in the presence of a gaseous phase, if the liquid does not completely cover the solid, the liquid will form with the solid a contact angle θ, as depicted in Figure 4.
The equilibrium value of θ, which is used to define the wetting behavior of a fluid, can be written in terms of the balance of three energies present in the interfacial vapor-liquid-solid system, according to the well-known Young equation (1805):
where σ_{ij} is the interfacial energy between the phases i and j (s, l, and v represent the solid, liquid, and vapor phases, respectively). A contact angle of less than 90° identifies a liquid that wets the substrate, while larger angles identify liquids that do not wet it. An ideal wetting liquid will have an assigned value of θ equal to 0°.
The interfacial energy σ_{lv}, also called surface tension (
If now this equation is inserted in the Young’s equation, the so-called Young-Dupré equation is obtained:
All infiltration processes are driven by a differential pressure that allows the liquid matrix to penetrate into the porous reinforcement. This difference of pressure between the two fluids (liquid matrix and gas phase filling porosity) is defined as the capillary pressure (
Where
The energy shift that occurs when immersing the capillary into a liquid can be described in terms of work as:
where
where
Since the measurement of the surface area of the particles becomes rather inaccurate when the particle diameter surpasses a certain value, and the surface roughness is very difficult to estimate, it is possible to assume a particular geometry for the porous media, e.g., cylindrical pores with a radius r as an approximation, in which case the surface area can be written as
Where λ is a geometrical factor introduced to describe deviation from the assumed geometry of the pores and the surface roughness of the material (it usually oscillates between 2 and 4 [10, 21-23, 25, 29],
4. The simplest case: Obtaining Darcy’s law
As previously mentioned, Darcy’s law predicts how a fluid moves through a porous system. This empirical law is based on the study of water movement in sand beds, and represents an ideal model in the manufacture of a composite material. Darcy’s law can be derived from the resolution of the Navier-Stokes equation, considering the ideality of the system. In this section, the basis for describing the process of manufacturing of composite materials will be provided, and the most important factors that affect the process will be identified.
As aforementioned, the flow of viscous liquids is governed by the Navier-Stokes equation. This differential equation does not have a general solution. However, in some cases, it is possible to solve, e.g., for an incompressible fluid, in laminar flow (low Reynolds number), assuming that the time derivative of the fluid velocity is much smaller than its spatial derivatives, and that the effect of gravity is negligible. In this case, the equation can be reduced to:
where
where z is the flow direction,
Combining both equations and integrating, the Darcy’s law, which shows the relationship between infiltration distance (h), time (t) and the pressure drop in the infiltrating liquid (ΔP=P-Po) is obtained.
It must be pointed out that the permeability of the porous media is proportional to the square of the average pore radius:
Darcy’s law has been widely used in basic studies of manufacture of metal matrix composites. Despite its simplicity, it has proven to be very useful in metallurgical engineering. As an example, the results of infiltration of liquid aluminium into a porous bed of particles, obtained by Narciso et al., are shown in Figure 6 [10, 25]. Figure 6a provides the relationship between infiltration distance and applied pressure for liquid aluminium into a porous bed of alumina particles, while Figure 6b shows the results for a liquid infiltration of an aluminium alloy into preforms of graphite particles. As predicted by Darcy’s law, the square of infiltration distance is proportional to the pressure drop (h^{2} vs ΔP) in these systems.
As discussed in the preceding paragraphs, to solve the Navier-Stokes equation a full saturation and a laminar flow (Re < 2100) have been assumed. Besides, the effect of gravity has been considered to be negligible. Figure 6b shows a good correlation of the experimental data with Darcy’s law, and thus, the assumed conditions should be valid. Reynolds number (Re) can be obtained from the following expression:
Where
And considering that the liquid column is always less than 10 cm, in the case of aluminium for example, the pressure would be about 2 kPa, which is almost 3 orders of magnitude lower than the capillary forces, and therefore negligible.
In the following sections different real situations that can occur are considered and will be analyzed separately. The combination of all these particular cases simultaneously, although more realistic, is much more complex and goes beyond the scope of this chapter.
5. The effect of saturation
On systems where the liquid does not wet the solid, which is most common in MMCs, it must be taken into account that the permeability of the porous system is not constant [37-39] and that it depends on the saturation and this latter, and at the same time, on the applied pressure [21, 40-42].
The study of unsaturated flow is well known in geology and in 1931 it was described mathematically by Richards [34], whose partial differential equations are difficult to solve and usually solved numerically for each particular case. Numerous researchers have worked on this particular hydrodynamics and mathematical problem [44,45], and the books by Prof. Bear collect these study cases [46-48], which are essential reading for researchers working in this field. It is also noteworthy that Prof. Bear highlights in his books the goodness of Darcy’s law, and that its application range is wider (more useful) than it was expected at first (i.e., unsaturated flow, high Reynolds number, etc.).
These studies were not limited to the case study of hydrology, but have crossed over to the field of metal matrix composites (MMC), performed especially by the group of Andreas Mortensen [40-42], but as mentioned above without analytical solution of the differential equations describing the phenomenon. The highlight to be found in the recent literature are two articles that show that close to the threshold pressure, infiltration is actually governed by percolation phenomena (instead of slug flow) [49,50].
Another aspect to consider is that pressure also depends on the position along the flow path, and is usually assumed that the pressure drop between the front (P_{0}) and head (applied pressure) is linear, which is reasonably sane. Nevertheless, in order for the differential equations to have an analytic solution, it will be assumed in this chapter that practically all the preform is at the same pressure, which in principle is something reasonable as working pressures oscillate between 1 and 5% of the threshold pressure (i.e., head=1.05P_{0}; front=P_{0}).
In the literature, several models showing the relationship between pressure and saturation have been described, but the most famous are the Brooks-Corey (BC) [19] and van Genuchten (VG) [20] models. In this section, the Navier-Stokes equation has been solved considering that the permeability is not constant, relating it to the saturation through the BC and VG models. Different authors have used these models to describe the infiltration process in MMCs [51]. However, previous works have shown that for low pressures (close to the threshold pressure) these models, especially the BC, do not describe the infiltration process correctly [52, 53].
Starting from reduced Navier-Stokes equation:
Where z is the infiltration distance, K the permeability, t the time of infiltration, η the dynamic viscosity of the fluid,
Relative permeability, depends on the degree of saturation according to:
Where n is a parameter that ranges from 1 to 3. In the case of composites fabricated by infiltration, it has been found that 1 is the most accurate value for n [51]. There are several models that relate saturation to the applied pressure (P) and the threshold pressure (P_{0}) of the system. One of the most used models of this kind is the Brooks-Corey equation:
Where λ is a parameter that depends on the solid material. Replacing K in eq. 16 with eqs. 17, 18, and 19 yields a modified reduced Navier-Stokes equation:
Solving this equation provides a model to estimate depth of infiltration as a function of applied pressure. In this form, there is no general analytical solution to the equation, so, to obtain an analytical model, a particular case must be considered. As discussed above, experimental data fit better when
Where A is the integration constant, and its value can be determined using initial conditions, i.e., for P = P_{0} and h = 0
Then the final expression for h(P) is:
It must be noted that if
If, instead of using the Brooks-Corey equation, the Van Genuchten model is applied:
where
As in the case of the Brooks-Corey model, a general analytical solution is not available, so the case of n = 1 will be considered again:
.The second integral does not have an analytical solution either, so in order to obtain a solution as general as possible, the value of n1=2 is introduced [42, 51]:
Applying the substitution
Replacing back
To estimate the value of C, initial conditions must be applied
Then, the general solution for n1 = 2 is:
In order to analyze the accuracy of the results, the evolution of the square of the infiltrated height as a function of pressure drop is plotted in Figures 7a and 7b. As shown, experimental results are perfectly framed between the line defined by Darcy’s law (assuming that the liquid permeability is invariant and therefore permeability measured with polyethylene glycol is used) and the curve modified with the Van Genuchten equation.
The results shown here clearly demonstrate that it is not necessary to solve the Richards equation (modified Darcy’s law to take into account saturation), which has only numerical solutions to predict infiltration in composite materials unlike what was suggested by Mortensen et al. [54]. It is noteworthy that many authors in soil mechanics indicate that Darcy’s law is more universal than it was originally supposed.
6. The effect of variable radius
Although the variation of the radius is not a very common phenomenon in the manufacture of composite materials, it is a process that can occur. The variation of the radius is a fundamental issue for example in the manufacture of SiC by the RBSC process or reactive infiltration and SiC derived CMCs. This is a rather complex situation since the variation of the radius affects the permeability of the porous media.
To solve the Navier-Stokes equation it has been assumed that the reduction of the radius, caused by a chemical reaction, is controlled by diffusion. The equations presented so far are valid for infiltration in nonreactive systems. However, they are not valid for reactive infiltration because the pore radius, in this kind of systems, is not constant but decreases with the generation of new solid. As discussed above, the threshold pressure and the permeability depend on the radius, so they are not invariant over time.
So for reactive systems, it is necessary to re-solve the reduced Navier-Stokes equation considering the variation of the radius:
After integration, the following expression is obtained:
Usually in reactive systems, infiltration takes place spontaneously (without applying an external pressure), so
Combining this expression with eq. 6 yields:
Simplifying and rearranging:
In principle, porosity of the system is not constant, but for a first approximation it can be considered so. To proceed with the solution of the equation, a relationship between radius and time is necessary. Assuming that the controlling step of the reaction is the diffusion of C through the SiC layer, then the model of decreasing nucleus developed in the decade of the 1950s can be used [55, 56]:
Where
Introducing these relations into eq. 46 and defining:
Results in the following expression:
Solving the first integral yields:
Where A is an integration constant. Solving the second integral results in:
where B is an integration constant. In order to obtain the final solution, boundary conditions must be used to find the values of A and B:
Where
By applying condition a), the next expression is obtained:
Applying condition b) leads to:
and
Solving eq. 55 leads to:
Replacing the value of A in eq. 53 results in:
Replacing the value of B in eq. 54 provides the expression for h_{MAX}:
And after replacing the values of A and B, and rearranging the final expression h(r) is obtained:
It must be taken into consideration that this equation is valid only as long as the modification of the porosity is negligible. For evaluating h(t), eq. 59 and eq. 47 must be combined.
Figure 8 shows the calculated infiltration distance with time for a carbon preform with an initial pore radius of 30 micron and 50% open porosity, infiltrated with pure Si, using Darcy’s law and the variable radius model for cylindrical and spherical pores. The variable radius models predict a slower rate of infiltration, which is closer to the experimental data obtained by several authors [12, 57-61].
Another approach is to consider that the pore closure is limited by the chemical reaction that originates it. In this case, the equations of the decreasing nucleus model to be used are:
Introducing this equation into the reduced Navier-Stokes equation and solving the integrals results in:
The integration constants A and B are calculated using contour conditions:
For condition b):
(65) |
and replacing in condition a):
and
Replacing the values of A, B and h_{MAX}:
Figure 9 shows the square of the infiltration distance as a function of time for a carbon preform with an initial pore radius of 30 micron and 50% open porosity infiltrated with pure Si, using Darcy’s law and the variable radius model for a chemically controlled pore closure. As seen, infiltration rate decreases with time, just as in the previous case (diffusion controlled pore closure). However, the experimental data [12, 42, 59-66] show that the infiltration rate is lower than the one predicted by these models; furthermore, infiltrated height is proportional to t instead of to t^{(1/2)}. Experiments indicate that infiltration is controlled by the reaction at the triple line and not by the viscous phenomena described here. However, both interpretations are compatible. It is probable that in the first stage of the process, infiltration is controlled by the reaction at the triple line, and after the SiC interface has been formed, the viscous flow becomes the controlling process. A completely different approach to the process of infiltration is performed by Sergi et al. [67, 68], in which they use Lattice-Bolztman models to describe infiltration in a single capillary and predict that infiltration distance is proportional to time when considering pore narrowing of the capillary due to the reaction. However, their predictions have not been validated experimentally up to date. It is possible to conclude that a deeper understanding of the reactive infiltration phenomenon should be gained in order to discern what is happening in the overall process.
7. Effect of variable viscosity: Curing polymeric resins
The properties of PMCs depend strongly on the type of matrix used. Among the possible binders, thermoset resins present a singularity: their viscosity is time dependent [69, 70]. These kinds of polymers present an initial low viscosity that grows exponentially with time. During this process, the material undergoes a series of chemical reactions that cause the cross-linking of polymeric chains yielding a 3D structure with excellent thermal stability and solvent resistance. The process of cross-linking that takes place in these resins is known as curing. In polymer matrix composites in which the matrix cures with time, thus hardening and losing fluidity, viscosity becomes a limiting factor in infiltration.
The evolution of viscosity is usually expressed with the following equation:
Where
After integration, the following expression is obtained:
Introducing eq. 71 and solving the integral yields a model for predicting infiltration distance as a function of time for curing polymeric resins:
To calculate the values of integration constants A and B, initial conditions are applied:
For condition a):
for condition b):
Then the final equation results in:
It must be highlighted that the negative sign is only indicative of the direction of flow. Figure 10 shows the evolution of the square of infiltrated height for a preform infiltrated with curing resins with an initial viscosity of 100 Pa s at different temperatures. The drastic effect of temperature on infiltration distance becomes evident, since increasing temperature to 50°C reduces infiltrated distance by one order of magnitude.
8. Concluding remarks
The present chapter has introduced the reader to the processes involved in the manufacture of different types of composite materials (PMCs, MMCs, CMCs) via infiltration. Over 90% of the world production of composite materials is carried out by infiltration, especially in the case of PMCs where it is almost the only available processing technique, while in case of metal matrix, production via infiltration represents only 70% [7]. In order to accurately control the properties of the final material and to be able to optimize the production, it is key to identify and understand the different processes involved in the manufacture of these materials. For this purpose, in this chapter, different models have been elucidated, starting from the differential equation of Navier-Stokes and taking into account the basic laws of wetting, by applying different real conditions often found in industry. For PMCs, the main discussed issue was the change in viscosity that thermoset matrixes experiment during infiltration. By introducing this phenomenon in the infiltration equations, it becomes clear that it can be a limiting factor during production since it shows a great deviation from Darcy’s law. In the case of MMCs, the main phenomenon analyzed was the effect of saturation on infiltration. It became clear that variable saturation not only affects in a huge manner infiltration kinetics, but also changes the dominating phenomena. The last publications in this area have identified percolation as probably the leading phenomenon on infiltration in unsaturated systems at very low overpressure. Finally, for CMCs, the role of pore narrowing, an issue that has been widely discussed in the literature over the past 10 years was revised. This phenomenon is especially relevant in processing of ceramics and ceramic matrix composites via reactive techniques like the RBSC method used for SiC production. Even though the research on this area is vast, there is still not a full understanding of the effect of pore reduction on infiltration. Most models proposed still yield higher infiltration rates than the ones observed experimentally, because they do not consider the effect of chemical reactivity.
Appendix
Table 1 shows the main characteristic of the carbon particles used in this work.
Particle | D(90) | D(50) | D(10) | D(4,3) | D(3,2) | Span | V_{p} |
P1 | 25.8 | 13.4 | 7.0 | 15.1 | 11.6 | 1.39 | 0.511 |
P2 | 41.0 | 25.8 | 15.5 | 27.2 | 23.0 | 0.99 | 0.515 |
P3 | 95.1 | 62.0 | 36.8 | 64.0 | 53.0 | 0.94 | 0.518 |
P4 | 173.6 | 117.6 | 56.6 | 124.1 | 110.4 | 0.99 | 0.482 |
Table 2 shows the value of the parameters required to assess h(t) for the several cases exposed in this work. All parameters are in SI units.
Symbol | Magnitude | Value | SI units | Ref. |
λ | Geometrical factor | 2-4 | - | [10, 21-23,25, 29] |
a | Proportionality constant | 4·10^{-4} | [61] | |
γ_{lv} | Liquid-vapor surface tension of Si (at melting point) | 0.750 | N/m | [34-36] |
θ | Contact angle of Si/C (at Si melting point) | 40 | ° | [31] |
η | Dynamic viscosity of molten Si (at melting point) | 0.605·10^{-3} | Pa·s | [34] |
ρ_{SIC} | Density of SiC | 3210 | kg/m^{3} | [9, 12] |
M_{SiC} | Molar weight of SiC | 40·10^{-3} | kg/mol | [9, 12] |
M_{Si} | Molar weight of Si | 28·10^{-3} | kg/mol | [9, 12] |
| Diffusion coefficient of C through SiC | 2·10^{-14} | m^{2}/s | [62] |
ρ_{SI} | Density of Silicon (at melting point) | 2570 | kg/m^{3} | [65] |
k_{s} | kinetic constant in the SiC formation reaction | 4·10^{-8} | | [71, 72] |