二维流体力学有限差分和有限元程序的耦合计算.pdfVIP

第 1期 高 能 量 密 度 物 理 NO．1 2009年 3月 HIGH ENERGY DENSITY PHYSICS M ar．，2009 二维流体力学有限差分和 有限元程序的耦合计算 李 蕾，柏劲松，刘 坤，李 平 (中国工程物理研究院流体物理研究所 105室，四川绵阳 621900) 摘 要：根据 压 力和法 向速度 连续 准则，将 Euler方 法为基础 的 MFPPM (Piecewise—ParabolicMethod)程序和 Lagrange方法为基础 的 DEFEL(2一D Finite ElementsCode，二维流体弹塑性动力有限元)程序进行耦合 ，发展 了基于 levelset 的GEL(Ghost—fluidEuler-Lagrange)方法。该方法在处理大变形流场与小变形结 构以及复杂流动与多物体相互作用等问题具有优越性。通过二维算例的计算结果与 文献比较 ，检验 了GEL方法和耦合程序的正确性，并对水下爆炸形成的流场对多物 体作用过程进行 了数值模拟 。。 关键词 ：GEL方法；Euler—Lagrange程序耦合；数值模拟 1 引 言 由于数值模拟中遇到的很多实际问题涉及到多介质流或流 固之间的相互作用，采用纯 Euler 数值方法或纯 Lagrange数值方法都难 以得到满意的计算结果，因此产生了Euler和 Lagrange混 合或耦合方法。在计算物理和计算数学的研究课题 中，该方法的研究 日益受到关注，已成为爆炸力 学数值计算研究的热点 问题之一。在爆炸力学数值计算中，通常一部分介质高速流动或剧烈变形 ， 甚至还出现涡旋流动，另一部分介质变形相对较小，而我们关心的则是这些小变形介质(材料)的结 构响应或动力学响应(如水下爆炸对船体的影响，低速侵彻 问题和一些爆轰驱动问题等)。对于这 类问题的处理，往往希望在整体计算区域 中流体或大变形的部分采用 Euler方法计算，对小变形 以 及固体材料的部分采用 Lagrange方法计算，以期获得对整体计算区域的不同区域部分更细致的描 述 ，而这正是 GEL(Ghost—fluidEuler—Lagrange)Ll方法的优势和特点所在。 本文研究的GEL方法 ，借助 GFM (GhostFluidMethod)-2和 levelset函数思想处理因移动 界面出现的混合网格，即对界面处的Lagrange域用 ghost网格覆盖 ，ghost网格的压力，密度和切 向速度由Euler域外插赋值 ，界面的法向速度由Lagrange计算给出。Lagrange域在边界受到来 自 Euler域真实流体作用力，由此给出界面上 Lagrange节点受到的力。这个过程体现了Euler域和 Lagrange域的相互作用。GEL方法主要有两个特点：一是对 自由边界的处理过程和程序的编制比 较简单 ，计算量小；二是可以自由地将独立的Euler和 Lagrange计算程序连接起来，不必对它们进 行任何实质上的改动，连接过程比较简单易行。本文中Euelr域采用多流体 PPM(Piecewise—Para— bolicMethod)方法求解[3]，Lagrange域采用 DEFEL(2一DFiniteElementsCode，二维流体弹塑性 动力有限元)程序求解 。 * 基金项 目：国家 自然基金项 目，中国工程物理研究院科学技术发展基金(2008B0201011)资助 。 第 1期 李 蕾等：二维流体力学有限差分和有限元程序的耦合计算 11 2 基本方程 2．1 Euler方程 理想流体的Euler动力学方程为 f+ ·()一o l+(1f·)ll+-~-Vp=F (1) I箦+(1I·V)E+lpV·一R 其 中，E===pe+ 。／2；状态