Table of Contents
Fetching ...

A Family of Even-Order Central-Upwind WENO Schemes with Averaged Downwind and Novel Global Smoothness Indicators

Jiaxi Gu, Bao-Shan Wang, Wai Sun Don, Jae-Hun Jung

Abstract

We propose a simple yet effective local smoothness indicator for the downwind stencil in central-upwind weighted essentially non-oscillatory (WENO) schemes of even order for hyperbolic conservation laws. Starting from an odd-order upwind WENO scheme, we construct an even-number-of-points stencil by incorporating a downwind substencil whose smoothness indicator is the arithmetic mean of all local smoothness indicators. This straightforward averaging approach incorporates regularity information from the entire stencil without requiring additional tuning parameters or complex formulations. Combined with affine-invariant Z-type nonlinear weights and a carefully designed global smoothness indicator, the resulting scheme, termed WENO-ZA6 for the sixth-order case, achieves optimal convergence rates at critical points up to second order, exhibits favorable dispersion and dissipation properties as confirmed by approximate dispersion relation analysis, and provides sharp, essentially non-oscillatory resolution of discontinuities. Numerical experiments on scalar problems and the one- and two-dimensional Euler equations demonstrate that WENO-ZA6 achieves accuracy comparable to or better than existing sixth-order central-upwind schemes (WENO-CU6, WENO-S6) and the seventh-order WENO-Z7, while requiring approximately 15\%--21\% less computational time. The framework extends naturally to fourth-, eighth-, and tenth-order schemes.

A Family of Even-Order Central-Upwind WENO Schemes with Averaged Downwind and Novel Global Smoothness Indicators

Abstract

We propose a simple yet effective local smoothness indicator for the downwind stencil in central-upwind weighted essentially non-oscillatory (WENO) schemes of even order for hyperbolic conservation laws. Starting from an odd-order upwind WENO scheme, we construct an even-number-of-points stencil by incorporating a downwind substencil whose smoothness indicator is the arithmetic mean of all local smoothness indicators. This straightforward averaging approach incorporates regularity information from the entire stencil without requiring additional tuning parameters or complex formulations. Combined with affine-invariant Z-type nonlinear weights and a carefully designed global smoothness indicator, the resulting scheme, termed WENO-ZA6 for the sixth-order case, achieves optimal convergence rates at critical points up to second order, exhibits favorable dispersion and dissipation properties as confirmed by approximate dispersion relation analysis, and provides sharp, essentially non-oscillatory resolution of discontinuities. Numerical experiments on scalar problems and the one- and two-dimensional Euler equations demonstrate that WENO-ZA6 achieves accuracy comparable to or better than existing sixth-order central-upwind schemes (WENO-CU6, WENO-S6) and the seventh-order WENO-Z7, while requiring approximately 15\%--21\% less computational time. The framework extends naturally to fourth-, eighth-, and tenth-order schemes.

Paper Structure

This paper contains 17 sections, 70 equations, 15 figures, 7 tables.

Figures (15)

  • Figure 1: The stencil $S^8$ together with its five $4$-point substencils $\{S_0,S_1,S_2,S_3\}$ (upwind) and $S_4$ (downwind) used in the reconstruction of $\hat{f}^+_{i+\frac{1}{2}}$ for the eighth-order central-upwind WENO scheme. The corresponding configuration for $\hat{f}^-_{i+\frac{1}{2}}$ is obtained by mirror symmetry.
  • Figure 2: The stencil $S^6$ together with its four $3$-point substencils $\{S_0,S_1,S_2\}$ (upwind) and $S_3$ (downwind) used in the reconstruction of $\hat{f}^+_{i+\frac{1}{2}}$ for the sixth-order central-upwind WENO scheme. The corresponding configuration for $\hat{f}^-_{i+\frac{1}{2}}$ is obtained by mirror symmetry.
  • Figure 3: ADR comparison: $\text{Re}(\Phi(\varphi))$ (phase) versus $\varphi$ (Left) and $\text{Im}(\Phi(\varphi))$ versus $\varphi$ (Right). Curves: WENO-CU6 (purple), WENO-S6 (green), WENO-ZA6 (red), WENO-Z7 (blue), CT6 (grey solid), spectral reference (black dashed).
  • Figure 4: ADR comparison: $\text{Re}(\Phi(\varphi))$ (phase) versus $\varphi$ (Left) and $\text{Im}(\Phi(\varphi))$ versus $\varphi$ (Right). Curves: WENO-ZA6 with $\tau$\ref{['eq:GSI']} (red solid circle), WENO-ZA6 with $\tau'$\ref{['eq:GSI_alt']} (azure hollow square), CT6 (grey solid), spectral reference (black dashed).
  • Figure 5: Example \ref{['ex:multiwave']}, (Left) numerical solution, (Middle) close-up view of the square wave, and (Right) log-scale absolute pointwise error of multi-waves advection problem computed by the WENO-CU6 (purple), WENO-S6 (green), WENO-ZA6 (red), and WENO-Z7 (blue) schemes with $N = 400$ at $T = 20$.
  • ...and 10 more figures

Theorems & Definitions (11)

  • Remark 1
  • Example 1
  • Example 2: Scalar conservation laws with convex and nonconvex fluxes
  • Example 3
  • Example 4
  • Example 5
  • Example 6
  • Example 7
  • Example 8
  • Example 9
  • ...and 1 more