Experimental Study on the Transition Process to the Oscillatory Thermocapillary Convections in a Floating Half Zone

The transition process of the thermocapillary convection from a steady and axisymmetric mode to the oscillatory mode in a liquid bridge with a fixed aspect ratio and varied volume ratio was studied experimentally. To ensure the surface tension to play an important role in the ground-based experiment, the geometrical configuration of the liquid bridge was so designed that the associated dynamic Bond number Bd ≈ 1. The velocity fields were measured by Particle Image Velocimetry (PIV) technique to effectively distinguish the different flow modes during the transition period in the experiments. Our experiments showed that as the temperature difference increased the slender and fat bridges presented quite different features on the evolution in their flow feature: for the former the thermocapillary convection transformed from a steady and axisymmetric pattern directly into an oscillatory one; but for the latter a transition flow status, characterized by an axial asymmetric steady convection, appeared before reaching the oscillatory mode. Experimental observations agree with the results of numerical simulations and it is obvious that the volume of liquid bridge is a sensitive geometric parameter. In addition, at the initial stage of the oscillation, for the former a rotating oscillatory convection with azimuthal wave number m = 1 was observed while for the latter a pulsating oscillatory pattern with azimuthal wave number m = 2 emerged, and then with further increase of the temperature difference, the pulsating oscillatory convection with azimuthal wave number m = 2 evolved into a rotating oscillatory pattern with azimuthal wave number m = 2.

A simulation model of a floating half zone

A simulation model of a floating half zone was suggested by steady numerical simulation and experiment respectively, in the previous papers [Q.S. Chen, W.R. Hu, Int. J. Mass Heat Transfer 40 (1997) 757; J.H. Han, Y. Ar, R. Zhou, W.R. Hue, Int. J. Mass Heat Transfer 40 (1997) 2671]. In the present paper, the simulation model is studied by using the method of unsteady and three-dimensional numerical simulation, and the transient process from steady convection to oscillatory convection is especially analyzed. Comparison of onsets of oscillation for both simulation model and the usual model were obtained, and the results show that the critical Marangoni number of the simulation model is obviously smaller than that of the usual model for the same slender liquid bridge. This implies that the usual model of a floating half zone gives a lower estimation on the onset of oscillation for floating zone convection.

Magnetic field design for floating zone crystal growth

The magnetic fields produced by electrical coils are designed for P-doped Si crystal growth in a floating full zone in microgravity environment. The fields are designed specially to reduce the how near the free surface and then in the melt zone by adjusting the coil positions near the melt zone. The effects of the designed magnetic fields on reducing the Row velocity and the non-uniformity of the concentration distribution in the melt zone are better than those of the case of a uniform longitudinal magnetic field, obtained by numerical simulation. It is expected to improve the radial macro-segregation and reduce the convection in the crystal growth at the same time by using the designed magnetic field.

Influence of Forced Convection on Columnar Microstructure During Directional Solidification of Al-Ni Alloys

This paper presents a summary of cellular and dendritic morphologies resulting from the upward directional solidification of Al - Ni alloys in a cylindrical crucible. We analysed the coupling of solid-liquid interface morphology with natural and forced convection. The influence of natural convection was first analyzed as a function of growth parameters (solute concentration, growth rate and thermal gradient). In a second step, the influence of axial vibrations on solidification microstructure was investigated by varying vibration parameters (amplitude and frequency). Experimental results were compared to preliminary numerical simulations and a good agreement is found for natural convection. In this study, the critical role of the mushy zone in the interaction between fluid flow and solidification microstructure is pointed out.

Thermocapillary Convection in Floating Zones

This paper provides an overview of ongoing studies in the area of thermocapillary convection driven by a surface tension gradient parallel to the free surface in a floating zone. Here, research interests are focused around the onset of oscillatory thermocapillary convection, also known as the transition from quasisteady convection to oscillatory convection. The onset of oscillation depends on a set of critical parameters, and the margin relationship can be represented by a complex function of the critical parameters. The experimental results indicate that the velocity deviation of an oscillatory flow has the same order of magnitude as that of an average flow, and the deviations of other quantities, such as temperature and free surface radii fluctuations, are much smaller when compared with their normal counterparts. Therefore, the onset of oscillation should be a result of the dynamic process in a fluid, and the problem is a strongly nonlinear one. In the past few decades, several theoretical models have been introduced to tackle the problem using analytical methods, linear instability analysis methods, energy instability methods, and unsteady 3D numerical methods. The last of the above mentioned methods is known to be the most suitable for a thorough analysis of strong nonlinear processes, which generally leads to a better comparison with the experimental results. The transition from oscillatory thermocapillary convection to turbulence falls under the studies of chaotic behavior in a new system, which opens a fascinating new frontier in nonlinear science, a hot research area drawing many recent works. This paper reviews theoretical models and analysis, and also experimental research, on thermocapillary connection in floating zones. It cites 93 references.

Hydrothermal Instability Of Thermocapillary Convection In Large-Prandtl-Number Liquid Bridges Under Microgravity

