Fabrication of GaN Films in a Chemical Vapor Deposition Reactor

Author(s):  
J. Meng ◽  
Y. Jaluria ◽  
S. Wong

A three-dimensional numerical study has been carried out on the rotating disk GaN MOCVD process, and it is also coupled with an experimental study on the flow and thermal transport processes in the system. An impingement type reactor, with a rotating base, is considered. The dependence of the thin film growth rate and uniformity on operating conditions such as inflow velocity, rotational speed, and susceptor temperature are investigated in detail. Similarly, the effect of the geometry and configuration of the reactor are studied. The study also considers the effect of thermal and solutal buoyancy on the resulting flow. The flow and the associated transport processes are discussed in detail on the basis of the results obtained to suggest approaches to improve the uniformity of the film, minimize fluid loss and reduce flow recirculation that could affect growth rate and uniformity.

2008 ◽  
Vol 15 (01n02) ◽  
pp. 111-116 ◽  
Author(s):  
JAE-SANG BAEK ◽  
JIN-HYO BOO ◽  
YOUN-JEA KIM

A numerical study is needed to gain insight into the growth mechanism and improve the reactor design or optimize the deposition condition in chemical vapor deposition (CVD). In this study, we have performed a numerical analysis of the deposition of gallium arsenide ( GaAs ) from trimethyl gallium (TMG) and arsine in a vertical CVD reactor. The effects of operating parameters, such as the rotation velocity of susceptor, inlet velocity, and inlet TMG fraction, are investigated and presented. The three-dimensional model which is used in this investigation includes complete coupling between the thermal-fluid transport and species transport with chemical reaction.


Author(s):  
J. Meng ◽  
S. Wong ◽  
Y. Jaluria

A numerical study has been carried out on the metalorganic chemical vapor deposition (MOCVD) process for the fabrication of gallium nitride (GaN) thin films, which range from a few nanometers to micrometers in thickness. The numerical study is also coupled with an experimental study on the flow and thermal transport processes in the system. Of particular interest in this study is the dependence of the growth rate of GaN and of the uniformity of the film on the flow, resulting from the choice of various design and operating parameters involved in the MOCVD process. Based on an impingement type rotating-disk reactor, three-dimensional simulations have been preformed to indicate the deposition rate increases with reactor pressure, inlet velocity, and wafer rotating speed, while decreases with the precursor concentration ratio. Additionally, a better film uniformity is caused by reducing the reactor pressure, inlet velocity and wafer rotating speed, and increasing precursor concentration ratio. With the impact of wafer temperature included in this study as well, these results are expected to provide a quantitative basis for the prediction, design, and optimization of the process for the fabrication of GaN devices. The flow and the associated transport processes are discussed in detail on the basis of the results obtained to suggest approaches to improve the uniformity of thin film, minimize fluid loss, and reduce flow recirculation that could affect growth rate and uniformity.


Author(s):  
Jiandong Meng ◽  
Yogesh Jaluria

A numerical study has been carried out to characterize the metalorganic chemical vapor deposition (MOCVD) growth of Gallium Nitride (GaN) in a rotating-disk reactor. The major objective of this work is to examine the dependence of the growth rate and thin film uniformity on the primary parameters. First of all, for a rotating-disk system, the governing equations involved are obtained. Then, with the effect of thermal buoyancy included and based on the detailed mathematical model and chemical reaction mechanisms, the 3D simulation study is conducted for a rotating reactor. A comparison between the predicted growth rate and experimental data is presented. In addition, the effect of various primary operating and design parameters on the growth rate of GaN and thin-film uniformity is also examined. This provides further insight into the reactor performance and the characteristics of the entire process. The results obtained can also form the basis for the future design and optimization of this system.


Author(s):  
Jiandong Meng ◽  
Yogesh Jaluria

