应用五阶WENO格式和MGFM计算多介质可压缩流.pdfVIP

472 I革等：应甬五阶WENO格式和MGFM计算多介质司医缩漉 应用五阶WENO格式和MGFM计算 多介质可压缩流 王 革，谢昌坦 (哈尔滨工程大学航天与建筑工程学院，哈尔滨150001) 摘要：多介质可压缩流计算的主要困难之一在于如何合理的定义界面边界条件，在单介质流中广泛应 用的高阶WENO格式并不能直接用来模拟多介质可压缩流问题。本文采用Levelset方法捕获物质界面， 求解Euler方程，并在此基础上计算了几个经典的一维和二维数值算例。结果表明，改进虚拟流体方法在 计算多介质可压缩流中具有较强的通用性，五阶WENO格式对激波等间断具有较好的捕捉效果和较高的分 辨率． 关键词：五阶WENO格式：改进虚拟流体方法；Levelset方法；Euler方程 0引言 面的压力、速度及其两侧的熵，再将熵拓展到 Ghost区域，由此解决了传统GFM不易解决的强 多介质可压缩流的数值模拟是当今计算流体 间断问题，而且守恒型误差也得到有效控制。 力学领域研究的热点与难点。在模拟单介质流中 关于多介质可压缩流问题，作者已进行过初 广泛应用的高精度格式如TVD、ENO和WENO 格式等并不能直接用于求解多介质可压缩流问 题，这是因为在多介质流中物质界面两侧的流体 程序对激波和膨胀波的捕获效果一般，且二维问 ～般具有不同的状态方程，如果不加处理的在整 题的对称性也不够好。本文则在前面工作的基础 个流场中采用同一种计算格式，常常会在界面附 上进行了改进，应用Levelset方法追踪物质界 近出现明显的非物理振荡甚至使计算无法进行。 面，采用MGFM定义虚拟节点参数，采用五阶 为了解决这一难题，各种改进方法相继被提出 来，其中最简单且最有价值的是Fedl【iw等于1999 Euler方程的空间项和时间项，并对几个典型的一 年提出的虚拟流体方法(GFM)tq和Liu等于2003 维和二维多介质可压缩流问题进行了数值模拟， 年提出的改进虚拟流体方法(MGFM)【21。GFM的基 获得了较为满意的计算结果。 本思想是利用Levelset方法捕获物质界面，在界 1控制方程 面附近定义3．5个虚拟节点，虚拟节点上的压力和 法向速度取真实流体的压力和法向速度，密度由 1．1二维Euler方程 等压装配法单侧熵插值求得。该方法在处理弱激 波和中等强度激波与物质界面作用的问题上取得 不考虑传热、粘性及化学反应等的影响，描 了非常好的效果，但在处理强间断问题时仍然会 述多介质可压缩流的方程为二维Euler方程 在界面附近产生明显的非物理振荡。MGFM利用 (1) ～ iOU+_OF(U)+掣：o 新 锄 加 特征关系和激波跳跃条件，通过迭代方法求出界