Improvement of Z-Weighted Function Based on Fifth-Order Nonlinear Multi-Order Weighted Method for Shock Capturing of Hyperbolic Conservation Laws

Based on a 5-point stencil and three 3-point stencils, a nonlinear multi-order weighted method adaptive to 5-3-3-3 stencils for shock capturing is presented in this paper. The form of the weighting function is the same as JS (Jiang–Shu) weighting; however, the smoothness indicator of the 5-point stencil adopts a special design with a higher-order leading term similar to the τ in Z weighting. The design maintains that the nonlinear weights satisfy sufficient conditions for the scheme to avoid degradation even near extreme points. By adjusting the linear weights to a specific value and using the τ in Z weighting, the method can be degraded to Z weighting. Analysis of linear weights shows that they do not affect the accuracy in the smooth region, and they can also adjust the resolution and discontinuity-capturing capability. Numerical tests of different hyperbolic conservation laws are conducted to test the performance of the newly designed nonlinear weights based on the weighted compact nonlinear scheme. The numerical results show that there are no obvious oscillations near the discontinuity, and the resolution of both the discontinuity and smooth regions is better than that of Z weights.


Introduction
Low-order schemes have been widely adopted to solve hyperbolic equations to avoid non-physical oscillations near discontinuities [1].However, the relatively large dispersion and dissipation errors of low-order schemes will smear out a lot of flow details in smooth regions [2], which has promoted the development of high-order nonlinear schemes.Representative achievements of high-order nonlinear schemes include ENO (essentially non-oscillatory) [3], WENO (weighted ENO) [4], TENO (targeted ENO) [5], CNS (compact nonlinear scheme) [6], WCNS (weighted CNS) [7], TCNS (targeted CNS) [8], etc. Nonlinear weighting has played an important role in developing some of these schemes, which provide largely improved scheme accuracy and computational efficiency.
The research on nonlinear weighted schemes has lasted for nearly 30 years.Many kinds of nonlinear weighting functions based on the WENO scheme have been proposed for better accuracy in the smooth region and better non-oscillatory properties near discontinuities.They mainly focus on improving the three main components of a nonlinear weighting scheme, which are the candidate interpolation choice, smoothness indicator, and nonlinear weighting function.Research focusing on improving candidate interpolation choices includes the bandwidth optimization of Martin et al. [9] and the adaptive centralupwind weighted scheme of Hu et al. [10].Meanwhile, Ha [11] modified the calculation of smoothness indicators, the special contribution of which lies in the parameter that adjusts the contribution of the first-order difference and the second-order difference to the smoothness indicators.Henrick [12] improved the relationship between the smoothness indicator and the nonlinear weight function by introducing an extra mapping function.Borges [13] reconstructed the smoothness indicators of each sub-stencil by introducing a higher-order smoothness indicator, thereby altering the expression form of the weight function.This approach has subsequently been extensively studied and applied by numerous researchers.The review papers [14][15][16] give a relatively comprehensive overview of the development of the WENO schemes.
Another branch of nonlinear weighted schemes is the weighted compact nonlinear scheme (WCNS) [7,[17][18][19], which is the scheme the current research is based on.It was summarized to have a series of advantages over the WENO schemes [2].A key property that distinguishes WCNS from WENO is that it has some advantages in the simulation of complex geometries [20] because of satisfying the geometric conservation law [21][22][23][24].Note that WCNS and WENO share similar nonlinear weighting processes; therefore, nonlinear weights for one of them can be easily generalized to another.
The improvements in the weighting functions share a common goal, which is to maintain accuracy in smooth regions and capture discontinuities without introducing obvious oscillations.After the proposal of the original fifth-order WENO-JS [25] scheme, an issue of order reduction near extreme points was observed.As a result, several researchers have made a series of modifications [9][10][11][12][13]26] to further improve the performance of the weighting functions.Among these improved weighting functions, Zhu et al. [27] proposed the WENO-ZQ scheme by adopting sub-stencils of different sizes, which is different from the WENO-JS scheme using same-size sub-stencils.The sub-stencils of WENO-ZQ consist of one 5-point sub-stencil along with two upwind and central 2-point sub-stencils, which is represented by the notation (5-2-2) for simplicity.The WENO-ZQ scheme can achieve optimal order even at extreme points, but it still exhibits non-physical oscillations near strong discontinuities.Subsequently, the WENO-MR (5-3-1) scheme [28,29] was proposed based on the combination of (5-3-1) sub-stencils, which reduced non-physical oscillations near discontinuities and enhanced robustness.However, it is important to note that the introduction of the first-order sub-stencil may potentially result in decreased resolution.Zhang [30] applied a similar strategy of [28] to WCNS, resulting in the WCNS-MR (5-3-1) scheme, which exhibits better resolution in smooth regions compared to the WENO-MR (5-3-1) scheme.However, there are still significant non-physical oscillations near strong discontinuities.Wang [31] further improved the efficiency of the WCNS-MR scheme.In the above research, except for the 5-point sub-stencil, the interpolation accuracy of the remaining sub-stencils is very low.To ensure accuracy, not only will the design complexity of the nonlinear weight function increase, but the linear weight of the 5-point stencil in actual examples also needs to be close to 1.This linear weight deteriorates the ability to capture discontinuities.Since its proposal, the Z function has been widely used and provides good stability and resolution near discontinuities.However, as the complexity of the problem increases, it is necessary to obtain more flow details in smooth regions, and the resolution of Z weights can no longer meet practical needs.
The goal of this article is to propose a new nonlinear weighted format based on the WCNS-JS and 5-3-3-3 stencils.Its characteristic feature is that the order of smoothness indicator on the stencil is related to the highest interpolation accuracy of the stencil, which maintains that nonlinear weights meet sufficient conditions to ensure that the scheme does not degrade in accuracy near extreme points.The smoothness indicator of the 5-point sub-stencil is obtained by combining the smoothness indicators of two 3-point sub-stencils and with an order higher than that of the 3-point sub-stencil.In addition, the ratio of linear weights on the three 3-point stencils meets the requirements of convex combination to fifth-order accuracy.A theoretical derivation and numerical tests are conducted to compare the resolution of the weighted method with Z nonlinear weights.These tests demonstrate the new weighted method's ability to achieve the desired accuracy in the smooth region and near the first-and second-order extreme points as well as its capability to accurately capture discontinuities.
The organization of the paper is as follows: In Section 2, we briefly introduce the fifth-order WCNS scheme.In Section 3, we establish a multi-order nonlinear weighted interpolation method based on the 5-3-3-3 stencils.In Section 4, some traditional numerical examples are used to demonstrate the characteristics of the multi-order weighting strategy in terms of smooth region resolution and discontinuous capture.In Section 5, we provide the conclusions obtained in this article.

