基于浸入边界方法的二维流场数值模拟.docVIP

第１０卷第１６期２０１０年６月１６７１—１８１５（２０１０）１６－３８２３—０５ 科学技术与工程 ＳｃｉｅｎｃｅＴｅｃｈｎｏｌｏｇｙａｎｄＥｎｇｉｎｅｅｒｉｎｇ Ｖｏｋ１０Ｎｏ．１６ Ｊｕｎｅ２０ｌＯ ④２０１０ Ｓｃｉ．Ｔｅｅｈ．Ｅｎｇｎｇｒ 力学 基于浸人边界方法的二维流场数值模拟 冯春娟１ 刘锋１，２ （西北工业大学翼型、叶栅空气动力学国防科技重点实验室１．西安７１００７２；美国加州大学尔湾分校机械与宇航工程系２，加利福尼亚州，尔湾９２６９７－－３９７５） 摘要采用浸入边界方法进行了圆柱绕流的数值模拟。浸入边界方法是通过在Ｎａｖｉｅｒ—Ｓｔｏｋｅｓ方程中加入一项体积力，以满足物面边界的无滑移条件。整个流场在笛卡尔网格上求解，物面边界用均匀分布的拉格朗日点表示。拉格朗日点上的力通过离散了的Ｎａｖｉｅｒ—Ｓｔｏｋｅｓ方程求解后再分布到周围的网格点上。在有限差分法的基础上利用浸入边界方法分别计算了Ｒｅ＝４０和Ｒｅ＝１００时圆柱绕流，并将升力系数，阻力系数，斯特鲁哈尔数等结果与文献中的数值和试验结果进行了比较。关键词浸入式边界方法中图法分类号０３５１．３； 圆柱绕流 拉格朗日点 Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程 文献标志码Ａ 浸入边界方法最初是由Ｐｅｓｋｉｎ…提出并用于模拟人类心脏中的血液流动。这种方法最主要的特征是：整个计算在笛卡尔网格上进行。因为物面与网格线相交，要正确施加边界条件，就需要修改物面边界周围的运动方程。浸入边界方法是通过在Ｎａｖｉｅｒ．Ｓｔｏｋｅｓ方程右端加入一项体积力来体现物面边界对流场的影响。根据体积力的两种不同的处理方式，浸入边界方法可以分为连续力方法（ｃｏｎｔｉｎ－ ＵＯＵＳ 廊（ｔ）。其中ａ，卢为两个经验常数，通常会对计算 时间步长造成严格的限制。Ｓａｉｋｉ＆Ｂｉｆｉｎｇｅｎ【３１利用这种反馈力的浸入边界方法计算了低雷诺数（Ｒｅ、＜４００）时的静止，旋转和做横向振动的圆柱，说明了这种方法可以处理固体边界甚至运动边界问题。 反馈力方法有严格的稳定性限制，难以应用到三维复杂的几何外型，Ｍｏｈｄ．Ｙｕｓｏｆ［４１提出了一种直接力的浸入边界方法如下：Ｎａｖｉｅｒ—Ｓｔｏｋｅｓ方程在时 ｆｏｒｃｉｎｇ ａｐｐｒｏａｃｈ）和离散力方法（ｄｉｓｃｒｅｔｅ ｆｏｒｃｉｎｇ ａｐｐｒｏａｃｈ）。连续力方法是先求出体积力并加入Ｎａｖｉｅｒ．Ｓｔｏｋｅｓ方程中再将方程离散，而离散力方法则是先将Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程离散，再修改物面边界附近离散后的Ｎａｖｉｅｒ．Ｓｔｏｋｅｓ方程以满足无滑移边界条件。因此，连续力方法体积力的求解不依赖与方程的空间离散格式，而离散力方法则与方程的离散格式紧密相关。 所要处理的物体是固体边界，对于这种边界的处理Ｇｏｌｄｓｔｅｉｎ ｅｔ 间方向上被离散为生｝＝ＲＨＳｈ寺＋／”，其中 ｆ＋Ｉ ｆ 。 ， ４‘ ＲＩｆＳ“彳包括对流项，黏性项和压力梯度项。为了满足物面边界条件扯“１＝矿“，则物面上的体积力应 ． ． Ｔ矗＋１ｌ 为／峙＝一ＲＨＳ“寺＋Ｌ—ｆ旦。采用这种方法计算 ／．Ｉｔ 不再受到稳定性的限制，也不需要经验常数去构造体积力。直接作用力浸入边界法已经成功模拟了一些三维复杂流场问题¨】和运动边界问题［６】。 本文的浸入边界方法是基于Ｍｏｈｄ—Ｖｕｓｏｆ＿［４ｏ的方法，并进行了一定的改进：（１）物面用均匀分布 ａｌ【２１曾提出一种反馈力的浸入边界 方法，这种方法中物面边界对流场的影响通过一个 一 反馈力来体现，求解公式为声（ｔ）＝ａＩ．五（丁）打＋ ２０１０年３月２０日收到 的点表示，称之为拉格朗日点。将体积力直接构造与拉格朗日点，因此物面边界条件可以直接得 到满足。而Ｍｏｈｄ．Ｙｕｓｏｆｔ４１则是将体积力设置在物 万方数据 科学技术与工程１０卷 体内部最靠近物面的网格结点上，这些点上的速度边界条件需要通过相应的物体外部流场内的点的速度和物面上的点的速度插值得到。因此需要判断物体内部最靠近物面的结点及其相应的流场结点和边界结点的位置，这对于复杂的物面边界很难做到。（２）为了得到拉格朗日点上的体积力，需要根据网格点计算的速度插值出物面点的速度，求出物面结点上的体积力后又需要将其外推 至周围结点，此处采用的插值／Ｊ＇ｂ推公式是Ｓａｉｋｉ＆Ｂｉｒｉｎｇｅｎ【３１中提出的双线性插值（即面积加权插 值）。 １数值方法 笛卡尔坐标系下非定常不可压Ｎａｖｉｅｒ．Ｓｔｏｋｅｓ方程可以表示为 鲁＋警＝一老＋瓦１面０２Ｉｇ虿ｉ＋ｚ （１）垫：０ （２） ｄ菇‘ 式（３）中气是笛卡尔网格坐标，％是对应的速度向量，Ｐ是压力Ｚ是加入的体积力。上式中对流项采用三阶Ｒｕｎｇｅ—Ｋｕｔｔａ（ＲＫ３）格式，黏性项采用Ｃｒａｎｋ—Ｎｉｃｏｌｓｏｎ半隐式格式；其中空间离散均采用二阶中心差分式。则离散后的Ｎａｖｉｅｒ．Ｓｔｏｋｅ