求解流动与传热问题的一种高效稳定的分离式算法—IDEAL上.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083382 求解流动与传热问题的一种高效稳定的 分离式算法—IDEAL(上) 孙东亮 屈治国 何雅玲 陶文铨 （西安交通大学动力工程多相流国家重点实验室,陕西 西安 710049） （联系电话：02982669106 E-mail：wqtao@mail.xjtu.edu.cn） 摘要：本文提出了一种求解流动与传热问题的高效稳定的分离式算法—IDEAL(Inner Doubly-iterative Efficient Algorithm for Linked-equations) 。在IDEAL 算法中每个迭代层次上对压力方程进行两次内迭 代计算，第一次内迭代过程用于克服 SIMPLE 算法的第一个假设，第二次内迭代过程用于克服 SIMPLE 算法的第二个假设。这样在每个迭代层次上充分满足了速度和压力之间的耦合，从而大大 提高了计算的收敛速度和计算过程的稳定性。本文详细阐述了 IDEAL 算法在三维坐标下的计算步 骤，在本文下篇中通过 3 个三维算例对 IDEAL 算法和 SIMPLER、SIMPLEC、PISO 算法进行了比较 分析，证明了该算法的可行性。 关键词：分离式算法、IDEAL 、内迭代、收敛性、健壮性 0 前言 数值求解流动和传热问题的方法分为以压力为基本求解变量和以密度为基本求解 变量两大类。由于以压力为求解变量的方法可以同时适用于不可压缩及可压缩流动， 所以该方法的得到了普遍的应用。在以压力为求解变量的方法中，压力修正算法是最 普遍被使用的方法。压力修正算法于 1972 年由 Patankar 和 Spalding[1]提出，并被命名 为 SIMPLE 算法。SIMPLE 算法有两个基本假设，(1)迭代初始压力和速度场单独进行 设定，压力和速度的内在联系没有得到体现；(2)引入了压力修正，并在压力修正方程 的推导中，忽略了邻点速度修正的影响。以上两个假设不影响计算的最终结果，但会 影响计算过程的稳定性和收敛性，为了克服或者减轻这两点假设而引起的缺陷近三十 [2] [3] [4] 年来提出了多种改进方案，如 SIMPLER 算法 ，SIMPLEC 算法 ，PISO 算法 等。 以上算法均采用压力修正值p ′来改进压力和速度，但是一旦引入了压力修正值概念， 为了获得压力修正方程，在速度修正方程中不可避免的忽略邻点速度修正的影响，因 此为了克服 SIMPLE 算法中的第二个假设，我们需要直接改进压力，而不应引进压力 修正值，以上思想是CLEAR 算法[5] 的成功所在。但是直接求解压力方程会导致 CLEAR 算法稳定性较差，为了克服这一缺点本文在 CLEAR 算法的基础上提出了 IDEAL 算法 (Inner Doubly-Iterative Efficient Algorithm for Linked-Equations) 。 1 IDEAL算法的提出 SIMPLE 算法的两个基本假设是影响计算收敛速度和计算过程稳定性的根本原 因，所以一个理想的算法应该同时克服以上两个基本假设，也就是同时具有以下两个 基本特征：(1) 初始压力应尽可能接近当前迭代层次上的最终压力，该条件可以通过 资助项目: 国家自然科学基金重点项目和国家重点基础研究发展计划(2007CB206902) 克服 SIMPLE 算法的第一个假设获得。这样在该初始压力下计算动量方程所得的中间 速度基本同时满足了连续性方程，从而保证了算法的收敛性和健壮性；(2) 每个迭代 层次的最终速度应该尽可能同时满足动量和连续性方程，该条件可以通过克服 SIMPLE 算法的第二个假设获得。这是另外一个提高算法收敛性和健壮性的条件。通 过以上分析我们对压力方程进行两次内迭代计算，从而使