Hybrid numerical methods for exponential models of growth

A series of novel numerical methods for the exponential models of growth are proposed. Based on these methods, hybrid predictor-corrector methods are constructed. The hybrid numerical methods can increase the accuracy and the computing speed obviously, as well as enlarge the stability domain greatly. (c) 2005 Published by Elsevier Inc.

Modeling of Ammonothermal Growth of Nitrides

Bulk single crystals of GaN and AlN can be grown from supercritical fluids using the ammonothermal method, which utilizes ammonia as fluid rather than water as in the hydrothermal process. In this process, a mineralizer such as amide, imide or nitride is used to attack a bulk nitride feedstock at temperatures from 200°C to 500°C and pressures from 1 to 4 kbar. Ammonothermal systems have been modeled here using fluid dynamics, thermodynamics and heat transfer models. The nutrient is considered as a porous media bed and the fluid flow is simulated using the Darcy-Brinkman-Forchheimer model. The resulting governing equations are solved using the finite volume method. The effects of particle size on flow pattern and temperature distribution in an autoclave are analyzed.

Numerical Analysis of LEC Growth of GaAs with an Axial Magnetic Field

A set of numerical analyses for momentum and heat transfer For a 3 in. (0.075 m) diameter Liquid Encapsulant Czochralski (LEC) growth of single-crystal GaAs with or without all axial magnetic field was carried Out using the finite-element method. The analyses assume a pseudosteady axisymmetric state with laminar floats. Convective and conductive heat transfers. radiative heat transfer between diffuse surfaces and the Navier-Stokes equations for both melt and encapsulant and electric current stream function equations Cor melt and crystal Lire considered together and solved simultaneously. The effect of the thickness of encapsulant. the imposed magnetic field strength as well as the rotation rate of crystal and crucible on the flow and heat transfer were investigated. (C) 2002 Published by Elsevier Science Ltd.

Effects of Baffle Design on Fluid Flow and Heat Transfer in Ammonothermal Growth System of Nitrides

Efforts have been made in growing bulk single crystals of GaN front supercritical fluids using the ammonothermal method, which utilizes ammonia as fluid rather than water as in the hydrothermal process. Different mineralizers such as amide or azide and temperatures in the range of 200-600degreesC have been used to increase the solubility. The pressure is from 1 to 4 kbar. Modeling of the ammonothermal growth process has been used to identify factors which may affect the temperature distribution, fluid flow and nutrient transport. The GaN charge is considered as a porous media bed and the flow in the charge is simulated using the Darcy-Brinkman-Forchheimer model. The resulting governing equations are solved using the finite volume method. The effects of baffle design and opening on flow pattern and temperature distribution in an autoclave are analyzed. Two cases are considered with baffle openings of 15% and 20% in cross-sectional area, respectively.

The effect of strain softening in the matrix material during void growth

Furthermore, the compressed flow driven by the piston is discussed. The consistent solution of gasdynamical equations including solar gravity is obtained for the unsteady and two-dimensional configuration, which is applied to the region between the piston and shock wave. This solution may satisfy the jump conditions of shock wave, which separates the region of compressed flow and quiet corona.

Numerical Simulation of Heat Transfer and Kinetics in the Bulk Growth of Silicon Carbide

The growth process of 2-inch silicon carbide (SiC) single crystals by the physical vapor transport method (or modified Lely method) has been modeled and simulated. The comprehensive process model incorporates the calculations of radio frequency (RF) induction heating, heat and mass transfer and growth kinetics. The transport equations for electromagnetic field, heat transfer, and species transport are solved using a finite volume-based numerical scheme called MASTRAPP (Multizone Adaptive Scheme for Transport and Phase Change Process). Temperature distribution for a 2-inch growth system is calculated, and the effects of induction heating frequency and current on the temperature distribution and growth rate are investigated. The predicted results have been compared with the experimental data.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

Effects of Thermocapillary Convection on the Growth Process in Industrial 8-inch Silicon Czochralski Growth

