高分辨率有界的有限差分方法在铸件充型过程计算中的应用研究.pdfVIP

2005中国铸造活动周论文集 高分辨率有界的有限差分方法 在铸件充型过程计算中的应用 孙逊，关洋，李海兰，张士彦，苏仕方，王君卿 (沈阳铸造研究所，辽宁沈阳110022) 【摘要I引入规正变量定义，在非均匀网格内采用五点格式的高分辨率有界的MINMOD格式离散动 量守恒方程和能量守恒方程，用改进的VOF法处理流体自由表面变化，采用SOLA技术计算流体流动过 程。数值实验和水力模拟实验结果表明，上述处理不仅保证了计算速度，而且数值计算过程稳定，数值模 拟结果与水力模拟实验结果近似。 I关键词J充型过程；数值模拟；高分辨率；有界性；有限差分格式 ^_●■●_-o●一 O 刖嚣 数值模拟精度和计算效率是互相矛盾的。铸件充型过程本身是一个复杂的带有自由表面的流体湍流 流动和传热过程，但如果将湍流流动和传热过程数学模型用于铸件充型过程计算，虽然精度上可能更接 近于实际过程，但与层流流动和传热过程数值计算相比，要花费成倍的计算时间【l】。目前用于铸造充型 过程数值模拟的软件中大部分采用的是层流流动和传热模型。即便是采用层流模型，对于不同的算法，计 算精度和计算速度也有很大的差别，这一点已经被有关的研究成果证实【2】【3】。SOLA—VOF是一种计算 效率很高的算法，但其存在速度和压力迭代计算不容易收敛的缺陷。不仅如此，近十年的研究成果表明， 以前用于流体流动计算的有限差分格式如上风格式、混合格式以及幂函数(乘方)格式等一阶格式存在假 扩散、物理上不真实(有界性问题)和计算稳定性问题。为了克服这些同题，出现了多种高分辨率有界的有 限差分格式。本研究的目的，就是在保持充型过程计算效率的前提下，从提高动量守恒方程和能量守恒方 程的离散格式精度人手，研究高分辨率有界的有限差分格式对铸件充型过程模拟计算精度的影响，以期 达到在不过分损失计算效率的前提下改进计算精度和计算稳定性的目的。 1 数学模型 假设液态金属为常密度不可压缩的粘性流体，采用层流流动和传热过程数学模型描述铸件充型过 程，流体自由表面变化采用VOF方法处理。控制方程以张量形式描述如下。 议胂)。秽(肛J卿 a伊鼢、。a … —■o一十——；一=iU毒-J十a口 【1) 0l oX{ Ox{ 锨i 上式中，垂一变量。当①为l，Ui)F，T时，(1)式分别代表质量守恒、动量守恒、体积函数和能量守恒方 程，s由为源项。旷液态金属密度，g／em3。R一广义扩散系数，对于动量守恒方程，其量纲为Pa·s；对于 函数，0≤F≤l。T一温度，℃。 2 计算方法 方程(1)为对流扩散方程。对于对流扩散方程中的扩散项，一般均采用具有二阶截差的离散格式(中 心差分)离散，对于大多数实际问题，这种离散格式已完全能满足需要。离散格式的准确性主要受到对流 项的影响。数值扩散(假扩散false 2005中国铸造活动周论文集 (overshoot／undershoot)等使数值模拟的精度下降。 所谓对流项离散格式的稳定性，指的是数值计算过程中不会出现解的振荡现象。铸件充型过程数值 模拟中经常采用的迎风格式是绝对稳定的。 对于某些类型问题，即当物理量的在计算域中的某一位置或某一微小区域发生急剧变化时，对流项 的离散格式就会显示出另外一种重要特性，称为有界性(bound—ness)。 铸件形状千差万别，液态金属充型过程中经常遇到速度和压力梯度突变的情况，因此存在有界性问 题。 有界性问题是由对流项引起，与扩散项无关。所谓对流项离散格式的有界性，是指当采用某种格式对 物理量的场进行纯对流数值模拟时，数值计算的结果不会超出物理问题本身所规定的物理量的上下限， 即不会“越界”。越界现象不仅使物理量不真实，而且也有可能使数值求解过程发散。 凡是具备有界性的格式一定是绝对稳定的。格式具有有界性的条件是要求插值函数有界，即由插值 函数所规定的界面函数值永远不会超过获得这一插值的有关节点的最大值或最小值。 基于上述观点，本研究选择了计算效率较高的并适用于非均分网格的高分辨率有界的有限差分格式 MINMODt．j151离散动量守恒方程和能量守恒方程，为简化差分格式