可压缩不可压缩流体交界面高精度数值方法研究.pdf

摘 要 Meshkov不稳定性以及表面张力诱导的液滴动力学等．此外，本文还考虑了高效并行 计算方法在界面数值模拟中的应用，并讨论了浅水自由面问题的数值模拟方法． 计算模型和经典的PPM方法数值模拟高密度比的多介质可压缩流动问题．同时本 AC／WENO—AC)相结合应用于一般状态方程的可压缩多介质流动模拟．算例表明两 种方法可以有效应用于可压缩多介质交界面的数值模拟，相对于多介质PPM方法， Level—Set界面跟踪方法能更好地无数值耗散地确定交界面的位置，适合于高分辨率交 界面计算，但是后者的健壮性不如前者，对于状态方程差异较大且高密度比的情况还 不能够很好地处理．本文还探讨了计算一般状态方程多介质可压缩流动问题的数值方 法，本文应用Level—Set技术捕捉界面位置，在界面附近采用守恒数值离散，并用双波 近似求解一般状态方程Riemann问题．一维算例表明，该方法能准确地模拟交界面的 位置，可处理任意状态方程的多介质问题． G脯st，Fluid方法的并行算法，以进行大规模界面不稳定性的数值模拟．由于两种方法 均为显式格式，本文利用粗粒度数据并行，可达到很高的并行效率．本文还将两种并 界面的两种典型情况：重流体或轻流体中激波加速受正弦扰动的交界面，并清晰地得 到了明显的钉状结构和泡状结构，以及后期发展的蘑菇状交界面．三维结果同二维情 况演化趋势极为相似，区别在于定量的交界面扰动强度和扰动增长率的不同．此外。 本文还考虑了周向和轴向内聚爆炸波诱导的柱形R_M不稳定性． 本文还介绍了两种高精度的不可压缩流动界面模拟方法： VOF方法和Level—Set 方法，以及它们的相关改进，并将它们归结到一个统一的特征相函数框架下；通过不 可压动量方程粘性项的半隐处理，提高了计算方法的稳定性；同时借用了多重网格方 法，可以准确高效地处理变密度问题的压力Poisson方程，因而可以模拟高密度比和 高粘度比的实际问题．数值算例表明两种方法都能够正确地模拟交界面的演化以及拓 扑结构变化；Level—Set方法在处理表面张力问题的时候较VOF方法精度更高一些。 也更加准确一些；但是VOF方法的体积守恒保持得要比Level—Set方法好得多． 最后，本文将计算气体动力学中的混合有限体积／有限元方法推广到一般浅水流 动的数值模拟，通过基于水面高度的浅水方程来解决由于底坡和摩阻源项带来的通量 梯度和源项不平衡问题，使用简单的中心差源项离散得到满足精确守恒C性质的流场 解．同时本文将非结构网格自适应方法应用到浅水方程的模拟中来，并仔细讨论了一 种基于截断误差的误差指示器．数值算例表明，本非结构网格自适应混合有限体积／ 有限元方法是一种高效算法，能够很好地捕捉流场的细节，适用于复杂几何区域的一 般浅水自由面流动·y Abstract This thenumerical ofthe resolutioninterf靠 thesis high presents investigation methodsand Volume—Of-FluidLevel—Set cialsimulation methods， methods．especially andtheir in and flows．suchas incompressible applicationscompressible Richtmyer- inducedsurface soon．Fur- Meshkov(R-M)instability，dropdynamicsby tension，and