Numerical simulation of dendritic solidification with convection: Two-dimensional geometry

Nabeel Al-Rawahi, Gretar Tryggvason

Research output: Contribution to journalArticle

124 Citations (Scopus)

Abstract

A front tracking method is presented for simulations of dendritic growth of pure substances in the presence of flow. The liquid-solid interface is explicitly tracked and the latent heat released during solidification is calculated using the normal temperature gradient near the interface. A projection method is used to solve the Navier-Stokes equations. The no-slip condition on the interface is enforced by setting the velocities in the solid phase to zero. The method is validated through a comparison with an exact solution for a Stefan problem, a grid refinement test, and a comparison with a solution obtained by a boundary integral method. Three sets of two-dimensional simulations are presented: a comparison with the simulations of Beckermann et al. (J. Comput. Phys. 154, 468, 1999); a study of the effect of different flow velocities; and a study of the effect of the Prandtl number on the growth of a group of dendrites growing together. The simulations show that on the upstream side the dendrite tip velocity is increased due to the increase in the temperature gradient and the formation of side branches is promoted. The flow has the opposite effect on the downstream side. The results are in good qualitative agreement with published experimental results, even though only the two-dimensional aspects are examined here.

Original languageEnglish
Pages (from-to)471-496
Number of pages26
JournalJournal of Computational Physics
Volume180
Issue number2
DOIs
Publication statusPublished - Aug 10 2002

Fingerprint

Thermal gradients
solidification
Solidification
convection
Geometry
Latent heat
Prandtl number
Computer simulation
dendrites
geometry
Flow velocity
Navier Stokes equations
temperature gradients
simulation
boundary integral method
latent heat
Liquids
liquid-solid interfaces
Navier-Stokes equation
upstream

ASJC Scopus subject areas

  • Physics and Astronomy (miscellaneous)
  • Computer Science Applications

Cite this

Numerical simulation of dendritic solidification with convection : Two-dimensional geometry. / Al-Rawahi, Nabeel; Tryggvason, Gretar.

In: Journal of Computational Physics, Vol. 180, No. 2, 10.08.2002, p. 471-496.

Research output: Contribution to journalArticle

@article{466c405f594545e1a42964d7bdd150ff,
title = "Numerical simulation of dendritic solidification with convection: Two-dimensional geometry",
abstract = "A front tracking method is presented for simulations of dendritic growth of pure substances in the presence of flow. The liquid-solid interface is explicitly tracked and the latent heat released during solidification is calculated using the normal temperature gradient near the interface. A projection method is used to solve the Navier-Stokes equations. The no-slip condition on the interface is enforced by setting the velocities in the solid phase to zero. The method is validated through a comparison with an exact solution for a Stefan problem, a grid refinement test, and a comparison with a solution obtained by a boundary integral method. Three sets of two-dimensional simulations are presented: a comparison with the simulations of Beckermann et al. (J. Comput. Phys. 154, 468, 1999); a study of the effect of different flow velocities; and a study of the effect of the Prandtl number on the growth of a group of dendrites growing together. The simulations show that on the upstream side the dendrite tip velocity is increased due to the increase in the temperature gradient and the formation of side branches is promoted. The flow has the opposite effect on the downstream side. The results are in good qualitative agreement with published experimental results, even though only the two-dimensional aspects are examined here.",
author = "Nabeel Al-Rawahi and Gretar Tryggvason",
year = "2002",
month = "8",
day = "10",
doi = "10.1006/jcph.2002.7092",
language = "English",
volume = "180",
pages = "471--496",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "2",

}

TY - JOUR

T1 - Numerical simulation of dendritic solidification with convection

T2 - Two-dimensional geometry

AU - Al-Rawahi, Nabeel

AU - Tryggvason, Gretar

PY - 2002/8/10

Y1 - 2002/8/10

N2 - A front tracking method is presented for simulations of dendritic growth of pure substances in the presence of flow. The liquid-solid interface is explicitly tracked and the latent heat released during solidification is calculated using the normal temperature gradient near the interface. A projection method is used to solve the Navier-Stokes equations. The no-slip condition on the interface is enforced by setting the velocities in the solid phase to zero. The method is validated through a comparison with an exact solution for a Stefan problem, a grid refinement test, and a comparison with a solution obtained by a boundary integral method. Three sets of two-dimensional simulations are presented: a comparison with the simulations of Beckermann et al. (J. Comput. Phys. 154, 468, 1999); a study of the effect of different flow velocities; and a study of the effect of the Prandtl number on the growth of a group of dendrites growing together. The simulations show that on the upstream side the dendrite tip velocity is increased due to the increase in the temperature gradient and the formation of side branches is promoted. The flow has the opposite effect on the downstream side. The results are in good qualitative agreement with published experimental results, even though only the two-dimensional aspects are examined here.

AB - A front tracking method is presented for simulations of dendritic growth of pure substances in the presence of flow. The liquid-solid interface is explicitly tracked and the latent heat released during solidification is calculated using the normal temperature gradient near the interface. A projection method is used to solve the Navier-Stokes equations. The no-slip condition on the interface is enforced by setting the velocities in the solid phase to zero. The method is validated through a comparison with an exact solution for a Stefan problem, a grid refinement test, and a comparison with a solution obtained by a boundary integral method. Three sets of two-dimensional simulations are presented: a comparison with the simulations of Beckermann et al. (J. Comput. Phys. 154, 468, 1999); a study of the effect of different flow velocities; and a study of the effect of the Prandtl number on the growth of a group of dendrites growing together. The simulations show that on the upstream side the dendrite tip velocity is increased due to the increase in the temperature gradient and the formation of side branches is promoted. The flow has the opposite effect on the downstream side. The results are in good qualitative agreement with published experimental results, even though only the two-dimensional aspects are examined here.

UR - http://www.scopus.com/inward/record.url?scp=0037055751&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0037055751&partnerID=8YFLogxK

U2 - 10.1006/jcph.2002.7092

DO - 10.1006/jcph.2002.7092

M3 - Article

AN - SCOPUS:0037055751

VL - 180

SP - 471

EP - 496

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 2

ER -