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

Vok10 No．16 20lO 第10卷第16期2010年6月科学技术与工程 June Science and TechnologyEngineering ④2010 1671—1815(2010)16-3823—05 Sci．Teeh．Engngr 力 学 基于浸人边界方法的二维流场数值模拟 冯春娟1 刘锋1，2 (西北工业大学翼型、叶栅空气动力学国防科技重点实验室1．西安710072； 美国加州大学尔湾分校机械与宇航工程系2，加利福尼亚州，尔湾92697--3975) 足物面边界的无滑移条件。整个流场在笛卡尔网格上求解，物面边界用均匀分布的拉格朗日点表示。拉格朗日点上的力通 40和Re=100时圆柱绕流，并将升力系数，阻力系数，斯特鲁哈尔数等结果与文献中的数值和试验结果进行了比较。 关键词浸入式边界方法 圆柱绕流 拉格朗日点 Navier-Stokes方程 中图法分类号0351．3； 文献标志码A 浸入边界方法最初是由Peskin…提出并用于模 廊(t)。其中a，卢为两个经验常数，通常会对计算 拟人类心脏中的血液流动。这种方法最主要的特 征是：整个计算在笛卡尔网格上进行。因为物面与 用这种反馈力的浸入边界方法计算了低雷诺数 网格线相交，要正确施加边界条件，就需要修改物 (Re、400)时的静止，旋转和做横向振动的圆柱， 面边界周围的运动方程。浸入边界方法是通过在 说明了这种方法可以处理固体边界甚至运动边界 问题。 Navier．Stokes方程右端加入一项体积力来体现物面 边界对流场的影响。根据体积力的两种不同的处 反馈力方法有严格的稳定性限制，难以应用到 理方式，浸入边界方法可以分为连续力方法(contin- UOUS forcingapproach)和离散力方法(discreteforcing f+I f 。 ， approach)。连续力方法是先求出体积力并加入 间方向上被离散为生}=RHSh寺+／”，其中 4‘ Navier．Stokes方程中再将方程离散，而离散力方法 则是先将Navier-Stokes方程离散，再修改物面边界RIfS“彳包括对流项，黏性项和压力梯度项。为了满 附近离散后的Navier．Stokes方程以满足无滑移边界足物面边界条件扯“1=矿“，则物面上的体积力应 条件。因此，连续力方法体积力的求解不依赖与方 ． ． T矗+1 l 为／峙=一RHS“寺+L—f旦。采用这种方法计算 程的空间离散格式，而离散力方法则与方程的离散 ／．It 不再受到稳定性的限制，也不需要经验常数去构造 格式紧密相关。 体积力。直接作用力浸入边界法已经成功模拟了 所要处理的物体是固体边界，对于这种边界的 一些三维复杂流场问题¨】和运动边界问题[6】。 et 处理Goldsteinal【21曾提出一种反馈力的浸入边界 方法，这种方法中物面边界对流场的影响通过一个 本文的浸入边界方法