可压缩多介质粘性流体的数值计算  

Computation of compressible multi-viscosity-fluid flows

作　　者：柏劲松[1] 李平[1] 王涛[2] 谢彬[2] 钟敏[2] 陈森华[2] 

机构地区：[1]中国工程物理研究院流体物理研究所冲击波物理与爆轰物理实验室,四川绵阳621900 [2]中国工程物理研究院流体物理研究所,四川绵阳621900

出　　处：《爆炸与冲击》2007年第6期515-521,共7页Explosion and Shock Waves

基　　金：国家自然科学基金项目(10672151);冲击波物理与爆轰物理实验室基金项目(9140C6704010704)

摘　　要：将考虑热传导和粘性情况下的Navier Stokes方程描述的物理过程分解成3个子过程进行数值计算,即把整个流量计算分解成无粘性流量、粘性流量和热流量3部分,采用多介质流体高精度parabolic piece-wise method(PPM)方法、二阶空间中心差方法和两步Rung-Kutta时间推进方法相结合进行数值计算。给出了激波管中Riemann问题和二维、三维Richtmyer-Meshkov界面不稳定性的Navier Stokes方程和Euler方程对比计算结果,显示了粘性对界面不稳定性的影响。By using a splitting technique, the flux of the Navier Stokes governing equations including heat exchange and viscosity are divided into three parts, so called inviscid flux, viscous flux and heat flux, to calculate. The inviscid part of flux is evaluated using a high resolution multi-fluid parabolic piecewise method, the viscous part of flux is computed with a second order central difference in space and two-step Rung-Kutta scheme in time, and the heat flux part is not considered. Influences of fluid viscosity on instability of the interface between two fluids are revealed by computation results of the examples of Riemann problem in shock tube, two and three dimensional Richtmyer-Meshkov instabilities using the Navier Stokes governing equations and Euler governing equations.

关 键 词：流体力学 NAVIER STOKES方程 多流体PPM方法 粘性流体 

分 类 号：O357[理学—流体力学]