Numerical Simulation of Bridgman Solidification of Binary Alloys

S Battaglioli, S McFadden, A Robinson

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

A transient 2D axisymmetric numerical model for the Bridgman solidification process for a stationary furnace and moving sample is presented. The model is able to predict the evolution of temperature and solid fraction of binary alloys in cases where buoyancy induced convection is negligible, such as in microgravity conditions. A dimensionless form of the governing equations was derived in order to identify the dimensionless parameters that characterise the process, those being the Stefan, Péclet, and Biot numbers. The problem was solved using a finite volume method and an explicit time stepping scheme. To test the efficacy of the model, simulated results were compared with experimental data from the literature and acceptable agreement was obtained. Finally, a parametric analysis was performed for understanding the influence of the process parameters on solidification. One key feature of this study was the inclusion of a term describing the advection of latent heat due to the translation of the mushy zone with varying solid fraction. This thermal transport mechanism was shown to be significant, since its magnitude was comparable to the advection of sensible heat. It was also found that when small Biot numbers were due to low values of the heat transfer coefficients at the surface of the sample, rather than to small sample radii, advective mechanisms were enhanced resulting in more convex shapes of the liquidus isotherm. This highlighted the importance of considering both axial and radial heat fluxes when describing the process.
LanguageEnglish
Pages199-211
JournalInternational Journal of Heat and Mass Transfer
Volume104
Early online date21 Aug 2016
DOIs
Publication statusPublished - Jan 2017

Fingerprint

Binary alloys
Advection
binary alloys
solidification
Biot number
Solidification
advection
Latent heat
Microgravity
Computer simulation
Finite volume method
Buoyancy
mushy zones
Heat transfer coefficients
Isotherms
Heat flux
Numerical models
Furnaces
Peclet number
simulation

Keywords

  • Materials
  • thermal engineering
  • modelling

Cite this

Battaglioli, S ; McFadden, S ; Robinson, A. / Numerical Simulation of Bridgman Solidification of Binary Alloys. 2017 ; Vol. 104. pp. 199-211.
@article{93d4b93405cc46e88b9be45fb618c37f,
title = "Numerical Simulation of Bridgman Solidification of Binary Alloys",
abstract = "A transient 2D axisymmetric numerical model for the Bridgman solidification process for a stationary furnace and moving sample is presented. The model is able to predict the evolution of temperature and solid fraction of binary alloys in cases where buoyancy induced convection is negligible, such as in microgravity conditions. A dimensionless form of the governing equations was derived in order to identify the dimensionless parameters that characterise the process, those being the Stefan, P{\'e}clet, and Biot numbers. The problem was solved using a finite volume method and an explicit time stepping scheme. To test the efficacy of the model, simulated results were compared with experimental data from the literature and acceptable agreement was obtained. Finally, a parametric analysis was performed for understanding the influence of the process parameters on solidification. One key feature of this study was the inclusion of a term describing the advection of latent heat due to the translation of the mushy zone with varying solid fraction. This thermal transport mechanism was shown to be significant, since its magnitude was comparable to the advection of sensible heat. It was also found that when small Biot numbers were due to low values of the heat transfer coefficients at the surface of the sample, rather than to small sample radii, advective mechanisms were enhanced resulting in more convex shapes of the liquidus isotherm. This highlighted the importance of considering both axial and radial heat fluxes when describing the process.",
keywords = "Materials, thermal engineering, modelling",
author = "S Battaglioli and S McFadden and A Robinson",
note = "Compliant in UIR; evidence uploaded to 'Other files'",
year = "2017",
month = "1",
doi = "10.1016/j.ijheatmasstransfer.2016.08.030",
language = "English",
volume = "104",
pages = "199--211",

}

Numerical Simulation of Bridgman Solidification of Binary Alloys. / Battaglioli, S; McFadden, S; Robinson, A.

Vol. 104, 01.2017, p. 199-211.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Numerical Simulation of Bridgman Solidification of Binary Alloys

AU - Battaglioli, S

AU - McFadden, S

AU - Robinson, A

N1 - Compliant in UIR; evidence uploaded to 'Other files'

PY - 2017/1

Y1 - 2017/1

N2 - A transient 2D axisymmetric numerical model for the Bridgman solidification process for a stationary furnace and moving sample is presented. The model is able to predict the evolution of temperature and solid fraction of binary alloys in cases where buoyancy induced convection is negligible, such as in microgravity conditions. A dimensionless form of the governing equations was derived in order to identify the dimensionless parameters that characterise the process, those being the Stefan, Péclet, and Biot numbers. The problem was solved using a finite volume method and an explicit time stepping scheme. To test the efficacy of the model, simulated results were compared with experimental data from the literature and acceptable agreement was obtained. Finally, a parametric analysis was performed for understanding the influence of the process parameters on solidification. One key feature of this study was the inclusion of a term describing the advection of latent heat due to the translation of the mushy zone with varying solid fraction. This thermal transport mechanism was shown to be significant, since its magnitude was comparable to the advection of sensible heat. It was also found that when small Biot numbers were due to low values of the heat transfer coefficients at the surface of the sample, rather than to small sample radii, advective mechanisms were enhanced resulting in more convex shapes of the liquidus isotherm. This highlighted the importance of considering both axial and radial heat fluxes when describing the process.

AB - A transient 2D axisymmetric numerical model for the Bridgman solidification process for a stationary furnace and moving sample is presented. The model is able to predict the evolution of temperature and solid fraction of binary alloys in cases where buoyancy induced convection is negligible, such as in microgravity conditions. A dimensionless form of the governing equations was derived in order to identify the dimensionless parameters that characterise the process, those being the Stefan, Péclet, and Biot numbers. The problem was solved using a finite volume method and an explicit time stepping scheme. To test the efficacy of the model, simulated results were compared with experimental data from the literature and acceptable agreement was obtained. Finally, a parametric analysis was performed for understanding the influence of the process parameters on solidification. One key feature of this study was the inclusion of a term describing the advection of latent heat due to the translation of the mushy zone with varying solid fraction. This thermal transport mechanism was shown to be significant, since its magnitude was comparable to the advection of sensible heat. It was also found that when small Biot numbers were due to low values of the heat transfer coefficients at the surface of the sample, rather than to small sample radii, advective mechanisms were enhanced resulting in more convex shapes of the liquidus isotherm. This highlighted the importance of considering both axial and radial heat fluxes when describing the process.

KW - Materials

KW - thermal engineering

KW - modelling

U2 - 10.1016/j.ijheatmasstransfer.2016.08.030

DO - 10.1016/j.ijheatmasstransfer.2016.08.030

M3 - Article

VL - 104

SP - 199

EP - 211

ER -