R-M 不稳定性数值模拟方法NUMERICAL SIMULATION METHODS .PDF

CNIC-01683 IAPCM-0029 R-M 不稳定性数值模拟方法 NUMERICAL SIMULATION METHODS TO RICHTMYER-MESHKOV INSTABILITIES (In Chinese) 中 国 核 情 报 中 心 China Nuclear Information Centre 165 CNIC-01683 IAPCM-0029 R-M 不稳定性数值模拟方法 赵 宁 余 彦 唐维军 （北京应用物理与计算数学研究所，100088 ） 摘 要 通常，界面追踪方法都是由界面将计算区域分割为多个部分来 分别处理，同时用求解黎曼问题来追踪全部波系的发展。然而在多 维情况下，波系间相互作用的黎曼问题将非常复杂且目前尚处于分 析阶段，故在实际应用中只能追踪界面本身。提出了一种新的界面 追踪方法，该方法只追踪界面的发展而不再追踪波系，使计算量和 复杂性大大减少并保持了原界面追踪方法精确捕捉界面的特点，数 值结果令人满意。 关键词：界面追踪 黎曼解 高分辨 MmB 格式 强波探测开关 分类号：TL632 ; TL612.2 166 Numerical Simulation Methods to Richtmyer-Meshkov Instabilities (In Chinese) ZHAO Ning YU Yan TANG Weijun (Institute of Applied Physics and Computational Mathematics, Beijing, 100088) ABSTRACT Front tracking algorithms have generally assumed that the computational medium is divided into piece-wise smooth subdomains bounded by interfaces and that strong wave interactions are solved via Riemann solutions. However, in multi-dimensional cases, the Riemann solution of multiple shock wave interactions are far more complicated and still subject to analytical study. For this reason, it is very desirable to be able to track contact discontinuities only. A new numerical algorithm to couple a tracked contact surface and an untracked strong shock wave are described. The new tracking algorithm reduces the complication of computation, and maintains the sharp resolution of the contact surfac