研究论文

铸造充型过程中液固转变影响流动行为的数值计算

  • 曹流 ,
  • 孙飞 ,
  • 陈涛 ,
  • 滕子浩 ,
  • 唐玉龙 ,
  • 廖敦明
展开
  • 华中科技大学材料成形及模具技术国家重点实验室 武汉 430074

作者简介 曹 流,男,1991年生,博士生

收稿日期: 2017-03-15

  网络出版日期: 2017-07-13

基金资助

教育部新世纪优秀人才支持计划项目No.NCET-13-0229和国家数控重大专项项目No.2012ZX04010-031

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
Expand
  • State Key Laboratory of Materials Processing and Die & Mould Technology, Huazhong University of Science and Technology, Wuhan 430074, China

Received date: 2017-03-15

  Online published: 2017-07-13

Supported by

Supported by New Century Excellent Talents in University (No.NCET-13-0229) and National Science & Technology Key Projects of Numerical Control (No.2012ZX04010-031)

摘要

为准确预测浇不足及冷隔,在已有的处理液固转变方法基础上,提出基于固相率变化的糊状区流动行为计算模型,该模型可以有效地处理液固转变过程中糊状区不同阶段的流动行为,即高固相率糊状区采用临界固相率方法,低固相率糊状区采用变黏度方法,中等固相率糊状区采用多孔介质拖拽模型。模拟了S型铸型水模拟实验,模拟结果与实验结果吻合很好,验证了不考虑液固转变时所采用模型的准确性。针对简单形状的底注式铸造工艺,对比分析了处理液固转变过程中采用不同控制参数的计算效果,证明了糊状区流动行为计算模型的合理性。

本文引用格式

曹流 , 孙飞 , 陈涛 , 滕子浩 , 唐玉龙 , 廖敦明 . 铸造充型过程中液固转变影响流动行为的数值计算[J]. 金属学报, 2017 , 53(11) : 1521 -1531 . DOI: 10.11900/0412.1961.2017.00083

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.

参考文献

[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
文章导航

/