Numerical modelling of microwave plasma-assisted chemical vapour deposition reactors for diamond synthesis
Publication Date
July 24, 2020
Creator
Abstract
This thesis studies the numerical modelling of microwave plasma-assisted chemical vapour deposition (MPA-CVD) reactors for synthetic diamond synthesis. Specifically we model the formation of a plasma within these reactors, which is responsible for enabling the deposition of carbon onto a substrate for diamond growth. This involves first introducing a chemically reacting binary gas flow model for a hydrogen gas mixture, which in turn is coupled to an applied electric field operating at resonance. This electric field acts as the energy source for the ionising collisions that characterise plasma formation.
Following this, we introduce numerical schemes for general hyperbolic and elliptic operators based on employing the discontinuous Galerkin finite element method (DGFEM). Notably, in the case of Maxwell's equations, we employ a non-stabilised formulation of the symmetric interior penalty method. These schemes form the basis of the underlying numerical approximation of the MPA-CVD model.
A large portion of this thesis is dedicated to improving the predictive capabilities of the existing MPA-CVD model. Noteworthy changes include a more physically-accurate description of ionisation processes, the introduction of a secondary heat transfer mechanism, namely, enthalpic diffusion, and the use of an electron energy conservation equation to better approximate the average electron temperature. Additionally, the success of the model is demonstrated by its application to three existing MPA-CVD reactor geometries.
The majority of MPA-CVD reactors are largely axisymmetric, and thus, the majority of MPA-CVD simulation in the literature is performed in 2D. This disregards the possibility of non-axisymmetric perturbations to the plasma distribution. Here, we perform a 3D Maxwell eigenvalue analysis of two MPA-CVD reactor geometries to investigate the possibility of the electric field distribution being polluted by neighbouring non-axisymmetric resonant modes that operate close to the applied frequency. Additionally, we uncover a curious occurrence whereby a spurious eigenvalue spectrum is obtained when employing a higher order boundary representation in conjunction with a standard DGFEM for Maxwell's equations; as far as we are aware, this is not covered in the literature, but was overcome by application of a Piola covariant transformation.
Item Type
ethesis
Thesis Type
PhD
Supervisors
Keywords
Subjects (LC)
Associated Schools / Departments
School of Mathematical Sciences (UK)
eprints ID
60919
UoN Repository URI
Except where otherwise noted, this item's license is described as
File(s)![Thumbnail Image]()
Name
thesis.pdf
Type
Full-text
Description
Examined
Size
21.2 MB
Format
Adobe PDF
Checksum (MD5)
550461c222d077b8441b6ece65f62756