Search

The influence of thermosolutal convection on solute segregation in crystals grown by vertical directional solidification of binary metallic alloys or semiconductors is studied. Finite differences are used in a two-dimensional time-dependent model which assumes a planar crystal-melt interface to obtain numerical results. It is assumed that the configuration is periodic in the horizontal direction. Consideration is given to the possibility of multiple flow states sharing the same period. The results are represented in bifurcation diagrams of the nonlinear states associated with the critical points of linear theory. Variations of the solutal Rayleigh number can lead to the occurrence of multiple steady states, time-periodic states, and quasi-periodic states. This case is compared to that of thermosolutal convection with linear vertical gradients and stress-free boundaries.

Stommel et al. (1956) have first described an instability, known as thermosolutal convection, thermohaline convection, or double-diffusive convection. This instability may occur in the case of a fluid in a gravitational field with two diffusing components present. The present study is concerned with the effect of sidewalls on flow in the fingering regime in the absence of applied horizontal gradients. The work was motivated by numerical results obtained on the basis of a simulation of thermosolutal convection occurring during the unidirectional solidification of a binary alloy. In this case, the unperturbed solute field in the liquid ahead of the solidifying planar interface has an exponential vertical profile because of the rejection or preferential incorporation of solute by the solid phase.

Linear stability of the one-dimensional flow between infinite vertical coaxial cylinders induced by heating the inner cylinder is considered for various ratios (kappa) of the inner radius to the outer radius, and for Prandtl numbers (P) appropriate to air and water. For air with P = 0.71 the least stable disturbance is nonaxisymmetric for kappa less than 0.44 and is axisymmetric for kappa greater than 0.44, and in either case the instability is due to the action of the shear forces. For P = 3.5, the situation is similar, except that the asymmetric shear mode is superseded by an axisymmetric instability driven by buoyancy forces for kappa = 0.03-0.16. Wave speeds, wavenumbers, and critical Grashof numbers for these cases and for the case of zero Prandtl number are given.