计算流体驱动方腔程序

驱动方腔（粘性不可压流），方腔示意图如左下：物理背景：如左图所示,装满黏性不可压流体的方腔上底面以速度U运动,其它壁面固定,其内部流体将作类似湍流的运动.流体运动复杂度将依赖于以下要素: (1)方腔形