高等流体力学第六章..ppt

6.7 线性代数方程的解 改写成： 与式（2）相比，有 6.7 线性代数方程的解 可以由左端点的离散方程来确定： 其中 ， 所以， 6.7 线性代数方程的解 当消元到最后一行时，有： 而 ， 所以， 从上式出发，逐个回代，得出Ti（N-1，…1）。 例题 k＝1；源项S＝－T，SC＝0，Sp＝－1，T1＝0 设有一导热型方程 ，边界条件为x＝0， T＝0；x＝1， 。将该区域三等分，求该问题的解。 解：由导热方程知 ， 6.7 线性代数方程的解 6.7 线性代数方程的解 a1=a3=3,a2=19/3,b2=0 a2=a4=3,a3=19/3,b3=0 a3=3,a4=19/6,b4=-1 ， 6.7 线性代数方程的解 即 19/3T2=3T3+3T1 19/3T3=3T4+3T2 19/6T4=3T3-1 T1=0 T2=-243/1121 T3=-27/59 T4=-840/1121 6.7 线性代数方程的解 高斯—塞始尔(Gauss—5eldel)逐点计算法 其中Tnb*代表在计算机存贮器中所存在的相邻点的温度值．对于那些在本次迭代过程中已经被访问过的相邻点是 新鲜的计值，而对于那些尚待访问的相邻点是由前一次达代所得到的值． 6.7 线性代数方程的解 高斯—塞始尔(Gauss—5eldel)逐点计算法 迭代序列号 0 1 2 3 4 T1 0 -1 -4 -11.5 -30.25 T2 0 -3 -10.5 -29.5 -76.13 6.7 线性代数方程的解 高斯—塞始尔(Gauss—5eldel)逐点计算法 迭代的序号 0 1 2 3 4 5 ··· T1 0 0．2 0．68 0．872 0．949 0．920 ··· 1．0 T2 0 1．2 1．68 1．872 1．949 1．980 ··· 2．0 6.7 线性代数方程的解 逐行计算法 6.7 线性代数方程的解 松弛 我们在方程的右侧加上 ，再减它，我们就有： 括号内的部分代表由本次迭代所产生的 的变化。这一变化可以通过引进一个 6.7 线性代数方程的解 松弛 括号内的部分代表由本次迭代所产生的 的变化。这一变化可以通过引进一个 松弛因子加以修改，所以 6.8 超松弛与欠松弛 在代数方程的迭代求解过程中，或是用于处理非线性问题的整体迭代模式中，往往希望加快或是减慢前后二次迭代之间因变量值的变化。 因变量的变化被加速——超松弛；因变量的变化被减慢——欠松弛。 对于离散方程 取 作为前一次迭代所得TP值。 采用一个松弛因子 离散方程可以 写成 如果我们在方程的右侧加上 变换 可以得到 括号内部分代表本次迭代所产生的 TP的变化。可以通过引进一个松弛因子加以修改 , 为松弛因子 当迭代收敛时， TP和 相等 6.8 超松弛与欠松弛 6.8 超松弛与欠松弛 松弛因子 在0与1之间时，他的作用是欠松弛的；对一很小的 值，TP的变化很慢。当 大于1时，就产生超松弛。 的最佳值与许多因素有关。如研究问题本身的特性，离散网格点的数目，网格点之间的间距以及采用的迭代方法。 在不同的迭代次数之间可以改变 的，对每一个网格点也可选用各自不同的 值。 6.8 超松弛与欠松弛 通过惯量进行松弛 用下面的公式来代替上述离散方程 式中i时惯量，i为正值时，上述方程具有欠松弛的作用，i为负值时具有超松弛的作用。 i应当于ap不相上下；同时i的数值愈大，松弛的作用也愈烈。 第六章 导热问题的数值解 6.1 一维稳态导热 6.2 边界条件与源项的处理 一维稳态导热 源项线性化 6.2 边界条件与源项的处理 一维稳态导热 其中 6.2 边界条件与源项的处理 一维稳态导热 式中下标nb表示一个相邻结点， 表示对所有的相邻结点求和 6.3 边界条件 在热传导问题中有三类典型的边界条化： 1．已知边界温度 2．已知边界热流密度 3．通过放热系数和周围流体的 温度来规定边界的热流密度． 6.3 边界条件 6.3 边界条件 对控制容积积分，考虑热流与温度关系 6.3 边界条件 6.3 边界条件（第三类） 如果热流密度 系由放热系数h以及环境流体温度T 那么 于是对于 的方程变为： 式中 6.4