传热与流体流动的数值计算(45章.ppt

5-6假扩散关于假扩散的一般观点中心差分格式具有二阶的精度，而上风方案只具有一阶的精度；上风方案引起严重的假扩散由于在对流-扩散问题中所产生的φ~x变化是指数的，除了极小的Δx值外，任何情况下泰勒级数不再是一个好的表达式。较大Δx值时，泰勒级数的分析给人误解；此时上风方案给出比中心差分格式更为合理的结果。把中心差分格式同上风方案的系数比较可以看到，上风方案似乎把真正的扩散系数扩大了一个虚拟的扩散系数ρuδx/2。但需要注意的是，争辩的麻烦建立在把中心差分格式假设为精确而又标准的基础上，并且用这个参考来观察上风方案。按这种做法，就会发现本身就是精确解的指数方案也有某种假扩散。实际上所谓的假扩散系数是在大的贝克列数条件下的一种理想补充。*有关假扩散的正确看法假扩散是一种多维的现象，在稳态的一维问题中绝对不会有相应形式的假扩散。为了具体化，讨论图中所示状态。两股速度相等而温度不同的平行流相遇。如果扩散系数不等于0，就形成一温度逐渐由高温到低温变化的混合层；如果扩散系数不为0，就不会形成混合层，并且在流线方向将维持温度得不连续性。如果对Γ=0情况的数值解产生一个逐渐变化的温度分布，我们可以得出结论：该数值解方案引起假扩散。*对两种不同网格方向用上风方案求解5.13b中的问题。在x方向的均匀流动。Φp=φW，结果，每一条水平线上上游的已知值将构成在该线上所有点上的值。没有假扩散现象发生。在与网格线成45。角的方向上的均匀流动。有设左边界温度为100，底边界温度为0来表示不连续。在内点上所得的解写在每个网格点的旁边。所得到的实际值却代表着一个逐渐变化的温度分布。*要点流场方向与网格点线成一倾斜角，并在与流动方向相垂直的方向上存在有非零的因变量梯度时，就会有假扩散出现。对二维情况的假扩散系数的一个近似表达式由公式给定：其中，Γ假是假扩散系数，U是合速度，θ是速度向量与x方向间夹角。可以看出：合流动方向与其中一组网格线相重时，不存在假扩散；流动方向与网格线间夹角45。时，假扩散最严重。减少Δx和Δy可以减少假扩散大小，只要可能，应当把网格线布置在接近于流动的方向。应使假扩散与实际值相比足够小。采用中心差分格式不是解决假扩散的灵丹妙药。产生假扩散的基本原因是由于把流过每一控制容积面的流动处理成局部一维流动这样的一种做法。那些可能给出较小假扩散的方案应当考虑流动的多维性质。*逐行法把TDMA和高斯-赛德尔法结合起来。选择一条网格行（设在y方向选取这样的网格行），假定沿相邻的行上的T值批最新值构成。用TDMA法求得所选行上的T值。将在同一方向的所有行进行这种计算。如果想做的话，再按相同的方法在其他方向重复上述程序。以二维为例，如图所示的情况需要注意：其它一些迭代方法ADI（方向交替的隐式）的逐行求解法；解多维离散化方程的强隐式法（SIP）*依前后二次迭代之间因变量的变化究竟是被加速还是被减慢的过程称为超松弛或欠松弛。超松弛常用于和高斯-赛德尔法相结合，叫做持续超松弛（SOR）；欠松弛在强烈非线性方程组的迭代求解中用来避免发散。取T*p作为前一次迭代所得Tp值。引进松弛因子，得到：可以根据经验以及对所给定的问题所作的试探性计算求得一个合适的?值。4-5超松弛和欠松弛(overrelaxationandunderrelaxation)*么么么么方面Sds绝对是假的通用惯量进行松弛。用下面公式代替离散化方程：式中i是所谓的惯量。对于正的i值，方程具有欠松弛作用；对于负的i则产生超松弛。*控制容积面的位置讨论控制容积面构成的两种不同的替代形式，并讨论它们各自有关的优点。为方便起见，描述针对二维问题。方法A：控制容积面放在两个网格之间的中点。4-6某些几何上的考虑结果是：一个典型的网格P并不落在包围该点得控制容积的几何中心上。*方法B：网格点放在控制容积的中心：克服了A的缺点。具有方便性。*我们所提出的这种方法不只限于直角坐标系，还可以用于任意一种正交坐标系。以二维极坐标问题为例，与方程对应的r?形式是：其中的网格与控制容积如图示：设控制容积在z方向厚度为1，方程两边乘以r，并在整个控制容积范围内对r和进行积分，得到下面的离散化方程：4-6其它坐标系*由一个新的坐标系引入的补充特征主要是几何上的特征。*5-1任务在通用微分方程中将对流项考虑进去，只要对流项的加入不改变离散化的形式，同样的处理方法仍然适用。本章任务是：在已知的流场（即速度分量和密度）的情