Stabilized seventh-order dissipative compact scheme for two-dimensional Euler equations
Qin Jia-Xian, Chen Ya-Ming, Deng Xiao-Gang
       

Density contours of the vortex of Eqs. (15)–(18) on the grid consists of two subdomains with 41 × 21 and 21 × 21 points. Here the time step Δ t = 0.1 h x 1 h y 1 is chosen to implement the fourth-order Runge–Kutta scheme. (a) At the moment t = 0. (b) At the moment t = 10.