Czochralski (CZ) crystal growth process is a widely used technique in manufacturing of silicon crystals and other semiconductor materials. The ultimate goal of the IC industry is to have the highest quality substrates, which are free of point defect, impurities and micro defect clusters. The scale up of silicon wafer size from 200 mm to 300 mm requires large crucible size and more heat power. Transport phenomena in crystal growth processes are quite complex due to melt and gas flows that may be oscillatory and/or turbulent, coupled convection and radiation, impurities and dopant distributions, unsteady kinetics of the growth process, melt crystal interface dynamics, free surface and meniscus, stoichiometry in the case of compound materials. A global model has been developed to simulate the temperature distribution and melt flow in an 8-inch system. The present program features the fluid convection, magnetohydrodynamics, and radiation models. A multi-zone method is used to divide the Cz system into different zones, e.g., the melt, the crystal and the hot zone. For calculation of temperature distribution, the whole system inside the stainless chamber is considered. For the convective flow, only the melt is considered. The widely used zonal method divides the surface of the radiation enclosure into a number of zones, which has a uniform distribution of temperature, radiative properties and composition. The integro-differential equations for the radiative heat transfer are solved using the matrix inversion technique. The zonal method for radiative heat transfer is used in the growth chamber, which is confined by crystal surface, melt surface, heat shield, and pull chamber. Free surface and crystal/melt interface are tracked using adaptive grid generation. The competition between the thermocapillary convection induced by non-uniform temperature distributions on the free surface and the forced convection by the rotation of the crystal determines the interface shape, dopant distribution, and striation pattern. The temperature gradients on the free surface are influenced by the effects of the thermocapillary force on the free surface and the rotation of the crystal and the crucible.

Modeling of Ammonothermal Growth of Gallium Nitride Single Crystals

Ammonothermal growth of GaN crystals with a retrograde solubility has been modeled and simulated here using fluid dynamics, thermodynamics and heat transfer models. The nutrient is considered as a porous media bed and the flow in the porous charge is simulated using the Darcy-Brinkman-Forchheimer model. The resulting governing equations are solved using the finite volume method. For the case of retrograde solubility, the charge is put above the baffle. The temperature difference between the dissolving zone and growth zone is found smaller than that applied on the sidewall of autoclave. The baffle opening has a strong effect on the nutrient transport and supersaturation of GaN species in the growth zone.

Transport Phenomena in Ammonothermal Growth of Gallium Nitride

GaN can be used to fabricate blue/green/UV LEDs and high temperature, high power electronic devices. Ideal substrates are needed for high quality III-nitride epitaxy, which is an essential step for the manufacture of LEDs. GaN substrates are ideal to be lattice matched and isomorphic to nitride-based films. Bulk single crystals of GaN can be grown from supercritical fluids using the ammonothermal method, which utilizes ammonia as fluid rather than water as in the hydrothermal process. In this process, a mineralizer such as amide, imide or azide is used to attack a bulk nitride feedstock at temperatures from 200 - 500癈 and pressures from 1 - 4 kbar. Baffle design is essential for successful growth of GaN crystals. Baffle is used to separate the dissolving zone from the growth zone, and to maintain a temperature difference between the two zones. For solubility curve with a positive coefficient with respect to temperature, the growth zone is maintained at a lower temperature than that in the dissolving zone, thus the nutrient becomes supersaturated in the growth zone. The baffle opening is used to control the mixing of nutrients in the two zones, thus the transfer of nutrient from the lower part to the upper part. Ammonothermal systems have been modeled here using fluid dynamics, thermodynamics and heat transfer models. The nutrient is considered as a porous media bed and the flow is simulated using the Darcy-Brinkman-Forchheimer model. The resulting governing equations are solved using the finite volume method. We investigated the effects of baffle opening and position on the transport phenomena of nutrient from dissolving zone to the growth zone. Simulation data have been compared qualitatively with experimental data.

Temperature-Dependent Growth and Life Cycle of Nemoura sichuanensis (Plecoptera: Nemouridae) in a Chinese Mountain Stream

Plecoptera constitute a numerically and ecologically significant component in mountain streams all over the world, but little is known of their life cycles in Asia. The life cycle of Nemoura sichuanensis and its relationship to water temperature was investigated during a 4-year study in a headwater stream (known as the Jiuchong torrent) of the Xiangxi River in Central China. Size structure histograms suggest that the life cycle was univoltine, and the relationships between the growth of Nemoura sichuanensis, physiological time, and effective accumulated water temperature were described using logistic regressions. The growth pattern was generally similar within year classes but growth rates did vary between year-classes. Our field data suggest a critical thermal threshold for emergence in Nemoura sichuanensis, that was close to 9 degrees C. The total number of physiological days required for completing larval development was 250 days. The effective accumulated water temperature was 2500 degree-days in the field. Development during the life cycle increased somewhat linearly with the physiological time and the effective accumulated water temperature, but some non-linear relationships were best developed by logistic equations.

