Numerical simulation of compressible flow with shocks using the OUCS2 upwind compact scheme

Abhishek Kundu

Authors

Abhishek Kundu Department of Mechanical Engineering, Swami Vivekananda Institute of Science and Technology, West Bengal, India

Keywords:

AUSM , Compressible flow, Euler equations, Flux limiter, High Mach number

Abstract

The OUCS2 upwind compact scheme for calculation of first derivatives in the Euler and NavierStokes equations is the focus of the present paper. Derived and analyzed by Sengupta et al. and primarily meant for incompressible flow problems, the scheme has recently been applied to compressible flow with shocks. To handle shocks effectively, a combination of second and fourth derivative artificial dissipation terms were used by Sengupta et al., instead of the inbuilt sixth derivative dissipation of the OUCS2. For unsteady flow problems with shocks and vortices, use of high order numerical dissipation is desirable for proper resolution of small scale structures. In a present paper we show that when used with flux limiters, one can use the inbuilt high order dissipation in smooth regions and still bring down the scheme to small stencil formulae with low order dissipation near shocks. The basic solver is built upon the AUSM+ algorithm. OUCS2 comes in while computing the left and right states of the primitive variables at the cell faces, requiring little effort to modify an existing AUSM+ based solver to a high resolution version with the help of OUCS2. We demonstrate the procedure for a number of problems in one and two dimensions solving the Euler equations.

References

Artebrant, R. and Schroll, H. J. (2006). Limiter-free third-order logarithmic reconstruction. SIAM J. Sci. Comput. 28(1):359–381.
De, S. and Murugan, T. (2011). Numerical simulation of shock tube generated vortex: effect of numerics. Int. J. Comput. Fluid Dyn. 25(6):345–354.
Dora, C. L., Murugan, T., De, S. and Das, D. (2014). Role of slipstream instability in formation of counterrotating vortex rings ahead of a compressible vortex ring. J. Fluid Mech. 753:29–48.
Hahn, M. and Drikakis, D. (2009). Assessment of largeeddy simulation of internal separated flow. J. Fluids Eng. 131(071201).
Halder, P., De, S., Sinhamahapatra, K. P. and Singh, N. (2013). Numerical simulation of shock-vortex interaction in schardin’s problem. Shock Waves. 23:495–504.
Jiang, G. S. and Shu, C. W. (1996). Efficient implementation of weighted eno schemes. J. Comput. Phys. 126:202–228.
Kawai, S., Shankar, S. K. and Lele, S. K. (2010). Assessment of localized artificial-difusivity scheme for compressible turbulent flows. J. Comput. Phys. 229:1739–1762.
Kundu, A. and De, S. (2014). High resolution euler simulation of a cylindrical blast wave in an enclosure. J. Vis. 18:733–738.
Liou, M. S. (1996). A sequel to the ausm: Ausm+. J. Comput. Phys. 129:364–382.
Meinke, M., Schroder, W., Krause, E. and Rister, T. (2002). A comparison of second- and sixth-order methods for large-eddy simulations. Comput. Fluids. 31:695–718.
Murugan, T., De, S., Dora, C. L. and Das, D. (2012). Numerical simulation and piv study of compressible vortex ring evolution. Shock Waves. 22:69–83.
Pirozzoli, S. (2002). Conservative-hybrid compact-weno schemes for shock-turbulence interaction. J. Comput. Phys. 178:81–117.
Ravichandran, K. S. (1994). Flux limited nonoscillatory cud-3 schemes for 1-d euler equations. Int. J. Comput. Fluid. Dyn. 3:141–152.
Ravichandran, K. S. (1997). Explicit third order compact upwind difference schemes for compressible flow calculations. Int. J. Comput. Fluid. Dyn. 8:311–316.
Rizzetta, D. P. and Morgan, P. E. (2008). A high-order compact finite-difference scheme for large eddy simulation of active flow control. Prog. Aerosp. Sci. 44:397–426.
Schulz-Rinne, C. W., Collins, J. P. and Glaz, H. M. (1993). Numerical solution of the riemann problem for two-dimensional gas dynamics. SIAM J. Sci. Comput. 14(6):1394–1414.
Sengupta, T. K., Bhole, A. and Sreejith, N. A. (2013). Direct numerical simulation of 2d transonic flows around airfoils. Comput. Fluids. 88:19–37.
Sengupta, T. K., Ganeriwal, G. and De, S. (2003). Analysis of central and upwind compact schemes. J. Comput. Phys. 192:677–694.
Sengupta, T. K., Ganerwal, G. and Dipankar, A. (2004). High accuracy compact schemes and gibbs’ phenomenon. J. Sci. Comput. 21(3):253–268.
Serna, S. (2006). A class of extended limiters applied to piecewise hyperbolic methods. SIAM J. Sci. Comput. 28(1):123–140.
Shu, C. W. and Osher, S. (1989). Efficient implementation of essentially non-oscillatory shock-capturing schemes, ii. J. Comput. Phys. 83:32–78.
Tolstykh, A. I. (1973). On a method for the numerical solution of navier-stokes equations of a compressible gas over a wide range of reynolds numbers. Dokl. Akad. Nauk. 210(1):48–51.
Tolstykh, A. I. (1994a). High-accuracy non-centered compact difference schemes for fluid dynamics applications. World Scientific., Singapore. Tolstykh, A. I. and Lipavskii, M. V. (1998). On performance of methods with third- and fifth-order compact upwind differencing. J. Comput. Phys. 140:205– 232.
Tolstykh, M. A. (1994b). Application of fifth-order compact upwind differencing to moisture transport in the atmosphere. J. Comput. Phys. 112:394–403.