传热与流体流动的数值计算(4-5章).ppt

逐行法 把TDMA和高斯-赛德尔法结合起来。选择一条网格行（设在y方向选取这样的网格行），假定沿相邻的行上的T值批最新值构成。用TDMA法求得所选行上的T值。将在同一方向的所有行进行这种计算。如果想做的话，再按相同的方法在其他方向重复上述程序。 以二维为例，如图所示的情况需要注意： 其它一些迭代方法 ADI（方向交替的隐式）的逐行求解法；解多维离散化方程的强隐式法（SIP） 依前后二次迭代之间因变量的变化究竟是被加速还是被减慢的过程称为超松弛或欠松弛。超松弛常用于和高斯-赛德尔法相结合，叫做持续超松弛（SOR）；欠松弛在强烈非线性方程组的迭代求解中用来避免发散。 取T*p作为前一次迭代所得Tp值。 引进松弛因子，得到： 可以根据经验以及对所给定的问题所作的试探性计算求得一个合适的?值。 4-5 超松弛和欠松弛(overrelaxation and underrelaxation) 通用惯量进行松弛。用下面公式代替离散化方程： 式中i是所谓的惯量。对于正的i值，方程具有欠松弛作用；对于负的i则产生超松弛。 控制容积面的位置 讨论控制容积面构成的两种不同的替代形式，并讨论它们各自有关的优点。为方便起见，描述针对二维问题。 方法A：控制容积面放在两个网格之间的中点。 4-6 某些几何上的考虑 结果是：一个典型的网格P并不落在包围该点得控制容积的几何中心上。 方法B： 网格点放在控制容积的中心： 克服了A的缺点。 具有方便性。 我们所提出的这种方法不只限于直角坐标系，还可以用于任意一种正交坐标系。以二维极坐标问题为例，与方程 对应的r?形式是： 其中的网格与控制容积如图示： 设控制容积在z方向厚度为1，方程 两边乘以r，并在整个控制容积范围 内对r和进行积分，得到下面的离散 化方程： 4-6 其它坐标系 由一个新的坐标系引入的补充特征主要是几何上的特征。 5-1 任务 在通用微分方程中将对流项考虑进去，只要对流项的加入不改变离散化的形式，同样的处理方法仍然适用。 本章任务是：在已知的流场（即速度分量和密度）的情况下，求得对φ的解。 已知流场必须满足连续性方程： 第五章 对流与扩散 通用微分方程 也可以改写为： 对于已知的?、uj、?以及S的分布，任何解φ及其变体（φ 加一常量）将同时满足方程，关于系数和的基本原则仍然适用。 讨论只有对流与扩散这两项存在的情况下的一维稳态问题。 控制微分方程： 应用图示三网点群： 5-2 一维稳态对流与扩散 P W E x w e Dx (dx)w (dx)e - 预备性的推导 在整个控制容积内对方程（5.4）积分： 由对φ的一个分段线性分布表示项Γdφ/dx。结果是： 因子1/2出自界面位于中点的假设；对不同的界面位置要采用其它内插因子。则方程（5.6）写成： F x P W E w e Fp Fw FE 定义两个新的符号： 两者具有相同因次， F表示对流或流动的强度；D是扩散传导性。（注意，D永远为正，而F不同） 离散化方程变为: 讨 论 由于连续性Fe=Fw，得到ap=aE+aW 上述离散化方程隐含着φ分段线性分布——中心差分 假定De=Dw=1及Fe=Fw＝4，则 若φE＝200及φW＝100，φP＝50 若φE＝100及φW＝200，φP＝250 方程（5.11)表明系数可能出现负值。 当|F|小于2D时，系数才可能始终为正，即中心差分只能限于低Reynold数流动。 若扩散项为零，则中心差分格式导致ap=0，无法使用逐点法求解 - 上风方案 亦称为上风差分格式、迎风格式、上游差分格式以及供体（施主）室法等。 方案认为预备性公式的弱点在于假设：界面上的对流性质φe 是φE和φP的平均值。提出：保留扩散项的公式不变，而对流项则按下列假设计算： 界面上φ的值等于界面上风侧网格点上的φ值。于是 类似方法可确定φw值。 定义 代表A，B中大者。则上风方案意味： 离散化方程可写为： 不会产生负的系数； 可以把这个方案说成是建立在“槽与管”的模型基础上，管内的流体不会“知道”将要流入那个槽内的任何情况，但它却携带了它所来自那个槽内的全部信息。这就是上风方案的本质。 精确解(Exact solution) 如果Γ取作常数，且边界条件为： 则 其中 为贝克列数（ Peclet Number）， 是对流与扩散强度之比。 不同的贝克列数时的φ~x变化如图 P为0的极限条件下，问题成了纯扩散（或热传导）问题。 除非lP