Melt convection and solid phase movement play an important role in solidification processes, which significantly influence the formation of grain structures and solute segregations. In general, the melt convection and grain movement are a result of buoyancy forces. The densities within melt are different due to the variation of temperature and concentration, leading to thermally and solutally driven melt convection. Similarly, the density differences between the grains and the bulk melt cause the grain movement, leading to solid sedimentation or grain floating, as the case may be. Free, unattached solid grains are produced by partial remelting and fragmentation of dendrites, by mechanical disturbances such as stirring or vibration and by heterogeneous nucleation of grains in solidification of grain-refined alloys. In this way, movement of solid crystals during solidification can be ascertained in the following two cases. In the first case, during columnar solidification of non-grain-refined alloys, solid movement is possible in the form of dendrite fragments detached from the columnar stalks by the process of remelting and fragmentation. Movement of grains during columnar solidification gives rise to altogether different microstructure from columnar to equiaxed. In the second case, during equiaxed solidification of grain-refined alloys, the movement of solid crystals is possible in the form of equiaxed dendrite crystals nucleated due to presence of grain refiners. The rate and manner by which the free solids settle (or float) will influence macrosegregation in metal castings. Control of the solidification process is possible through an understanding of the solid movement and its effect on macrosegregation and microstructure. With this viewpoint, the overall objective of the present thesis is to study, experimentally and numerically, the phenomenon of solid phase movement during solidification. Through this study, deeper insights of the role of solid phase movement in solidification are developed which can be used for possible control of quality in castings. Both columnar and equiaxed solidification are considered. Models for transport phenomena associated with columnar solidification with solid phase movement are rarely found in the literature, because of inherent difficulty associated with consideration of microscopic features such as remelting and fragmentation. To tackle this problem, solidification modules for remelting and fragmentation are developed first, followed by integration of these molecules in a macroscopic solidification model. A Rayleigh number based fragmentation criterion is developed for detachment of dendrite fragments from the developing mushy zone, which determines the conditions favorable for fragmentation of dendrites. The criterion developed is a function of net concentration difference, liquid fraction, permeability, growth rate of mushy layer, and thermophysical properties of the material. The effect of various solidification parameters on fragmentation is highlighted. The integrated continuum model developed is applied to stimulate the solidification of aqua-ammonia system in a side-cooled rectangular cavity. The numerical results are in good qualitative agreement with those of experiments reported in literature. A gentle ramp of the mushy zone due to settling of solid crystals, as also noticed in experimental literature, is observed towards the bottom of the cavity. The influence of various modeling parameters on solid phase movement and resulting macrosegregation is investigated through a parametric study. Movement of grains during columnar solidification gives rise to altogether different microstructure and sometimes may initiate a morphological transition of the microstructure from columnar to equiaxed if the number and size of equiaxed grains ahead of the columnar front become sufficient to arrest the columnar growth. The generalised model developed, considering solid phase movement during columnar solidification is used to predict columnar-to-equiaxed transition (CET) based on a prescribed cooling rate criterion. It is found that presence of convection significantly affects the solidification behaviour. Moreover, the movement of dendrite fragments and their accumulation at the columnar front further trigger the occurrence of CET. Cooling configuration, too significantly affects the nature of CET. In unidirectional solidification cases, the locations of CET are found to be in a plane parallel to the chill face. However, for the case of the non-unidirectional solidification (as in side-cooled cavity), the locations of CET need not be in a plane parallel to the chill face. In contrast to fixed columnar solidification, equiaxed solidification is poorly understood; in particular, the phenomena associated with solid crystal movement. Movement of unattached solid crystals, formed due to heterogeneous nucleation on grain-refiners, is induced by the convective currents as well as by buoyancy effects, causing the solid to sediment or to float, depending on density of solid compared to that of the bulk melt. While moving in the bulk melt these crystals can also remelt or grow. A series of casting experiments with AI-based alloys are performed to investigate the role and influence of movement of solid crystals on macrosegregation and microstructure evolution during equiaxed solidification. Controlled experiments are designed for studying, separately, settling and floatation of equiaxed crystals for different cooling conditions and configurations. Further, these experiments are carried out in convective and non-convective cases to understand the effect of convection on solid phase movement. Temperature measurements are performed at various locations in the mould during the experiments. After the cavity is solidified, microstructural and chemical analyses of the experimental samples are carried out, several notable features are observed in temperature histories, macrosegregation pattern, and microstructures due to settling/flotation phenomenon of solid crystals. It is found that the flow behavior of solid grains has a profound influence on the progress of solidification (in terms of grain size distribution and fraction eutectic) and macrosegregation distribution. In some cases, the induced flow due to solid phase movement can cause a flow reversal. The observations and quantitative data obtained from experiments, with the help of detailed solidification conditions provided, can be used for future validations of models for equiaxed solidification. Subsequently, numerical studies are carried out, using a modified version of the macroscopic model developed for columnar solidification with motion of solid crystals, to predict the transport phenomena during equiaxed solidification. The model is applied to simulate the solidification processes corresponding to each of the experimental cases performed in this study. For a better understanding of the phenomenon of movement of solid crystals, the following two special cases of solidification are also presented: 1) without movement of solid crystals and 2) movement of solid crystals without any relative velocity between solid and liquid phases. The numerical predictions showing nature of flow field and progress of solidification are substantiated by the experimental data for the thermal analysis, qualitative microstructural Images and quantitative microstructural analysis. It is concluded, with the help of various experiments and simulations, that movement of solid crystals influences the casting quality appreciably, in terms of macrosegregation and microstructures. It is expected that the improved understanding of the role and influence of solid phase movement during solidification processes (both columnar and equiaxed) obtained through this thesis will be useful for possible control of quality of as-cast products.

