模拟任意曲线边界型腔充填过程的sola-mac 方法 - 东南大学学报.pdf

第28 卷第4 期 东 南 大 学 学 报 Vol28 No4 1998 年7 月 JOURNAL OF SOUTHEAST UNIVERSITY July 1998 SOLA-MAC 王云华 高志强 苏华钦 ( 东南大学机械工程系, 南京21009 ) 针对压铸、挤压等铸型型 具有曲线、曲面这一情况, 对具有任意曲线边界的 型 充填过程进行了研究. 对任意曲线边界进行了精确处理, 并通过判断自由表面位 置, 较好地满足了自由表面边界条件, 从而提高了数值模拟的精度. 采用了易于实现的 SOLA 方法进行压力-速度场迭代, 采用改进的MAC 方法跟踪了自由表面, 并给出了一 个算例. 自由表面流动; 充填过程; 铸造CAE T G214 采用数值方法模拟型腔的 填过程, 可以预测气体的卷入、夹渣、冷隔、浇不足等缺陷以及 熔接线位置, 还可得到初始温度场以提高缩孔、缩松的预测精度, 因此近年来受到了国内外学 者的广泛重视[ 1~ 5] . 型腔 填属于具有自由表面的粘性不可压缩流体的流动问题, 目前主要 有拉格朗日和欧拉两种解法. 前者采用移动单元, 可以精确地满足自由表面边界条件, 但难以 [ ~ 9] 解决单元扭曲、相交等问题 . 因此这种方法在商业软件中很少使用. 人们用得最多的是欧 拉方法或欧拉与拉格朗日相结合的方法, 如MAC、SOLA-VOF、SMAC、ABMAC、PIC、VOF 等. 其中得到广泛应用的是MAC 和SOLA-VOF 两种方法. MAC ( M arker and Cell) 法通过在单元中布置大量粒子( Marker) 来跟踪自由表面, 假定粒 子没有质量且随流体的运动而运动. 而为了标识这些粒子需要额外占用大量计算机内存与计 [ 10] 算时间. 为了解决这一问题, Nichols 等人 又提出了SOLA-VOF ( SOLution Algorithm-Vol- ume of Fluid) 算法. 它只用一个称为体积函数( Volume of Fluid Function) 的标量来标识自由表 面, 所以提高了计算效率. 但对自由表面的近似程度有所下降. 为了解决计算效率与计算精度 之间的矛盾, 作者首先对MAC 算法进行了改进, 只在自由表面单元布置粒子, 从而大大减少 了粒子数目. 其次将MAC 方法与SOLA-VOF 方法相结合, 采用SOLA 方法求解压力与速度 场, 采用改进的MAC 方法跟踪自由表面. 在实际生产中, 模具型腔往往包含很多曲线、曲而. 传统的控制体积法只是简单地以规则 网格来逼近型腔形状, 曲线、曲面处的边界条件处理得较为粗糙, 降低了分析精度. 针对这一情 况作者重点对任意曲线边界的处理进行了研究. 国家自然科学基金资助项目. 收稿日期: 1997- 10- 27, 修改稿收到日期: 1998- 03- 25. 114 东 南 大 学 学 报 第28 卷 1 1. 1 型腔中往往具有液体和气体两相. 当同时考虑液相和气相时, 控制方程为: 1) 不可压缩液体 v = 0 ( 1) v 2 = - ( v ) v - p + b + v ( 2) l t l l 2) 气体 air t +