热流问题数值计算.pptVIP

关于源项线性化的说明 1. 当源项为未知量的函数时，线性化处理比假定源项为常数更合理。因为如果S = f(T)，则把各控制容积内的S作为常数处理就是以上一次迭代计算所得的T*来计算S，这样源项对于T永远有一个滞后； 2. 线性化处理是建立线性代数方程所必须的。如果采用二阶或高阶的 多项式，则所形成的离散方程就不是代数方程； 3.为保证代数方程迭代求解的收敛，要求 。 因为离散方程都可表示为 线性代数方程迭代求解收敛的一个充分必 要条件是对角占优，即： 要求 4. 由代数方程迭代求解的公式 SP绝对值的大小影响到迭代过程中温度的变化速度，SP的绝对值越大，系统的惯性越大，相邻两次迭代之间TP的变化越小，收敛速度下降，但有利于克服迭代过程的发散；SP的绝对值小，可使变化率加快，但容易引起发散。 例1： ，取SC = 3，SP = -5。 例2： ，可取， 例3： 所以 当计算区域的边界为第二、三类边界条件，可采用补充节点法和附加源项法 1.补充节点法： 对于区域离散法A，以无限大平板的第二类边界条件为例，图4-9边界条件的表达式： 4.3.2 边界条件的处理 用差分表达式表示为： （2） 为得到二阶截差的公式，采用虚拟点法，在 右边设一个虚拟点M1+1，M1即可视为内节点，其一阶导数的中心差分为： （3） 为消去TM1+1，由一维、稳态、含内热源的控制 方程可得在M1点的离散形式为： (4) 由式（3）和（4）消去TM1+1得 (5) 其中 对于第三类边界条件，用 代入式 （2）、（5），得相应于一阶与二阶截差的节点离散方程为： （6） （7） 用控制容积平衡法，对边界节点M1的控制容积 作能量平衡，得 求解上式，得 （8） 控制容积平衡法在边界节点离散方程的建立中得到广泛的应用。 二阶精度 当采用区域离散法B时，边界节点可看作是第一种区域离散法中当边界节点所代表的控制容积趋近于零时的极限，对右端点，由式（5）和（7）得第二类边界条件：