WCNS
We focus on hyperbolic conservation laws of the form ∂u ∂t The WCNS procedure consists of three components: a high-order flux difference scheme, numerical flux evaluation, and weighted nonlinear interpolation.WCNS-E5 [32] is adopted for spatial discretizations, in which explicit flux differencing is used to calculate the ∂ f /∂x term in Equation ( 1), which is where a 0 = 75/64, a 1 = −25/384, a 2 = 3/640, fj+l+1/2 is the numerical flux obtained at cell-edge, and f ′ j represents a numerical approximation to the spatial derivative of the numerical flux f at the j cell-node.Equation ( 2) is sixth-order in accuracy with the parameters adopted.
In the WCNS, the flux at the half nodes is calculated using the flow variables û+ j+l+1/2 and û− j+l+1/2 , which are interpolated based on variables on the cell nodes.Numerical dissipation and the capture of discontinuities can be achieved through the use of either linear or nonlinear upwind interpolation and a suitable choice of the numerical flux.The interpolations of û+ j+l+1/2 and û− j+l+1/2 are symmetric to each other.For clarity, the superscripts of + and − are dropped, and only the formulas for û+ j+l+1/2 are presented.To ensure fifth-order accuracy of the overall scheme in smooth regions, the interpolation of ûj+l+1/2 also requires fifth-order accuracy.The WCNS-E5 scheme designs a nonlinear interpolation method with an ideal accuracy of fifth-order on the stencil x j+l−2 , x j+l−1 , x j+l , x j+l+1 , x j+l+2 .The specific implementation of the interpolation process is shown in Section 3.

