采用多组分混合物质量分数模型和最小色散可控耗散格式的高分辨率有限体积方法,数值模拟了弱激波冲击不同角度的"V"形空气/SF_6界面的Richtmyer-Meshkov不稳定性问题.激波冲击界面后,在界面附近沉积涡量,形成沿界面规则排列的旋涡结构,同时界面扰动发展形成气泡和尖钉结构.本文统计了界面左端移动速度和界面混合宽度增长率等特征量的演化规律,并与已有的实验结果进行了对比,两者符合较好.讨论了物质界面处的流体向湍流混合发展的过程,随着界面旋涡结构的演化,涡结构之间开始发生相互诱导、并对等现象,并逐渐聚集在几个区域,而多尺度结构也因旋涡的诱导作用在这些区域中产生.通过对由雷诺数定义的惯性尺度进行分析,发现了具有上下边界的惯性尺度区域的形成,对动能能谱的分析发现了-5/3对数率的出现,这同样说明了惯性尺度区域的形成.由于湍流混合转捩与惯性尺度区域的形成是一致的,界面附近流场将发展为湍流.
Based on the mass fraction model of multicomponent mixture,the interactions between weak shock wave and"V" shaped air/SF6 interface with different vertex angles are numerical simulated.The numerical scheme used in the simulation is the high-resolution finite volume method with minimized dispersion and controllable dissipation scheme,in which the dissipation can be adjusted without affecting the already optimized dispersion property of the scheme.The grid sensitivity study is performed to guarantee that the resolution is sufficient in the numerical simulation.After the shock wave interacts with the interface,the baroclinic vorticity is deposited near the interface due to the misalignment of the density and pressure gradient,which is the manifestation of the Richtmyer-Meshkov instability,leading to the vortical structures forming along the interface.The interface perturbations lead to the bubbles and spikes appearing.The predicted leftmost interface displacement and interface width growth rate in the early stage of interface evolution agree well with the experimental results.The process of transition to turbulence at the material interface is studied in detail.The numerical results indicate that with the evolution of the interfacial vortical structure due to KelvinHelmholtz instability,the array of vortices begins to merge.As a result,the vortices accumulate in several distinct regions.It is in these regions that the multi-scale structures are generated because of the interaction between vortices.It is shown clearly that in the regions where vortices are accumulated,the fluctuation energy spectrum has many large and smallscale elements,which indicates there may be turbulent structures in these regions.To further examine if there is mixing transition in these regions,the characteristic length scales of the flow fields are calculated.The separation between the Lipemann-Taylor scale and inner viscous scale is observed based on the circulation-based Reynolds number,leading to the appearance of an uncoupled inertial