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 obtained by the scheme of Eqs. (6)–(14) with penalty coefficients τ1=τ2=τ3=τ4=5. Here N×M=41×81, the time step Δt=0.1hx1hy1 is chosen to implement the fourth-order Runge–Kutta scheme. (a) At the moment t=4.2. (b) At the moment t=6.3.