Linear stability analysis was performed to study the mechanism of transition of thermocapillary convection in liquid bridges with liquid volume ratios ranging from 0.4 to 1.2, aspect ratio of 0.75 and Prandtl number of 100. 2-D governing equations were solved to obtain the steady axi-symmetric basic flow and temperature distributions. 3-D perturbation equations were discretized at the collocation grid points using the Chebyshev-collocation method. Eigenvalues and eigenfunctions were obtained by using the Q-R. method. The predicted critical Marangoni numbers and critical frequencies were compared with data from space experiments. The disturbance of the temperature distribution on the free surface causes the onset of oscillatory convection. It is shown that the origin of instability is related to the hydrothermal origin for convections in large-Prandtl-number liquid bridges. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

Proper orthogonal decomposition of oscillatory Marangoni flow in half-zone liquid bridges of low-Pr fluids

Proper orthogonal decomposition (POD) using method of snapshots was performed on three different types of oscillatory Marangoni flows in half-zone liquid bridges of low-Pr fluid (Pr = 0.01). For each oscillation type, a series of characteristic modes (eigenfunctions) have been extracted from the velocity and temperature disturbances, and the POD provided spatial structures of the eigenfunctions, their oscillation frequencies, amplitudes, and phase shifts between them. The present analyses revealed the common features of the characteristic modes for different oscillation modes: four major velocity eigenfunctions captured more than 99% of the velocity fluctuation energy form two pairs, one of which is the most energetic. Different from the velocity disturbance, one of the major temperature eigenfunctions makes the dominant contribution to the temperature fluctuation energy. On the other hand, within the most energetic velocity eigenfuction pair, the two eigenfunctions have similar spatial structures and were tightly coupled to oscillate with the same frequency, and it was determined that the spatial structures and phase shifts of the eigenfunctions produced the different oscillatory disturbances. The interaction of other major modes only enriches the secondary spatio-temporal structures of the oscillatory disturbances. Moreover, the present analyses imply that the oscillatory disturbance, which is hydrodynamic in nature, primarily originates from the interior of the liquid bridge. (C) 2007 Elsevier B.V. All rights reserved.

Transition from steady to oscillatory convection with chaotic feature in thermocapillary convection

The transition process from steady to turbulent convection via subharmonic bifurcation in thermocapillary convection of half floating zone was studied by numerical simulation and experimental test. Both approaches gave structure of period doubling bifurcations in the present paper, and the Feigenbaum universal law was checked for the system of thermocapillary convection.

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.

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.

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.

Experimental study on oscillatory thermocapillary convection

The transition processes from steady flow into oscillatory flow in a liquid bridge of the half floating zone are studied experimentally. Two methods of noncontacted diagnoses are developed to measure the distribution of critical Marangoni numbers described by the onset of the oscillation st the free surface of the liquid bridge.The experimental results obtained for both cases of the upper rod heated and the lower rod heated agree with the prediction by Rayleigh's instability theory.The sensitive relations between the relatively fat or slender liquid bridge and the onset of oscillatory convection are also discussed to reveal the insight of the pressure distribution near the free surface. The experiments have been performed in a small liquid bridge, where the Bond number is much smaller than 1, and the results can be used to simulate the experiment in the microgravity environment.

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.

Effect of interfacial heat transfer on the onset of oscillatory convection in liquid bridge

In present study, effect of interfacial heat transfer with ambient gas on the onset of oscillatory convection in a liquid bridge of large Prandtl number on the ground is systematically investigated by the method of linear stability analyses. With both the constant and linear ambient air temperature distributions, the numerical results show that the interfacial heat transfer modifies the free-surface temperature distribution directly and then induces a steeper temperature gradient on the middle part of the free surface, which may destabilize the convection. On the other hand, the interfacial heat transfer restrains the temperature disturbances on the free surface, which may stabilize the convection. The two coupling effects result in a complex dependence of the stability property on the Biot number. Effects of melt free-surface deformation on the critical conditions of the oscillatory convection were also investigated. Moreover, to better understand the mechanism of the instabilities, rates of kinetic energy change and "thermal" energy change of the critical disturbances were investigated (C) 2009 Elsevier Ltd. All rights reserved.

两层流体热毛细对流空间实验研究=Space Experiments of Thermo Capillary Convection in Two-liquid Layers

1999年,在我国实践5号卫星上完成了两层流体空间实验,实验研究两层不相混合流体的纯Marangoni对流（温度梯度与界面垂直）与热毛细对流（温度梯度方向与流体界面平行）.前者存在发生Marangoni对流的最小临界温差值△Tc,低于该值流体系统处于静止状态;后者中只要存在沿界面的温度梯度便会产生热毛细对流.空间实验采用石蜡和氟化液两层流体新体系,实现了平整的液-液交界面,并从卫星上传回上万幅数字图像.通过多幅图像叠加处理得到了定量的流速场.数值模拟计算分析了相应工况时对流流动的速度场,两者的流场结构和速度大小基本一致,实验验证了理论模型.