激波驱动下柱壳混合的高分辨数值模拟的研究.pdfVIP

激波驱动下柱壳混合的高分辨数值模拟研究 *1) * * 王 丽 丽 ， 张 树 道 ， 杭 义 洪 * （北京应用物理与计算数学研究所，北京 100094 ） 摘要：本文采用高阶WENO 方法结合具有TVD 性质的高阶Runge-Kutta 方法，对汇聚激波驱动下柱壳结构分层 气体混合过程进行了一系列数值模拟，研究了初始扰动尺度对混合过程的影响。数值模拟结果表明：对于本文计 算中采用的小尺度条带形谱分布初始扰动，混合区的发展对扰动尺度较为敏感。扰动尺度越小，混合区宽度增长 越慢，混合区内混合沿展向更为充分。 关键词：汇聚激波；界面不稳定性；混合；柱壳；扰动 引 言 Richtmyer-Meshkov (RM)、Rayleigh-Taylor (RT)和Kelvin-Helmholtz (KH)三种界面不稳定性在自 然界、天体物理和工程应用领域广泛存在，是诱发湍流混合的一个重要机制，对武器物理和惯性约束聚变 （ICF）中的有关过程具有重要影响。ICF靶丸一般为多层球壳结构，其内爆压缩过程同时包含RT、RM和 KH 三种界面不稳定性，易造成物质的混合。本文利用高阶WENO 方法处理对流项差分，采用中心差分方法 处理空间二阶导数项，结合高阶Runge-Kutta方法进行时间积分，对汇聚激波驱动下柱壳结构分层流体的 混合过程进行数值模拟研究。该内爆模型在减速压缩阶段存在RT 失稳机制；汇聚激波扫过界面时存在RM 失稳机制，此外激波的多次聚心反射以及在各层界面之间的来回反射，使整个内爆过程中发生多次RM不稳 定性增长；RT、RM 不稳定性增长的同时又伴随发生KH 不稳定性增长。以上三种不稳定性的发展使混合区 动力学状态十分复杂，同时收缩几何效应进一步加剧了不稳定性增长和混合的发展。目前对此类汇聚几何 中界面不稳定性和混合的发展规律还缺乏足够的认识，本文通过高分辨数值模拟重点研究了初始扰动尺度 对混合过程的影响。 1 控制方程 采用如下多组分Navier-Stokes方程组描述流体运动：  ui  0 t x i u (u u p ) d i  i j i j ij F t xj xj i (1) E (E p )u  T d u  j ( )  ji i Fu t x x x x i i j j j j  ( u   m  m j (D m ) , m 1, N 1 t x x x j j j 1）E-mail ： wanglili@iapcm.ac.cn RT 其中 表示第 m 种组分的质量份额。采用多方气体状态方程确定混合气体的压力：P ， m