← Back to Portfolio

Microwave Plasma Deposition Resonator

Consulting & Reimplementation Available: This model was originally implemented in PLASIMO during my PhD. We can discuss how to implement it in another plasma simulation software upon request, or adapt your own model, add new chemistry, or create verification workflows to test it. Contact me to discuss.

Disclaimer: PLASIMO and COMSOL Multiphysics® are not affiliated with Deep Why and require the client to hold their own licenses.

Microwave Induced Plasmas (MIPs) are widely used in industrial applications such as the Plasma-activated Chemical Vapour Deposition (PCVD) process, which is the foundational step in the fabrication of high-quality optical fibers (such as those produced at Draka).

During PCVD, a moving microwave resonator surrounds a quartz tube, creating a non-local thermal equilibrium plasma. By feeding gases like $\mathrm{SiCl_4}$ and $\mathrm{O_2}$ into the tube, thin glass layers are deposited on the inner wall, forming the core refractive index profile of the preform. To prevent electromagnetic radiation leakage to the environment and maximize power coupling efficiency, metal chokes are designed into the cavity walls.

Optical Fiber Fabrication Steps
Figure 1: The three steps of optical fiber production: 1) Deposition (PCVD), 2) Collapsing of the hollow quartz tube into a preform rod, and 3) Stress-guided fiber drawing in a furnace tower.

Computational Domain & Resonator Design

A self-consistent 2D axisymmetric multiphysics model was built using the PLASIMO simulation platform to analyze the interaction between the electromagnetics, transport processes, and argon chemistry inside the PCVD reactor.

The electromagnetic fields are solved as TM modes. Discontinuities in the cavity walls (such as quarter-wave short-circuit chokes) are modeled to maximize reflections and prevent waves from leaking out of the reactor region into the environment.

2D Axisymmetric Computational Domain
Figure 2: 2D axisymmetric computational domain showing the dielectric quartz tube, air gaps, reactor metallic boundary, and the short-circuit chokes placed at the boundaries.

2D Spatial Plasma Distributions

The model self-consistently computes the 2D spatial distribution of the main plasma quantities: the electron density ($n_e$), electron temperature ($T_e$), and gas temperature ($T_h$). These profiles explain how the power coupled through the waveguide gap is absorbed and transported axially.

2D Electron Density distribution
Figure 3: 2D spatial distribution of the electron density ($n_e$) showing strong axial confinement within the reactor region.
2D Electron Temperature distribution
2D electron temperature ($T_e$) distribution.
2D Gas Temperature distribution
2D heavy-particle gas temperature ($T_h$) distribution.

Confinement & Choke Optimization

The position and depth ($d_{ch}$) of the chokes play a critical role in wave reflection. The resonance coupling of the microwave power into the plasma shifts depending on the choke depth because the plasma is a highly dispersive medium. The figure below shows the reflected/absorbed power as a function of the choke depth, highlighting the confinement optimum.

Power Absorption vs Choke Depth
Figure 4: Influence of the quarter-wave short-circuit choke depth ($d_{ch}$) on the power absorbed in the plasma, showing distinct resonance modes.

Radial Densities

Detailed radial species distribution profiles show how charged particles (electrons, atomic and molecular argon ions) and excited states behave inside the quartz boundaries, which dictates the radial plasma conductivity and chemical reactivity.

Radial Density Profiles
Figure 5: Radial distribution of heavy species and charged particle densities inside the quartz tube.

Bibliography

  • Jimenez-Diaz M. (2011). Modelling of microwave induced plasmas: the interplay between electromagnetism, plasma chemistry and transport. PhD Thesis, Eindhoven University of Technology (TU/e). [ read thesis pdf ]