高精度有限差分法模拟Kelvin-Helmholtz不稳定性  

Simulation of Kelvin-Helmholtz instability with high-order finite difference method

作　　者：王立锋[1,2] 叶文华[3,4,2] 范征锋[2] 李英骏[1] 

机构地区：[1]中国矿业大学力学与建筑工程学院,北京100083 [2]北京应用物理与计算数学研究所计算物理实验室,北京100088 [3]浙江大学物理系,杭州310028 [4]北京大学应用物理与技术研究中心,北京100871

出　　处：《强激光与粒子束》2009年第3期381-385,共5页High Power Laser and Particle Beams

基　　金：国家重点基础研究发展计划(2007CB815100);国家自然科学基金项目(10775020,10874242);高等学校博士学科点专项科研基金资助课题(20070290008)

摘　　要：WENO有限差分格式有较高的分辨精度,适合复杂流场的计算,在国际上被广泛采用。本文利用WENO有限差分格式求解2维守恒型欧拉方程,实现了对无粘流体中Kelvin-Helmholtz不稳定性的数值模拟。速度剪切方向采用周期边界条件;扰动增长方向采用嵌边出流边界条件,一个不稳定波长分布64个网格。数值模拟给出的扰动幅值线性增长率与线性稳定性分析给出的结果很好符合,显示了该格式的有效性和精度。数值模拟给出了清晰的密度等值线,表明该方法还具有较好的界面变形捕捉能力。The Kelvin-Helmholtz instability in the inviscid fluid was simulated by applying high-order weighted essentially non-oscillatory（WENO） scheme to solve the two-dimensional conservative Euler equations. The simulation adopted the periodic boundary condition for the velocity shear direction and the fringe method as the outflow boundary condition for the disturbance growth direction, where an unstable wavelength contained 64 grids. The linear growth rate of the disturbance amplitude derived from this numerical simulation agrees well with that from the linear stability analysis, which proves the validity and accuracy of the simulation method. The method can provide clear density contours, showing a good capturing ability of interface deformation.

关 键 词：KELVIN-HELMHOLTZ不稳定性 高精度有限差分法 WENO格式 欧拉方程 

分 类 号：O534.2[理学—等离子体物理] O539[理学—物理]