NONLINEAR-ANALYSIS OF BUOYANT CONVECTION IN BINARY SOLIDIFICATION WITH APPLICATION TO CHANNEL FORMATION

We consider the problem of nonlinear thermal-solutal convection in the mushy zone accompanying unstable directional solidification of binary systems. Attention is focused on possible nonlinear mechanisms of chimney formation leading to the occurrence of freckles in solid castings, and in particular the coupling between the convection and the resulting porosity of the mush. We make analytical progress by considering the case of small growth Peclet number, delta, small departures from the eutectic point, and infinite Lewis number. Our linear stability results indicate a small O(delta) shift in the critical Darcy-Rayleigh number, in accord with previous analyses. We find that nonlinear two-dimensional rolls may be either sub- or supercritical, depending upon a single parameter combining the magnitude of the dependence of mush permeability on solids fraction and the variations in solids fraction owing to melting or freezing. A critical value of this combined parameter is given for the transition from supercritical to subcritical rolls. Three-dimensional hexagons are found to be transcritical, with branches corresponding to upflow and lower porosity in either the centres or boundaries of the cells. These general results are discussed in relation to experimental observations and are found to be in general qualitative agreement with them.