热流问题的数值计算.pptVIP

要计算Pi, Qi，需知道Pi-1, Qi-1, 最终需知道P1及Q1。P1及Q1可由左端点的离散方 程来确定： 其中 所以 当消元到最后一行时，有： ,而 所以 将TM1回代即可得出Ti。 该算法又称为三对角阵算法（Tridiagonal matrix method ，TDMA算法） 4.4.2 交替方向隐式算法 ADI 原理：在某一方向用隐式求解而使y，z方向成为显式，然后再类似地在y, z方向上各自用隐式求解依次。 1.Peaceman-Rachford ADI格式 在这一格式中，对于三维非稳态导热，在每一方向都以?t/3的时间间隔用隐式方法求解一次，三次连续求解就完成了一个时层方向上的推进。 若将二阶导数的中心差分表示为 ， 则这种ADI算法可以表示为： ? 其中 及 分别为 及 时刻上 的中间温度。 该格式稳定的条件为： 对于二维问题P-R ADI方法是绝对稳定的， 没有时间步长的限制。 2.? Brian ADI格式 在此格式中，所用时间间隔为 。 该格式可以表示为： 该格式无条件稳定。 特点：是利用相邻两层的平均值来计算空间导数，用?t计算时间导数。 该格式无条件稳定。 3. Douglas ADI格式 谢谢指导！ 汇报完毕 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）