Open Access

ARTICLE

A Meshless Approach Towards Solution of Macrosegregation Phenomena

by Gregor Kosec¹, Miha Založnik², Božidar Šarler¹, Hervé Combeau²

1 University of Nova Gorica, Vipavska 13, SI-5000 Nova Gorica, Slovenia, gregor.kosec@ung.si,bozidar.sarler@ung.si, www.ung.si
2 Institut Jean Lamour, CNRS – Nancy-Université – UPV Metz, Ecole des Mines de Nancy,Parc de Saurupt CS14234, F-54042 Nancy cedex, France. miha.zaloznik@mines.inpl-nancy.fr,herve.combeau@ijl.nancy-universite.fr, www.ijl.nancy-universite.fr

Computers, Materials & Continua 2011, 22(2), 169-196. https://doi.org/10.3970/cmc.2011.022.169

Abstract

The simulation of macrosegregation as a consequence of solidification of a binary Al-4.5%Cu alloy in a 2-dimensional rectangular enclosure is tackled in the present paper. Coupled volume-averaged governing equations for mass, energy, momentum and species transfer are considered. The phase properties are resolved from the Lever solidification rule, the mushy zone is modeled by the Darcy law and the liquid phase is assumed to behave like an incompressible Newtonian fluid. Double diffusive effects in the melt are modeled by the thermal and solutal Boussinesq hypothesis. The physical model is solved by the novel Local Radial Basis Function Collocation Method (LRBFCM). The involved physical relevant fields are represented on overlapping 5-noded sub-domains through collocation by using multiquadrics Radial Basis Functions (RBF). The involved first and second derivatives of the fields are calculated from the respective derivatives of the RBFs. The fields are solved through explicit time stepping. The pressure-velocity coupling is calculated through a local pressure correction scheme. The evolution of the solidification process is presented through temperature, velocity, liquid fraction and species concentration histories in four sampling points. The fully solidified state is analyzed through final macrosegregation map in three vertical and three horizontal cross-sections. The results are compared with the classical Finite Volume Method (FVM). A surprisingly good agreement of the numerical solution of both methods is shown and therefore the results can be used as a reference for future verification studies. The advantages of the represented meshless approach are its simplicity, accuracy, similar coding in 2D and 3D, and straightforward applicability in non-uniform node arrangements. The paper probably for the first time shows an application of a meshless method in such a highly non-linear and multi-physics problem.

---

Keywords

---