通过在界面标识方程中引入人工压缩项,采用基于非结构网格的高分辨率有界格式,编程实现了一类适合模拟复杂界面流问题的VOF(Volume of Fluid)方法.这类VOF方法可以自动处理界面的拓扑演化,无须进行实时的界面重构,提高了计算效率.应用这种基于有界压缩思想的界面流模拟方法,对三维Rayleigh-Taylor不稳定性问题进行了并行模拟.计算结果与相关文献的数据吻合较好,证明了算法的可靠性.研究了密度差和Reynolds数对界面演化的影响,分析了在不同Atwood数下鞍点结构演化规律的异同,对于Reynolds数小于282的界面流问题黏性起着很明显的作用,而Atwood数的影响限于低密度差的界面流动问题,即Atwood数小于0.90的情形.
