Localisation of convection in mushy layers by weak background flow

2011 ◽  
Vol 675 ◽  
pp. 518-528 ◽  
Author(s):  
S. M. ROPER ◽  
S. H. DAVIS ◽  
P. W. VOORHEES

It is known that freckles form at the sidewalls of directionally solidified materials. We present a weakly nonlinear analysis of the effects of a weak and slowly varying background flow formed by non-axial thermal gradients on convection near onset in a mushy layer. We find that in the two-dimensional case, the onset of mush convection occurs away from the walls. However if three-dimensional disturbances are allowed, the onset occurs near the walls of the container confining the mush. We derive amplitude equations governing this behaviour and simulate their evolution numerically.

2000 ◽  
Vol 408 ◽  
pp. 53-82 ◽  
Author(s):  
C. A. CHUNG ◽  
FALIN CHEN

A weakly nonlinear analysis is employed to investigate the onset of plume convection in the mushy layer of a binary solution directionally solidified from below. An improved mathematical model including the constant pressure condition at the melt/mush interface is applied to analytically analyse the nonlinear behaviour of the convection. Results show that, due to the consideration of the constant pressure condition at the interface, the present analytical results are much better in comparison with the experiments than those obtained by previous studies, in which the no-vertical-flow condition at the interface was considered. It is also shown that the bifurcation to three-dimensional hexagon convection, corresponding to the onset of plume convection, is subcritical. For the case of small concentration ratio [Cscr ] (equation (5b)) the channel of the plume is generated at the top of the mush and grows downwards into the mush. For the case of large [Cscr ], the channel may be generated at the interior of the mush and grow upwards to the top of the mush. The possible parameter regime in which the flow is of a stable down-centre hexagon is discussed.


2010 ◽  
Vol 645 ◽  
pp. 411-434 ◽  
Author(s):  
PETER GUBA ◽  
M. GRAE WORSTER

We study nonlinear, two-dimensional convection in a mushy layer during solidification of a binary mixture. We consider a particular limit in which the onset of oscillatory convection just precedes the onset of steady overturning convection, at a prescribed aspect ratio of convection patterns. This asymptotic limit allows us to determine nonlinear solutions analytically. The results provide a complete description of the stability of and transitions between steady and oscillatory convection as functions of the Rayleigh number and the compositional ratio. Of particular focus are the effects of the basic-state asymmetries and non-uniformity in the permeability of the mushy layer, which give rise to abrupt (hysteretic) transitions in the system. We find that the transition between travelling and standing waves, as well as that between standing waves and steady convection, can be hysteretic. The relevance of our theoretical predictions to recent experiments on directionally solidifying mushy layers is also discussed.


1995 ◽  
Vol 291 ◽  
pp. 57-81 ◽  
Author(s):  
S. M. Churilov ◽  
I. G. Shukhman

We consider the nonlinear spatial evolution in the streamwise direction of slightly three-dimensional disturbances in the form of oblique travelling waves (with spanwise wavenumber kz much less than the streamwise one kx) in a mixing layer vx = u(y) at large Reynolds numbers. A study is made of the transition (with the growth of amplitude) to the regime of a nonlinear critical layer (CL) from regimes of a viscous CL and an unsteady CL, which we have investigated earlier (Churilov & Shukhman 1994). We have found a new type of transition to the nonlinear CL regime that has no analogy in the two-dimensional case, namely the transition from a stage of ‘explosive’ development. A nonlinear evolution equation is obtained which describes the development of disturbances in a regime of a quasi-steady nonlinear CL. We show that unlike the two-dimensional case there are two stages of disturbance growth after transition. In the first stage (immediately after transition) the amplitude A increases as x. Later, at the second stage, the ‘classical’ law A ∼ x2/3 is reached, which is usual for two-dimensional disturbances. It is demonstrated that with the growth of kz the region of three-dimensional behaviour is expanded, in particular the amplitude threshold of transition to the nonlinear CL regime from a stage of ‘explosive’ development rises and therefore in the ‘strongly three-dimensional’ limit kz = O(kx) such a transition cannot be realized in the framework of weakly nonlinear theory.


1999 ◽  
Vol 395 ◽  
pp. 253-270 ◽  
Author(s):  
Y.-J. CHEN ◽  
S. H. DAVIS

