Please wait a minute...
Acta Metall Sin  2017, Vol. 53 Issue (11): 1521-1531    DOI: 10.11900/0412.1961.2017.00083
Orginal Article Current Issue | Archive | Adv Search |
Numerical Simulation of Liquid-Solid Conversion Affecting Flow Behavior During Casting Filling Process
Liu CAO, Fei SUN(), Tao CHEN, Zihao TENG, Yulong TANG, Dunming LIAO
State Key Laboratory of Materials Processing and Die & Mould Technology, Huazhong University of Science and Technology, Wuhan 430074, China
Cite this article: 

Liu CAO, Fei SUN, Tao CHEN, Zihao TENG, Yulong TANG, Dunming LIAO. Numerical Simulation of Liquid-Solid Conversion Affecting Flow Behavior During Casting Filling Process. Acta Metall Sin, 2017, 53(11): 1521-1531.

Download:  HTML  PDF(6335KB) 
Export:  BibTeX | EndNote (RIS)      
Abstract  

Misrun and cold shut are common defects in casting productions, which could make surface accuracy of castings poorer, even leading to cracking and casting scraps in them. The formation process of misrun and cold shut is hard to be observed directly only by experiment measures, since casting filling process is in a state of high temperature flow inside mold. The key to predict the defects accurately is the way to handle the effect of liquid-solid conversion on flow behavior. On the basis of existing methods for treating liquid-solid conversion, a calculation model of mushy region flow behavior through measurement of solid-fraction is developed, which can effectively investigate the flow behavior of mushy region in different stages. Generally, the critical solid-fraction method is adopted for mushy region with high solid-fraction, in consideration of that only the speed of high solid-fraction region is supposed to be zero during casting filling process. The variable viscosity method is applied for mushy region with low solid-fraction, due to casting filling process being unlikely to form toothpaste-like flow. However, the porous medium drag-based model is used for mushy region with middle solid-fraction, because only the middle solid-fraction region can be equivalent to porous medium. Combining the above three methods, a flow-field calculation program considering the effect of liquid-solid conversion on flow behavior during casting filling process is developed, in which finite volume method (FVM) is included for discretization equations; the pressure implicit with splitting of operator (PISO) algorithm is added for coupling pressure and velocity; the volume of fluid (VOF) algorithm is also combined for interface tracking. An numerical simulation of water-filled S-shaped channel is performed in the case of taking no account of liquid-solid conversion, and the simulated results coincide better with the experimental results, which certifies for its accuracy as an adopted model. Since the bottom filling casting craft is commonly used in single-shape casting, a comparison between the calculated results obtained using other single models and those using this model at different control parameters, is needed. The better agreement between them indicates that this new model is appropriate for calculating the flow behavior in mushy region.

Key words:  liquid-solid conversion      mushy region      porous medium      variable viscosity      finite volume method      numerical simulation     
Received:  15 March 2017     
ZTFLH:  TG245  
Fund: Supported by New Century Excellent Talents in University (No.NCET-13-0229) and National Science & Technology Key Projects of Numerical Control (No.2012ZX04010-031)

URL: 

https://www.ams.org.cn/EN/10.11900/0412.1961.2017.00083     OR     https://www.ams.org.cn/EN/Y2017/V53/I11/1521

