Influence of convections on dopant segregation in floating-zone crystal-growth system

The steady and axisymmetric crystal growth process of floating zone model was studied numerically to concern with the influence of convection and phase change on effective segregation. An iteration method of numerical simulation considering both thermocapillary and buoyancy effects for GaAs crystal growth gave the effective segregation coefficient, which was compared with the space experiment of GaAs on board the Chinese recoverable satellite. The calculated segregation coefficient of a two-dimensional model was found to be smaller than the one suggested by space experiment with the simplified assumption of an one-dimensional model.

Numerical-simulation of ac plasma-arc thermodynamics

A mathematical model and approximate analysis for the energy distribution of an ac plasma arc with a moving boundary is developed. A simplified electrical conductivity function is assumed so that the dynamic behavior of the arc may be determined, independent of the gas type. The model leads to a reduced set of non-linear partial differential equations which governs the quasi-steady ac arc. This system is solved numerically and it is found that convection plays an important role, not only in the temperature distribution, but also in arc disruptions. Moreover, disruptions are found to be influenced by convection only for a limited frequency range. The results of the present studies are applicable to the frequency range of 10-10(2) Hz which includes most industry ac arc frequencies. (C) 1994 Academic Press, Inc.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

Influence of low-gravity level on crystal-growth in floating-zone

The influence of low gravity level on crystal growth in the floating zone, which involves thermocapillary convection, phase change convection, thermal and solutal diffusion, is investigated numerically by a finite element method for the silicon crystal growth process. The velocity, temperature, concentration fields and phase change interfaces depending on heating temperature and growth rates are analyzed. The influence of low gravity level on the concentration is studied especially. The results show that the non-uniformities of concentration are about 10(-3) for growth rate nu(p) = 5.12 x 10(-8) m/s, 10(-2) for nu(p) = 5.12 x 10(-7) m/s and relatively larger for larger growth rate in the gravity level g = 0-9.8 m/s2. The thermocapillary effect is strong in comparison with the Bridgman system, and the level of low gravity is relatively insensitive for lower growth rates.

A study on stratification of concentration in crystal-growth of alpha-liio3 by holographic-interferometry

A new thermoplastic-photoconductor laser holographic recording system has been used for real-time and in situ observation of alpha-LiIO3 crystal growth. The influence of crystallization-driven convection on the concentration stratification in solution has been studied under gravity field. It is found that the stratification is closely related to the seed orientation of alpha-LiIO3 crystal. When the optical axis of crystal seed C is parallel to the gravity vector g, the velocity of the concentration stratification is two times larger than that in the case of C perpendicular-to g. It needs 40 h for the crystalline system of alpha-LiIO3 to reach stable concentration distribution (expressed as tau) at 47.6-degrees-C. The time tau is not sensitive to the seed orientation. Our results provide valuable data for designing the crystal growth experiments ia space.

H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.

Crystal-growth in floating zone with phase-change and thermosolutal convections

Floating zone crystal growth in microgravity environment is investigated numerically by a finite element method for semiconductor growth processing, which involves thermocapillary convection, phase change convection, thermal diffusion and solutal diffusion. The configurations of phase change interfaces and distributions of velocity, temperature and concentration fields are analyzed for typical conditions of pulling rates and segregation coefficients. The influence of phase change convection on the distribution of concentration is studied in detail. The results show that the thermocapillary convection plays an important role in mixing up the melt with dopant. The deformations of phase change interfaces by thermal convection-diffusion and pulling rods make larger variation of concentration field in comparison with the case of plane interfaces.

Some theoretical aspects of the simplified Navier-Stokes(SNS) equations

This study deals with the formulation, mathematical property and physical meaning of the simplified Navier-Stokes (SNS) equations. The tensorial SNS equations proposed is the simplest in form and is applicable to flow fields with arbitrary body boundaries. The zones of influence and dependence of the SNS equations, which are of primary importance to numerical solutions, are expounded for the first time from the viewpoint of subcharacteristics. Besides, a detailed analysis of the diffusion process in flow fields shows that the diffusion effect has an influence zone globally windward and an upwind propagation greatly depressed by convection. The maximum upwind influential distance of the viscous effect and the relative importance of the viscous effect in the flow direction to that in the direction normal to the flow are represented by the Reynolds number, which illustrates the conversion of the complete Navier-Stokes (NS) equations to the SNS equations for flows with large Reynolds number.

