将考虑热传导和粘性情况下的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.