This paper describes a model for the growth of gallium nitride in a vertical impinging metalorganic chemical vapor deposition (MOCVD) reactor. With trimethylgallium (TMGa) and ammonia (NH3) carried by hydrogen (H2) as precursors, the flow, temperature and concentration profiles are predicted by numerical modeling, which is performed using a commercial CFD software package CFD-ACE+. The growth rate is predicted based on detailed reaction mechanisms given in the literature, and related studies are carried out to verify the reliability and adaptability of the chosen chemical kinetics. A detailed mathematical model is developed first, and the complete chemical mechanisms are introduced. Then, the dependence of the growth rate and uniformity of the deposited layers on operating conditions, such as reactor operating pressure, susceptor temperature, inlet velocity and concentration of the precursors, is investigated to gain greater insight into the reactor performance and characteristics. Based on the simulation results, discussion is presented in this paper to offer the possibility of better control of the GaN film growth process, and to ultimately lead to an optimization of the process, with respect to production rate and film quality.


Author(s):  
H. X. Liang ◽  
Q. W. Wang ◽  
L. Q. Luo ◽  
Z. P. Feng

Three-dimensional numerical simulation was conducted to investigate the flow field and heat transfer performance of the Cross-Wavy Primary Surface (CWPS) recuperators for microturbines. Using high-effective compact recuperators to achieve high thermal efficiency is one of the key techniques in the development of microturbine in recent years. Recuperators need to have minimum volume and weight, high reliability and durability. Most important of all, they need to have high thermal-effectiveness and low pressure-losses so that the gas turbine system can achieve high thermal performances. These requirements have attracted some research efforts in designing and implementing low-cost and compact recuperators for gas turbine engines recently. One of the promising techniques to achieve this goal is the so-called primary surface channels with small hydraulic dimensions. In this paper, we conducted a three-dimensional numerical study of flow and heat transfer for the Cross-Wavy Primary Surface (CWPS) channels with two different geometries. In the CWPS configurations the secondary flow is created by means of curved and interrupted surfaces, which may disturb the thermal boundary layers and thus improve the thermal performances of the channels. To facilitate comparison, we chose the identical hydraulic diameters for the above four CWPS channels. Since our experiments on real recuperators showed that the Reynolds number ranges from 150 to 500 under the operating conditions, we implemented all the simulations under laminar flow situations. By analyzing the correlations of Nusselt numbers and friction factors vs. Reynolds numbers of the four CWPS channels, we found that the CWPS channels have superior and comprehensive thermal performance with high compactness, i.e., high heat transfer area to volume ratio, indicating excellent commercialized application in the compact recuperators.


2014 ◽  
Vol 755 ◽  
pp. 705-731 ◽  
Author(s):  
Sasan Sarmast ◽  
Reza Dadfar ◽  
Robert F. Mikkelsen ◽  
Philipp Schlatter ◽  
Stefan Ivanell ◽  
...  

AbstractTwo modal decomposition techniques are employed to analyse the stability of wind turbine wakes. A numerical study on a single wind turbine wake is carried out focusing on the instability onset of the trailing tip vortices shed from the turbine blades. The numerical model is based on large-eddy simulations (LES) of the Navier–Stokes equations using the actuator line (ACL) method to simulate the wake behind the Tjæreborg wind turbine. The wake is perturbed by low-amplitude excitation sources located in the neighbourhood of the tip spirals. The amplification of the waves travelling along the spiral triggers instabilities, leading to breakdown of the wake. Based on the grid configurations and the type of excitations, two basic flow cases, symmetric and asymmetric, are identified. In the symmetric setup, we impose a 120° symmetry condition in the dynamics of the flow and in the asymmetric setup we calculate the full 360° wake. Different cases are subsequently analysed using dynamic mode decomposition (DMD) and proper orthogonal decomposition (POD). The results reveal that the main instability mechanism is dispersive and that the modal growth in the symmetric setup arises only for some specific frequencies and spatial structures, e.g. two dominant groups of modes with positive growth (spatial structures) are identified, while breaking the symmetry reveals that almost all the modes have positive growth rate. In both setups, the most unstable modes have a non-dimensional spatial growth rate close to $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\pi /2$ and they are characterized by an out-of-phase displacement of successive helix turns leading to local vortex pairing. The present results indicate that the asymmetric case is crucial to study, as the stability characteristics of the flow change significantly compared to the symmetric configurations. Based on the constant non-dimensional growth rate of disturbances, we derive a new analytical relationship between the length of the wake up to the turbulent breakdown and the operating conditions of a wind turbine.


