二维铸造充型过程数值模拟特征分数步长法.pdfVIP

二维铸造充型过程数值模拟的特征分数步长法∗ 1 2 鲁统超 ，葛亮 1 山东大学数学与系统科学学院, （250100 ） 2 山东大学数学与系统科学学院, （250100 ） E-mail ：lutc@ 摘 要：铸造充型过程的数学模型是包括连续性方程和动量方程的偏微分方程组。本文利用 分数步长法将动量方程分裂成两部分，对第一个方程采用特征差分法进行处理，对第二个方 程结合连续性方程进行处理后得到压力的 泊松方程，用迭代法进行求解，给出了收敛性分 析和稳定性条件。 关键词：分数步长；特征差分；收敛性；迭代。 1. 引 言 铸造生产的实质就是直接将液态金属浇入铸型并在铸型中凝固和冷却，进而得到铸件。 液态金属的充型过程是铸件形成的第一个阶段。许多铸造缺陷(如卷气、夹渣、浇不足、冷 隔及砂眼等)都是在充型不利的情况下产生的。因此，了解并控制充型过程是获得优质铸件 的重要条件。但是，由于充型过程非常复杂，长期以来人们对充型过程的把握和控制主要是 建立在大量实验基础上的经验准则。随着计算机的发展，铸件充型过程数值模拟才得到广泛 应用。 充型过程流场数值模拟的主控方程均为非线性方程。其计算使用有限差分或有限元等数 值方法求解质量守恒方程(连续性方程)和动量守恒 方程即Navier-Stokes 方程，以得出流体 运动规律。在以前的研究中，Chorin(1968)和Temam(1969)分别独立的提出投影法。1972 年 由Minnesota 大学的Patankar 与Spalding 提出了simple 算法，这是一个压力修正算法，在 以后的研究中又有simplec 方法,Raithby 提出的simplex 方法, Sheng 等提出的simplet 算法。 本文中利用分数步长法的思想将动量方程分裂成两部分，对第一个方程采用特征差分法 求解，对第二个方程结合连续性方程进行处理后得到压力的 泊松方程，我们用迭代法进行 求解，给出了收敛性分析和稳定性条件。 2. 问题的数学模型 铸造充型过程的模型主要由连续性方程和动量方程组成。 (a) 流体的动量方程 ∂u r 1 ∂p µ − ⋅ − + + ∆ (2.1) V u gx u ∂t ρ ∂x ρ ∗ 本课题得到教育部高等学校博士点基金资助，编号：20030422049 - 1 - ∂v r 1 ∂p µ − ⋅ − + + ∆ (2.2) V v gy v ∂t ρ ∂y ρ r r r 其中(x ,y ) ∈Ω, t ∈J (0,T ] ，V ui +vj 是速度向量，p 为压力，gx