统一坐标系下多介质流体力学计算的一种方法.pdfVIP

维普资讯 计 算 物 理 第23卷 第 1期 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Vo1．23．No．1 2006年 1月 Jan．，2006 [文章编号] 1001．246X(2006)01．0019．06 统一坐标系下多介质流体力学计算的一种方法 贾鹏彦 (中国工程物理研究院北京研究生部，北京 100088) [摘 要] 采用一种基于网格节点之间的吸引和排斥重新分布节点的移动网格方法确定网格移动的速度，方法 简单，易于实现 ．数值结果表明，本方法不仅可以显著地提高计算效率 ，而且仍然具有较高的流场分辨率 ． [关键词] 统一坐标系；移动网格方法；多介质流体 [中图分类号] 0241 [文献标识码] A O 引言 随着计算机技术的飞速发展，计算流体力学已经成为研究流体运动规律，解决工程实际问题的重要手段 之一．同时，用于流体力学计算的数值方法也得到迅猛发展，目前已经取得 了很多重要的研究成果． 二维非定常的流体运动是很复杂的运动 ，在流场上可以发生扭曲、涡流与滑移等现象 ．在有多种介质时， 这些机制会显出复杂性来 ．不同介质的界面运动的不稳定性使不同介质混淆起来，而一维流体运动因为流体 介质前后按序不能相互超越而简单得多．这些情况使得数值求解多介质流体力学运动方程时，二维问题 比一 维问题复杂很多．在实际计算中，很难有统一的格式计算各种问题，因此高维多介质流体力学数值方法仍然 是重要的研究课题 ． 近年来 ，Hui等n 提出了统一坐标系的思想，并对单介质流体力学问题进行数值模拟，取得了比较理想 的结果 ．在统一坐标系下，流体动力学 的各物理量看成是时间和伪粒子 (pseudo．particles)的某种 固有特征函 数，伪粒子运动速度为 幻 ，蟹是流体质点速度 ．统一坐标包括 Euler坐标和 Lagrange坐标两种特殊情况，当 h = 0时为Euler坐标，h=1时为 Lagrange坐标 ．在二维无粘流动中，自由函数 h的选择为保持网格角时，可以 避免 Euler坐标系下接触间断处的过度数值耗散 以及 Lagrange坐标系下严重的网格扭 曲，同时在滑移线附近 仍然具有很高的分辨率 ． 由于统一坐标系相对于传统的Euler坐标和 Lagrange坐标有较明显的优势，近来我们将二维单介质的流 体计算扩展到多介质的计算，取得 了比较好的结果 ]．在文[4]中，对 自由函数 h的选择仍采用保持网格角 的方法 ．数值实验表明，虽然方法可行，但由于在求解确定 h的偏微分方程时，要采用迭代法，所 以在部分算 例的计算 中，迭代收敛很慢，计算耗时长，计算效率低 ． 本文用一种基于网格节点之间的吸引和排斥重新分布节点的移动网格方法确定 h的值，不需要迭代求 解确定 h的偏微分方程 ，计算效率显著提高，而且仍然具有较高的流场分辨率 ． 1 统一坐标系下的扩展 Euler方程 二维可压缩多介质 7模型(考虑理想气体情况) D “ Dv D D “ pu +P D “ D “ pv +P a Dv a a +