The primary focus of the present work is the development of macro-models for numerical simulation of binary alloy solidification processes, consistent with microscopic phase-change considerations, with a particular emphasis on capturing the effects of non-equilibrium species redistribution on overall macrosegregation behaviour. As a first step, a generalised macroscopic framework is developed for mathematical modelling of the process. The complete set of equivalent single-phase governing equations (mass, momentum, energy and species conservation) are solved following a pressure-based Finite Volume Method according to the SIMPLER algorithm. An algorithm is also developed for the prescription of the coupling between temperature and the melt-fraction. Based on the above unified approach of solidification modelling, a macroscopic numerical model is devised that is capable of capturing the interaction between the double-diffusive convective field and a localised fluid flow on account of solutal undercooling during non-equilibrium solidification of binary alloys. Numerical simulations are performed for the case of two-dimensional transient solidification of Pb-Sn alloys, and the simulation results are also compared with the corresponding experimental results quoted in the literature. It is observed that non-equilibrium effects on account of solutal undercooling result in an enhanced macrosegregation. Next, the model is extended to capture the effects of dendritic arm coarsening on the macroscopic transport phenomena occurring during a binary alloy solidification process. The numerical results are first tested against experimental results quoted in the literature, corresponding to the solidification of an Al-Cu alloy in a bottom-cooled cavity. It is concluded that dendritic arm coarsening leads to an increased effective permeability of the mushy region as well as an enhanced eutectic fraction of the solidified ingot. Consequently, an enhanced macrosegregation can be predicted as compared to that dictated by shrinkage-induced fluid flow alone. For an order-of-magnitude assessment of predictions from the numerical models, a systematic approach is subsequently developed for scaling analysis of momentum, heat and species conservation equations pertaining to the case of solidification of a binary mixture. A characteristic velocity scale inside the mushy region is derived, in terms of the morphological parameters of the two-phase region. A subsequent analysis of the energy equation results in an estimation of the solid layer thickness. It is also shown from scaling principles that non-equilibrium effects result in an enhanced macro-segregation compared to the case of an equilibrium model For the sake of assessment of the scaling analysis, the predictions are validated against computational results corresponding to the simulation of a full set of governing equations, thus confirming the trends suggested by the scale analysis. In order to analytically investigate certain limiting cases of unidirectional alloy solidification, a fully analytical solution technique is established for the solution of unidirectional, conduction-dominated, alloy solidification problems. The results are tested for the problem of solidification of an ammonium chloride-water solution, and are compared with those from existing analytical models as well as with the corresponding results from a fully numerical simulation. The effects of different microscopic models on solidification behaviour are illustrated, and transients in temperature and heat flux distribution are also analysed. An excellent agreement between the present solutions and results from the computational simulation can be observed. The generalised numerical model is subsequently utilised to investigate the effects of laminar double-diffusive Rayleigh-Benard convection on directional solidification of binary fluids, when cooled and solidified from the top. A series of experiments is also performed with ammonium chloride-water solutions of hypoeutectic and hypereutectic composition, so as to facilitate comparisons with numerical predictions. While excellent agreements can be obtained for the first case, the second case results in a peculiar situation, where crystals nucleated on the inner roof of the cavity start descending through the bulk fluid, and finally settle down at the bottom of the cavity in the form of a sedimented solid layer. An eutectic solidification front subsequently progresses from the top surface vertically downwards, and eventually meets the heap of solid crystals collected on the floor of the cavity. However, comparison of experimental observations with corresponding numerical results from the present model is not possible under this situation, since the associated transport process involves a complex combination of a number of closely interconnected physical mechanisms, many of which are yet to be resolved. Subsequent to the development of the mathematical model and experimental arrangements for macroscopic transport processes during an alloy solidification process, some of the important modes of double-diffusive instability are analytically investigated, as a binary alloy of any specified initial composition is directionally solidified from the top. By employing a close-formed solution technique, the critical liquid layer heights corresponding to the onset of direct mode of instability are identified, corresponding two a binary alloy with three different initial compositions. In order to simulate turbulent transport during non-equilibrium solidification processes of binary alloys, a modified k-8 model is subsequently developed. Particular emphasis is given for appropriate modelling of turbulence parameters, so that the model merges with single-phase turbulence closure equations in the pure liquid region in a smooth manner. Laboratory experiments are performed using an ammonium chloride-water solution that is solidified by cooling from the top of a rectangular cavity. A good agreement between numerical and experimental results is observed. Finally, in order to study the effects of three-dimensionality in fluid flow on overall macrosegregation behaviour, the interaction between double-diffusive convection and non-equilibrium solidification of a binary mixture in a cubic enclosure (cooled from a side) is numerically investigated using a three-dimensional transient mathematical model. Investigations are carried out for two separate model systems, one corresponding to a typical metal-ally analogue system and other corresponding to an actual metal-alloy system. As a result of three-dimensional convective flow-patterns, a significant solute macrosegregation is observed in the transverse sections of the cavity, which cannot be captured by two-dimensional simulations.