针对大平板铸件的反重力高速充填.DOC

针对大平板铸件的反重力高速充填 流场数值模拟的对称性研究 申发田，李日，李菲，张建强 （河北工业大学材料科学与工程学院，天津 300130） 摘要:采用Simple结合守恒标量法进行大平板反重力充填的数值模拟计算，已经能够精确地模拟充型速度小于2m/s情况下的流动过程，但当速度大于该速度时流场就会发生偏转，产生流场不对称现象。本文对Simple算法的流量计算方法用调和平均值法进行修正，使流场不对称现象得到了极大的改善，充型速度提高到40m/s。 关键词：数值模拟；高速流场对称性；调和平均值；薄壁铸件 The Study on Numerical Simulation of Anti-Gravity Mould Filling in High-speed of Thin Wall Plate with Large area SHEN Fa-tian ,LI Ri ,LI Fei ,ZHANG Jian-qiang (Materials Science and Engineering School of Hebei University of Technology ,Tianjin 300130) ABSTRACT：The antigravity mould filling of thin wall plate with large area was simulated by SIMPLE method, but the maximum filling rate for simulation was only 2m/s, and if the speed was greater than 2m/s, the flow field would be asymmetric. In this paper，the algorithm for flow was amended with harmonic mean method, the simulation results showed that the symmetry of the flow field in high speed could greatly be improved to 40m/s. Key Words: Numerical Simulation，the symmetry of flow field in high-speed，harmonic mean method， thin wall plate 薄壁铸件在充型过程中易产生浇不足、冷隔等缺陷，模拟其流动充型过程对提高铸件质量有重要意义。 论文[1、2]采用Simple法结合守恒标量法的Van-leer格式对大平板薄壁铸件进行了充型过程数值模拟，当速度小于2m/s时，可以获得与实验吻合较好的计算结果，但当速度大于2m/s时，流场对称性会出现较大偏差，文献[3]采用各种方法如Simpler算法(Simple的修订算法)、Simple的DATE改进方案和多时间步长技术，对解决该问题作出了努力，但效果不理想。本文采用流量调和平均值法显著地提高了流场的对称性。 1 计算原理 1.1流量计算控制方程 流场计算控制方程的张量形式为： （1）N-S方程 　　　　 （），为xi方向的速度分量和重力加速度分量；，时间；，压力；，流体密度；，运动粘性系数。 （2）连续性方程 守恒标量法是用单相的控制方程同时求解计算域中的气液两相流动，对自由表面格子而言，流入流体 密度与流出流体密度差别很大，因此，相应的连续性方程采用如式（2）所示形式。 作者简介：申发田（1986-），男，河北邯郸人，硕士，主要研究方向：材料加工CAD/CAE；电话: E-mail: shenfatian@126.com； 通讯作者：李日，E-mail:sdzllr@163.com （2） 在本文的研究中，因排气通畅，空气可视为不可压缩流体，所以连续性方程为： （3）　　　　　　 　　　　　　　　　　　　　　　　　　　 式（2）、（3）两种形式的本质是体积守恒。 （3）自由表面控制方程： （） 式中为满格子；为自由表面格子；为空气格子。上述各控制方程均用Simple法