计算流动与传热格式稳定性分析.docVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：123363 计算流动与传热格式稳定性分析? ?基金项目：国家自然科学基金（NO 张晓晖 仲启刚 张明杰 方攀 ( 苏州大学物理科学与技术学院，能源学院，苏州 215006) (通信作者：电话：0512 E-mail:zhangxhem@126.com) 摘要:对SGSD格式在一维、线性、无源、非稳态对流扩散方程离散时，应用von Neumann分析方法给出显示格式和隐式格式的时间（类时间）与空间步长匹配结果，以指导实践应用。在此基础上应用SGSD格式对具体数值问题进行模拟试验，即使在隐式格式的情况下，对于固定网格系统，数值解会随着时间步长的不同而产生不同结果，研究认为，较强列的非线性问题产生的数值解不收敛情况，是研究问题物理属性所决定的，产生不同时间步长导致不同结果的原因在于离散方程求解过程中不同时层系数变化较大导致不同结果。 关键词: 离散格式；稳定性；非线性；收敛 0 前言 热物理数值计算中，会遇到三种不同的稳定性问题[1]：一是代数方程迭代求解过程的不稳定性。它是指迭代方法选择不当，使迭代收敛条件不满足，而导致迭代发散的情形。二是初值问题离散格式的不稳定性。它是指采用离散格式计算时，由于时间（类时间）步长取得过大或者时间（类时间）与空间步长匹配不当，导致解的振荡发散现象。第三则是对流格式的不稳定性。它是指在采用某些离散格式求解对流扩散方程时，即使稳态情况，由于空间步长过大或者流速过高而形成网格peclet数过大，也会使解振荡发散。 所谓初值问题离散格式的稳定性，是指该格式在时间或类时间的推进计算过程中，对在某个时层由于某种原因引入的误差，都有抑制能力，即它能确保引入的各种误差不产生实质性增长，以致变得无界。 为了考察计算中误差的演化特性，OBrien[2]等人提出了这样的稳定性定义：如果计算中包含了舍入误差在内的全部误差是增长的，则格式是强不稳定的，反之是强稳定的；若只是舍入误差是增长的，则格式是弱不稳定的，反之是弱稳定的。若稳定仅用舍入误差的增长和衰减来定义，这就为实际进行稳定性分析提供了方便，可以采用Fourier分析方法来解决，又称之为von Neumann分析法。同时还假定，应用von Neumann分析方法证明是弱稳定的就暗示强稳定。 von Neumann稳定性分析方法是将一个时层上在网格点的误差函数分布展成为由不同谐波组成的有限Fourier级数，再考察其在时间演化过程中这些谐波分量的变化。如果各谐波分量的振幅都衰减或保持有界，则误差也衰减或保持有界，求解函数不致因误差引入无限制增长，离散格式稳定；反之不稳定。 对流扩散方程中对流项的离散格式一直是计算传热学的一个基础研究课题。国内研究者已提出一些具有一定通用性的格式，其中一种是稳定性可控的二阶差分(SCSD)[3]格式．该格式以CD和SUD格式为基础，在SCSD格式的基础上，国内研究者提出了一种保证稳定性的二阶差分(SGSD)格式[4]。 为了克服SIMPLE算法的假设，国内学者提出了IDEAL算法[5]，在 IDEAL算法中每个迭代层次上对压力方程进行两次内迭代计算，第一次内迭代过程用于克服SIMPLE算法的第一个假设，第二次内迭代过程用于克服 SIMPLE算法的第二个假设．这样在每个迭代层次上充 分满足了速度和压力之间的耦合，从而大大提高了计算的收敛速度和计算过程的稳定性。 本文首先本文对SGSD格式在一维、线性、无源、非稳态对流扩散方程离散时，应用von Neumann分析方法给出显示格式和隐式格式的时间（类时间）与空间步长匹配结果，应用SGSD格式对具体数值问题进行模拟试验，发现，即使在隐式格式的情况下，对于固定网格系统，数值解会随着时间步长的不同而产生不同结果，初步研究认为，除了常见的三种不同的稳定性问题外，还有一类稳定性问题即较强列的非线性问题产生的数值解不收敛情况，是研究问题物理属性所决定的。 1 一维对流扩散方程数值稳定性分析 1.1物理模型 无源项，速度u、密度、扩散系数均为常数的一维非稳态对流换热通用控制方程为的 其中为常数。 界面插值格式图案见图1所示。 图1 界面插值格式图案 常用格式界面插值的定义见表1.表1中为之间的一个参数，表示稳定性可控的（SCSD）格式界面插值中中心差分（CD)与二阶迎风（SUD）格式的相对分量。稳定性（SGSD