Thermosolutal convection during directional solidification

During solidification of a binary alloy at constant velocity vertically upward, thermosolutal convection can occur if the solute rejected at the crystal-melt interface decreases the density of the melt. We assume that the crystal-melt interface remains planar and that the flow field is periodic in the horizontal direction. The time-dependent nonlinear differential equations for fluid flow, concentration, and temperature are solved numerically in two spatial dimensions for small Prandtl numbers and moderately large Schmidt numbers. For slow solidification velocities, the thermal field has an important stabilizing influence: near the onset of instability the flow is confined to the vicinity of the crystal-melt interface. Further, for slow velocities, as the concentration increases, the horizontal wavelength of the flow decreases rapidly — a phenomenon also indicated by linear stability analysis. The lateral in-homogeneity in solute concentration due to convection is obtained from the calculations. For a narrow range of solutal Rayleigh numbers and wavelengths, the flow is periodic in time. Formerly with the Mathematical Analysis Division, Center for Applied Mathematics, National Bureau of Standards, Washington, DC 20234. This paper is based on a presentation made at the symposium “Fluid Flow at Solid-Liquid Interfaces” held at the fall meeting of the TMS-AIME in Philadelphia, PA on October 5, 1983 under the TMS-AIME Solidification Committee.