附录F-有限体积法计算方腔流(F).docVIP

附录F二维不可压缩黏性流体方腔流动问题的

有限体积算法与计算程序

二维方腔流动问题是一个不可压缩黏性流动中典型流动。虽然目前尚不能求得它的解析解，但是它常被用来作为检验各种数值算法计算精度和可靠性的算例。文献中几乎大多数算法都对它进行过计算。在本算例中采用有限体积算法三阶迎风型离散格式对它进行数值求解。同时，为了初学者入门和练习方便,这里给出了用语言和语言编写的计算二维不可压缩黏性方腔流动问题计算程序，供大家学习参考。

F-1利用有限体积算法三阶迎风型离散格式求解

二维不可压缩黏性流体方腔流动问题

1.二维不可压缩黏性流体方腔流动问题

二维不可压缩黏性流体方腔流动（cavityflow图F.1二维不可压缩黏性方腔流问题示意图）：有一正方形腔室，其量纲为一的宽度为，里面充满静止的不可压缩黏性流体，方腔内初始时刻压力和密度为它周围壁面（左右壁面和底面）固定不动，上壁面以量纲为一的速度沿着上壁面方向自左向右运动（图F.1）。

图F.1二维不可压缩黏性

方腔流问题示意图

2.基本方程组、初始条件和边界条件

图F.1二维不可压缩黏性方腔流动问题示意图设流体是黏性流体。二维方腔流动问题在数学上可以由二维不可压缩黏性流动N-S方程组来表示，把它写成通用变量的微分方程组形式，有：

图F.1二维不可压缩黏性方腔流动问题示意图

（F.1）