-边界层方程【DOC精选】.docVIP

§9. 边界层方程 考虑不可压粘性流体的平板边界层流动。在边界层流动中， 的作用类似于非定常流动中的时间坐标 ，因此可以沿着 方向推进。这种方法称为空间推进方法。空间推进方法将从 线上的来流开始，沿 方向，向下游推进，逐步计算出 线、 线、… 等各条线上的流动参数。 平板边界层流动的控制方程 这里 是边界层外缘处的速度分布。通过求解边界层外的无粘流动，可事先求得 。因此在边界层流动的求解中， 这一项是已知的。 边界层流动 线），给定来流速度 ， 而边界条件是 沿壁面（ 线），速度 ， 取特征速度 和特征长度 ，将方程无量纲化，平板边界层流动的控制方程成为 方程中 、 、 ， 、 而 是（来流）雷诺数。 下面为了简洁起见，省略所有无量纲量上方的横杠。 按照空间推进的观点，流向坐标 是类时间坐标。因此边界层方程中的第一个方程，就是 的对流方程，可用迎风格式 离散。这里， 方向网格的指标 写成了上标。 这一格式的精度较低， 方向和 方向都是一阶的。稳定性条件是 。作为改进方案，可考虑类似于交错网格的处理，在网格中点 处列出差分格式，并改用隐式格式，即 这样可将 方向的精度提高到二阶，并改善格式的稳定性。同时还可以抑制“棋盘”模式的振荡。 类似地，控制方程中的第二个方程，就是 、 的对流扩散方程。可采用的差分格式有 显式中心差分格式 这一格式的精度较高， 方向一阶， 方向可达二阶。但稳定性条件较为苛刻，推进步长的限制为 实际计算中，考虑到速度剖面的单调性， 可取其最小值 。再取 ，则有 显式迎风格式 这一格式的精度较低， 方向和 方向都是一阶的。但是稳定性条件放宽了，推进步长的限制是 即 平均隐式格式 这一格式的精度最好， 方向和 方向都是二阶的，而且是无条件（线性）稳定的。但是需要求解方程组（如果不用冻结系数法，方程组还是非线性的）。 完全隐式格式 这一格式的精度为 方向一阶， 方向二阶，无条件（线性）稳定的。也需要求解方程组。 （注）上面谈到的稳定性，是根据冻结系数法所做的线性化稳定性分析。网格比 ， 这里不考虑边界层分离等复杂的流动现象，因此可以认为速度分量 。如果出现分离，整个空间推进的做法都不成立了，因为类时间坐标是不可逆的。 显式格式的优点在于计算简单，计算步骤是一个递推过程。以显式中心差分格式 为例，具体的计算步骤如下。 在 线给定初始 ， ，（） 但是在前缘（，）处，需按边界条件给定 ， 给定推进计算的最大步数 。对 ： 先给定 线的边界条件 ， 利用 线上的计算结果，根据取定的CFL数和格式的稳定性条件确定 方向上的网格尺寸 （也就是空间推进的步长） 其中 。 动量方程的差分格式是 的显式格式 （对于无压力梯度的平板边界层流动， ， 上式右边方括号内的最后一项为零。） 对 ，由此式依次计算出 ，直到在某一个 值 处，计算出的 满足条件 此时可认为已经到达边界层外缘，停止 线上 的计算。同时也确定了 线边界层外缘的位置 ， 在求出 之后，连续性方程的差分格式成为 的显式格式 可对 ，依次计算出 。 如果 ，将 加 ，返回；否则，结束计算。 考虑雷诺数为400的具体算例（无压力梯度）。 建议计算参数按下述方式取值： ， 方向取300个网格，即 ， ； CFL数取0.8； 推进1000000步（）。 下面的两张图给出了计算结果。一张是平板边界层的外缘，另一张是五个位置的速度剖面，并与Blasius解进行了比较。 8