EulerLagrange耦合计算的GEL方法研究.pdfVIP

姚阳，李平，柏劲松，姜洋，裴晓阳，钟敏 (中国工程物理研究院流体物理研究所四川绵阳) 蕾要=本文研究一种Euler-Lagrange耦合数值方法一Ghost-fluid 二维动力有限元程序。通过一维黎曼问题的计算结果与高精度PPM方法进行的比较，以及二维移动边界 li丘-offproblem’的计算结果与文献【l】的对比，验证了GEL方法和本文程序的正确性。 “cylinder FluidMemod内边界条件 关键词：Euler-Lagrange耦合GEL方法Ghost 1引言 在爆炸力学问题的数值计算中，有这样一类问题：一部分介质高速流动或剧烈变形，甚至还出 现涡旋流动，另一部分介质变形相对较小，同时还需关心这些小变形介质(材料)的结构响应或动力 学响应，如水下爆炸对船体的影响，低速侵彻问题和一些爆轰驱动问题等等。对于这类问题进行数 方法在流体力学计算以及强激波和大变形流场有很强的优势，而Lagrange数值方法却在很多固体力 学计算、以及小变形流场的材料动力学计算中有较好的表现，因此，如何有效利用这两种计算方法 的优越性，使它们优势互补，是爆炸力学数值计算研究的热点和难点问题之一。 在Euler法中吸收了Lagrange法的一些做法，如跟踪界面，这种结合两者优点方法的代表是质 in 点网格法PIC(Particle 物理量，对计算机的速度和存贮量的要求都比较高【4J，而示踪粒子法难于推广至三维计算，且逻辑 和程序结构复杂。另一种是在Lagrange法中吸收Euler法的一些做法，如重分网格，以任意 则构造新网格，而在给定边界的区域上如何构造合适的网格是这种方法的难点之一14J。另外一些方 法如PISCESt61，CELl7 格而带来的特殊处理方法。但它们普遍存在的难题之一是对混合网格的处理，如CEL方法需计算 混合网格中的各力学量及它们在边界上的输运量，而在程序中实现这样的计算是一项十分复杂的工 作141。 Fluid Method【9J的思想处理因移动界 本文研究一种Euler-Lagrange耦合计算方法，借助Ghost 面出现的混合网格，文献[1】中称之为Ghost．fluid 一是对自由边界的处理过程和程序的编制要简单，计算量小；二是可以自由地将独立的Euler和 Lagrange计算程序连接起来，不必对它们进行任何实质上的改动，连接过程简单易行【1】。 2基本方程 2．1Euler方程 理想可压缩流体Euler控制方程是： V 7 等+V·(卅=o 西 昙(删+V·(p沥+IP)=o (2) c3， 昙fp(P+譬)]+V·[p面(P+万P+季)]=。 其中，p为密度，历为速度矢量，P为压力，e为单位质量内能，f为时间，I是单位矩阵。在本文 中状态方程形式取为： P=尸0，p) (4) Euler计算采用二阶SCB格式㈣编制的计算程序。 2．2Lagrange方程 连续介质力学给出的基本方程有： p+矽·疗=0 (5) p历=V·矛+6 (6) 庄=矛：V厅 (7)