Nonlinear Multi-Order Weighted Interpolation
The nonlinear weighted interpolation procedure in WCNS is very similar to the weighted reconstruction in popular schemes such as WENO, TENO, etc.These weighted nonlinear procedures commonly incorporate smoothness indicators on sub-stencils to effectively detect and capture discontinuities.In theory, the research results obtained can be generalized to other nonlinear weighting schemes as well.
After first introducing the general nonlinear weighted interpolation procedure, the classical JS nonlinear interpolation and the Z nonlinear interpolation will be presented in Sections 3.1 and 3.2, respectively.After that, the newly proposed weighting method is presented in Section 3.3, followed by its spectral analysis in Section 3.4.
The Lagrangian interpolation expression for WCNS-E5 based on 5-point stencil S = x j−2 , x j−1 , x j , x j+1 , x j+2 is The above interpolation is decomposed into a combination of three third-order interpolations based on 3-point sub-stencils S 1 = x j−2 , x j−1 , x j , S 2 = x j−1 , x j , x j+1 , and S 3 = x j , x j+1 , x j+2 .The specific form of the combination and the third-order interpolations are where the ideal weights are (5)

JS Nonlinear Interpolation
The famous JS nonlinear weights from Jiang and Shu [25] are defined as where the constant ε in the above equation is a small quantity to avoid the denominator becoming 0, and β k is the smoothness indicator.Usually, q takes a value of 2. The term β JS k [14] is defined as Previous studies have shown that such nonlinear weights can lead to a decrease in interpolation accuracy at the extreme points in smooth regions, resulting in a decrease in the overall scheme accuracy.To address this issue, numerous scholars have conducted extensive research, proposing the mapped WENO (M) [12], Z nonlinear weights [13], and ZQ nonlinear weights [27], etc.

Z Nonlinear Interpolation
Z nonlinear weights proposes [13] a high-order smoothing indicator and uses it to provide a construction method of smoothing indicators that differs from JS, theoretically solving the problem of extreme point reduction.Although Z weights are proposed based on the WENO format, this idea can be ported to nonlinear interpolation methods in the WCNS scheme.A brief introduction is as follows: The calculation method for the smoothness indicator of each stencil is Equation (7), and the high-order smoothness indicator τ z is defined as follows: The redefined smoothness indicator for each stencil is represented as The famous Z nonlinear weights from Jiang and Shu [25] are defined as The order of the nonlinear weighting function at non-extreme points is This can naturally meet the requirement of fifth-order accuracy for q = 1, and the order of the nonlinear weighting function at the extreme point is when q ≥ 2 can meet the accuracy requirements.

A New Fifth-Order Multi-Order Z Nonlinear Interpolation
Although Z weighting theoretically overcomes the problem of JS weighting decreasing at extreme points, in practical example testing, the nonlinear weight values of the substencil often deviate from the linear weights on smooth regions, leading to smooth regions not achieving fifth-order accuracy.Inspired by the Z weighting high-order smoothing indicator, this article designs a nonlinear weighting method for the 5-3-3-3 stencils.The weighting method is completely consistent with JS, and the interpolation expression for the 5-point stencil adopts Equation (3).The smoothing indicators of each sub-stencil are The nonlinear weighting is called MOZ (multi-order-Z) candidate interpolation weighting since the order of the smoothness indicators for candidate interpolations of different orders is also different.Their orders of magnitude are The term β MOZ 0 depends on all the nodes in the whole stencil.Near discontinuities, 3 is a small quantity related to the grid scale.This design improves the order of the smoothness indicator of the 5-point sub-stencil, which increases the nonlinear weight of the 5-point sub-stencil in smooth regions to improve the scheme's resolution and to avoid loss of extreme point accuracy.
The nonlinear weights of MOZ are specifically expressed as where , and γ 0 , γ k are the linear weights, which satisfy Based on order analysis using Taylor series expansions, the leading terms of the nonlinear weights are The final nonlinear interpolation expression at the half node is in which u k j+1/2 is given by Equation (4).To ensure that the denominator is not zero, all subsequent examples in this article use ε = 10 −20 .Through accuracy analysis, it is easy to obtain that MOZ weights avoid degrading the accuracy at both first-and second-order extreme points.To improve the resolution near discontinuities, the linear weights of the three 3-point stencils are set to γ k = d k (1 − γ 0 ), k = 1, 2, 3, which can achieve fifth-order accuracy.In addition, it should be pointed out that based on the design concept of MOZ, it is easy to prove that if β MOZ 0 = τ z , p = 2, q = 2 and the linear weights satisfy MOZ can degenerate into Z weighting; When the linear weights satisfy Equation ( 19) and p = 1, q = 2, it degenerates into JS weighting.
To further illustrate the influence of the linear weight γ 0 on the nonlinear weight values of the 5-point sub-stencil in the MOZ weighting function, the variation of the nonlinear weight values at the first-order extreme points and discontinuities is given based on the following equation: From Figure 1, it can be seen that although γ 0 is much smaller than the grid scale, the nonlinear weight ω MOZ 0 quickly approaches 1.0.In the vicinity of discontinuities, as γ 0 approaches 1.0, ω MOZ 0 is also much smaller than the grid scale.It can be seen that the MOZ designed in this paper achieves fifth-order accuracy and is relatively lenient with regard to the selection of linear weights.

