Influence of Local Magnetic Fields on P-Doped Si Floating Zone Melting Crystal Growth in Microgravity

A two-dimensional axisymmetric numerical model is presented to study the influence of local magnetic fields on P-doped Si floating zone melting crystal growth in microgravity. The model is developed based on the finite difference method in a boundary-fitted curvilinear coordinate system. Extensive numerical simulations are carried out, and parameters studied include the curved growth interface shape and the magnetic field configurations. Computed results show that the local magnetic field is more effective in reducing the impurity concentration nonuniformity at the growth interface in comparison with the longitudinal magnetic field. Moreover, the curved growth interface causes more serious impurity concentration nonuniformity at the growth interface than the case with a planar growth interface.

Effect of non-uniform magnetic field on crystal growth by floating-zone method in microgravity

The magnetic damping effect of the non-uniform magnetic field on the floating-zone crystal growth process in microgravity is studied by numerical simulation. The results show that the non-uniform magnetic field with designed configuration can effectively reduce the flow near the free surface and then in the melt zone. At the same time, the designed magnetic field can improve the impurity concentration non-uniformity along the solidification interface. The primary principles of the magnetic field configuration design are also discussed.

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 liquid bridge volume on the onset of oscillation in floating-zone convection III. Three-dimensional model

An unsteady and three-dimensional model of the floating-half-zone convection on the ground is studied by the direct numerical simulation for the medium of 10 cSt silicon oil, and the influence of the liquid bridge volume on the critical applied temperature difference is especially discussed. The marginal curves for the onset of oscillation are separated into two branches related, respectively, to the slender liquid bridge and the fat liquid bridge. The oscillatory features of the floating-half-zone convection are also discussed.

Influence of Liquid Bridge Volume on the Floating Zone Convection

The liquid bridge volume is a critical geometrical parameter in addition to the aspect ratio for onset of oscillation in the floating zone convection. The oscillatory features are generally divided into two characteristic regions: slender liquid bridge region and fat liquid bridge region. The oscillatory modes in two regions are discussed in the present paper.

Numerical investigation on a simulation model of floating zone convection

A simulation model with adiabatic condition at the upper rod and constant temperature at the lower rod is studied numerically in this paper. The temperature distribution in a simulation model is closer to the one in the half part of a floating full zone in comparison with the one in a usual floating half zone model with constant temperature at both rods, because the temperature distribution of a floating full zone is symmetric for the middle plane in a microgravity environment. The results of the simulation model show that the temperature profiles and the how patterns are different from those of the usual floating half zone model. Another type of half zone model, with a special non-uniform temperature distribution at the upper rod and constant temperature at the lower rod, has been suggested by recent experiments. The temperature boundary condition of the upper rod has a maximum value in the center and a lower value near the free surface. This modified simulation model is also simulated numerically in the present paper. Copyright (C)1996 Elsevier Science Ltd.

Numerical simulation of g-jitter effects on half floating zone convection under microgravity environment

The g-jitter effects on the thermocapillary convection in liquid bridge of floating half zone were studied by numerical simulation for unsteady and axi-symmetric model in the cylindrical coordinate system. The g-jitter field was given by a steady microgravity field in addition to an oscillatory low-gravity field, and the effects on the flow field, temperature distribution and free surface deformation were analyzed numerically.

How to determine critical Marangoni number in half floating zone convection

In the present paper, the coordinated measurements of the temperature profile inside the liquid bridge and the boundary variation of Free surface, in addition to other quantities, were obtained in the same time for the half floating zone convection. The results show that the onset of free surface oscillation is earlier than the one of temperature oscillation during the increasing of applied temperature difference, and the critical Marangoni numbers, defined usually by temperature measurement, are larger than the one defined by free surface measurement, and the difference depends on the volume of liquid bridge. These results induce the question, ''How to determine experimentally the critical Marangoni number?'' Copyright (C) 1996 Elsevier Science Ltd.

Fractal features of oscillatory convection in the half-floating zone

Unsteady and two-dimensional numerical simulation is applied to study the transition process from steady convection to turbulence via subharmonic bifurcation in thermocapillary convection of a liquid bridge in the half-floating zone. The results of numerical tests show clearly the fractal structure of period-doubling bifurcations, and frequency-locking at f/4, f/8, f/16 with basic frequency f is observed with increasing temperature difference. The Feigenbaum universal constant is given by the present paper as delta(4) = 4.853, which can be compared with the theoretical value 4.6642016.

Hardware development and ground-based studies for drop shaft experiment on oscillatory convection in half floating-zone

Coordinated measurement of temperature, velocity and free surface oscillation were obtained by using the drop shaft facility for microgravity experiments of half floating zone convection. The ground-based studies gave transition from steady to oscillatory convection for multi-quantities measurement.

G-jitter effects on half floating-zone convection in intermediate-frequency range

The g-jitter influence on thermocapillary convection and critical Marangoni number in a liquid bridge of half-floating rone was discussed in the low frequency range of 0.4 to 1.5 Hz in a previous paper. This paper extended the experiments to the intermediate frequency range of 2 to 18 Hz, which htrs often been recorded as vibration environment of spacecrafts. The experiment was completed on the deck of a vibration machine, which gave a periodical applied acceleration to simulate the effects of g-jitter. The experimental results in the intermediate frequency range are different from that in the low frequency range. The velocity field and the shape of the free surface have periodical fluctuations in response to g-jitter. The amplitude of the periodical varying part of the temperature response decreases obviously with increasing frequency of g-jitter and vanishes almost when the frequency of g-jitter is high enough. The critical Marangoni number is defined to describe the transition from a periodical convection in response to g-jitter to an oscillatory convection due to internal instability, and will increase with increasing g-jitter frequency. According to the spectral analysis, it can be found that the oscillatory part of temperature is a superposition of two harmonic waves if the Marangoni number is larger than a critical value.

Influence of liquid bridge volume on the onset of oscillation in floating-zone convection .1. Experiments

The onset of oscillation in the floating zone convection driven by the gradient of surface tension was experimentally studied, and discussions were concentrated on the influence of liquid bridge volume on the onset of oscillation. Distributions of critical applied temperature difference and frequency depending on the volume of the liquid bridge were obtained, and there was a gap range of liquid volume which separated the curve of marginal stability into two parts for fixed rod diameter and aspect ratio. The results imply that the volume of the liquid bridge is a sensitive critical parameter for the onset of oscillation. The implication on the instability is also discussed in the present paper.

Influence of liquid bridge volume on the onset of oscillation in floating-zone convection .2. Numerical-simulation

The onset of oscillation in the floating zone convection driven by the gradient of surface tension was studied numerically for an unsteady and two-dimensional model, and studies were concentrated on the influence of liquid bridge volume on the onset of oscillation in comparison with the experimental results in the Paper I. The numerical results agree with the experimental ones presented in the previous paper, in which the distributions of critical applied temperature difference depending on the volume of liquid bridge and a gap range of liquid volume in marginal stability curve were obtained.

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 buoyancy on floating-zone convection of small bond number

Gravity may influence the velocity and temperature distributions, bouyancy may induce Rayleigh instability and the instability may be excited due to the change of free surface shape associating with gravity in the thermocapillary convection. These effects have been studied in the present paper. The results show that gravity may have an important effect in thermocapillary oscillatory convection even for the cases of small Bond number experiments either on the ground or in space.