Minimizing errors from linear and nonlinear weights of WENO scheme for broadband applications with shock waves

Abstract Improvements in the numerical algorithm for the dynamics of flows that involve discontinuities and broadband fluctuations simultaneously are proposed. These two flow features suggest numerical strategies of a paradoxical nature because the discontinuities demand dissipation, and the small-scale smooth features require the opposite. There may be several ways to approach such a complicated issue, but the natural choice is a numerical technique that can adjust adaptively with flow regimes. The weighted essentially non-oscillatory (WENO) scheme may be this choice. However, there are two sources of dissipation associated with the WENO procedure: the upwind optimal stencil and the nonlinear adaption mechanism. The current work suggests a robust and comprehensive treatment for the minimization of dissipation error from these two sources. The optimization technique, which is guided by restriction on the linear optimal weights derived from stability and consistency requirement, is used to delay the dissipation of the upwind optimal stencil to those wavenumbers for which the dispersion error is large. The parallel advantage of this technique is the improvement of the dispersion property. Nevertheless, optimization decreases the formal order of accuracy of the optimal stencil from fifth order to third order. This loss of accuracy is derived by Taylor series expansion. Using Taylor-series expansion and WENO procedure, the third-order accuracy is verified in the smooth region, except at the critical point of order two, where the order of accuracy reduces to at least second order. This possible loss of accuracy at the second-order critical point is restored in an attempt to reduce the dissipation induced by the nonlinear adaptive weights. Modification of the nonlinear weights to reduce the dissipation is introduced by redefining them with an additional smoothness indicator. Other suggestions to minimize the dissipation of the nonlinear weights are also reviewed. The numerical approximation of the spatial derivative is performed by means of a conservative and consistent finite difference method based on monotone local Lax–Friedrichs Riemann solver. The resulting scheme is then integrated by the optimal third-order TVD Runge–Kutta method to ensure the nonlinear stability of the overall numerical method. A variety of benchmark problems, ranging from non-broadband to broadband, are solved using the proposed schemes and compared with the existing ones. Most test problems are validated against exact or reference data. The numerical results with bandwidth optimization and modification of the nonlinear weights are consistently superior.

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

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

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

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

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

[6]  S. A. Orsag,et al.  Small-scale structure of the Taylor-Green vortex , 1984 .

[7]  San-Yih Lin,et al.  Parametric Study of Weighted Essentially Nonoscillatory Schemes for Computational Aeroacoustics , 2001 .

[8]  H. Lomax,et al.  Computation of shock wave reflection by circular cylinders , 1987 .

[9]  P. Lax,et al.  Systems of conservation laws , 1960 .

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

[11]  David P. Lockard,et al.  High accuracy algorithms for computational aeroacoustics , 1994 .

[12]  G. Taylor,et al.  Mechanism of the production of small eddies from large ones , 1937 .

[13]  Yuxin Ren,et al.  A characteristic-wise hybrid compact-WENO scheme for solving hyperbolic conservation laws , 2003 .

[14]  Parviz Moin,et al.  Assessment of high-resolution methods for numerical simulations of compressible turbulence with shock waves , 2010, J. Comput. Phys..

[15]  V. Gregory Weirs,et al.  Optimization of weighted ENO schemes for DNS of compressible turbulence , 1997 .

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

[17]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[18]  Yong-Tao Zhang,et al.  Resolution of high order WENO schemes for complicated flow structures , 2003 .

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

[20]  Wai-Sun Don,et al.  An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws , 2008, J. Comput. Phys..

[21]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[22]  Chi-Wang Shu,et al.  Anti-diffusive flux corrections for high order finite difference WENO schemes , 2005 .

[23]  R. Rogallo.  Numerical experiments in homogeneous turbulence , 1981 .

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

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

[26]  D. Pullin,et al.  Hybrid tuned center-difference-WENO method for large eddy simulations in the presence of strong shocks , 2004 .

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

[28]  Jang-Hyuk Kwon,et al.  A high-order accurate hybrid scheme using a central flux scheme and a WENO scheme for compressible flowfield analysis , 2005 .