Spectral and Accuracy Analysis
In this section, a simple analysis is conducted for comparing MOZ weighting with different linear weights and Z weighting in terms of accuracy and resolution.Using the ADR (approximate dispersion relation) [33] method, the spectral properties of fifth-order WCNS using ideal weights, Z weights, and MOZ weights with different linear weights are presented.Figure 2 shows the modified wavenumber between the ideal linear scheme and nonlinear weighted schemes.Figure 2a shows the real part, and Figure 2b shows the imaginary part.MOZ weighting with different linear weights has better dispersion and dissipation spectral characteristics than Z weighting.The convergence accuracy of Z weighting and MOZ weighting with different linear weights are tested at extreme points using the function u = e 0.75(x−1) x p [34].The point x = 0 is a first-order extreme point for p = 2, while a second extreme point is p = 3. Grid spacing of h = 2 −m h 0 (h 0 = 0.04) is adopted.The accuracy test results are presented in Tables 1 and 2. The MOZ weights successfully achieve the ideal order at both the first-and second-order extreme points.The Z function still degrades in order at the second-order extreme point.The actual errors of MOZ weighting are better than those of Z weighting, which are the same up to the third significant digits in Tables 1 and 2.
In addition, the MOZ weights constructed in this article have the following two characteristics: (1) Based on the original three 3-point stencils, a 5-point stencil is added and matched with a high-resolution high-order smoothness indicator.The choice of the smoothness indicator is consistent with the sub-stencil width.This choice helps the overall nonlinear interpolation to achieve the ideal order of accuracy.This choice ensures that despite the use of the JS weighting method, the format accuracy can still be maintained.In addition, by changing the linear weights, MOZ weighting can degrade to Z weighting and JS weighting, resulting in high design flexibility.(2) The linear weights of three 3-point stencils are designed in an ideal allocation ratio of a convex combination to fifth-order accuracy, which increases the utilization of the sub-stencil information and improves the resolution near discontinuities.
Table 1.Errors and accuracy of nonlinear weighted interpolation method at first-order extreme point x = 0 for the case u = e 0.75(x−1) x 2 for grids of h = 2 −m h 0 (h 0 = 0.04).Errors and accuracy of nonlinear weighted interpolation method at second-order extreme point x = 0 for the case u = e 0.75(x−1) x 3 for grids of h = 2 −m h 0 (h 0 = 0.04).

Numerical Tests
The 4-stage fourth-order Runge-Kutta method in reference [35] is used for temporal discretization.To maintain the same temporal truncation error as the fifth-order spatial discretization, the time-marching step calculation method is where (a, CFL, h) are the maximum eigenvalue, the Courant number, and the minimum grid scale, respectively.