地幔对流的数值实验研究

数值实验是地幔对流研究的基本方法。本文概述了近二十年地幔对流研究的主要进展,介绍并评价了每个研究阶段所提出的新概念和新思想。最后,着重讨论了上地幔与全地幔对流的论争。

Gain saturation in gas flow and chemical lasers

A theoretical model for gain saturation in gas flow and chemical lasers is presented. The theory is applicable to all possible numerical values of τ/τc, where τ is the characteristie flow time for the flowing gas to move across the laser action region and τc is the characteristic collision relaxation time. The saturation effects of the convection and the "source flow" of the inverted population are revealed. A general relation of gain coefficient and some new gain saturation laws are obtained. For the special case of τ/τc1, the present theoretical results agree with the experimental results on the "anomalous" saturation phenomena in the supersonic diffusion HF chemical laser determined recently by Gross and Coffer[8]. The theory also agrees with the measured results of saturation intensity varying with τ/τc in gas flow CO2 lasers[7]. For the special case of τ/τc1, the present theory is consistent with both the standard theory[1] for gas lasers where the gas has no macroscopic motion and the known gain saturation theory[2-5] for gas flow and chemical lasers.

Steady axi-symmetrical thermocapillary motion of a short melting column

This paper deals with the steady axi-symmetric thermo-capillary motion in a short meltingcolumn.with the assumptions that the Marangoni number M<<1, the Reynolds number Re<<1 andthe capillary number C<<1, at the leading order, the solutions of the problem are obtained inthe form of series. For two kinds of typical cases, symmetric and anti-symmetric distributionof air temperature, the numerical calculations are made. The results describe the effect ofendwalls on thermo-capillary flow.

A model experiment of the plate movement

Detailed analysis of some difficult aspects has been made from modeling the platemovement. A new method of using differential density of material (instead of differentialtemperature) has been developed in the experiments. The effect of convection of mantleon the plate movement has been studied using a centrifugal technique, and a patternshowing the recurrence of the plate movement has been successfully obtained. In this paper, a criterion De=Dm is presented for the similarity of the model to thecounterpart of the original mantle. According to the criterion, what happens in the modelin a span of ten minutes suggests a process of the "original model" going on in geologi-cal time of three million years.

Generalized rayleigh principle and its applications

In this paper, we mainly deal with cigenvalue problems of non-self-adjoint operator. To begin with, the generalized Rayleigh variational principle, the idea of which was due to Morse and Feshbach, is examined in detail and proved more strictly in mathematics. Then, other three equivalent formulations of it are presented. While applying them to approximate calculation we find the condition under which the above variational method can be identified as the same with Galerkin's one. After that we illustrate the generalized variational principle by considering the hydrodynamic stability of plane Poiseuille flow and Bénard convection. Finally, the Rayleigh quotient method is extended to the cases of non-self-adjoint matrix in order to determine its strong eigenvalne in linear algebra.

池沸腾现象中热毛细对流的成因

针对微重力条件下单组份液体池沸腾现象，探讨了有关热毛细效应作用及其成因的不同观点，指出正确描述气液相变界面两侧的温度关系是解决争议的关键．详细评述了现有的气液相变界面温度模型，发现气液相变的非平衡特征会导致气液界面上产生温度梯度，引起表面张力的显著差异，从而驱动热毛细对流的形成．但现有气液相变界面温度模型机理差异迥然，预测结果也差别很大，因此，需要更深入研究．

短时微重力环境下的浮力-热毛细对流

在NMLC落塔中,以大Pr数硅油作为实验介质,对短时微重力环境下矩形液池中的液面变形和浮力热毛细对流进行了初步的实验研究.观测了气液自由界面的变形过程,确定了重定位时间,并验证了毛细爬升的两个特征时间;通过实验发现楔形结构和挡板对于液面控制有一定的作用,但是要想获得更好的控制效果,还需要采取其它措施使内角处不满足Concus-Finn条件;运用PIV方法得到了液体内部对流速度场,观察到了在重力水平突变和液面大变形的条件下浮力热毛细对流的转变过程.