2014 ◽  
Vol 136 (9) ◽  
Author(s):  
Mike Probyn ◽  
Ben Thornber ◽  
Dimitris Drikakis ◽  
David Youngs ◽  
Robin Williams

This paper presents an investigation into the use of a moving mesh algorithm for solving unsteady turbulent mixing problems. The growth of a shock induced mixing zone following reshock, using an initial setup comparable to that of existing experimental work, is used to evaluate the behavior of the numerical scheme for single-mode Richtmyer–Meshkov instability (SM-RMI). Subsequently the code is used to evaluate the growth rate for a range of different initial conditions. The initial growth rate for three-dimensional (3D) SM Richtmyer–Meshkov is also presented for a number of different initial conditions. This numerical study details the development of the mixing layer width both prior to and after reshock. The numerical scheme used includes an arbitrary Lagrangian–Eulerian grid motion which is successfully used to reduce the mesh size and computational time while retaining the accuracy of the simulation results. Varying initial conditions shows that the growth rate after reshock is independent of the initial conditions for a SM provided that the initial growth remains in the linear regime.


Author(s):  
Anton Pylypenko ◽  
Yevgenii Rastigejev ◽  
Lijun Wang ◽  
Abolghasem Shahbazi

The objective of this work is to analyze the dynamics and regimes of cold gas-solid flow in a biomass gasifier that is built at North Carolina Agricultural and Technical State University and to identify its corresponding ranges of operating conditions. The value of the minimum fluidization velocity Umf ≈ 8 cm/s has been found experimentally in a series of measurements of a pressure drop in the fluidized bed filled with Gledart type-B silica sand for the range of superficial gas velocities between 0 and 40 cm/s. To complement the experimental results, a set of three-dimensional numerical simulations of the isothermal gas-solid flow based on Eulerian-Eulerian approach has been performed. The analysis of the fluidization characteristics such as axial void fraction distributions has allowed us to evaluate the dependence of the bed expansion ratios from the flow superficial velocity. Good agreement between experimental and numerical results for the considered operating conditions of the gasifier has been observed.


2019 ◽  
Vol 141 (8) ◽  
Author(s):  
Omar Jumaah ◽  
Yogesh Jaluria

Gallium nitride (GaN) is an attractive material for manufacturing light emitting diodes and other electronic devices due to its wide band-gap and superb optoelectronic performance. The quality of GaN thin film determines the reliability and durability of these devices. Metal-organic chemical vapor deposition (MOCVD) is a common technique used to fabricate high-quality GaN thin films. In this paper, GaN growth rate and uniformity in a vertical rotating disk MOCVD reactor are investigated on the basis of a three-dimensional computational fluid dynamics (CFD) model. GaN growth rate is investigated under the influence of reactor pressure, precursor concentration ratio, and composition of the carrier gas mixture. The numerical simulation shows that the carrier gas mixture and the reactor pressure have significant effects on growth rate and uniformity of GaN thin films. It is also found that an appropriate mixture of N2 and H2 may be employed as the carrier gas to improve the flow field characteristic in the reactor. This results in an improved crystal growth of GaN thin films.


1999 ◽  
Vol 14 (6) ◽  
pp. 2377-2380 ◽  
Author(s):  
Gyeong Soon Hwang ◽  
Sang Heup Moon ◽  
Suk Woo Nam ◽  
Chee Burm Shin

Profile evolution simulations during chemical vapor deposition based on a 2D continuum model reveal that the type of surface kinetics plays an important role in determining step coverage of films deposited in high aspect ratio trenches and vias. Linear surface kinetics, resulting from an adsorption rate limited process, is found to cause difficulty in bringing about conformal step coverage in deep narrow trenches without reducing the growth rate considerably. Under such condition, void-free filling cannot be achieved while maintaining a growth rate acceptable to integrated circuit (IC) manufacturing. The numerical study also suggests that the high tendency of the precursor for chemical equilibrium on a surface, resulting in nonlinear kinetics by a surface reaction limited process, is crucial to achieve a uniform step coverage as typically observed in SiO2 deposition from tetraethylorthosilicate (TEOS).


Sign in / Sign up

Export Citation Format

Share Document