Linear Scalar Equation
∂u ∂t Example 1: This example is calculated in x ∈ [0, 1] using periodic boundary conditions with initial conditions [30] of Based on the linear convection equation, the influence of MOZ weighting on the accuracy of the scheme is tested.Different from the accuracy tests, especially only for the extreme points in Section 3.4, the accuracy of the whole spatial and temporal discretization is tested here.The computational efficiency is compared, as shown in Table 3.The CPU time for MOZ is very close to that of Z. Firstly, we test the impact of the linear weights γ 0 = {h min , 0.01, 0.5, 0.8} in MOZ weighting on the calculation results.Taking u(x, 0) = e −300(x−0.5) 2 as the initial value, the solution is advanced to the non-dimensional time of 1 using a grid spacing of h = 0.1 × 2 −m for the m-th level of the grid.Table 4 shows the errors and convergence orders of the results using MOZ weighting.L 1 is the average of the sum of errors of all grid points relative to the analytical solution.L 2 is the average of the sum of squares of errors of all grid points relative to the analytical solution.L 3 is the maximum value among the absolute values of errors of all grid points relative to the analytical solution.The different choices of the linear weights γ 0 do not affect the convergence order of the schemes.Figure 3 shows the error curves concerning grid resolutions using the three nonlinear weights.The resulting error of MOZ weighting with different linear weights is better than that of Z weighting.Example 2: This example is calculated in x ∈ [−1, 1] using periodic boundary condi- tions with initial conditions [13] where z = −0.7,δ = 0.005, β = log 2/(36δ 2 ), a = 0.5, α = 10.The simulations are run to t = 8.From Figure 4, it can be seen that Z weighting and MOZ weighting with different linear weights have almost no oscillations near discontinuities.Meanwhile, MOZ weighting exhibits higher resolution near extreme points and discontinuities than Z weighting.

Burgers' Equation
The Burgers' equation is a nonlinear scalar equation with the form ∂u ∂t This example is calculated in x ∈ [0, 2] using periodic boundary conditions with initial conditions [12] u(x, 0) = sin(πx).( 26) Table 5 and Figure 5 show the errors and convergence orders using four types of weights at t = 1 when discontinuities of the solution have not formed.From the test results, it can be seen that the convergence orders using Z and MOZ weights can converge to the fifth order.The errors of MOZ weights are similar to those of Z weights.Discontinuity forms in the solution after t = 1.5/π.From Figure 6, it can be seen that MOZ weights with different linear weights and Z weights effectively capture the discontinuity.Next, the basic examples of one-dimensional and two-dimensional Euler equations are used to test the resolution and discontinuous capture ability of MOZ weighting and Z weighting.Based on the theoretical analysis and simple testing of scalar equations mentioned above, it is known that the linear weight of MOZ weighting will not affect the accuracy of the scheme, but the actual error will further decrease with the increase in linear weights γ 0 .The following examples in this article were tested using linear coefficients of γ 0 = {h, 0.01, 0.5, 0.8}, and the test results of γ 0 with the smallest results are still better than those of Z weighting.The maximum γ 0 can also achieve stable discontinuous capture.Recommended values for linear weights are h ⩽ γ 0 ⩽ 0.8.By adjusting γ 0 , the resolution and discontinuous capture ability of specific examples can be flexibly controlled and have a default value of γ 0 = 0.5.Due to limited space, we only provide test results for γ 0 = 0.5.

