TY - JOUR
T1 - An analysis of convection in a mushy layer with a deformable permeable interface
AU - Roper, S. M.
AU - Davis, S. H.
AU - Voorhees, P. W.
PY - 2008/1/25
Y1 - 2008/1/25
N2 - We study the dynamics of a mushy layer in directional solidification for the case of a thin near-eutectic mush with a deformable and permeable mush - liquid interface. We examine the onset of convection using linear stability analysis, and the weakly nonlinear growth of liquid inclusions that signal the onset of chimneys. This analysis is compared to past analyses in which the mush - liquid interface is replaced by a rigid impermeable lid. We find qualitative agreement between the two models, but the rigid-lid approximation gives substantially different quantitative behaviour. In linear theory, the rigid-lid approximation leads to an over-estimate of the critical Rayleigh number and wavenumber of the instability. The condition for the onset of oscillatory instability is also changed by a factor of about 5 in composition number C. In the weakly nonlinear theory, the location of the onset of liquid inclusions is near the undisturbed front for the free-boundary analysis, whereas it lies at the centre of the mushy layer when the rigid-lid approximation is used. For hexagonal patterns, the boundary between regions of parameter space in which up and down hexagons are stable, shifts as a result of coupling between the liquid and mush regions.
AB - We study the dynamics of a mushy layer in directional solidification for the case of a thin near-eutectic mush with a deformable and permeable mush - liquid interface. We examine the onset of convection using linear stability analysis, and the weakly nonlinear growth of liquid inclusions that signal the onset of chimneys. This analysis is compared to past analyses in which the mush - liquid interface is replaced by a rigid impermeable lid. We find qualitative agreement between the two models, but the rigid-lid approximation gives substantially different quantitative behaviour. In linear theory, the rigid-lid approximation leads to an over-estimate of the critical Rayleigh number and wavenumber of the instability. The condition for the onset of oscillatory instability is also changed by a factor of about 5 in composition number C. In the weakly nonlinear theory, the location of the onset of liquid inclusions is near the undisturbed front for the free-boundary analysis, whereas it lies at the centre of the mushy layer when the rigid-lid approximation is used. For hexagonal patterns, the boundary between regions of parameter space in which up and down hexagons are stable, shifts as a result of coupling between the liquid and mush regions.
UR - http://www.scopus.com/inward/record.url?scp=38349109421&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=38349109421&partnerID=8YFLogxK
U2 - 10.1017/S0022112007009469
DO - 10.1017/S0022112007009469
M3 - Article
AN - SCOPUS:38349109421
SN - 0022-1120
VL - 596
SP - 333
EP - 352
JO - Journal of Fluid Mechanics
JF - Journal of Fluid Mechanics
ER -