基于自适应有限元的浸入边界法.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：123597 基于自适应有限元的浸入边界法 1,2 1 1,2 1,2 1,2 王志成 ，唐大伟 ，温娟 李宴君 张璟 , , ⒈ 中国科学院工程热物理研究所 北京 ⒉中国科学院研究生院 北京 ( , 100190; , 100039 ) (Tel: 010 Email: wangzhicheng@iet.cn) 摘 要： 本文发展了一种基于自适应网格的浸入边界法(Immersed Boundary Method,IBM)， 其特点是用一套简单、固定的网格处理多个具有复杂几何机构的计算，同时也用于移动界面 模拟。用有限元离散不可压缩N‐S 方程，压力与速度间用投影法耦合。通过引入虚拟网格 (Ghost Cell)满足浸入边界上的无滑移无穿透边界条件。消减虚拟网格内压力额外的速度散度 实现浸入边界上压力梯度为 的条件。最后，我们用浸入边界法研究了 数为 ，振荡 0 Re 100 频率0.8 1.2 范围内，圆柱的绕流特性。 关键词：浸入边界法；虚拟网格；圆柱绕流；振荡 0前言 在传统的计算流体力学中，模拟带有运动边界流动的方法有动网格法或参考坐标系法。前者 需要对网格做一些复杂的处理，后者难以模拟多个固体在流场中的运动。一些学者提出在正 交网格上离散Navier‐Stokes 方程的来计算这类流动问题。其中最常用的一种是由Peskin [1] 首先提出的浸入边界法，并由Lai 等人[2]将之改进，使其具有二阶精度。更多关于IBM 的 细节，可以参考Mittal 等人[3]的综述。为了提高计算精度，浸入边界法往往要求在固体边 界附近分布足够多的计算点，即局部加密网格。开源有限元程序库Deal.ii[4]不仅为N‐S 方程 求解器的开发者提供了各种函数，还具有网格自适应的功能。我们的工作是在这个程序库的 框架内实现浸入边界法，从而得以数值研究振荡圆柱的绕流特性。 1数值方法 在我们的浸入边界法中，计算区域中的固体采用一个指示函数 来表征。当这个函数的 值等于 时，表明所在区域为固体，函数值等于 时，表明处在流体区域。 1 0 无量纲后的二维流动控制方程为： (1) (2) 在这里， 是雷诺数。方程 中的 是用来表明流体与固体之间相互作用的体积力密度。 Re (1) 由于我们选用了投影法求解上述不可压缩 方程，通过插值的方法满足浸入边界内固体的 N‐S 迁移速度，因此不需要显示地计算该体积力密度。具体来说，浸入边界法的计算步骤为： 第一步，外推速度和压力： (3) 在上式中， k 是计算推进的时间步。 第二步，预测速度： (4)