论文部分内容阅读
本文提出了一种新的数值计算方法来模拟不可压粘性流动,这种方法是针对传统浸入边界法的改进。传统浸入边界法采用的是直角网格,浸入边界对周围流场的影响是通过力的形式来体现的,通常的做法是把浸入边界上的集中力分布到其周围的流场中,在分布力的作用下,重新求解控制方程得到新的流场。通过不断地迭代,试图满足边界无滑移条件。传统浸入边界法的主要优势在于网格生成简单,效率高。同时它也存在两大缺点。一是有流线穿透的现象,二是边界附近精度很低。造成这两大问题的主要原因是在传统浸入边界法中无滑移边界条件只得到近似满足。受到当地DFD(Domain-Free Discretization)方法的启发,我们把当地DFD方法中边界处理的方式与传统的浸入边界法结合起来,这就是本文提出的新方法。这种方法的主要步骤是,先在规则的网格上求解流场控制方程,再去修正边界周围的速度场使得边界能够严格满足无滑移条件。最后,应用改进后的浸入边界法做了一定量的数值计算,模拟了多种雷诺数条件下的静态圆柱绕流和翼型绕流,也模拟了较为复杂的动边界问题,如不同频率条件下的圆柱振荡。上述算例的计算结果与前人的实验或数值结果相符合,从而验证了该方法的有效性。