One-Dimensional Euler Equation
The Euler' equation is a nonlinear scalar equation with the form Example 1: The Sod and Lax problems [13] are both shock tube problems.The initial conditions of the Sod problem are while the initial conditions of the Lax problem are (ρ, u, p, γ) T = (0.445, 0.698, 3.528, 1.4) T , x ∈ [−0.5, 0), (0.5, 0, 0.571, 1.4 All the boundaries adopt the Dirichlet boundary condition using the left and right state values.
Figures 7 and 8 show the density results of nonlinear weighted functions at t = 0.25 on the Sod problem and t = 0.15 on the Lax problem, respectively.It can be seen that there are almost no non-physical oscillations near the discontinuity using Z weighting and MOZ weighting.The resolution of MOZ weighting near the discontinuity is higher than that of Z weighting.Example 2: The Shu-Osher problem [25] has initial conditions of (ρ, u, p, γ) T = (3.857,2.629, 10.333, 1.4 All the boundaries adopt the Dirichlet boundary condition.Figure 9 shows the density results of nonlinear weighted functions at t = 1.8 on the Shu-Osher problem.From the enlarged density, it can be seen that MOZ weighting has better resolution for high-frequency waves compared to Z weighting, and it is closer to the exact solution.
The upper and lower boundaries use Dirichlet boundary conditions according to Equation (36), while a symmetric boundary condition is adopted for the left and right boundaries.Using a Courant number of 0.3, the simulation is run until t = 1.95.Figure 12 shows that MOZ weighting has already captured the vortices near y = 0.5 on a 481 × 121 grid.However, these vortices have not yet fully formed when using Z weighting.With refinement of the grid, MOZ weighting can distinguish smaller-scale vortex structures near y = 0.6 compared to Z weighting.The results of MOZ weighting exhibit better resolution compared to Z weighting.Example 5: A Richtmyer-Meshkova instability problem [31] has initial conditions of (3.33333, 2.07063, 0.0, 7.125, 1.4), 0 ⩽ x < 0.06, (1.0, 0.0, 0.0, 1.0, 1.4), 0.06 ⩽ x < 0.1 + 0.008 • cos(nπy), (0.138, 0.0, 0.0, 1.0, 1.4), 0.1 + 0.008 • cos(nπy) ⩽ x ⩽ 0.6.(37) The calculation domain is [0, 0.6] × [0, 0.1].The upper and lower boundary conditions are periodic, and the left and right boundary conditions are fixed values that corresponds to the initial field.Using a Courant number of 0.3, the simulation is run until t = 0.16.Three scenarios are calculated: namely, 1, 2, and 3 shock bubbles, corresponding to n = 20, 40, and 60, respectively, in the initial conditions.From the results in Figure 13, it can be seen that the calculation results for all three cases show that MOZ weighting can capture more flow details about the "wake vortex" and "rod-shaped" positions of the shock bubble compared to Z-weighting, and its resolution is significantly better than that of Z weighting.Example 6: Shock vortex interaction problem [37].This example describes the interaction between a moving vortex and a steady shock wave.The calculation domain is [0, 2] × [0, 1].A shock wave with a Mach number of 1.1 is vertically located at x = 0.5 and has values of (ρ, u, v, p) = (1, 1.1 √ γ, 0, 1) on the left-hand side, and the values on the right-hand side are calculated according to the shock relationship.The left and right boundaries are fixed values, and the top and bottom are symmetric boundaries.The vortex is initially located at x = 0.5, and the vortex is described by ũ = ετ v e α(1−τ v 2 ) y, ṽ = −ετ v e α(1−τ v 2 ) x, T = − (γ−1)ε 2 e 2α ( 1−τv 2 ) 4αγ , S = 0, where r = (x − x c ) 2 + (y − y c ) 2 is the distance from the vortex center, r c = 0.05 is a parameter controlling the vortex size.ε = 0.3 controls the strength of the vortex, and α = 0.204 controls the attenuation rate of the vortices.The simulations are run until t = 0.8. Figure 14 shows that MOZ weighting can effectively simulate the shock.In addition, the results of Z weighting and MOZ weighting are compared at t = 0.6.Figure 15 shows the distribution of pressure at the vortex center y = 0.5.Compared with the finer-grid result of JS weights, it can be seen that MOZ weights have lower dissipation, and no obvious overshoot occurs post the shock wave.MOZ weighting has higher resolution before and after shock waves than Z weighting.zoomed in for the pre-shock position.The reference solution is calculated using JS weights and 1600 × 800 grid points.(The red square symbol denotes Z weighting, the blue upper triangle denotes MOZ weighting, the blue long dash line denotes MOZ weighting (1600 × 800 grid points), the black solid line denotes the reference solution.)

Conclusions
Based on 5-3-3-3 stencil interpolation, a new high-order smoothness indicator for the fifth-order sub-stencil is proposed; the resulting multi-order weighting can maintain the ideal order of the (5-3-3-3) nonlinear weighting.The linear weight selection of three 3-point stencils satisfies the convex combination condition for the optimal accuracy of fifth-order interpolation.These strategies form multi-order-Z nonlinear weighting, which we denote as MOZ weighting.The main conclusions are as follows: (1) Maintains accuracy: MOZ weighting ensures the ideal fifth-order accuracy of the scheme at the first-and second-order extremum points.(2) High aspect ratio: Numerical validation tests show that the format using MOZ weighting has better resolution than Z weighting in smooth regions, shear layers, and discontinuities, which is consistent with the conclusions of spectral analysis, and there are no obvious non-physical false fluctuations at discontinuities.From a theoretical analysis perspective, the number of calculations for MOZ is equivalent to that required by Z weighting, so increasing resolution does not decrease computational efficiency, resulting in a significant improvement in the ratio of resolution to efficiency.

Figure 2 .
Figure 2. Results of spectral characteristics of nonlinear weights using MOZ and Z: (a) real part; (b) imaginary part.(Red square symbol denotes the result based on Z weights; Yellow right triangle, blue left triangle, green top triangle, and purple circle have MOZ weights of γ 0 = {h, 0.01, 0.5, 0.8}, respectively; Solid lines indicate linear format).

Figure 3 .
Figure 3. Convergence orders of fifth-order schemes based on Z and MOZ weighting with different linear weights at t = 1 for the case of u(x, 0) = e −300(x−0.5) 2 based on the linear advection equation.(The horizontal axis represents the number of grid points (N) used in the calculation.The error of (a), (b), and (c) has already been provided in the previous text respectively.The circular symbol denotes the result of Z weighting; values of 1, 2, 3, and 4 denote the results of different linear weights of MOZ weighting).

Figure 4 .
Figure 4. Test of three types of nonlinear weights based on piecewise discontinuous solutions using 200 grid points at t = 8: (a) variable u; (b) variable u zoomed in.(Red square symbol denotes the result based on Z weighting; yellow right triangle, blue left triangle, green top triangle, and purple circle have MOZ weightings of γ 0 = {h, 0.01, 0.5, 0.8}, respectively; solid lines denotes the exact solution).

Figure 5 .
Figure 5. Convergence orders of four nonlinear fifth-order schemes at t = (0.5/π) of the onedimensional inviscid Burgers' equation with u(x, 0) = sin(πx), t = (0.5/π) (The horizontal axis represents the number of grid points (N) used in the calculation.The error of (a), (b), and (c) has already been provided in the previous text respectively.The circular symbol denotes the result of Z weighting; values of 1, 2, 3, and 4 denote the results of different linear weights of MOZ weighting).

Figure 6 .
Figure 6.The results of four nonlinear fifth-order schemes of the one-dimensional inviscid Burgers' equation with (x, 0) = sin(πx) at t = 1.5/π, grid N:100.(Red square symbol denotes the result based on Z weighting; yellow top triangle, blue right triangle, and green left triangle have MOZ weighting of γ 0 = {0.01, 0.5, 0.8}, respectively; solid lines denote the exact solution.)

Figure 7 .Figure 8 .
Figure 7. Numerical and exact solutions of the Sod problem at t = 0.25 using 200 grid points: (a) density; (b) density zoomed in (red square symbol denotes Z weighting, blue upper triangle denotes MOZ weighting, and solid line denotes the exact solution).

Figure 9 .
Figure 9. Numerical and exact solutions of the Shu-Osher problem at t = 1.8 using 200 grid points: (a) density; (b) density zoom in (red square symbol denotes Z weighting, blue upper triangle denotes MOZ weighting, and solid line denotes the exact solution).

Figure 14 .
Figure 14.Pressure contour lines of shock vortex interaction problem using 300 × 150 grid points.(a,b) Present the pressure distribution of Z weighting at t = 0.6 and t = 0.8, respectively; (c,d) present the pressure distribution of MOZ weighting at t = 0.6 and t = 0.8, respectively.There are 60 pressure contour lines ranging from 1.19∼13.7.CFL = 0.3.

Figure 15 .
Figure 15.Pressure distribution for the line y = 0.5 across the center of the vortex at t = 0.6 using 300 × 150 grid points.(a) Pressure distribution using different weights.(b)Pressure zoomed in for the post-shock position.(c) Pressure zoomed in near the center of the vortex.(d) Pressurezoomed in for the pre-shock position.The reference solution is calculated using JS weights and 1600 × 800 grid points.(The red square symbol denotes Z weighting, the blue upper triangle denotes MOZ weighting, the blue long dash line denotes MOZ weighting (1600 × 800 grid points), the black solid line denotes the reference solution.)

Table 3 .
A comparison of CPU times obtained by solving the scalar equation with different nonlinear weights (the grid N is 1600).

Table 6 .
Errors and convergence orders based on JS, Z, and MOZ weights for two-dimensional Euler vortex simulation, at t = 1.