Fig.1  Schematics of the evolution process of dendrite morphology in fluid liquid
(a) starting stage of dendrite growth (b) middle stage of dendrite growth (c) final stage of dendrite growth
Fig.2  Calculation flow chart used for the pressure implicit with splitting of operator (PISO) algorithm
Parameter Value Unit
Water density 1000 kgm-3
Air density 1 kgm-3
Water dynamic viscosity 1×10-3 Pas
Air dynamic viscosity 1×10-5 Pas
Water-air surface tension coefficient 0.07275 Nm-1
Acceleration of gravity {0, 0, -9.8} ms-2
Inlet velocity {0, 8.7, 0} ms-1
Outlet pressure 0 Pa
Table 1  Parameters used in simulation of S-shaped channel
Fig.3  Comparisons between experimental results (a~c) and simulated results of volume fraction of water (d~f) and velocity of liquid-gas phase (g~i) in S-shaped channel filled with water at 7.15 ms (a, d, g), 25.03 ms (b, e, h) and 39.34 ms (c, f, i)
Fig.4  Geometric and mesh model of bottom filling casting craft with simple shape
Parameter Value Unit
Aluminum alloy density 2385 kgm-3
Air density 1 kgm-3
Aluminum alloy dynamic viscosity 0.003 Pas
Air dynamic viscosity 1×10-5 Pas
Aluminum alloy-air surface tension coefficient 0.871 Nm-1
Acceleration of gravity {0, 0, -9.8} ms-2
Inlet velocity {0, 0, 0.1} ms-1
Outlet pressure 0 Pa
Liquidus temperature of aluminum alloy 660
Solidus temperature of aluminum alloy 560
Latent heat of aluminum alloy 300 Jg-1
Inlet temperature 670
Mold temperature 30
Heat transfer coefficient between lower surface and mold 6000 Wm-2K-1
Heat transfer coefficient between other surfaces and mold 100 Wm-2K-1
Table 2  Parameters used in simulation of bottom filling casting craft
Fig.5  Comparisons of simulated results of solid-fraction of aluminum phase (a~c), velocity of liquid-gas phase (d~f) and volume fraction of aluminum phase (g~i) under critical solid-fractions of 0.5 (a, d, g), 0.75 (b, e, h) and 0.95 (c, f, i) at 5 s Color online
Fig.6  Comparisons of simulated results of solid-fraction of aluminum phase (a~c), velocity of liquid-gas phase (d~f) and volume fraction of aluminum phase (g~i) with incremental multiples as 1 (a, d, g), 10 (b, e, h) and 100 (c, f, i) at 5 s (The incremental multiple means the viscosity difference between solid and liquid phase) Color online
Fig.7  Comparisons of simulated results of solid-fraction of aluminum phase (a~c), velocity of liquid-gas phase (d~f) and volume fraction of aluminum phase (g~i) with porous medium drag coefficients of 1 (a, d, g), 10 (b, e, h) and 100 (c, f, i) at 5 s
Fig.8  Simulated results of solid-fraction of aluminum phase (a~c), velocity of liquid-gas phase (d~f) and volume fraction of aluminum phase (g~i) at different times of 1.5 s (a, d, g), 3 s (b, e, h) and 5 s (c, f, i) by using the calculation model for mushy region flow behavior through measurement of solid-fraction
[1] Rajkolhe R, Khan J G.Defects, causes and their remedies in casting process: A review[J]. Int. J. Res. Advent. Technol., 2014; 2: 375
[2] Dabade U A, Bhedasgaonkar R C.Casting defect analysis using design of experiments (DoE) and computer aided casting simulation technique[J]. Procedia CIRP, 2013, 7: 616
[3] Liu D R, Yang Z P, Wang L P, et al.Development of simulation of mould filling during casting: A review[J]. J. Harbin Univ. Sci. Technol., 2016, 21(3): 96(刘东戎, 杨智鹏, 王丽萍等. 铸造充型过程数值模拟技术的发展及现状评述[J]. 哈尔滨理工大学学报, 2016, 21(3): 96)
[4] Vazquez V, Juarez-Hernandez A, Mascarenas A, et al.Cold shut formation analysis on a free lead yellow brass tap[J]. Eng. Fail. Anal., 2010, 17: 1285
[5] Lee J H, Won C W, Cho S S, et al.Effects of melt flow and temperature on the macro and microstructure of scroll compressor in direct squeeze casting[J]. Mater. Sci. Eng., 2000, A281: 8
[6] Lewis R W, Ravindran K.Finite element simulation of metal casting[J]. Int. J. Numer. Meth. Eng., 2000, 47: 29
[7] Jin Z M, He J C, Xu G J.Numerical simulation of flow, temperature and thermal stress fields during twin-roll casting process[J]. Acta Metall. Sin., 2000, 36: 391(金珠梅, 赫冀成, 徐广?. 双辊连续铸轧工艺中流场、温度场和热应力场的数值计算[J]. 金属学报, 2000, 36: 391)
[8] Yoshizawa A, Nisizima S.A nonequilibrium representation of the turbulent viscosity based on a two-scale turbulence theory[J]. Phys. Fluids Fluid Dynam., 1993, 5A: 3302
[9] Neuman S P.Theoretical derivation of Darcy's law[J]. Acta Mech., 1977, 25: 153
[10] Whitaker S.Flow in porous media I: A theoretical derivation of Darcy's law[J]. Transport Porous Med., 1986, 1: 3
[11] Jana S, K?ttlitz O, Hediger F, et al.Predictions of misruns using three-phase coupled mold-filling and solidification simulations in low pressure turbine (LPT) blades [A]. IOP Conference Series: Materials Science and Engineering [C]. Bristol: IOP Publishing, 2012: 1
[12] Wang C H, Hu H J, Luo J.Computer simulation of investment casting based on procast software[J]. Foundry Technol., 2007, 28: 1360(王春欢, 胡红军, 罗静. 基于Procast软件的熔模铸造计算机模拟[J]. 铸造技术, 2007, 28: 1360)
[13] Frehse J, Málek J, Steinhauer M.On analysis of steady flows of fluids with shear-dependent viscosity based on the Lipschitz truncation method[J]. SIAM J. Math. Anal., 2003, 34: 1064
[14] Arnberg L, B?ckerud L, Chai G.Solidification Characteristics of Aluminum Alloys: Dendrite Coherency[M]. Schaumburg: American Foundrymen's Society, 1996: 30
[15] Carman P C.Fluid flow through granular beds[J]. Chem. Eng. Res. Des., 1937, 75: S32
[16] Voile V R, Prakash C.A fixed grid numerical modeling methodology for convection diffusion mushy region phase-change problems[J]. Int. J. Heat Mass Transfer, 1987, 30: 1709
[17] Mitchell A R, Griffiths D F.The Finite Difference Method in Partial Differential Equations[M]. Hoboken: John Wiley & Sons, 1980: 35
[18] Dhatt G, Lefran?ois E, Touzot G.Finite Element Method[M]. Hoboken: John Wiley & Sons, 2013: 40
[19] Versteeg H K, Malalasekera W.An Introduction to Computational Fluid Dynamics: The Finite Volume Method[M]. New York: Pearson Education, 2007: 45
[20] Vladimir G.Finite-difference methods for simulating the solidification of castings[J]. Mater. Tech., 2009, 43: 233
[21] Cao L, Liao D M, Lu Y Z, et al.Heat transfer model of directional solidification by LMC process for superalloy casting based on finite element method[J]. Metall. Mater. Trans., 2016, 47A: 4640
[22] Kim J, Kim D, Choi H.An immersed-boundary finite-volume method for simulations of flow in complex geometries[J]. J. Comput. Phys., 2001, 171: 132
[23] Patankar S V, Spalding D B.A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows[J]. Int. J. Heat Mass Transfer, 1972, 15: 1787
[24] Issa R I, Gosman A D, Watkins A P.The computation of compressible and incompressible recirculating flows by a non-iterative implicit scheme[J]. J. Comput. Phys., 1986, 62: 66
[25] Bridson R, Houriham J, Nordenstam M.Curl-noise for procedural fluid flow[J]. ACM Trans. Graphic., 2007, 26: 46
[26] Ménard T, Tanguy S, Berlemont A.Coupling level set/VOF/ghost fluid methods: Validation and application to 3D simulation of the primary break-up of a liquid jet[J]. Int. J. Multiphas. Flow, 2007, 33: 510
[27] Wu S S, Liu Y Q.Principle of Material Forming [M]. 2nd Ed., Beijing: China Machine Press, 2008: 29(吴树森, 柳玉起. 材料成形原理 [M]. 第2版. 北京: 机械工业出版社, 2008: 29)
[28] Hirt C W, Nichols B D.Volume of fluid (VOF) method for the dynamics of free boundaries[J]. J. Comput. Phys., 1981, 39(1): 201
[29] Constantin P, Foias C.Navier-Stokes Equations [M]. Chicago: University of Chicago Press, 1988: 50
[30] Mathur S R, Murthy J Y.Pressure boundary conditions for incompressible flow using unstructured meshes[J]. Numer. Heat Transfer, 1997, 32B: 283
[31] Park T S.Effects of time-integration method in a large-eddy simulation using the PISO Algorithm: Part I—flow field[J]. Numer. Heat Transfer, 2006, 50A: 229
[32] Tavakoli R, Babaei R, Varahram N, et al.Numerical simulation of liquid/gas phase flow during mold filling[J]. Comput. Methods Appl. Mech. Eng., 2006, 196: 697
[33] Pang S Y, Chen L L, Zhang M C, et al.Numerical simulation two phase flows of casting filling process using SOLA particle level set method[J]. Appl. Math. Modell., 2010, 34: 4106
[1] BI Zhongnan, QIN Hailong, LIU Pei, SHI Songyi, XIE Jinli, ZHANG Ji. Research Progress Regarding Quantitative Characterization and Control Technology of Residual Stress in Superalloy Forgings[J]. 金属学报, 2023, 59(9): 1144-1158.
[2] WANG Chongyang, HAN Shiwei, XIE Feng, HU Long, DENG Dean. Influence of Solid-State Phase Transformation and Softening Effect on Welding Residual Stress of Ultra-High Strength Steel[J]. 金属学报, 2023, 59(12): 1613-1623.
[3] ZHANG Kaiyuan, DONG Wenchao, ZHAO Dong, LI Shijian, LU Shanping. Effect of Solid-State Phase Transformation on Stress and Distortion for Fe-Co-Ni Ultra-High Strength Steel Components During Welding and Vacuum Gas Quenching Processes[J]. 金属学报, 2023, 59(12): 1633-1643.
[4] XIA Dahai, DENG Chengman, CHEN Ziguang, LI Tianshu, HU Wenbin. Modeling Localized Corrosion Propagation of Metallic Materials by Peridynamics: Progresses and Challenges[J]. 金属学报, 2022, 58(9): 1093-1107.
[5] HU Long, WANG Yifeng, LI Suo, ZHANG Chaohua, DENG Dean. Study on Computational Prediction About Microstructure and Hardness of Q345 Steel Welded Joint Based on SH-CCT Diagram[J]. 金属学报, 2021, 57(8): 1073-1086.
[6] LI Zihan, XIN Jianwen, XIAO Xiao, WANG Huan, HUA Xueming, WU Dongsheng. The Arc Physical Characteristics and Molten Pool Dynamic Behaviors in Conduction Plasma Arc Welding[J]. 金属学报, 2021, 57(5): 693-702.
[7] WANG Fuqiang, LIU Wei, WANG Zhaowen. Effect of Local Cathode Current Increasing on Bath-Metal Two-Phase Flow Field in Aluminum Reduction Cells[J]. 金属学报, 2020, 56(7): 1047-1056.
[8] LIU Jizhao, HUANG Hefei, ZHU Zhenbo, LIU Awen, LI Yan. Numerical Simulation of Nanohardness in Hastelloy N Alloy After Xenon Ion Irradiation[J]. 金属学报, 2020, 56(5): 753-759.
[9] WANG Bo,SHEN Shiyi,RUAN Yanwei,CHENG Shuyong,PENG Wangjun,ZHANG Jieyu. Simulation of Gas-Liquid Two-Phase Flow in Metallurgical Process[J]. 金属学报, 2020, 56(4): 619-632.
[10] XU Qingyan,YANG Cong,YAN Xuewei,LIU Baicheng. Development of Numerical Simulation in Nickel-Based Superalloy Turbine Blade Directional Solidification[J]. 金属学报, 2019, 55(9): 1175-1184.
[11] Peiyuan DAI,Xing HU,Shijie LU,Yifeng WANG,Dean DENG. Influence of Size Factor on Calculation Accuracy of Welding Residual Stress of Stainless Steel Pipe by 2D Axisymmetric Model[J]. 金属学报, 2019, 55(8): 1058-1066.
[12] ZHANG Qingdong, LIN Xiao, LIU Jiyang, HU Shushan. Modelling of Q&P Steel Heat Treatment Process Based on Finite Element Method[J]. 金属学报, 2019, 55(12): 1569-1580.
[13] LU Shijie, WANG Hu, DAI Peiyuan, DENG Dean. Effect of Creep on Prediction Accuracy and Calculating Efficiency of Residual Stress in Post Weld Heat Treatment[J]. 金属学报, 2019, 55(12): 1581-1592.
[14] Jun LI, Mingxu XIA, Qiaodan HU, Jianguo LI. Solutions in Improving Homogeneities of Heavy Ingots[J]. 金属学报, 2018, 54(5): 773-788.
[15] Xinhua LIU, Huadong FU, Xingqun HE, Xintong FU, Yanqing JIANG, Jianxin XIE. Numerical Simulation Analysis of Continuous Casting Cladding Forming for Cu-Al Composites[J]. 金属学报, 2018, 54(3): 470-484.
No Suggested Reading articles found!