Numerical Study on Flow Field and Temperature Distribution in Growth Process of 200 mm Czochralski Silicon Crystals

The melt flow and temperature distribution in a 200 mm silicon Czochralski furnace with a cusp magnetic field was modeled and simulated by using a finite-volume based FLUTRAPP ( Fluid Flow and Transport Phenomena Program) code. The melt flow in the crucible was focused, which is a result of the competition of buoyancy, the centrifugal forces caused by the rotations of the crucible and crystal, the thermocapillary force on the free surfaces and the Lorentz force induced by the cusp magnetic field. The zonal method for radiative heat transfer was used in the growth chamber, which was confined by the crystal surface, melt surface, crucible, heat shield, and pull chamber. It was found that the cusp magnetic field could strength the dominant counter-rotating swirling flow cell in the crucible and reduce the flow oscillation and the pulling-rate fluctuation. The fluctuation of dopant and oxygen concentration in the growing crystal could thus be smoothed.

Periodic vortex breakdown in wide spherical gaps

The flow field with vortex breakdown in wide spherical gaps was studied numerically by a finite difference method under the axisymmetric condition. The result shows that the flow bifurcates to periodic motion as the Reynolds number or the eccentricity of the spheres increases. (C) 1997 American Institute of Physics.

Drag on an ellipsoid particle in free-molecular plasma flow

Based on the analysis of molecular gas dynamics, the drag and moment acting on an ellipsoid particle of revolution X-2/a(2) + Y-2/a(2) + Z(2)/c(2) = 1, as an example of nonspherical particles, are studied under the condition of free-molecular plasma flow with thin plasma sheaths. A nonzero moment which causes nonspherical particle self-oscillation and self-rotation around its own axis in the plasma flow-similar to the pitching moment in aerodynamics-is discovered for the first time. When the ratio of axis length c/a is unity, the moment is zero and the drag formula are reduced to the well-known results of spherical particles. The effects of the particle-plasma relative velocity, the plasma temperature, and the particle materials on the drag and moment are also investigated.

Numerical study of bifurcation solutions of spherical Taylor-Couette flow

The steady bifurcation flows in a spherical gap (gap ratio sigma=0.18) with rotating inner and stationary outer spheres are simulated numerically for Re(c1)less than or equal to Re less than or equal to 1 500 by solving steady axisymmetric incompressible Navier-Stokes equations using a finite difference method. The simulation shows that there exist two steady stable flows with 1 or 2 vortices per hemisphere for 775 less than or equal to Re less than or equal to 1 220 and three steady stable flows with 0, 1, or 2 vortices for 1 220

Numerical simulation of axisymmetric unsteady incompressible flow by a vorticity-velocity method

A new numerical method for solving the axisymmetric unsteady incompressible Navier-Stokes equations using vorticity-velocity variables and a staggered grid is presented. The solution is advanced in time with an explicit two-stage Runge-Kutta method. At each stage a vector Poisson equation for velocity is solved. Some important aspects of staggering of the variable location, divergence-free correction to the velocity held by means of a suitably chosen scalar potential and numerical treatment of the vorticity boundary condition are examined. The axisymmetric spherical Couette flow between two concentric differentially rotating spheres is computed as an initial value problem. Comparison of the computational results using a staggered grid with those using a non-staggered grid shows that the staggered grid is superior to the non-staggered grid. The computed scenario of the transition from zero-vortex to two-vortex flow at moderate Reynolds number agrees with that simulated using a pseudospectral method, thus validating the temporal accuracy of our method.

A new type of boundary integral-equation for plane problems of elasticity including rotational forces

Based on the idea proposed by Hu [Scientia Sinica Series A XXX, 385-390 (1987)], a new type of boundary integral equation for plane problems of elasticity including rotational forces is derived and its boundary element formulation is presented. Numerical results for a rotating hollow disk are given to demonstrate the accuracy of the new type of boundary integral equation.

An application of adjoint variational method to the capillary instability in liquid

In this paper, applying the direct variational approach of first-order approximation to the capillary instability problem for the eases of rotating liquid column, toroid and films on both sides of cylinder, we have obtained the necessary and sufficient conditions for motion stability of the "cylindrical coreliquid-liquid-cylindrical shell" systems. The results obtained before are found to be special cases of the present investigation. At the same time, we have explained physical essence of rotating instability and settled a few disputes in previous investigations.

Non-linear theory of a positive-column in a magnetic-field

A nonlinear theory of an intermediate pressure discharge column in a magnetic field is presented. Motion of the neutral gas is considered. The continuity and momentum transfer equations for charged particles and neutral particles are solved by numerical methods. The main result obtained is that the rotating velocities of ionic gas and neutral gas are approximately equal. Bohm's criterion and potential inversion in the presence of neutral gas motion are also discussed.

复合材料飞轮破坏转速的算法和高速旋转破坏实验

