用户名: 密码: 验证码:
铸造充型过程气液两相流动数值模拟的研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
气液两相流动数值模拟技术是近几年来计算流体力学(CFD)领域研究的热点和前沿。本研究将计算流体力学最新方法,应用于铸造充型过程气液两相流动的数值模拟,对于深化铸造仿真技术和提高预测水平具有重要的理论意义和实际应用价值。
     精确描述相间运动界面的变化是模拟铸造充型过程气液两相流动的前提和基础。本研究采用新的界面追踪技术Level Set方法描述铸造气液两相界面,并利用高精度离散格式求解界面输运方程。在此基础上,发展了新的适用于铸造充型过程气液两相流动的SOLA-Level Set算法。该算法改进了传统的SOLA方法使其具备求解可压缩气体运动方程的能力,使用Ghost Fluid方法解决了气液两相界面处密度跨越较大易引起的计算发散问题。
     为了验证本研究算法模拟铸造气液两相流动的能力,采用了实验验证与软件验证相结合的方法。利用水力学相似原理设计了水模拟实验方案,实验与计算结果基本一致,表明该算法可以准确处理气液两相流动过程中的复杂界面变化,能够适用于密度跨越较大的气液两相流动计算。利用现有商业化软件FLUENT与自行开发系统进行模拟对比,验证了自行开发系统的正确性和可靠性。
     高温液态金属的充型过程总是伴随着热量的散失,欲准确模拟铸造气液两相流动过程,流动场与温度场耦合计算是必需的。本文建立了铸造气液两相流动场与温度场耦合计算的数学模型,编制了铸造气液两相流动传热耦合软件。
     为了检验本文算法在三维复杂边界条件下的稳定性和可靠性,并分析两相流计算与单相流计算的异同,采用国内知名铸造商业软件华铸CAE进行计算对比。分别采用华铸CAE和本文算法对实际铸件的生产工艺进行了模拟分析,对比两种计算的充型过程,表明气液两相流动算法能够准确模拟形状较为复杂的实际铸件的充型过程,该算法不仅反映了型腔背压对液体流动的影响,而且可以定量预测卷气形成过程,比单相流计算更能准确再现铸造充型过程的复杂现象。本文提出的气液两相流动数值模拟方法对铸造生产具有较强的理论意义和实用价值。
The gas-liquid two-phase flow (GLTF) simulation is an advancing research on Computational Fluid Dynamics (CFD). Based on the application of the latest technique of CFD, the simulation of GLTF during the mold filling process is developed in this research, which has a great meaning on improving casting simulation technique.
     Precisely describing the interface movement is the key and basis to the GLTF simulation during mold filling process. Level Set method, which adopts high order discrete schemes, is a new interface tracking method developed in recent years. In this study, Level Set is introduced to simulate the interface movement between the gas and the liquid. And the GLTF during mold filling process is modeled by the combination of SOLA and Level Set, in which the conventional algorithm SOLA is modified to calculate the compressible gas movement, and Ghost Fluid Method is adopted to deal with the problem of large jump in density across the interfaces.
     The validity of the presented method is verified by experiments and commercial CFD software. Firstly, water analog experiments based on similarly laws are carried out. Comparison between experimental and simulation results shows that the self-developed method is efficient to deal with complex interfaces movement, and also suitable for the calculation with large density jump. And then, the simulation results are compared with that by the commercial CFD software FLUENT, which confirms that this study has a good reliability.
     The casting process involves coupled effects of fluid flow and heat transfer. In order to simulate the mold filling process properly, heat transfer during the fluid flow must be taken into account. Therefore, a numerical simulation program coupling fluid field and temperature field of GLTF during casting process has been developed.
     Finally, using the commercial casting software InteCast CAE to simulate the complex 3-D castings, and compare the results with that of self-developed method, the reliability and validity of this study can be further tested. It can be found that the method in this study can correctly simulate the complex filling process of 3-D castings, and furthermore, the effect of backpressure on the liquid metal, and also the air entrapment during the filling process, can be quantificationally simulated in this study. Therefore, the presented method for GLTF has a profound meaning on casting production.
