用FEMLevel+Set方法数值模拟RayleighTaylor不稳定性问题.pdfVIP

第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 用FEM．LevelSet方法数值模拟 万德成 (』：海交通大学船舶海洋．与建筑丁程学院海洋T程国家霞点实验室200240) 摘要： 本研究采用FEM．LeveI set方法中的距离 性流动问题进行数值模拟。应用QuanapelIe—selmin四阶格式离散求解Level 函数纯对流方程；通过有限元的弱积分形式，把表面张力自动包含在N．S方程的动量方程中， 避免了人为引进参数问题。本研究分别考虑了界面有表面张力和无表面张力存在的情况，模 拟一个直立二维管道中上面比重大的流体落到下面比重较轻流体中的流动失稳情况。通过计 以看出本研究的数值方法能较好模拟出Raylei曲．Taylor不稳定性现象。 关键词： FEM；Level 1引言 问题，而且还要考虑界面的表面张力影响。冈此对数值方法和格式的要求就更高。比如传统 的VoF界面数值模拟方法就不能高精度地给山界面的形状，从而影响对界面表面张力的计算。 虽然LevelSet方法可以较为准确给出界面的形状，但由丁二距离函数的纯对流方程数值求解中 的数值黏性问题，使得Levelset方法给出的界面形状容易失真，而且质量不守恒。此外，在 考虑表面张力时，一般需要引进人Ij厚度来处理表面张力的计算I卜珥j。 本研究把Level LevelSet方法中的距离函数纯对流方程，保证距离函数求解的高精度，并同时进行距离函数 的重初始化修止和质量守恒校止，确保LevelSet的距离函数为真实的距离函数；应用二步算 子分裂法处理N．S方程的求解，简化速度与压力的耦合求解困难；通过有限元的弱积分形式， 可以把表面张力臼动包含在N．S方程的动量方程中，这样就避免了人为引进参数问题。本研 究采用上述的FEM．Level 了数值模拟，分别考虑了界面有表面张力存在和无表面张力存在的情况，模拟得到了上面比 重大的流体落剑-卜．面比重较轻流体中的流动火稳情况。通过计算发现表面张力能否被止确计 流动问题的计算。 ．264． 第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 2控制方程 的流体时出现的流动欠稳现象。我秆J这里考虑的两相流体都是不可压缩流体，每相不可压流 体流动的控制方程为 V．y=o (1) 以f：}+矿．Vyl—V．(2鸬s)+vp=pg，f，1．曰，f=l，2，… ＼仞 ／ 这里p和“分别是第f相流体的密度和动力黏性系数，p为压力，y为速度，g是重力加速 度，s是流体变形率张量，可表示为 1 r 一1 ，、l Vy+(Vy)1、 ， I (2)’’ s=寺I I 令力=曰u露u厂为整个计算区域，其中厂=露n瑶为两相流体的交界面。在交界 baIance 面厂上要满足下面两个条件：第一是质量平衡方程(m舾s equation)。由于没有质量穿透 交界面，因此交界面应随流体一起运动，即： y=K·n=％·n (3) 这里y是交界面上的速度，一是交界面的单位法向矢量，K和巧分别是两相流体的速度。第 balance 二是动量平衡方程(momentumequation)。在交界面上要满足应力平衡条件，即： I石一疋)·珂=舸，? (4) 这里Z=一Ⅳ+2∥，S，，是单位张量，盯是交界面表面