分别基于平面应力型全弹性模型和三维数值模型建立了计算复合材料飞轮破坏转速的二维和三维算法.这两个算法均采用了正交各向异性材料的最大拉应力(材料主方向)准则,其中三维算法还采用了两种强度判据,即基于轴向大部分区域每层应力的平均值判据和轴向边界区域每层应力的最大值判据.对张紧力缠绕的3个实验复合材料飞轮成功实施了高速旋转破坏实验,破坏均发生在径向强度最弱的飞轮与金属芯轴的界面处.实验结果表明,飞轮的实际破坏转速与理论破坏转速十分接近,证实本文中建立的二维和三维算法是可靠的;二维算法得到的理论破坏转速偏高,而飞轮的实际破坏转速落在三维算法分别按最大值判据和平均值判据得到的两个理论破坏转速之间,说明三维算法的精度更高.

复合材料转动壳体动力特性分析

根据转动壳体的动力学方程，在利用DMAP对大型有限元软件MSC／NASTRAN二次开发的基础上，呆用有限元方法对复合材料层合壳体转动状态下的动力特性进行了研究，分析了科氏力、离心力和横向剪力对固有频率的影响，并与实验结果进行了比较．结果显示由于科氏力的作用，即使在动坐标系下，层合壳体也会出现频率不等的前后行波，横向剪力使壳体的固有频率降低，且随着壳体厚径比的增加其影响显著增加．

用非结构网格研究旋转振荡圆柱绕流

研究了雷诺数Re=200,1000,线速度比α=0.5,2.0,4.0,强迫振荡频率fs=0.1～2.0情况下的旋转振荡圆柱绕流问题.通过基于非结构同位网格有限体积法对Navier-Stokes方程进行数值求解.对流项、扩散项和非恒定项的离散格式均具有二阶精度,利用SIMPLE算法处理压力-速度耦合.计算得到了作用力系数随不同控制参数的变化规律.通过对升力系数的频谱分析得到自然脱落频率和强迫振荡频率下的作用力振幅.通过对不同频率作用力幅值的分析,得到频率之间的竞争关系,进而定量地给出了不同尾迹涡脱落模式的分区图.

Distribution of Temperature in GlidArc Discharge Plasma

The GlidArc discharge is one of the main generation methods of non-equilibrium plasma near atmospheric pressures. In general, Gliding Arc discharge is driven by gas flow [1] in axial direction or by magnetic field in circumferential direction. [2] In this paper, a GlidArc discharge driven by rotating-gas-flow in circumferential direction is presented. The principle of the plasma generator is analyzed. The distribution of the temperature in axial direction is measured by a digital thermometer for three different gases. The experimental set-up of the GlidArc plasma is shown in Fig.1. It consists of a center electrode, an outside electrode, a power supply and a gas supply. The shortest distance between the electrodes is 2-3 mm. When a power supply with 10000 volts is attached to the electrodes, the arc will be ignited at the shortest distance. The small plasma column is rotated by the rotating gas flow in circumferential direction and then the rotating arc is driven towards the exit of the setup by the gas flow.

Research on Q-Switched Nd : YAG laser engraving

IN this paper, the engraving process with Q-Switched Nd:YAG laser is investigated. High power density is the pre- requisition to vapor materials, and high repetition rate makes the engraving process highly efficient. An acousto- optic Q-Switch is applied in the cavity of CW 200 W Nd:YAG laser to achieve the high peak power density and the high pulse repetition rate. Different shape craters are formed in a patterned structure on the material surface when the laser beam irradiates on it by controlling power density, pulse repetition rate, pulse quantity and pulse interval. In addition, assisting oxygen gas is used for not only improving combustion to deepen the craters but also removing the plasma that generated on the top of craters. Off-focus length classified as negative and positive has a substantial effect on crater diameters. According to the message of rotating angle positions from material to be engraved and the information of graph pixels from computer, a special graph is imparted to the material by integrating the Q- Switched Nd:YAG laser with the computer graph manipulation and the numerically controlled worktable. The crater diameter depends on laser beam divergence and laser focal length. The crater diameter changes from 50 micrometers to 300 micrometers , and the maximum of crater depth reaches one millimeter.

A Breakdown Criterion of Free Molecular Flows and an Optimum Analysis of EBPVD

Two important issues in electron beam physical vapor deposition (EBPVD) are addressed. The first issue is a validity condition of the classical cosine law widely used in the engineering context. This requires a breakdown criterion of the free molecular assumption on which the cosine law is established. Using the analytical solution of free molecular effusion flow, the number of collisions (N-c) for a particle moving from an evaporative source to a substrate is estimated that is proven inversely proportional to the local Knudsen number at the evaporation surface. N-c = 1 is adopted as a breakdown criterion of the free molecular assumption, and it is verified by experimental data and DSMC results. The second issue is how to realize the uniform distributions of thickness and component over a large-area thin film. Our analysis shows that at relatively low evaporation rates the goal is easy achieved through arranging the evaporative source positions properly and rotating the substrate.

固体力学进展及应用: 庆贺李敏华院士90华诞文集

《固体力学进展及应用：庆贺李敏华院士90华诞文集》收录了近代固体力学基础理论及其应用领域的重要科技成果和最新进展。作者是在同体力学领域工作多年的资深研究员，他们来自各行各业，有丰富的科研与丁作经验。他们提供的论文在相当程度上反映当前同体力学的发展现状与成就，并能看出发展趋势，对未来研究的课题选择有参考价值。《固体力学进展及应用：庆贺李敏华院士90华诞文集》还收集了李敏华院士的珍贵照片和纪念李敏华院士90华诞的庆贺和回忆文章，具有重要的史料价值。