第六章回流问题流动－传热耦合计算的数值方法
前章讲述了引入流动项后，相对纯粹的导热扩散方程，对流扩散方程离散求解的特点，
侧重讨论了对流项一阶导数的离散方法。但是并没有涉及至关重要的流动物理量——速度的
求解方法，也未论及流动和传热之间的耦合。本章将针对回流问题讨论此问题。不同于至少
有一个空间坐标是单通性质的边界层流动问题，回流问题的空间坐标均为双通性质，其控制
方程类型，对稳态问题为椭圆型，对非稳态问题则是以时间坐标作为行进坐标的抛物型。实
际工程中的流动和传热问题，回流问题更为普遍，其中不可压流体又居多数，因此，我们把
讨论范围限定于不可压流体的回流问题。
6.1不可压缩流体流－热耦合问题数值计算概述
流动—传热耦合问题，其控制方程是由连续、动量和能量方程在内的一组耦合的方程组
成，不能单独求解。即使物性都是常数，对流项的存在，动量方程是非线性的；若速度与温
度相关，如自然对流问题，则能量方程也是非线性的。控制方程的耦合性和非线性，决定了
流热耦合问题的求解必须使用迭代，而且包含多种性质的迭代，如各个关联求解量间的迭代，
解非线性问题的线化求解迭代，解线性代数方程组的迭代等。在使用迭代方法求解具体的流
热耦合问题时，要注意不同迭代的意义。
动量和能量方程中的一阶导数对流项的离散，会遇到前章讲到的诸多问题，必须予以特
殊的考虑，使其有必要的精度，又能保证解的稳定性，不致使解失去物理的真实性。
求解流热耦合问题的关键是流场的求解。求解流场，可以用原始方程中的速度、压力作
为基本变量，这称为原始变量法；也可用涡函数、流函数作为变量，这称为非原始变量的涡
流函数法。本章将分别介绍这两种方法，但重点放在原始变量法。
对离散方程的代数求解方法，可以分为耦合（或称联立）求解法和顺序（或称分离）求
解法两类。前者又包括所有变量的代数方程组全场联立求解、部分变量的代数方程组全场联
立求解以及局部区域所有变量联立求解等多种。但耦合求解法对计算机资源要求较高，编程
相对困难，发展缓慢。而顺序求解法则是先不考虑各变量的耦合，一个方程针对一个主要变
量，而把其它变量暂且作为已知，一个一个独立求解。完成一轮计算后，再按改进变量的计
算公式进行更新；返回至起始计算作迭代，再更新，直到迭代收敛，所有变量值都能满足全
部的方程。此法相对简单，可以统一编程，对计算机要求也较低，是当前工程计算的主要方
法。本章将针对以速度、压力为变数的原始变量法和以涡、流函数为变量的非原始变量法的
离散方程，分别介绍顺序求解法的压力修正方法和涡流函数法。
6.2原始变量法顺序求解流场所遇困难及解决途径
6.2.1简化条件下原始变量法求解流场的控制方程
为便于说明用原始变量法顺序求解流场所遇到的问题，先考察一个不计质量力影响的二
维不可压缩流动问题。此时，流动可以脱开能量方程先行解出，控制方程为
?u?v
+0（6-1a）
?x?y
120
?u?(uu)?(uv)??u??u1?p
++(ν)=+(ν)?(6-1b)
?t?x?y?x?x?y?yρ?x
?v?(uv)?(vv)??v??v1?p
++(ν)=+(ν)?(6-1c)
???????ρ?
txyxxyyy
该方程组有三个待求变量u,v,p，但只有关于速度u,v的两个输运方程，压力p则是以一
阶导数的形式作为源项出现在速度的方程之中，它没有自身独立的输运方程。控制方程中压
力的这些特征，使流场的求解面临一些需要解决的新问题。
6.2.2常规网格下离散压力导数可能导出不合理的解
所谓常规网格，是指控制方程中所有求解变量u,v,