A steady, two-dimensional cellular convection modifies the morphological instability of a binary alloy that undergoes directional solidification. When the convection wavelength is far longer than that of the morphological cells, the behaviour of the moving front is described by a slow, spatial–temporal dynamics obtained through a multiple-scale analysis. The resulting system has a parametric-excitation structure in space, with complex parameters characterizing the interactions between flow, solute diffusion, and rejection. The convection in general stabilizes two-dimensional disturbances, but destabilizes three-dimensional disturbances. When the flow is weak, the morphological instability is incommensurate with the flow wavelength, but as the flow gets stronger, the instability becomes quantized and forced to fit into the flow box. At large flow strength the instability is localized, confined in narrow envelopes. In this case the solutions are discrete eigenstates in an unbounded space. Their stability boundaries and asymptotics are obtained by a WKB analysis. The weakly nonlinear interaction is delivered through the Lyapunov–Schmidt method.


2001 ◽  
Vol 437 ◽  
pp. 337-365 ◽  
Author(s):  
PETER GUBA

This study concentrates on a relatively simple model of a mushy layer originally proposed by Amberg & Homsy (1993) and later studied in further detail by Anderson & Worster (1995). We extend this model to the case in which the system is in a state of uniform rotation about the vertical. Of particular interest is to determine how the rotation of the system controls the bifurcating convection with both the oblique-roll planform and the planform of hexagonal symmetry. We find that two-dimensional oblique rolls can be either subcritically or supercritically bifurcating, depending on a pair of parameters (K1/CS, [Tscr ]), where K1 measures how the permeability linearly varies with the local solid fraction, CS relates the compositional difference between the liquid and solid phases to the variation of composition throughout the mushy layer, and the Taylor number [Tscr ] gives a measure of the local Coriolis acceleration relative to the viscous dissipation in a porous medium. The three-dimensional convection with hexagonal symmetry is found to be transcritical. Furthermore, distorted hexagons with upflow at the centres can be either subcritical or supercritical, depending on the value of the Taylor number [Tscr ].


1998 ◽  
Vol 370 ◽  
pp. 249-270 ◽  
Author(s):  
A. M. METCALFE ◽  
T. J. PEDLEY

Complex bioconvection patterns are observed when a suspension of the oxytactic bacterium Bacillus subtilis is placed in a chamber with its upper surface open to the atmosphere. The patterns form because the bacteria are denser than water and swim upwards (up an oxygen gradient) on average. This results in an unstable density distribution and an overturning instability. The pattern formation is dependent on depth and experiments in a tilted chamber have shown that as the depth increases the first patterns formed are hexagons in which the fluid flows down in the centre.The linear stability of this system was analysed by Hillesdon & Pedley (1996) who found that the system is unstable if the Rayleigh number Γ exceeds a critical value, which depends on the wavenumber k of the disturbance as well as on the values of other parameters. Hillesdon & Pedley found that the critical wavenumber kc could be either zero or non-zero, depending on the parameter values.In this paper we carry out a weakly nonlinear analysis to determine the relative stability of hexagon and roll patterns formed at the onset of bioconvection. The analysis is different in the two cases kc≠0 and kc=0. For the kc≠0 case (which appears to be more relevant experimentally) the model does predict down hexagons, but only for a certain range of parameter values. Hence the analysis allows us to refine previous parameter estimates. For the kc=0 case we carry out a two-dimensional analysis and derive an equation describing the evolution of the horizontal planform function.


1995 ◽  
Vol 302 ◽  
pp. 307-331 ◽  
Author(s):  
D. M. Anderson ◽  
M. Grae Worster

We consider the solidification of a binary alloy in a mushy layer and analyse the system near the onset of buoyancy-driven convection in the layer. We employ a neareutectic approximation and consider the limit of large far-field temperature. These asymptotic limits allow us to examine the rich dynamics of the mushy layer in the form of small deviations from the classical case of convection in a horizontal porous layer of uniform permeability. Of particular interest are the effects of the asymmetries in the basic state and the non-uniform permeability in the mushy layer, which lead to transcritically bifurcating convection with hexagonal planform. We obtain a set of three coupled amplitude equations describing the evolution of small-amplitude convecting states in the mushy layer. These equations are analysed to determine the stability of and competition between two-dimensional roll and hexagonal convection patterns. We find that either rolls or hexagons can be stable. Furthermore, hexagons with either upflow or downflow at the centres can be stable, depending on the relative strengths of different physical mechanisms. We determine how to adjust the control parameters to minimize the degree of subcriticality of the bifurcation and hence render the system globally more stable. Finally, the amplitude equations reveal the presence of a new oscillatory instability.


Sign in / Sign up

Export Citation Format

Share Document