TY - JOUR
T1 - A fully coupled finite element formulation for liquid–solid–gas thermo-fluid flow with melting and solidification
AU - Yan, J.
AU - Yan, W.
AU - Lin, S.
AU - Wagner, G. J.
N1 - Funding Information:
The support of this research by the National Science Foundation under CPS-1646592 for J. Yan and G.J.Wagner and a Graduate Research Fellowship under DGE-1324585 for S.E. Lin is gratefully acknowledged.
Publisher Copyright:
© 2018 Elsevier B.V.
PY - 2018/7/1
Y1 - 2018/7/1
N2 - Many important industrial processes, such as additive manufacturing, involve rapid mass, flow and heat transport between gas, liquid and solid phases. Various associated challenges, such as the large density ratio between gas and condensed phases, make accurate, robust thermal multi-phase flow simulations of these processes very difficult. In order to address some of the associated challenges, a computational framework for thermal multi-phase flows is developed based on the finite element method (FEM). A unified model for thermal multi-phase flows similar to the models widely used in the manufacturing community is adopted. The combination of the level-set method and residual-based variational multi-scale formulation (RBVMS) is used to solve the governing equations of thermal multi-phase flows. Phase transitions between solid and liquid phases, i.e., melting and solidification, are considered. Interfacial forces, including surface tension and Marangoni stress, are taken into account and handled by a density-scaled continuum surface force model. A robust fully coupled solution strategy is adopted to handle various numerical difficulties associated with thermal multi-phase flow simulations, and implemented by means of a matrix-free technique using Flexible GMRES. The mathematical formulation and its algorithmic implementation are described in detail. Four numerical test cases are presented to demonstrate the capability of the proposed formulation. The first case is a benchmark example of solidification of aluminum in a graphite mold, the second case is a thermo-capillary droplet migration problem, the third case is a spot laser melting problem, and the fourth case is the melting of metal with an interior gas bubble. The computational results are compared with analytical, experimental and simulation data from other researchers, with good agreement in cases where such data is available.
AB - Many important industrial processes, such as additive manufacturing, involve rapid mass, flow and heat transport between gas, liquid and solid phases. Various associated challenges, such as the large density ratio between gas and condensed phases, make accurate, robust thermal multi-phase flow simulations of these processes very difficult. In order to address some of the associated challenges, a computational framework for thermal multi-phase flows is developed based on the finite element method (FEM). A unified model for thermal multi-phase flows similar to the models widely used in the manufacturing community is adopted. The combination of the level-set method and residual-based variational multi-scale formulation (RBVMS) is used to solve the governing equations of thermal multi-phase flows. Phase transitions between solid and liquid phases, i.e., melting and solidification, are considered. Interfacial forces, including surface tension and Marangoni stress, are taken into account and handled by a density-scaled continuum surface force model. A robust fully coupled solution strategy is adopted to handle various numerical difficulties associated with thermal multi-phase flow simulations, and implemented by means of a matrix-free technique using Flexible GMRES. The mathematical formulation and its algorithmic implementation are described in detail. Four numerical test cases are presented to demonstrate the capability of the proposed formulation. The first case is a benchmark example of solidification of aluminum in a graphite mold, the second case is a thermo-capillary droplet migration problem, the third case is a spot laser melting problem, and the fourth case is the melting of metal with an interior gas bubble. The computational results are compared with analytical, experimental and simulation data from other researchers, with good agreement in cases where such data is available.
KW - Gas–liquid–solid systems
KW - Manufacturing
KW - Marangoni stress
KW - Melting
KW - Surface tension
KW - Thermal multi-phase flows
UR - http://www.scopus.com/inward/record.url?scp=85044764549&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85044764549&partnerID=8YFLogxK
U2 - 10.1016/j.cma.2018.03.017
DO - 10.1016/j.cma.2018.03.017
M3 - Article
AN - SCOPUS:85044764549
VL - 336
SP - 444
EP - 470
JO - Computer Methods in Applied Mechanics and Engineering
JF - Computer Methods in Applied Mechanics and Engineering
SN - 0374-2830
ER -