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.

QE and Suns-V-oc study on the epitaxial CSiTF solar cells

In order to clarify the major factors having confined the efficiencies of as-prepared crystalline silicon thin film (CSiTF) solar cells on the SSP (silicon sheets from powder) ribbons, QE (quantum efficiency) and Suns-V-oc study were performed on the epitaxial CSiTF solar cells fabricated on the SSP ribbons, the SSP ribbons after surface being zone melting recrystallized (ZMR) and single crystalline silicon (sc-Si) substrates. The results show that the epi-layers deposited on the SSP ribbons have rough surfaces, which not only increases the diffusion reflectance on the surfaces but also makes the anti-reflection coatings become structure-loosened, both of which would deteriorate the light trapping effect; in addition, the epi-layers deposited on the SSP ribbons possess poor crystallographic quality, so the heavy grain boundary (GB) recombination limits the diffusion length of the minority carriers in the epi-layers, which makes the as-prepared CSiTF solar cells suffer the worse spectra response at long-wavelength range. Nearly all the dark characteristic parameters of the CSiTF solar cells are far away from the ideal values. The performances of the CSiTF solar cells are especially affected by too high I-02 (the dark saturation current of space charge region) values and too low R-sh (parallel resistance) values. The higher 102 values are mainly caused by the heavy GB recombination resulting from the poor crystallographic qualities of the silicon active layers in the space charge regions, while the lower R-sh values are attributed to the electrical leakage at the un-passivated PN junction or solar cell edges after the solar cells are cut by the laser scriber.

Coupling of melting-solidification phase-change convection with floating-zone convection under microgravity

A finite element algorithm is used to analyze the process of floating zone crystal growth under microgravity. The effect of phase change convection coupled with surface tension convection is considered. The results show that the rate of crystal growth is very important. The single-crystal-melt interface is steeper than the feed-melt interface during the process of crystal growth. When the rate exceeds a critical value, the Marangoni vortex near the feed-melt interface will become so large that a secondary vortex will exist.

Application of enthalpy method for moving boundary problem in laser surface melting

A numerical analysis was carried out to study the moving boundary problem in the physical process of pulsed Nd-YAG laser surface melting prior to vaporization. The enthalpy method was applied to solve this two-phase axisymmetrical melting problem Computational results of temperature fields were obtained, which provide useful information to practical laser treatment processing. The validity of enthalpy method in solving such problems is presented.

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.

Influence of grain boundary on melting

The temperature behaviour of an Al bicrystal with surfaces consisting of (110) and (111) crystals is simulated using molecular dynamics. The result shows that the (110) crystal losses its crystalline order at 820K, whereas the disorder does not propagate through the (111) crystal at this temperature. Instead, some disordered atoms are recrystallized into the (111) crystal and the initial grain boundary changes into a stable order-disorder interface. Thus, it was discovered that at a temperature near its melting point, the (111) crystal grew and obstructed the propagation of disorder. Such an obstruction is helpful for understanding melting.

Numerical Simulation of Transient Thermal Field and Residual Stress in Laser Melting Process

以激光熔凝表面强韧化处理为背景,应用空间弹塑性有限单元和高精度数值算法同时考虑材料组织性能的变化模拟工件的温度场及残余应力,研究激光熔凝加工中瞬时温度场及残余应力数值模拟,同时考虑相变潜热及相变塑性的影响,用算例验证了模型的正确性,给出了不同时刻温度场分布及残余应力分布。

Three-Dimensional Modeling of Heat Transfer and Fluid Flow in Laminar-Plasma Material Re-Melting Processing

Modeling study is performed concerning the heat transfer and fluid flow for a laminar argon plasma jet impinging normally upon a flat workpiece exposed to the ambient air. The diffusion of the air into the plasma jet is handled by using the combined-diffusion-coefficient approach. The heat flux density and jet shear stress distributions at the workpiece surface obtained from the plasma jet modeling are then used to study the re-melting process of a carbon steel workpiece. Besides the heat conduction within the workpiece, the effects of the plasma-jet inlet parameters (temperature and velocity), workpiece moving speed, Marangoni convection, natural convection etc. on the re-melting process are considered. The modeling results demonstrate that the shapes and sizes of the molten pool in the workpiece are influenced appreciably by the plasma-jet inlet parameters, workpiece moving speed and Marangoni convection. The jet shear stress manifests its effect at higher plasma-jet inlet velocities, while the natural convection effect can be ignored. The modeling results of the molten pool sizes agree reasonably with available experimental data.

Numerical simulation of magnetic field design for damping thermocapillary convection in a floating half-zone

The magnetic fields produced by electrical coils are designed for damping the the thermocapillary convection in a floating half-zone in microgravity. The fields are designed specially to reduce the flow near the free surface and then in the melt zone by adjusting the longitudinal coil positions close to the melt zone. The effects of the designed magnetic fields on reducing the flow velocity and temperature distribution non-uniformity in the melt zone are stronger than those of the case of an uniform longitudinal magnetic field obtained by numerical simulation, particularly at the melt-rod interface. It brings fundamental insights into the heat and mass transfer control at the solidification interface by the magnetic field design for crystal growth by the floating full-zone method.

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.

Study on Flow Patterns of Oscillatory Thermocapillary Convection In Floating Half Zone by PIV

The velocity fields of oscillatory convection have been measured using the techniques of Particle Image Velocimetry (PIV) in a liquid bridge of half floating zone with small typical scales of a few millimeters for emphasizing the thermocapillary effect in comparison with the effect of buoyancy. The flow patterns of the oscillatory flow have been studied experimentally in a liquid bridge. The flow patterns in the liquid bridge are classified with mode numbers according to oscillatory flow characteristics. Results of the experiment show that the mode depends on the aspect ratio as well as the volume ratio of the liquid bridge. The experimental results are helpful for studying the structure of flow at the onset of oscillatory thermocapillary convection in a liquid bridge.

Investigation on Depth of Heat Affected Zone of Discharge Spot by High Repetitive Rate YAG Laser-Induced Discharge Texturing

It is known that the press formability and the elongation of laser textured sheet are improved, and the service life of textured roll is longer than that of the un-textured roll due to hardening of the treated surface. One of the goals to develop high repetitive rate YAG laser-induced discharge texturing (LIDT) is to get deeper hardening zone. By observing and measuring cross-section of LIDT spots in different discharge conditions, it is found that the single-crater, which is formed by the discharge conditions of anode, which is covered by an oil film and with rectangular current waveform, has the most depth of heat affected zone (HAZ) comparing with other crater shapes when discharge energy is the same. The depth of HAZ is mainly depends on pulse duration when the discharge spot is single-crater. The results are analyzed.

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.