水珠滴落的最小二乘粒子有限元方法模拟.pdfVIP

第 57卷 第 11期 2008年 11月 物 理 学 报 V01．57，No．11，November，2008 1000．3290／2008／57(11)／6722—08 ACTA PHYSICA SINICA ⑥2008Chin．Phys．Soc． 水珠滴落的最小二乘粒子有限元方法模拟 汤 波 李俊峰 王天舒 (清华大学航天航空学院，北京 100084) (2007年 12月 17日收到 ；2008年4月 10日收到修改稿) 提出了一种最小二乘粒子有限元方法 ，用其模拟了二维水珠滴落水面并飞溅散开 的过程 ．该法基于拉格 朗 日 描述 ，在每个时间步上使用扩展的Delaunay划分更新计算网格 ，并应用 a形方法识别 自由面形状 ；用最小二乘有限 元方法离散流体运动的Navier．Stokes方程 ，并推导了一种 自适应时间步长方案 以提高计算效率和鲁棒性 ；引入 网格 拉伸技术修正减小流体质量误差 ．对水滴飞溅进行仿真，得到了与商用软件 Flow一3d比较符合的结果，且具有更清 晰锐利的 自由面 ． 关键词 ：滴落 ，网格划分，a形，最小二乘有限元 PACC：0340G，4710，4755K，9530L 但都不能用来处理 自由面的拓扑变化．无 网格方 1．引 言 法 由数据点直接计算插值函数 ，通过配点法 或 背景网格积分来离散方程得到数据结点的速度，通 液珠撞击水面，反弹并散开 ，溅起漂亮水花 ，这 过拉格朗 日描述直接进行结点更新 ．由于它没有网 种复杂而美丽的图景是 自然界以及工业进程 中的一 格也就不存在网格的畸变 ，它的难点是如何施加强 种基本现象 ¨，譬如雨滴滴落到池塘里 ，空气中的水 制边界条件 以及进行精确积分h ．Gonzalez等 使 滴溅落在结冰的机翼上．此 问题可以归结为带 自由 用无网格 自然元方法 ，通过 a形算法识别边界来引 面的不可压缩流体运动问题，它在计算上存在如下 入强制边界条件，成功计算了坝溃等 自由面问题，它 两个困难 ． 解决 了边 界条件 的问题，但仍不能精 确积分． 困难之一是不可压缩 Navier-Stokes方程 的求 Idelsohn和 Onate 提 出 的 粒 子 有 限 元 方 法 解 ．方程中出现的强非线性的对流项，和作为拉格朗 (PFEM) 在使用 a形算法识别边界的基础上 ， 日乘子引入的压力项 ，都给计算带来很大的挑战 ． 引入 了网格重新划分技术得到计算网格 ，以克服精 第二个困难是 自由面的跟踪或捕捉 ．界面捕捉在 确数值积分的困难 ．由于方法较新，现阶段对 PFEM 欧拉描述下的固定 网格内，通过 VOF或者 LevelSet 的研究以Idelsohn和 Onate的小组为主，他们用其计 等方法的辅助辨识 出界面位置 ¨，这是当前最流行 算了坝溃 、晃动等一大类 自由面问题 ．但 PFEM文献 的自由面求解方案 ．但辅助方程 的精确求解是 困难 中，时间步长的选取比较简单，可能造成流体穿透固 所在 ，需要额外的处理和反复初始化 ．界面跟踪针 壁，所以作者在固壁外定义了一层额外的网格来阻 对拉格朗 日描述下的移动网格 ，自由面与计算 网格 止液体的穿透--．本文推导了一种 自适应时间步长 边界重合，在初始时刻指定并在计算过程 中保持不 格式，不再需要这层额外的网格，同时也提高了计算 变 ，用它可以得到精确的界面位置．近些年来 ，跟踪 的效率和鲁棒性 ；针对现有 PFEM文献中没有关于