A new coupled model for alloy solidification

A new coupled model in the binary alloy solidification has been developed. The model is based on the cellular automaton (CA) technique to calculate the evolution of the interface governed by temperature, solute diffusion and Gibbs-Thomson effect. The diffusion equation of temperature with the release of latent heat on the solid/liquid (S/L) interface is valid in the entire domain. The temperature diffusion without the release of latent heat and solute diffusion are solved in the entire domain. In the interface cells, the

A new coupled model for alloy solidification

A new coupled model in the binary alloy solidification has been developed. The model is based on the cellular automaton (CA) technique to calculate the evolution of the interface governed by temperature, solute diffusion and Gibbs-Thomson effect. The diffusion equation of temperature with the release of latent heat on the solid/liquid (S/L) interface is valid in the entire domain. The temperature diffusion without the release of latent heat and solute diffusion are solved in the entire domain. In the interface cells, the energy and solute conservation, thermodynamic and chemical potential equilibrium are adopted to calculate the temperature, solid concentration, liquid concentration and the increment of solid fraction. Compared with other models where the release of latent heat is solved in implicit or explicit form according to the solid/liquid (S/L) interface velocity, the energy diffusion and the release of latent heat in this model are solved at different scales, i.e. the macro-scale and micro-scale. The variation of solid fraction in this model is solved using several algebraic relations coming from the chemical potential equilibrium and thermodynamic equilibrium which can be cheaply solved instead of the calculation of S/L interface velocity. With the assumption of the solute conservation and energy conservation, the solid fraction can be directly obtained according to the thermodynamic data. This model is natural to be applied to multiple (> 2) spatial dimension case and multiple (> 2) component alloy. The morphologies of equiaxed dendrite are obtained in numerical experiments.