基于非结构化网格和浸入边界的流固耦合数值模拟.docx

1、基于非结构化网格和浸入边界的流固耦合数值模拟 万方数据 万方数据明平剑,等:基于非结构化网格和浸入边界的流同耦含数值模拟(a原始浸入式边界方法图l不同网格方法步的值;甜;和嵋均为上一时间步的速度值;D为单元体积。浸入式边界对流体的影响采用vOs法计算。质量力按下式计算:一“,:2i1,一1,2i其中”和矿为速度二次修正值,可通过流体和浸入刚体边界速度线性插值得到11】:甜,=(1一口甜7+劬?“矿=(1一倪/+叫+1(b非结构化网格浸入式边界方法图2与浸入边界相交单元的固体体积分数325删I=Pl珀P2I挖肿2_-_l-_._Ic2一几_-_lCl钾2I鼢2I啪姗俐挖图3单元内部形成链表示意图

2、(8222计算固体体积分数其中口是单元体中的固体体积分数(单元中固体体积与单元总体积的比值,2.2.2小节将详细介绍它的计算方法。将材和,代入(5、(6得,“j 一甜六钳tF,.匕一,钳F(10压力速度耦合采用sIMPLEc算法,界面流量计算采用Rhie.Chow插值方法。各变量求解顺序安排见2.2.3。定义了一个数据结构,形成链表来表示边界,其方向是使得固体区域始终在边界的左手侧。为了计算每个单元内部固体体积分数,定义了指针型结构数组,在每个单元内形成一个多边形的链表,多边形的边界具有逆时针方向,如图2阴影部分多边形,图3给出了单元ICl和IC2内的链表,图中Q岬表示网格节点,IBP为浸入边

3、界点,IBCP表示浸入边界与网格面交点。将多边形剖分为一系列三角形计算其面积和,从而确定每个单元内固体体积分数。下面较详细的介绍多边形的形成。从图2中完成IBPl操作后,单元IC2内的链表包含节点GNP6.-IBCPl。完成IBP2后,链表为GNP6一IBCPl一IBP2,浸入边界IBP2一IBP3与单元IC2的边界a岬2.GNP3相交于IBCP2。找到交点后,IC2内链表增加节点IBCP2,同时确定位于 固体区域内的网格节点GNP2。此时单元IC2的链万方数据 万方数据明平剑,等:基于非结构化网格和浸入边界的流固耦合数值模拟图7浸入边界前沿示意图RP:型(13弛朋=坠(14fd其中圆柱运动最

4、大速度,d圆柱直径厂振荡频率。圆柱直线振荡可以用简谐振动描述:x(f=一以sin(2妒(15其中x(f为圆柱中心位置,以为振荡幅值。由位移对时间求导可以得到圆柱运动速度327计算得到的单个周期内阻力系数变化与文献【14】的结果对比见图8。图9给出了不同相位角时圆柱周围的涡量分布与文献【ll】中数值计算结果的对比。从图中可以看出本文计算结果与文献中cartesian网格下计算结果一致,表明本文方法可以很好地模拟振荡圆柱周围流场。图lO分别给出了180、210、330三个相位角情形,不同空问位置的水平和竖直方向速度分布与文献【14】中实验数据的对比。从图中可以看出,除局部位置有些偏差外,计算结果与

5、实验值吻合较好,如图lO(b右图中,数值计算结果具有较好的对称性,(乒0,0侧与实验值吻合较好,而下半部分与实验值不一致,原因可能是实验测量引起的误差。%(,=2碱,cos(2妒(163.2阀门动态工作工程模拟因此最大运动速度和位移满足以下关系:=2彤k(17计算区域为30出20d。文献【13】采用了同样的计算域研究了类似问题,坐标原点位于计算区域中心,d=o.olm,流体密度p=1kgm.3,流体动力黏度=1.Oe.5Pas。给定RP=loo和xC=5,由此可以得到=0.1m/s,=2Hz,4=O.007957m图8R萨100和柏5时,阻力系数时间历程阀门是管路系统中重要元件之一,对管路系统

6、有较重要的影响,如文献15】利用商业软件Fluent 研究了阀门的流场特性。本文对二维管道内的阀门动态工作过程进行模拟,各管段长度如图ll所示,其中桔0.2m,阀门为椭圆形,长短轴长分别为0.75 Jjl和O.25厅,长轴方向与水平方向夹角为450。进口流速为0.03111/s,|jI:O.2m。先对阀门静止时进行了模拟。图12给出了阀门后的瞬时涡量分布图,横纵坐标分别为舶和m。图13给出了阀门周围速度矢量分布和固体体积分数分布。从图中可以看到,利用本文方法可以计算出阀门的流动特征,尽管由于阀门绕流在管道内形成了非常复杂流动。图14给出了阀门由图11所示位置顺时针转过90度角的动态过程。从图中

7、可以明显地看到阀门引起的漩涡脱落过程,表明本文方法可以较好的模拟 复杂结构运动的流固耦合问题。万方数据328水动力学研究与进展A辑20lO年第3期1.00.6o.2逞h吨2吨61.O 图9一个周期中不同相位角时的涡最值与文献中计算结果对比(左为文献【l】的结果,右为本文计算结果o.015_0.0lo.5O n0o.Ol0.015u7r1.On6o.28h吨2吨6一1.0(计算值计算值计剪值(计算值(试验值试验值试验值(试验值o.015o.010.0050o-5o.Ol o.015r/Is-I 相位角l。明平剑,等:基于非结构化网格和浸入边界的流固耦合数值模拟L0o.6o.2丑、吨2吨6一1.0

8、1.Oo.60.2逞h吨2吨6一1.0吨015吨01吨6O o.0惦nOl o.015吨015吨01吨O畸O0.5nOl o.015血s-1.,/Is-lCb相位角2lO。1.O0.6o.2逞吨2吨6一1.O1.0o.6o.2逞、吨2吨6一1.Oj迭Z季lo、。:亲A、Z-.,一一-7-:j,书:-二二:耄:妒一一一。;,了一:冀.、食、会.A.、可_一二:司”“o时一:。五=lZ。.、刍o-015o.Ol0.50o.0o.0l0.015o.015_o.Ol0.0050O.005o.0l o-015 u/h-一y/h一(c相位角瑚。图lO不同相位角时计算的速度场与文献中实验值得对比图11阀门结

9、构示意图图12阀门静止时计算涡量分布图13阀门周围速度矢量和固体体积分数分布4结论329本文提出了一种基于SIMPLE类算法的模拟流固耦合作用的数值方法,空间采用非结构化网格有限体积方法离散,时间项采用欧拉隐式格式,通过浸入边界方法和固体体积法(Volume of Solid考虑流体作用在固体上的作用力,给出了非结构化网 格上流体体积分数计算方法。该方法可以很方便地330水动力学研究与进展A辑2010年第3期(a吼5s (d_2.O s(b=1.O s(e=2.5s(c=1.5s(ff:3.O s刮.O s图14阀门动态T作过程中的涡量分布图添加到已有的有限体积方法的cFD求解器。为了验证方法

10、的正确性,本文对于振荡例球绕流做了模拟,所得圆柱表面阻力系数随时问变化也与文献中实验数据吻合较好,不同相位角时的涡量分布与文献中数值计算结果一致,不同相位角和不同位黄的水平方向和竖直方向速度分布与文献中实验数据也吻合较好,表明了本文方法的合理性和正确性。此外本文还对阀门的工作进程进行了模拟,进步验证了本文方法的适用性。参考文献t【l】PEsKIN C S.Flow pattems a毋咖d heart valves:a舢merical method【J】.Jo哪l of Computatial Physics,1972,lO(2:252-271.【21PAN D.An immersed bou

11、ndary method斯income-prsible nows璐ing Volume of bodyction【J】.Intematio眦l Jmal for N啪erical Methods in Fluids.2006,50(6:733-750.r3】SHEN L,CHN E.Numecal simuJatjon of丑ujd-蛐mctIlm inte均dion u咖g a combined voIume of fluid锄d immcrsed boundary me也od【刀.oc伽Engin咖g,2008,35(8:939-952.【41王则力,罗坤,樊建人.修正内嵌边界法对颗粒运动

12、的直接模拟【J】.工程热物理学报,2008,29(1:99一101.WANG Zc.h。LUo KllllFAN Ji锄-砌1.Di删simulatiof pmticula捃f】ow using modi矗ed imma嚣。dbo岫da叮memod【J】.Jo哪al of Enginee抽g 基于非结构化网格和浸入边界的流固耦合数值模拟作者:明平剑, 张文平, 卢熙群, 朱明刚, MING Ping-jian, ZHANG Wen-ping, LU Xi-qun , ZHU Ming-gang作者单位:哈尔滨工程大学动力与能源工程学院,哈尔滨,150001刊名: 水动力学研究与进展A辑英文刊名

13、:CHINESE JOURNAL OF HYDRODYNAMICS年,卷(期:2010,25(3参考文献(15条1.UDAYKUMAR H;MITTAL R;RAMPUNGGOON P A sharp interface cartesian grid method for simulating flows with complex moving boundaries 2001(012.ROMAN F;NAPOLI E;MILICI B An improved immersed boundary method for curvilinear grids 2009(083.林景松;王平阳;高红液体

14、推进剂管路局部损失及瞬变流动数值模拟期刊论文-水动力学研究与进展A辑2008(014.DUTSCH H;DURST F;BECKER S Low Reynolds number flow around an oscillating circular cylinder at low Keulegan-Carpenter numbers 19985.SHEN L;CHAN E;LIN P Calculation of hydrodynamic forces acting on a submerged moving object using immersed boundary method 2009

15、(036.明平剑基于非结构化网格气液两相流数值方法及并行计算研究与软件开发 20087.NG K C A collocated finite volume embedding method for simulation of flow past stationary and moving body 2009(028.MING Ping-jian;ZHANG Wen-ping Numerical simulation of low Reynolds number fluid-structure interaction with immersed boundary method 20099.GHI

16、AS R;MITTAL R;DONG H A sharp interface immersed boundary method for compressible viscous flows 外文期刊 2007(0110.宫兆新;鲁传敬;黄华雄浸入边界法及其应用期刊论文-力学季刊 2007(0311.钟国华;梁岸;孙晓峰基于浸入式边界的流固耦合的非定常数值模拟研究期刊论文-工程热物理学报 2007(0312.王则力;罗坤;樊建人修正内嵌边界法对颗粒运动的直接模拟期刊论文-工程热物理学报 2008(0113.SHEN L;CHAN E Numerical simulation of fluid-structure interaction using a combined volume offluid and immersed boundary method 2008(0814.PAN D An immersed boundary method for income-pressible flows using volume of body function外文期刊 2006(0615.PESKIN C S Flow patterns around heart valves:a numerical method 1972(02本文链接: