Level+Set方法和多介质可压缩流研究.pdfVIP

第25卷增刊 计算力学学报 V01．25，Sup． Mechanics December 2008年12月 a面雠J0明词ofComputational 2008 Level Set方法和多介质可压缩流 王革‘， 张斌 (哈尔滨工程大学航天工程系，黑龙江哈尔滨150001) 摘要：多介质可压缩流问题计算的关键是如何精确的捕获不同时刻物质界面的位置，从而将多介质问题分解 成多个单介质f=l题去处理．LevelSet方法的优点是不用显示的追踪物质界面，而用距离函数就能精确定位界面． 同时，用Level Set方法追踪界面运动易于处理界面拓扑结构的变化、易于处理大变形问题．本文成功地将Level Set方法应用在二维多介质可压缩流计算． 关键词：EuIcr方程：Level Set方法：MGFM方法：Iscntropic修正 中图分类号：037 文献标识码：A 2．1 Euler方程 1 引 言 二维可压缩流的基本方程是Eulcr方程： 多介质可压缩流涉及的范围很广，如激波与 U+F(u)，+G(u)，=o(1) 物质界面的相互作用，液体中气泡的变形和坍塌 U=(p，pu，∥，E)。 等。传统的多组分单相流在处理这类问题时遇到 了困难。现在对于多介质可压缩流的模拟方法有 e(u)=(肛，pu2+p，胛，“(E+p))2 Set[I】 很多种，如拉格朗日法、阵面跟踪法、Level G(u)=(∥，肛n，，∥2+p，V(E+p))7 方法、基于Godunov格式的两相流方法，以及格子 BGK方法等。1992年，Muldcr等人【2】首次用Level Set方法计算了可压缩两介质气体界面问题。但是 式中E：肛+旦幽，表示单位体积的 2 在接触间断处产生了很大的数值耗散，因此后来 总能量。P，P，e分别表示密度、压力和单位质量 有人陆续提出了不同的改进方法。近年来，由 的内能。甜，v分别表示x，Y方向上的速度。气体 Fedidw等人提出的GhostFluid方法在计算多介质 采用理想气体状态方程：p=p／驴一1)；而液体 流体问题时取得了良好效果。Liup。5】在此基础上又 采用Tait方程：pe=(p 提出了一种改进的虚拟流动方法(MGFM)。 程采用HLLC方法【61求解。 本文首先编制了求解LevelSet方程的程序，并 2．2 Level Set方程 用经典算例验证了程序的正确性。随后将求解 Level 由于界面随流体质点一起运动，描述界面位 Set方程的程序和求解Euler方程的程序结合 置的运动方程即为levelset方程。 起