Effect of dietary cyanobacteria on growth and accumulation of microcystins in Nile tilapia (Oreochromis niloticus)

A 12-week growth trial was conducted in a flow-through system to investigate the chronic toxic effect of dietary intake of cyanobacteria on growth, feed utilization and microcystins accumulation in Nile tilapia (Oreochromis niloticus L.) (initial body weight: 5.6 g). Six isonitrogenous and isocaloric diets were formulated to include different contents of cyanobacteria with the dietary microcystins increasing from 0 to 5460.06 ng/g diet. The results showed that dietary intake of cyanobacteria could increase the growth of tilapia while there are no impacts on feed conversion efficiency or mortality. Feeding rate was higher for the diets containing highest cyanobacteria. Microcystins were mostly accumulated in fish liver. The relationship between microcystins contents in muscle, liver, spleen and dietary intake could be described by quadratic equations. Microcystins content in the muscle of Nile tilapia in present study exceeded the upper limit of the tolerable daily intake (TDI) of microcystins suggested by the WHO (0.04 mu g/kg body weight/d). It is suggested that Nile tilapia fed on toxic cyanobacteria is not suitable for human food. (c) 2006 Elsevier B.V. All rights reserved.

Food consumption and growth of two piscivorous fishes, the mandarin fish and the Chinese snakehead

Rates of maximum food consumption and growth were determined for immature mandarin fish Siniperca chuatsi (47.2-540.2 g) and Chinese snakehead Channa argus (45.0-546.2 g) at 10, 15, 20, 25, 30 and 35 degrees C. The relationship between maximum rate of food consumption (C-max), body weight (W) and temperature (T) was described by the multiple regression equations: lnC(max) = -4.880 + 0.597 lnW+0.284T - 0.0048T(2) for the mandarin fish, and lnC(max)= -6.718 + 0.522 lnW+0.440T-0.0077T(2) for the Chinese snakehead. The optimum temperature for consumption was 29.6 degrees C for the mandarin fish and 28.6 degrees C for the Chinese snakehead. The relationship between growth rate (G), body weight and temperature was ln(G+0.25)= - 0.439 - 0.500 lnW+0.270T - 0.0046T(2) for the mandarin fish, and ln(G+0.25)= - 6.150+ (0.175 - 0.026T) lnW+0.571T - 0.0078T(2) for the Chinese snakehead. The weight exponent in the growth-weight relationship was -0.83 for the mandarin fish, but decreased with increasing temperature for the Chinese snakehead. The optimum temperature for growth was 29.3 degrees C for the mandarin fish, but tended to decrease with increasing weight for the Chinese snakehead, being 30.3 degrees C for a 45-g fish, and 26.1 degrees C for a 550-g fish. (C) 1998 The Fisheries Society of the British Isles.

Numerical analysis of LEC growth of GaAs with an axial magnetic field

A set of numerical analyses for momentum and heat transfer For a 3 in. (0.075 m) diameter Liquid Encapsulant Czochralski (LEC) growth of single-crystal GaAs with or without all axial magnetic field was carried Out using the finite-element method. The analyses assume a pseudosteady axisymmetric state with laminar floats. Convective and conductive heat transfers. radiative heat transfer between diffuse surfaces and the Navier-Stokes equations for both melt and encapsulant and electric current stream function equations Cor melt and crystal Lire considered together and solved simultaneously. The effect of the thickness of encapsulant. the imposed magnetic field strength as well as the rotation rate of crystal and crucible on the flow and heat transfer were investigated. (C) 2002 Published by Elsevier Science Ltd.

THEORETICAL TREATMENT OF A3B5 CHLORIDE VAPOR-PHASE EPITAXY - GROWTH, DOPING, OPTIMIZATION

A theoretical description of chloride vapour-phase epitaxy (CVPE) has been proposed which contains two-dimensional (2D) gas-dynamic equations for transport of reactive components and kinetic equations for surface growth processes connected by nonlinear adiabatic boundary conditions. No one of these stages is supposed to be the limiting one. Calculated variations of growth rate and impurity concentrations along the growing layer fit experimental data well.