基于浸入式边界的流固耦合的非定常数值模拟研讨.pdfVIP

中困T程热物理学会 热机气动热力学 学术会议论文 编tj．：062044 基于浸入式边界的流固耦合 的非定常数值模拟研究 钟国华，梁岸，孙晓峰 北京航守航天人学能浮与动力T程学院100083 lt【话mail：sxf@buaa．edu．cn 摘要 本文针对T程中关心的流例榴合以及复杂边界物仆绕流问题，采用算了分裂的方法求解-二维小 口r压N--S方程，采用完令正交的矩形刚格，用浸入式边界模拟静止和运动物体的边界。使用这种方法， 本文分别对低雷诺数下的光滑1‘板绕流．网桂绕流．低KC数下的振荡同柱进行J，数值模拟。计算结果 与以往的数值和试验结果吻合的很好．hF明厂浸入式边界方法Jf几所发展的程序的町靠性。 关键词浸入式边界，流l剞耦合．1F定常 1．引言 边界的处理一直是计算流体力学的一个重要环仃。对丁．不规!l!lJ的i又：域，一种直接和 成熟的处理方法是生成贴体网格，通过坐标变换将不规则的物理域映射剑规则的计算 域。而贴体网格的生成一股要迭代求解微分方程，这需要花费一定的时间。如果在计算 过程中边界保持静I}：，那么生成网格所需要的时间可以不予考虑。对于运动边界问题． 一个自然的思路就是对每一个时间步都生成新的一套贴体网格，如果采用这种方法，人 部分的计算时间将花费在网格的生成上。从目前的研究成果和趋势来看，通过凋整网格 米适应边界运动的处理方法，都不可避免地带来巨人的时间消耗。浸入式边界方法由丁． 避免了重新生成网格，成为处理运动边界问题的趋势。由丁．不需要生成贴体网格。我们 就可以使用最简单的矩形网格，将复杂的边界“浸入”在网格中。这也止是“浸入式边 界”名称的由米。 浸入式边界方法最甲．是1972年由CharlesPeskin和DavidMcQueen13．4J等人 提出并应用T心脏内血液流动的模拟，随后铍人鼙的麻刖丁．各种流动问题。90年代之后， 浸入式边界方法得剑了极人的重视，一方面，方法本身得剑不断的改进和完善，浸入式 边界与有限著分、有限体积、有限元、谱方法的结合以及高精度边界插值方法都得到了 发展；另一方面，方法的应用领域得剑了很人的扩展，从最初的生物医学扩展到湍流的 直接数值模拟、多相流动等等。 nl+轮机nl。片气动弹性问题一直是r：程上的焦点，由丁．这是复杂边界流同耦合的非定 常问题，使刚传统的方法处理凼难比较人，目前住这方面进行的数值模拟I：作还比较少。 止如前面所描述的，浸入式边界方法最显蒋的一个优势就是对复杂和运动边界的处理。 要描述边界的运动，只需改变边界点的坐标和给出边界的运动速度。使川浸入式边界也 许是一个解决这类问题的一个契机，这止是作者进行尝试的ilj发点和立足点。 权平均的方式求解流体与同体相Ⅱ耦合的过袢。本文使刚的方法就是基丁．文献『20J。 2．数值模拟方法 一：维1卜定常不可FR,N-S方稃可以弓为： 712 塑+(u．v)u：一vp+上v2U+F(1) 1 Ot Re V．U=0 (2) 其中，是时间，U(u，1，)是速度．P是压力，Re是雷诺数．F是力源项． 的数值方法进行求解，这四步川时间上标表示如卜．： ／,／”--“肿174_U”+’72一“”374寸／,t”’ (3) 第一步，求解对流项 (“”叫甜川74) 詈一(u．V)u“) 第二步 求解力源项 (甜”+174叫甜肿172) 空：F ) @