Development of nonlinear weighted compact schemes with increasingly higher order accuracy

In this paper, we design a class of high order accurate nonlinear weighted compact schemes that are higher order extensions of the nonlinear weighted compact schemes proposed by Deng and Zhang X. Deng, H. Zhang, Developing high-order weighted compact nonlinear schemes, J. Comput. Phys. 165 (2000) 22-44] and the weighted essentially non-oscillatory schemes of Jiang and Shu G.-S. Jiang, C.-W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys. 126 (1996) 202-228] and Balsara and Shu D.S. Balsara, C.-W. Shu, Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy, J. Comput. Phys. 160 (2000) 405-452]. These nonlinear weighted compact schemes are proposed based on the cell-centered compact scheme of Lele S.K. Lele, Compact finite difference schemes with spectral-like resolution, J. Comput. Phys. 103 (1992) 16-42]. Instead of performing the nonlinear interpolation on the conservative variables as in Deng and Zhang (2000), we propose to directly interpolate the flux on its stencil. Using the Lax-Friedrichs flux splitting and characteristic-wise projection, the resulted interpolation formulae are similar to those of the regular WENO schemes. Hence, the detailed analysis and even many pieces of the code can be directly copied from those of the regular WENO schemes. Through systematic test and comparison with the regular WENO schemes, we observe that the nonlinear weighted compact schemes have the same ability to capture strong discontinuities, while the resolution of short waves is improved and numerical dissipation is reduced.

[1]  Y. Hattori,et al.  Sound generation by shock–vortex interactions , 1999, Journal of Fluid Mechanics.

[2]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .

[3]  Sergio Pirozzoli,et al.  Simulations and analysis of the coupling process of compressible vortex pairs: Free evolution and shock induced coupling , 2001 .

[4]  D. Gottlieb,et al.  Numerical analysis of spectral methods : theory and applications , 1977 .

[5]  Xiaolin Zhong,et al.  Derivation of high-order compact finite difference schemes for non-uniform grid using polynomial interpolation , 2005 .

[6]  Yong-Tao Zhang,et al.  Interaction of an oblique shock wave with a pair of parallel vortices: Shock dynamics and mechanism of sound generation , 2006 .

[7]  D. Gottlieb,et al.  The Stability of Numerical Boundary Treatments for Compact High-Order Finite-Difference Schemes , 1993 .

[8]  Chi-Wang Shu,et al.  A New Smoothness Indicator for the WENO Schemes and Its Effect on the Convergence to Steady State Solutions , 2007, J. Sci. Comput..

[9]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[10]  P. Lax.  Weak solutions of nonlinear hyperbolic equations and their numerical computation , 1954 .

[11]  P. Moin,et al.  A dynamic subgrid‐scale model for compressible turbulence and scalar transport , 1991 .

[12]  R. Hirsh,et al.  Higher order accurate difference solutions of fluid mechanics problems by a compact differencing technique , 1975 .

[13]  Joel H. Ferziger,et al.  A robust high-order compact method for large eddy simulation , 2003 .

[14]  Li Jiang,et al.  Weighted Compact Scheme for Shock Capturing , 2001 .

[15]  T. A. Zang,et al.  Spectral methods for fluid dynamics , 1987 .

[16]  D. Gottlieb,et al.  Numerical analysis of spectral methods , 1977 .

[17]  Sergio Pirozzoli,et al.  Development of optimized weighted‐ENO schemes for multiscale compressible flows , 2003 .

[18]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[19]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[20]  J. M. Powers,et al.  Mapped weighted essentially non-oscillatory schemes: Achieving optimal order near critical points , 2005 .

[21]  Xiaolin Zhong,et al.  High-Order Finite-Difference Schemes for Numerical Simulation of Hypersonic Boundary-Layer Transition , 1998 .

[22]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[23]  C. Tam,et al.  Dispersion-relation-preserving finite difference schemes for computational acoustics , 1993 .

[24]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

[25]  P. Moin,et al.  Turbulence statistics in fully developed channel flow at low Reynolds number , 1987, Journal of Fluid Mechanics.

[26]  Krishnan Mahesh,et al.  The influence of entropy fluctuations on the interaction of turbulence with a shock wave , 1997, Journal of Fluid Mechanics.

[27]  Hiroshi Maekawa,et al.  Compact High-Order Accurate Nonlinear Schemes , 1997 .

[28]  Sanjiva K. Lele,et al.  Large Eddy Simulation of a Near Sonic Turbulent Jet and Its Radiated Noise , 2001 .

[29]  Chi-Wang Shu,et al.  High Order Weighted Essentially Nonoscillatory Schemes for Convection Dominated Problems , 2009, SIAM Rev..

[30]  Michel A. Saad,et al.  Compressible Fluid Flow , 1985 .

[31]  Chi-Wang Shu,et al.  Numerical viscosity and resolution of high-order weighted essentially nonoscillatory schemes for compressible flows with high Reynolds numbers. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[32]  Li Jiang,et al.  Weighted Compact Scheme , 2001 .

[33]  Chi-Wang Shu,et al.  Multidomain WENO Finite Difference Method with Interpolation at Subdomain Interfaces , 2003, J. Sci. Comput..

[34]  Sergio Pirozzoli,et al.  Shock wave–thermal inhomogeneity interactions: Analysis and numerical simulations of sound generation , 2000 .

[35]  Sergio Pirozzoli,et al.  Shock-Wave–Vortex Interactions: Shock and Vortex Deformations, and Sound Production , 2000 .

[36]  B. J. Boersma,et al.  A staggered compact finite difference formulation for the compressible Navier-Stokes equations , 2005 .

[37]  Bernardo Cockburn,et al.  Nonlinearly stable compact schemes for shock calculations , 1994 .

[38]  Nikolaus A. Adams,et al.  A High-Resolution Hybrid Compact-ENO Scheme for Shock-Turbulence Interaction Problems , 1996 .

[39]  M Israeli,et al.  Numerical Simulation of Viscous Incompressible Flows , 1974 .

[40]  Xiaogang Deng,et al.  Developing high-order weighted compact nonlinear schemes , 2000 .

[41]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[42]  Zhi J. Wang,et al.  Optimized weighted essentially nonoscillatory schemes for linear waves with discontinuity: 381 , 2001 .

[43]  Krishnan Mahesh,et al.  High order finite difference schemes with good spectral resolution , 1997 .

[44]  Parviz Moin,et al.  Interaction of isotropic turbulence with shock waves: effect of shock strength , 1997, Journal of Fluid Mechanics.

[45]  Yong-Tao Zhang,et al.  Multistage interaction of a shock wave and a strong vortex , 2005 .

[46]  Guang-Shan Jiang,et al.  A High-Order WENO Finite Difference Scheme for the Equations of Ideal Magnetohydrodynamics , 1999 .

[47]  Sergio Pirozzoli,et al.  Conservative Hybrid Compact-WENO Schemes for Shock-Turbulence Interaction , 2002 .