基本方程组的数值求解3.pptVIP

4）出口截面上的法向速度与切向速度都按局部质量守恒确定 在这一方法中，对图12所示情形， 计算式同式（85a），但ui,M的更新则通过对图13所示半个控制容积的质量守恒来获得。以0.5(vi,M+vi,M-1)作为该半个控制容积南侧流速，则可以写出 （86a） 将边界上的切向流速取本迭代层次之值而其余速度取上一迭代层次的值，得 （86b） 由图11（b）可见，在计算区域右边界， 之值为零，由此可推得 ，这样依次计算可以得出出口边界切向速度的更新值。 * 5）对流边界条件 用以下公式来更新出口截面法向速度及其它标量： （87a） 其中，c为对流速度或相速度，一般取为通道平均流速。上式离散后，可得 （87b） 这里，vm为出口截面平均法向流速。相当于假定在出口边界上无扩散作用, 为一个纯对流问题。 * 6）内点插值方法 在文献[36]中，提出采用由内点的速度来决定出口流速的方法： （88a） 其中，为与出口边界相邻的控制容积的宽度，x为出口边界上游第四个截面位置起算的坐标，当网格均分时（如图14所示），有 （88b） * 图14 插值公式（88a）的图示 * 7）各种高阶导数为零的方法 （89） 其中，x为出口方向，上角标j及q表示求导的阶数，要求 ，并建议取j=3，q=2。这一类方法实际上相当于提出了由内点来插值的一个方式。导数的阶数越高，涉及的内点数越多，该阶导数取零对解的影响也越微弱。 * 出口边界上边界条件的设定有多种看来相差很大的方法。在从事工程流动与传热问题的数值计算中实用的处理出口边界条件的方法： 出口边界上无回流时，可采用局部单向化或法向一阶导数为零的方法； 当出口边界有回流时，可以采用法向流速按质量守恒、切向流速按齐次Neumann条件的方法； 无论采用哪一种方法，出口截面上的法向速度分布必须满足总体质量守恒的条件。 * 4．固体边界 1）对于温度一类的变量 用附加源项法来处理 2）对于速度 固体壁面是受限流动中最常见的边界，原则上应采用无滑移条件，即流体的速度等于壁面的速度（第一类边界条件），层流时无论是对于动量方程中的流速还是连续性方程中的流速都是这样（后者表现在SIMPLE算法压力修正方程中的剩余质量流量的计算上） 当采用标准 模型时，平行于壁面的流速及垂直于壁面的流速还存在不同处理方法。对于平行于壁面的流速（例如u），采用无滑移条件，同时壁面上的动力粘性系数由壁面函数法的公式来计算；但对垂直于该壁面的流速（v）则应区分动量方程还是连续方程。在连续性方程中仍然采用无滑移条件（静止壁面上取v=0），而对v动量方程则采用法向一阶导数为零的条件： * 3）对于k 采用低Reynolds数 模型时，可直接令壁面上k=0 采用标准 模型时一般取第一个内节点处 ，可通过令该点与壁面之间的k方程的扩散系数为零来实现。 4）对于 O在低Reynolds数 模型中，可以取壁面上为零，也有取法向一阶导数为零的 在采用标准模型及壁面函数法时，一般对第一个内节点赋予规定值，同时令第一个内节点与边界之间的扩散系数为零 5）组分的浓度 如果固体壁面对化学反应无催化作用，固体壁面上各组分的浓度决定于通过壁面的质量渗透率，以及渗透介质的化学成分。如壁面对组分是不可渗透的，则该组分在固体壁面法线方向的导数等于零。 * 5．对称边界 对于一切标量及与对称线平行的速度分量，对称边界上采取法向一阶导数为零的条件 对于垂直于对称线的速度分量，则取对称线上该分量为零。 * 6．壁面函数 用上述给出的边界条件，在用数值求解燃烧问题时，对壁面附近的流场要特别小心，才能保证计算有足够的精度。这是因为在壁面附近，输运系数变化很陡 在层流流动中，是由于通过壁面经常有很强的热通量，因此温度变化很陡造成的；在湍流流动中，是由于湍流输运系数在接近壁面时，迅速趋于零所致 如果象建立内部结点的差分方程那样，用线性插值办法计算非结点位置的输运系数，在这种变化很陡的区域，则需相当密的网格，并仔细安排，才能合理计算输运系数值，但所需的计算机容量和计算时间都大大增加 在湍流流动中，由于壁面附近湍流雷诺数比较小，在湍流模型中需考虑湍流雷诺数的影响，求解也比流场内部复杂 若我们主要关心的是流场内部的情况，对壁面