引文
[1]尚俊玲,陈维平,朱权利.中国铸造行业可持续化发展的思考.材料导报,2007,5:1-20
    [2]张光跃,熊守美,柳百成.铸件充型过程数值模拟的实验验证及工程应用.铸造,2000,49(10):775-779
    [3]赵志维,高永,张虎.低压铸造中充型过程的水模拟研究.特种铸造及有色合金,2007,27(5):351-353
    [4]J.H.Kuo,W.S.Hwang.Development of an interactive simulation system for die cavity filling and its application to the operation of a low-pressure casting process.Modelling and Simulation in Materials Science and Engineering,2000,8(4):583-602
    [5]W.B.Lee,H.Y.Lu.Modeling of air back pressure in die-casting dies.Journal of Matedals Processing Technology,1999,91:264-269
    [6]E.Attar,P.Homayonifar,R.Babaei.Modelling of a/r pressure effects in casting moulds.Modelling and Simulation in Materials Science and Engineering,2005,13:903-917
    [7]E.S.Duff,T.T.Nguyen.Modelling of die cavity filling and solidification in low pressure die casting.The 17th NADCA International Die Casting Congress,Australia,1993,10:113-121
    [8]G.Backer,M.Ranganathan,J.Heimsch.Simulation of flow-induced gas entrapment and its effect on porosity in aluminum die castings.North American Die Casting Association Transactions,Indianapolis,USA,2001,5:165-176
    [9]J.Prindiville,S.G.Lee,A.Gokhale.An application of trapped-air analysis to large complex high-pressure magnesium casting.Computer Methods in Applied Mechanics and Engineering,2001,124(3):677-693
    [10]P.Homayonifar,R.Babaei,E.Attar.Numerical modeling of splashing and air entrapment in high-pressure die casting.The International Journal of Advanced Manufacturing Technology,2007,17(7):546-559
    [11]T.S.David,C.W.Hirt.Modeling entrainment of air at turbulent free surfaces.World Water and Environmental Resources Congress 2004,Salt Lake City,USA,2004,9:559-574
    [12]T.Rohallah,B.Reza,V.Naser.Numerical simulation of liquid/gas phase flow during mold filling.Computer Methods in Applied Mechanics and Engineering,2006,196(1-3):697-713
    [13]李德元,徐国荣,水鸿寿等.二维非定常流体力学数值方法.北京:科学出版社,1998
    [14]B.J.Daly,W.E.Pracht.Numerical study of density-current surges.Physics of Fluid,1968,11:15-30
    [15]D.Jude,G.Tryggvason.Numerical simulations of phase change in microgravity.Heat Transfer in Microgravity Systems,1996,332:33-44
    [16]汪德瓘.高等计算水力学.南京:河海大学出版社,1999
    [17]Y.L.Zhu.Three dimensional nonlinear numerical model with inclined pressure for saltwater in intrusion at the Yangtze River estuary.Hydrodynamic,2000,1:57-66
    [18]L.Balas,E.Ozhan.An implicit three-dimensional numerical model to simulate transport process in coastal water bodies.International Journal for Numerical Methods inFluids,2000,34:307-339
    [19]F.T.HaPPer.Simulation of liquid drop breakup behavior in a flow field using diserete element techniques.Proceedings of the 5th Intenational Meeting on Reactor Thermal Hydraulics,Salt Lake City,USA,1992,9:605-611
    [20]F.H.Harlow,J.E.Welch.Numerical study of large-amplitude free-surface motions.Physics of Fluid,1966,9..842-851
    [21]刘儒勋,舒其望.计算流体力学的若干新方法.北京:科学出版社,2004
    [22]刘全,水鸿寿,张晓轶.数值模拟界面流方法进展.力学进展,2002,32(2):259-274
    [23]S.Osher,R.Fedkiw.Level set methods:An overview and some recent results.Journal of Computational Physics,2001,169:463-502
    [24]F.Moukalled,M.Darwish.A pressure-based algorithm for multi-phase flow at all Speeds.Journal of Computational Physics,2003,190(2):550-571
    [25]C.W.Hirt,B.D.Nichols.Volume of fluid(VOF) method for the dynamics of free boundaries.Journal of Computational Physics,1981,39(1):201-225
    [26]J.Strain.Semi-lagrangian methods for Level Set equations.Journal of Computational Physics,1999,151:498-533
    [27]J.Strain.Tree methods for moving interfaces.Journal of Computational Physics,1999,151:616-648
    [28]J.A.Sethian.Fast marching methods.SIAM Review,1999,41(2):199-235
    [29]O.Ubbink,R.I.Issa.A method for capturing sharp fluid interfaces on arbitrary meshes.Journal of Computational Physics,1999,153:26-50
    [30]S.Koshizuka,H.TamakoH.Moving-particle semi-implicit method.Nuclear Seicnce Engineering,1996,123:421-434
    [31]N.Ashgriz,J.Y.Poo.FLAIR:Flux line segment model for advection and interface reconstruction.Journal of Computational Physics,1991,39:449-468
    [32]S.Osher,J.A.Sethian.Fronts Propagating with curvature dependent speed:algorithms based on a Hamilton-Jacobi formulation.Journal of Computational Physics,1988,79(1):12-49
    [33]S.Soukane,F.Trochu.Application of the Level Set method to the simulation of resin transfer molding.Composites Science and Technology,2006,66(7-8):1067-1080
    [34]K.Fukagata,N.Kasagi.Simplex free adaptive tree fast sweeping and evolution methods for solving Level Set equations in arbitrary dimension.International Journal of Heat and Fluid Flow,2007,28(2):72-82
    [35]李会雄,杨冬.Level Set方法及其在两相流数值模拟研究中的应用.工程热物理学报,2001,22(2):17-21
    [36]陈菲,张梦萍,徐胜利.运动激波和气泡串相互作用的初步数值模拟.计算物理,2004,21(5):443-448.
    [37]袁德奎,陶建华.用Levd Set方法求解具有自由面的流动问题.力学学报,2000,32(3):264-271
    [38]封建湖,蔡力.基于Levi Set的多维双曲守恒律标量方程的激波追踪法.应用力学学报,2005,22(1):17-21
    [39]袁政强,祝家麟.边界元中三维重调和方程的精确积分法.重庆大学学报,2002,25f4):43-46
    [40]J.T.Batina.A Gridless Euler/Navier-Stokes solution algorithm for complex aircraft applications.International Jounal of Numerical Method in Fluids,2003,26:151-169
    [41]谷汉斌,李炎保,李邵武.界面追踪的Levd Sa和Partide Levd Sa方法.水动力学研究与进展A辑,2005,20(2):152-160
    [42]陶建华,谢伟松.用Level Set方法计算溃坝波的传播过程.水利学报,1999,10:17-22.
    [43]Y.H.Tsai.Rapid and accurate computation of the distance function using grids.Intenational Journal of Numerical Methods Engineering,2004,17..29"--35
    [44]李俊,杨新,施鹏飞.对水平集方法鲁棒初始化的双向快速步进法.系统仿真学报,2001,13:163-167
    [45]T.Belytsehko,Y.Y.Lu,L.Gu.Element-free Galerkin methods.Intenational Journal of Numerical Methods Engineering,1994,37:229-240
    [46]B.J.Daly.Numerical fluid dynamics using the particle-and-force method.Journal of Computational Physics,2001,76:215-241
    [47]R.Malladi,J.A.Sethian,B.C.Vemuri.Shape modeling with front propagation:A Level Set approach.IEEE Transactions on PAMI,1995,17(2):158-175
    [48]L.Alvarez,P.L.Lions,J.M.Morel.Image selective smoothing and edge detection by nonlinear diffusion.Journal on Numerical Analysis,1992,29(3):845-866
    [49]N.Paragios,R.Deriehe.Geodesic active contours and level set for the detection and tracking of moving objects.IEEE Transactions on Pattern Analysis and Machine Intelligence,2000,22(3):266-280
    [50]M.Bertalmio,G.Sapiro,G.Randall.Region tracking on level-set methods.IEEE Transactions on Medical Imaging,1999,18(5):448-451
    [51]A.R.Masouri,B.Sirivong,J.Kortrad.Multiple motion segmentation with level sets.Proceedings of the SPIE,2000,39(74):584-595
    [52]鲁钟琪.两相流与沸腾传热.北京:清华大学出版社,2004
    [53]S.Koshizuka,Y.Oka.Nmerical analysis of breaking waves using the moving particle semi-implicit method.International Jounal of Numerical Method in Fluids,1998,26:751-769
    [54]H.GotohH,T.Sakai.Lagrangian simulation of breaking waves using particle method.Coastal Engingeering Jounal,1999,41(3):303-326
    [55]郭烈锦.两相与多相流动力学.西安:西安交通大学出版社,2005
    [56]C.G.Speziale.Turbulence modeling for time dependent RANS and VLES:A review.AIAAJ.,1998,36(2):173-184
    [57]R.A.Gingord.Kenel estimates as a basis for general particle methods in hydrodynamics.Journal of Computational Physics,1982,46:429-441
    [58]S.Koshizuka,H.Tamako,Y.Oka.A Particle method for incompressible viscous flow with fluid fragmentation.Computational Fluid Dynamics,1995,4(1):29-46
    [59]J.F.Sini,I.Dekeyser.Numerical prediction of turbulent plane buoyant jets discharging in a stratified stagnant or flowing ocean.Numerical Heat Transfer Part A,1989,16:371-387
    [60]陶文铨.数值传热学.西安:西安交通大学出版社,2001
    [6l]刘儒勋,王志峰.数值模拟方法和运动界面追踪.北京:中国科技大学出版社,2003
    [62]Ferzige,M.Perie.Computational Methods for Fluids Dynamics.Springer,1999
    [63]S.R.Mathur,Y.J.Murthy.A pressure based method for unstructured meshes. Numerical HeatTransfer,1997,31:195-216
    [64]D.Pan,J.Cheng.A second-order upwind finite-volume method for the Euler solution on unstructed triangle meshes.International Journal of Numerical Fluids,1993,16:1079-1098
    [65]于普兵.二维浅水水流数值模拟技术研究.南京:南京水利科学研究院,2006
    [66][美]S.V帕坦卡.传热与流体流动的数值计算.张政译.北京:科学出版社,1984
    [67]F.H.Harlow,J.E.Welch.Numerical calculation of time-dependent viscous incompressible flow of fluid with flee surface.The Physics of Fluids,1965,12(8):2182-2189
    [68]A.H.Nikseresht,M.M.Alishahi,H.Emdad.Complete flow field computation around an ACV(air-cushion vehicle) using 3D VOF with Lagrangian propagation in computational domain.Computers and Structures,2008,86(7):627-641
    [69]Wang Chengming.Computer modeling of fluid flow and heat transfer in metal castings.Ph.D thesis,University of Pittsburgh,1990
    [70]柳百成,荆涛等编著.铸造工程的模拟仿真与质量控制.北京:机械工业出版社,2001
    [71]H.J.Lin,W.S.Hwang.Three-dimensional mold filling simulation for casting and its experimental verification.AFS Transactions,1989,97:855-862
    [72]S.Sulaiman,A.M.S.Hamouda.Simulation of metal filling progress during the casting process.Journal of Material Processing Technology,2000,100:224-229
    [73]K.Yamarnoto,S.Satake,H.Yamashita.Lattice boltzmann simulation on porous structure and soot accumulation.Mathematics and Computers in Simulation,2006,72:257-263
    [74]Z.Y.Shi,Y.H.Yan.A lattice boltzmann method for simulation of a three-dimensional drop impact on a liquid film.Journal of Hydrodynamics,2008,20(3):267-272
    [75]D.Munger,A.Vincent.A Level Set approach to simulate magnetohydrodynamic instabilities in aluminum reduction cells.Journal of Computational Physics,2006,217(2):295-311
    [76]R.Caiden,R.P.Fedkiw.A numerical method for two-phase consisting of separate compressible and incompressible regions.Journal of Computational Physics,2001,166(1):1-27
    [77]R.R.Nourgaliev,S.Wiri.On improving mass conservation of Level Set by reducing spatial discretization errors.International Journal of Multiphase Flow,2005,31(12):1329-1336
    [78]李会雄,邓晟.Level Set输运方程的求解方法及其对气.液两相流运动界面数值模拟的影响.核动力工程,2005(26):242-248
    [79]C.Thomas,S.Osher.Simplex free adaptive tree fast sweeping and evolution methods for solving Level Set equations in arbitrary dimension.Journal of Computational Physics,2006,213(2):458-473
    [80]M.Sussman,E.G.Puckett.A coupled Level Set and Volume-of-Fluid method for computing 3D and axisymmetric incompressible two-phase flows.Journal of Computational Physics,2000,162(2):301-337
    [81]F.Losasso,P.R.Fedkiw.Spatially adaptive techniques for Level Set methods and incompressible flow.Computers and Fluids,2006,35(10):995-1010
    [82]J.Yu,S.Sakal.A coupled quadrilateral grid Level Set projection method applied to ink jet simulation.Journal of Computational Physics,2005,206(1):227-251
    [83]S.E.Hieber,P.Koumoutsakos.A lagrangian Particle Level Set method.Journal of Computational Physics,2005,210(1):342-367
    [84]M.P.Zhang,C.W.Shu.A weighted essentially non-oscillatory numerical scheme for a multi-class Lighthill-Whitham-Richards traffic flow model.Journal of Computational Physics,2003,191(2):639-659
    [85]RaimundB(u|¨)rger,Kozakevicius Alice.Adaptive multiresolution WENO schemes for multi-species kinematic flow models.Journal of Computational Physics,2007,224(2):1190-1222
    [86] R. R. Nourgaliev, T. G. Theofanous. High-fidelity interface tracking in compressible flows : Unlimited anchored adaptive level set. Journal of Computational Physics, 2007, 224(2): 836~866
    [87] N. Cavallini, V. Caleffi, V. Coscia. Finite volume and WENO scheme in one-dimensional vascular system modelling. Computers and Mathematics with Applications, 2008, 56(9): 2382-2397
    [88] W. S. Wang, S. F. Li. Conditional contractivity of Runge-Kutta methods for nonlinear differential equations with many variable delays. Communications in Nonlinear Science and Numerical Simulation, 2009, 14(2): 399~408
    [89] T. A. Gardiner, J. M. Stone. An unsplit Godunov method for ideal MHD via constrained transport in three dimensions. Journal of Computational Physics, 2008,8(1): 4123-4141
    [90] Jinhai Chen. The convergence analysis of inexact Gauss-Newton methods for nonlinear problems. Computational Optimization and Applications, 2008, 40(1):97-118
    [91] Y. H. Brandon, J. B. Yang. Using interactive multiobjective methods to solve DEA problems with value judgements. Computers and Operations Research, 2009,36(2): 623-636
    [92] R. P. Fedkiw, T. Aslam. A non-oscillatory Eulerian approach to interfaces in multimaterial flows (the Ghost Fluid Method). Journal of Computational Physics,1999, 152: 457~492
    [93] Yang Xiaofeng, James Ashley J. A Hybrid Particle Level Set method improved interface capturing. Journal of Computational Physics, 2001,183: 83~116.
    [94] S. P. Pij , A. Segal. Mass-conserving Level-Set method for modelling of multi-phase flows. International Journal for Numerical Methods in Fluids, 2005,47: 339~361
    [95] C. Farhat, A. Rallu, S. Shankaran. A higher-order generalized ghost fluid method for the poor for the three-dimensional two-phase flow computation of underwater implosions.Journal of Computational Physics,2008,227(16):7674-7700
    [96]陈立亮.低压铸造连续生产过程数值模拟及其质量控制.武汉:华中科技大学工学博士学位论文,1997
    [97]J.C.Martin,W.J.Moyce.An experimental study of the collapse of liquid columns on a rigid horizontal plane.Journal of Computational Physics,1992,244:312-324
    [98]E.J.Finnemore,J.B.Franzini.Fluid Mechanics with Engineering Appications.Springer,2003
    [99]王献孚,熊鳌魁.高等流体力学(研究生教学用书).武汉:华中科技大学出版社,2003
    [100]M.Schmid,F.Klein.Fluid flow in die cavities-experimental and numerical simulation.The 18th International Die Casting Congress and Exposition,Indianapolis,Indiana,1995,10:93-97
    [101]O.S.Asfour,M.B.Gadi.Using CFD to investigate ventilation characteristics of vaults as wind-inducing devices in buildings.Applied Energy,2008,85(12):1126-1140
    [102]K.P.Gertzos,P.G.Nikolakopoulos,C.A.Papadopoulos.CFD analysis of journal beating hydrodynamic lubrication by Bingham lubricant.Tribology International,2008,41(12):1190-1204
    [103]B.Sirrell,M.Holliday,J.Campbell.The benchmark test 1995.Proceedings of the 7th Conference on Modeling of Casting,Welding and Advanced Solidification Processes,London,UK,1995,9:915-933
    [104]M.Profit,S.Kulaesegaram.Modelling of high pressure die casting using the CSPH method.Proceedings of the 3rd International Conference on Engineering Computational Technology,2002,5:103-104

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700