Thermosolutal convection during dendritic solidification of alloys: Part II. Nonlinear convection

A mathematical model of thermosolutal convection in directionally solidified dendritic alloys has been developed that includes a mushy zone underlying an all-liquid region. The model assumes a nonconvective initial state with planar and horizontal isotherms and isoconcentrates that move upward at a constant solidification velocity. The initial state is perturbed, nonlinear calculations are performed to model convection of the liquid when the system is unstable, and the results are compared with the predictions of a linear stability analysis. The mushy zone is modeled as a porous medium of variable porosity consistent with the volume fraction of, interdendritic liquid that satisfies the conservation equations for energy and solute concentrations. Results are presented for systems involving lead-tin alloys (Pb-10 wt pct Sn and Pb-20 wt pct Sn) and show significant differences with results of plane-front solidification. The calculations show that convection in the mushy zone is mainly driven by convection in the all-liquid region, and convection of the interdendritic liquid is only significant in the upper 20 pct of the mushy zone if it is significant at all. The calculated results also show that the systems are stable at reduced gravity levels of the order of 10^-4 g ₀ (g ₀=980 cm·s^-1) or when the lateral dimensions of the container are small enough, for stable temperature gradients between 2.5≤G _l≤100 K·cm^-1 at solidification velocities of 2 to 8 cm·h^-1.