Next Article in Journal
Design of Polynomial Observer-Based Control of Fractional-Order Power Systems
Previous Article in Journal
Robust DDoS Attack Detection Using Piecewise Harris Hawks Optimizer with Deep Learning for a Secure Internet of Things Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved WENO-Z Scheme for Hyperbolic Conservation Laws with New Global Smoothness Indicator

1
School of Mathematics and Computational Science, Xiangtan University, Xiangtan 411105, China
2
Hunan Key Laboratory for Computation and Simulation in Science and Engineering, Xiangtan University, Xiangtan 411105, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(21), 4449; https://doi.org/10.3390/math11214449
Submission received: 15 September 2023 / Revised: 21 October 2023 / Accepted: 25 October 2023 / Published: 27 October 2023
(This article belongs to the Topic Computational Fluid Dynamics (CFD) and Its Applications)

Abstract

:
The fifth-order WENO-Z scheme proposed by Borges et al., using a linear combination of low-order smoothness indicators, is designed to provide a low numerical dissipation to solve hyperbolic conservation laws, while the power q in the framework of WENO-Z plays a key role in its performance. In this paper, a novel global smoothness indicator with fifth-order accuracy, which is based on several lower-order smoothness indicators on two-point sub-stencils, is presented, and a new lower-dissipation WENO-Z scheme (WENO-NZ) is developed. The spectral properties of the WENO-NZ scheme are studied through the ADR method and show that this new scheme can exhibit better spectral results than WENO-Z no matter what the power value is. Accuracy tests confirm that the accuracy of WENO-Z with q = 1 would degrade to the fourth order at first-order critical points, while WENO-NZ can recover the optimal fifth-order convergence. Furthermore, numerical experiments with one- and two-dimensional benchmark problems demonstrate that the proposed WENO-NZ scheme can efficiently decrease the numerical dissipation and has a higher resolution compared to the WENO-Z scheme.

1. Introduction

As an important class of differential equation, the hyperbolic conservation equation has a wide application in scientific and engineering fields. One of the characteristics of this equation is that the solution may generate discontinuities even if its initial condition is smooth. For capturing the discontinuities, many numerical schemes, such as total variation diminishing schemes (TVD) [1,2,3], essentially non-oscillatory schemes (ENO) [4,5,6], weighted essentially non-oscillatory schemes (WENO) [7,8,9,10], etc., have been developed in recent decades. Among them, the WENO schemes have attracted great attention because they can not only obtain high-order solutions in smooth regions but also efficiently avoid spurious oscillatory near discontinuities.
The WENO scheme, first proposed by Liu et al. [7], provided one-order higher solutions than the ones of the rth-order ENO scheme [4,5] by using a convex combination of r candidate sub-stencils. Subsequently, Jiang and Shu [8] devised a new smoothness indicator for the WENO scheme by utilizing the Lagrange form of interpolation polynomials and proposed an efficient WENO-JS scheme. For the case r = 3, the resulting WENO-JS scheme is about twice as fast as the fourth-order ENO scheme on computational efficiency. Most importantly, this scheme can achieve optimal fifth-order accuracy. Henrick et al. [11] derived a sufficient condition for the fifth-order convergence and conducted a further analysis of the accuracy of WENO-JS. They found that the fifth-order WENO-JS scheme was only third-order accurate at first-order critical points. To increase the accuracy of the WENO scheme, different types of finite difference WENO schemes by reconstructing flux functions have been proposed [12,13,14,15,16], such as the RBF-WENO scheme, which was based on radial basis functions (RBFs) and the nonlinear weights of WENO-JS, which was first proposed by Guo and Jung in [12] and further developed by Abedian et al. [14,15]. On the other hand, some scholars developed various techniques to calculate the weights to improve the accuracy. For example, Henrick et al. [11] introduced a corrected mapping function to the nonlinear weights of WENO-JS and presented a WENO-M scheme which could restore the optimal convergence order even at critical points but was computationally expensive. Unlike the idea of WENO-M, Borges et al. [9] designed a fifth-order global smoothness indicator using the linear combination of low-order smoothness indicators and proposed an innovative WENO-Z scheme with lower numerical dissipation. Meanwhile, this scheme was able to obtain almost the same accurate numerical solution as WENO-M and reduced computation by 25%.
Castro et al. [17] extended the WENO-Z scheme to arbitrary (2r − 1) th odd-order accuracy by designing a general formula for the global smoothness indicator. However, it should be noted that the numerical dissipation of WENO-Z is dependent upon the power parameter q. For example, the scheme is less dissipative with q = 1 but will lose accuracy at critical points, such as the fifth-order WENO-Z is only fourth-order accurate at first-order critical points. In contrast, the scheme will decrease the correction of nonlinear weights and is more dissipative with q = 2. Hence, Borges et al. [9] suggested the power q taken to be 1 to maintain the low dissipation property. Recently, some studies [18,19,20,21,22,23] reconstructed a series of higher-order (even up to eighth-order) reference smoothness indicators to recover the optimal convergence for the low-dissipation WENO-Z scheme (means q = 1). Although numerical results confirmed that these modified schemes with higher-order smoothness indicators can achieve the desired accuracy in smooth regions including critical points, most of these works either used larger stencils or introduced user-tunable parameters, which led to the increased computation time or limited the application.
On the other hand, Acker et al. [24] concluded the important role of increasing the contribution of less-smooth candidate sub-stencils for increasing the resolution. They added an extra term to WENO-Z weights to increase the relevance of less-smooth sub-stencils and proposed a fifth-order WENO-Z+ scheme. In terms of the lrest progress, a new fifth-order WENO-Z scheme was developed by Tang et al. [25] to perform lower dissipation at discontinuities by constructing a selector that can identify the less-smooth sub-stencils. Unfortunately, the presented WENO schemes in [24,25] were still fourth-order accurate at critical points. Therefore, how to construct a low-dissipation WENO-Z scheme with better performances, which has overall fifth-order accuracy including critical points and a higher resolution near discontinuities, should be investigated thanks to its significant computational benefit.
In this study, we divide the five-point global stencil of the WENO scheme into four smaller two-point sub-stencils. By utilizing the local smoothness indicators of these sub-stencils, we devise a novel fifth-order global smoothness indicator and develop a WENO-NZ scheme. The numerical results demonstrate that the WENO-NZ scheme is less dissipative and can restore the optimal convergence order whatever the value of power q. Additionally, this scheme shows improved capability in capturing discontinuities and small scales and exhibits higher resolution compared to the WENO-Z scheme.
The rest of the paper is structured as follows: In Section 2, a brief overview of the fifth-order WENO-JS and WENO-Z schemes is discussed. In Section 3, an improved WENO-Z scheme is developed by using a new fifth-order global smoothness indicator. The performance comparison between WENO-Z and the proposed scheme for the one- and two-dimensional Euler equations is given in Section 4. Finally, the conclusions are presented in Section 5.

2. Numerical Method

The one-dimensional scalar hyperbolic conservation laws can be written as
u t + f ( u ) x = 0 , x [ a , b ] , t > 0 ,
here, u ( x , t ) is conservative variable, f ( u ) is flux. Consider the uniform grids x j = j Δ x ( j = 0 , , N ) and define the cell interfaces as x j + 1 / 2 = x j + Δ x / 2 , Δ x is the uniform grid spacing, the semi-discretization of Equation (1) can be given as
d u j d t = ( f x ) x j ,
here, u j is a numerical approximation to u ( x j , t ) . The conservation property of Equation (2) is obtained by implicitly defining f ( u ) = 1 Δ x x Δ x / 2 x + Δ x / 2 h ( ξ ) d ξ , and have
d u j d t = 1 Δ x ( h j + 1 / 2 h j 1 / 2 ) ,
here, h j ± 1 / 2 = h ( x j ± 1 / 2 ) . Then, Equation (3) can be further approximated as
d u j d t 1 Δ x ( f ^ j + 1 / 2 f ^ j 1 / 2 ) ,
here, the numerical fluxes f ^ j ± 1 / 2 are higher order polynomial interpolations to h j ± 1 / 2 and can be obtained by the WENO reconstruction.
In general, the flux f ( u ) is always split into two parts, f + and f , to ensure numerical stability, such as
f ( u ) = f + ( u ) + f ( u ) ,
where d f + ( u ) d u 0 and d f ( u ) d u 0 . In this paper, the global Lax–Friedrichs splitting is used, there are
f ± ( u ) = 1 2 ( f ( u ) ± α u ) ,
where α = max u | f ( u ) | . Hence, f ^ j + 1 / 2 ± can be computed as the positive and negative parts of f ( u ) , respectively, and have
f ^ j + 1 / 2 = f ^ j + 1 / 2 + + f ^ j + 1 / 2 .
In fact, f ^ j + 1 / 2 + and f ^ j + 1 / 2 are symmetric with respect to x j + 1 / 2 for the WENO scheme, so we only describe how f ^ j + 1 / 2 + is approximated and the “+” will be dropped in the following content for convenience.
According to the construction of the fifth-order WENO scheme, the numerical flux f ^ j + 1 / 2 can be given as
f ^ j + 1 / 2 = k = 0 2 ω k f ^ k , j + 1 / 2 ,
where the numerical fluxes f ^ k , j + 1 / 2 ( k = 0 , 1 , 2 ) are given by
f ^ 0 , j + 1 / 2 = 1 3 f j 2 7 6 f j 1 + 11 6 f j , f ^ 1 , j + 1 / 2 = 1 6 f j 1 + 5 6 f j + 1 3 f j + 1 , f ^ 2 , j + 1 / 2 = 1 3 f j + 5 6 f j + 1 1 6 f j + 2 .
Now, there are many different WENO construction techniques to compute the nonlinear weights ω k , while the WENO-JS and WENO-Z are the most popular schemes.

2.1. WENO-JS Scheme

The well-known fifth-order WENO-JS scheme was proposed by Jiang and Shu [8], whose nonlinear weights were defined as
ω k J S = α k J S k = 0 2 α k J S , α k J S = d k ( β k + ε ) q ,
where d k are ideal weights that resulted in the fifth-order central upwind scheme, there are d 0 = 1 / 10 , d 1 = 6 / 10 , and d 2 = 3 / 10 respectively. The power q 1 is used to control the weights of less-smooth sub-stencils. The larger the q is set, the smaller the weight is obtained, but the scheme is more dissipative. Hence q = 2 is used for the WENO-JS scheme. ε is a positive sensitivity parameter used to avoid division by 0 and ε = 10 40 is adopted in this study. The local smoothness indicator β k is defined by
β k = l = 1 r 1 Δ x 2 l 1 x j 1 / 2 x j + 1 / 2 ( d l d x l f ^ k ( x ) ) 2 d x .
For the fifth-order WENO scheme (r = 3), the β k can be further expressed as
β 0 = 13 12 ( f j 2 2 f j 1 + f j ) 2 + 1 4 ( f j 2 4 f j 1 + 3 f j ) 2 , β 1 = 13 12 ( f j 1 2 f j + f j + 1 ) 2 + 1 4 ( f j + 1 f j 1 ) 2 , β 2 = 13 12 ( f j 2 f j + 1 + f j + 2 ) 2 + 1 4 ( 3 f j 4 f j + 1 + f j + 2 ) 2 .
The Taylor expansions of Equation (12) at x j can be given as
β 0 = f j 2 Δ x 2 + ( 13 12 f j 2 2 3 f j f j ) Δ x 4 ( 13 6 f j f j 1 2 f j f j ( 4 ) ) Δ x 5 + O ( Δ x 6 ) , β 1 = f j 2 Δ x 2 + ( 13 12 f j 2 + 1 3 f j f j ) Δ x 4 + O ( Δ x 6 ) , β 2 = f j 2 Δ x 2 + ( 13 12 f j 2 2 3 f j f j ) Δ x 4 + ( 13 6 f j f j 1 2 f j f j ( 4 ) ) Δ x 5 + O ( Δ x 6 ) .
Substituting Equation (13) into Equation (10), there are
{ ω k J S = d k + O ( Δ x 2 ) , f j 0 , ω k J S = d k + O ( Δ x ) , f j = 0 , f j 0 .
Henrick et al. [11] concluded a detailed analysis of the accuracy of the fifth-order WENO scheme and derived a sufficient condition for the fifth-order convergence as follows
ω k d k = O ( Δ x 3 ) .
Equation (15) is further weakened as ω k d k = O ( Δ x 2 ) by Don et al. in [26]. Although it is not a necessary condition, it can be used as a simple criterion for designing nonlinear weights. Obviously, the weights of WENO-JS cannot meet the sufficient condition of Equation (15), especially at first-order critical points. It was reported that the WENO-JS scheme was only third-order accurate at critical points [11].

2.2. WENO-Z Scheme

Borges et al. [9] developed another fifth-order WENO scheme, called WENO-Z. They defined the nonlinear weights as
ω k Z = α k Z k = 0 2 α k Z , α k Z = d k ( 1 + ( τ 5 β k + ε ) q ) , q = 1 , 2 .
Here, τ 5 is a global smoothness indicator given as
τ 5 = | β 0 β 2 | ,
which with the following Taylor expansion
τ 5 = | 13 3 f j f j f j f j ( 4 ) | Δ x 5 .
Substituting Equations (13) and (18) into Equation (16), then have
{ ω k Z = d k + O ( Δ x 3 q ) , f j 0 , ω k Z = d k + O ( Δ x q ) , f j = 0 , f j 0 .
It can be found that the weights of WENO-Z with q = 1 do not meet the sufficient condition (15) at first-order critical points. Furthermore, Borges et al., confirmed that the WENO-Z scheme with q = 2 achieved fifth-order accuracy even at critical points. However, an increase in q resulted in more numerical dissipation of WENO-Z. Therefore, the power parameter q = 1 is suggested by Borges et al. and used by other scholars [18,19,20,21,22,23,24,25].

3. The New WENO-Z Scheme

Some works [20,21,22,23] demonstrated that constructing more higher-order smoothness indicators was useful in recovering the desired order for the low-dissipation WENO-Z scheme (means q = 1). Unlike these studies, we reconstructed a novel fifth-order global smoothness indicator with lower dissipation for the WENO-Z scheme to overcome its inherent limitations, such as its failure to balance the low-dissipation property and fifth-order convergence.

3.1. Construction of Novel Global Smoothness Indicator

Within the framework for the original global smoothness indicator τ 5 of WENO-Z, it is found that a discontinuity would affect the interpolation on the range of 2 Δ x at least. Suppose a discontinuity is at any position to the left (or right) side of the center point x j , it would affect the interpolation on the three-point sub-stencil S 0 (or S 2 ). However, if a discontinuity is just right at x j , it would affect the interpolation of the entire five-point global stencil. To reduce the impact of the discontinuities on interpolation, we subdivided the global stencil into four smaller sub-stencils { S 00 , S 01 , S 12 , S 22 } , each one only involving two points, as shown in Figure 1. In fact, these sub-stencils are generally used to construct the third-order WENO-Z scheme (r = 2).
Based on these new smaller sub-stencils, a discontinuity will affect the nearby 2 Δ x range at most, that is, two two-point sub-stencils. Furthermore, only a two-point sub-stencil where the discontinuity is located may be affected. From Equation (11), the smoothness indicators on sub-stencils { S 00 , S 01 , S 12 , S 22 } can be computed as
β 00 = ( f i 2 f i 1 ) 2 , β 01 = ( f i 1 f i ) 2 , β 12 = ( f i f i + 1 ) 2 , β 22 = ( f i + 1 f i + 2 ) 2 ,
with the following Taylor expansions
β 00 = f j 2 Δ x 2 3 f j f j Δ x 3 + ( 9 4 f j 2 + 7 3 f j f j ) Δ x 4 ( 7 2 f j f j + 5 4 f j f j ( 4 ) ) Δ x 5 + O ( Δ x 6 ) , β 01 = f j 2 Δ x 2 f j f j Δ x 3 + ( 1 4 f j 2 + 1 3 f j f j ) Δ x 4 ( 1 6 f j f j + 1 12 f j f j ( 4 ) ) Δ x 5 + O ( Δ x 6 ) , β 12 = f j 2 Δ x 2 + f j f j Δ x 3 + ( 1 4 f j 2 + 1 3 f j f j ) Δ x 4 + ( 1 6 f j f j + 1 12 f j f j ( 4 ) ) Δ x 5 + O ( Δ x 6 ) , β 22 = f j 2 Δ x 2 + 3 f j f j Δ x 3 + ( 9 4 f j 2 + 7 3 f j f j ) Δ x 4 + ( 7 2 f j f j + 5 4 f j f j ( 4 ) ) Δ x 5 + O ( Δ x 6 ) .
Consider the WENO scheme will lose accuracy at first-order critical points, in order to minimize the influence of the first-order derivative, meanwhile ensuring the global smoothness indicator has high-order accuracy, the form of a new fifth-order global smoothness indicator can be given as
τ N 5 = | β 00 3 β 01 + 3 β 12 β 22 | / 6 ,
the Taylor expansion of τ N 5 is as follows
τ N 5 = | f j f j + 1 3 f j f j ( 4 ) | Δ x 5 .

3.1.1. Accuracy

Based on the nonlinear weights for the WENO-Z scheme, substituting τ N 5 instead of τ 5 into Equation (16) can obtain a new weight, as follows
ω k N Z = α k N Z k = 0 2 α k N Z , α k N Z = d k ( 1 + ( τ N 5 β k + ε ) q ) , q = 1 , 2 .
Substituting Equation (24) into Equation (8) results in a new WENO scheme, which is named WENO-NZ in this paper. Since τ N 5 has the same fifth-order accuracy as τ 5 , similar to the accuracy analysis of WENO-Z, there are
ω k N Z = { d k + O ( Δ x 3 q ) , f j 0 , d k + O ( Δ x q ) , f j = 0 , f j 0 ,
which indicates that the WENO-NZ scheme has the same theoretical convergence accuracy as WENO-Z. A detailed analysis of the accuracy properties can be found in [9]. For simplicity, the WENO-Z and WENO-NZ schemes with different powers (q = 1 or 2) are, respectively, denoted as WENO-Z1, WENO-Z2, WENO-NZ1, and WENO-NZ2 in this paper.

3.1.2. Weights of Less-Smooth Sub-Stencils

In [9], numerical results made evident that the WENO-Z places higher weights on less-smooth sub-stencils than WENO-JS, this is the reason why WENO-Z achieves less dissipative results. The results in [24] also showed that it is important to increase the weights of less-smooth sub-stencils. Therefore, here is a numerical example to analyze the weights of less-smooth sub-stencils for the new WENO-NZ scheme.
Consider the linear case of Equation (1) as follows
u ( x , t = 0 ) = f ( x ) = { sin ( π x ) 1 2 x 3 , 1 x < 0 , sin ( π x ) 1 2 x 3 + 1 , 0 x 1 ,
with the exact solution u ( x , t ) = u ( x t , 0 ) . This numerical example contains a jump discontinuity at x = 0, which is usually used to study the behavior of nonlinear weights in the smooth and discontinuous regions for the WENO scheme. We compute this numerical example to the final time t = 2.0 using N = 200 with periodic boundary, and the time step is set as Δ t = Δ x / 2 . As shown in Figure 2a, the numerical results of WENO-NZ1 are the closest to that of the exact solution, and WENO-Z2 generates the most dissipation near discontinuities compared to others. Figure 2b gives the values of smoothness indicators at the initial step t = 0, i.e., β k ( k = 0 , 1 , 2 ) , τ 5 , and τ N 5 . We can see that τ N 5 is smaller than τ 5 whether in smooth or discontinuous regions, hence, the WENO-NZ scheme is more accurate compared to WENO-Z.
Figure 3 further provides the distributions of nonlinear weights for WENO-Z1, WENO-Z2, WENO-NZ1, and WENO-NZ2. As Figure 3 shows, at the point x = 0.02 , where the discontinuity occurs for the first time, the ω 2 is smaller than ω 0 and ω 1 , that is, the nonlinear weight of the less-smooth sub-stencil is smaller than that of smooth ones. It can be seen that the WENO-NZ assigns relatively larger weights to the less-smooth sub-stencils compared to the corresponding WENO-Z scheme. It can also be seen that the nonlinear weights of the WENO-NZ scheme are nearer to the ideal ones, resulting in the WENO-NZ scheme being closer to the central upwind scheme. As a result, this new scheme is less dissipative than the WENO-Z scheme.

3.2. The Spectral Properties of New Scheme

The spectral properties of the WENO-NZ scheme are studied using the ADR method proposed by Pirozzoli in [27]. As shown in Figure 4, it can be found that the dispersion and dissipation results of the WENO-NZ scheme are closer to the fifth-order upwind scheme (UP5) compared to WENO-Z, and the dispersion and dissipation curves of WENO-NZ are very similar to those of UP5. In addition, the spectral results of WENO-Z1 are better than those of WENO-Z2, but the results of WENO-NZ2 are nearly equivalent to that of WENO-NZ1, which implies that compared to WENO-Z, the WENO-NZ is a lower-dissipation scheme regardless of power q.

4. Numerical Results

In this section, a series of typical numerical examples involving discontinuities and complex scale waves are computed to examine the accuracy, resolution, and low-dissipation properties of the proposed WENO-NZ scheme. The third-order Runge–Kutta method [1] is applied for time discretization with a CFL number of 0.5, as
u ( 1 ) = u n + Δ t L [ u n ] , u ( 2 ) = 3 4 u n + 1 4 ( u ( 1 ) + Δ t L [ u ( 1 ) ] ) , u ( 3 ) = 1 3 u n + 2 3 ( u ( 2 ) + Δ t L [ u ( 2 ) ] ) ,
where L is the spatial operator.

4.1. Linear Advection Problems

Consider the following linear advection problems
{ u t + u x = 0 , x [ a , b ] , u ( x , t = 0 ) = u 0 ( x ) , periodic boundary ,
with the exact solution u ( x , t ) = u 0 ( x t ) .

4.1.1. Accuracy Test

Case 1. The smooth initial condition is given as
u 0 ( x ) = sin ( π x ) , 1 x 1 .
This solution is computed up to t = 2.0 using the grids N = 40, 80, 160, 320, 640, and 1280. The time step is set as Δ x 5 / 4 . In Table 1, the L 1 , L 2 , and L errors, convergence orders, and CPU times are exhibited for the WENO-Z and WENO-NZ schemes with different powers, respectively. Here, the norm of errors is computed by
L 1 = 1 N j = 1 N | u j ( u e x a c t ) j | , L 2 = 1 N j = 1 N ( u j ( u e x a c t ) j ) 2 , L = max 1 j N | u j ( u e x a c t ) j | .
It can be seen that all schemes can achieve fifth-order accuracy, but the CPU time of WENO-Z1 is the shortest compared to other schemes.
Case 2. The second initial condition for the accuracy test is given as
u 0 ( x ) = sin ( π x sin ( π x ) π ) , 1 x 1 .
This solution has two critical points. We compute the solution up to t = 2.0 using the grids N = 40, 80, 160, 320, 640, and 1280 with Δ t = Δ x 5 / 4 . The L 1 , L 2 , and L errors, convergence orders, and CPU times for different schemes are exhibited in Table 2. It can be found that the L errors of WENO-Z2, WENO-NZ1, and WENO-NZ2 are one order of magnitude lower than WENO-Z1. For ease of observation, Figure 5 gives the line graphs of L errors and orders. From Figure 5a, it can be seen that the L errors of WENO-Z2, WENO-NZ1, and WENO-NZ2 are almost the same and are smaller than those of WENO-Z1. Moreover, it can also be seen from Figure 5b that the L order of WENO-Z1 is only fourth-order accurate with the increase in grid numbers, while other schemes can achieve the desired fifth-order accuracy. Therefore, WENO-NZ proves to be a precise and effective scheme.

4.1.2. Linear Problem with Several Critical Points

The initial condition is given as
u 0 ( x ) = e ( x 90 ) 2 / 400 ( cos ( π 8 ( x 90 ) ) + cos ( π 4 ( x 90 ) ) ) , 50 x 130 .
This solution involves several critical points and is often used to test the numerical dissipation of WENO schemes. The solution is calculated by using N = 100 until the final time t = 400 with Δ t = Δ x / 2 . As Figure 6 shows, the numerical results of WENO-NZ are closer to the exact solution than WENO-Z, which again reveals that the WENO-NZ is a lower-dissipation scheme. In addition, the numerical results of WENO-NZ2 are almost comparable or even better than those of WENO-NZ1 at critical points, this is because the weights of WENO-NZ1 cannot satisfy the fifth-order convergence sufficient condition at critical points, but WENO-NZ2 can satisfy and also provide the lower dissipation.

4.2. One-Dimensional Euler Problems

The one-dimensional Euler equation can be written as
U t + F ( U ) x = 0 ,
with
U = ( ρ ρ u E ) , F ( U ) = ( ρ u ρ u 2 + p u ( E + p ) ) .
Here, ρ p, u , and E are the density, pressure, x-velocity, and total energy, respectively. The total energy is given as E = p γ 1 + 1 2 ρ u 2 with the ratio of specific heat γ = 1.4 . In this subsection, the extrapolation boundary conditions are adopted and the “exact” solution is computed using the fifth-order WENO-JS scheme with N = 2000.

4.2.1. SOD Problem

The initial conditions of the sod problem are as follows
( ρ , u , p ) = { ( 1 , 0 , 1 ) , 0 x < 0.5 , ( 0.125 , 0 , 0.1 ) , 0.5 x 1 .
This problem is run up to the final time t = 0.25 with N = 200, resulting in a shock wave, a contact discontinuity, and an expansion wave. From Figure 7a, all the schemes can both simulate the sod problem without oscillatory. The zoom-in results of contact discontinuity for different schemes are displayed in Figure 7b. It can be found that the WENO-NZ1 scheme can capture the contact discontinuities well and is slightly closer to the “exact” solution than the others.

4.2.2. Lax Problem

The initial conditions of the lax problem are as follows
( ρ , u , p ) = { ( 0.445 , 0.698 , 0.3528 ) , 5 x < 0 , ( 0.5 , 0 , 0.571 ) , 0 x 5 .
Figure 8 shows a comparison of different schemes with N = 200 at t = 1.3. From Figure 8a, we can see that all schemes can pass through the sound velocity point smoothly. However, it can be seen from Figure 8b that, in comparison to WENO-Z, the corresponding WENO-NZ scheme provides lower numerical dissipation.

4.2.3. Shu–Osher Problem

The initial conditions of Shu–Osher shock-entropy wave interaction are as follows
( ρ , u , p ) = { ( 3.857143 , 2.629369 , 10.33333 ) , 5 x < 4 , ( 1 + 0.2 s i n ( 5 x ) , 0 , 1.0 ) , 4 x 5 .
This problem involves the presence of shocklets and fine-scale structures, it is frequently utilized to study the stability of capturing shocks. We solve this problem until t = 1.8 using N = 200 and present the density results for different schemes in Figure 9. Examining the high-frequency wave region as depicted in Figure 9b, it can be seen that the WENO-NZ scheme performs a better job in approaching the “exact” solution compared to WENO-Z. Especially, in the area [1.3, 2.3], the numerical dissipation of WENO-NZ2 is even equivalent to WENO-NZ1. Based on the observations from Figure 7, Figure 8 and Figure 9, it can be concluded that the proposed WENO-NZ scheme effectively reduces numerical dissipation and is capable of stable shock capture.

4.3. Two-Dimensional Euler Problems

The two-dimensional Euler equation can be written as
U t + F ( U ) x + G ( U ) y = 0 ,
with
U = ( ρ ρ u ρ v E ) , F ( U ) = ( ρ u ρ u 2 + p ρ u v u ( E + p ) ) , G ( U ) = ( ρ v ρ u v ρ v 2 + p v ( E + p ) ) .
Here, ρ p, u , and v are the density, pressure, x- and y-velocity, respectively. The total energy E is given as E = p γ 1 + 1 2 ρ ( u 2 + v 2 ) with the ratio of specific heat γ = 1.4 , unless otherwise indicated.

4.3.1. Riemann Problem

The initial conditions of the two-dimensional Riemann problem are as follows
( ρ , u , v , p ) = { ( 1.5 , 0 , 0 , 1.5 ) , 0.8 x 1 , 0.8 y 1 , ( 0.5323 , 1.206 , 0 , 0.3 ) , 0 x < 0.8 , 0.8 y 1 , ( 0.138 , 1.206 , 1.206 , 0.029 ) , 0 x < 0.8 , 0 y < 0.8 , ( 0.5323 , 0 , 1.206 , 0.3 ) , 0.8 x 1 , 0 y < 0.8 .
With the computational domain [ 0 , 1 ] × [ 0 , 1 ] , this problem involves some vortex structures generated by Kelvin–Helmholtz instability and is often used to investigate numerical dissipation. Figure 10 gives the density profiles computed by different schemes at t = 0.3 with uniform meshes of 800 × 800 and zero-order extrapolation boundary conditions. As depicted in the results, the WENO-NZ scheme captures reflected shock waves more accurately and clearly than WENO-Z, which indicates that the numerical dissipation of WENO-NZ is smaller than that of WENO-Z. Among these schemes, WENO-NZ1 exhibits the richest vortex structures. However, it is important to note that when the time is pushed forward to longer, WENO-NZ1 will break the symmetrical roll-up of the Kelvin–Helmholtz instability because its numerical dissipation is too small. Fortunately, WENO-NZ2 not only performs well but also can maintain symmetry to solve this problem.

4.3.2. Double Mach Reflection Problem

The double Mach reflection problem depicts a right-moving Mach 10 shock wave with an angle of 60° hitting a reflecting wall at x = 1/6, which lies at the bottom of the computational domain [ 0 , 4 ] × [ 0 , 1 ] . This results in the formation of complex flow patterns, including shock wave reflections and vortex structures (see Figure 11, which is computed by the WENO-Z1 scheme using grids 960 × 240 to t = 0.25). The initial conditions are as follows
( ρ , u , v , p ) = { ( 1.4 , 0 , 0 , 1 ) , y < 3 ( x 1 / 6 ) , ( 8 , 7.145 , 4.125 , 116.833 ) , y 3 ( x 1 / 6 ) .
The exact post shock condition is imposed from x = 0 to x = 1/6 and the reflective boundary condition is used for the rest of the bottom. The boundary condition at the top corresponds to the exact motion of a Mach 10 shock, while the left and right boundaries are subject to inflow and outflow boundary conditions, respectively. Figure 12 shows the density results near the Mach stem computed by different schemes with uniform meshes of 960 × 240 . It can be seen that the WENO-NZ scheme depicts richer vortex structures and has a higher resolution compared with the WENO-Z scheme. In addition, Figure 13 further shows the local density profiles computed by using denser meshes of 1600 × 400 , it can be also seen that the WENO-NZ scheme still performs better results than WENO-Z. However, Figure 12 and Figure 13 both show that WENO-NZ1 exhibits Gibbs oscillation due to its numerical dissipation being too small, while WENO-NZ2 does not. Based on the numerical results of the Riemann problem and double Mach reflection problem, it is not difficult to find that the numerical dissipation is too small and would not be a good choice for some two-dimensional problems involving small vortex structures, so we suggest taking the power to be q = 2 for WENO-NZ to calculate these problems.

4.3.3. Rayleigh–Taylor Instability Problem

The Rayleigh–Taylor instability problem describes the interface instability when a heavy fluid accelerates into a light one, whose initial conditions are given as
( ρ , u , v , p ) = { ( 2 , 0 , 0.025 γ p / ρ c o s ( 8 π x ) , 2 y + 1 ) , 0 y < 0.5 , ( 1 , 0 , 0.025 γ p / ρ c o s ( 8 π x ) , y + 1.5 ) , 0.5 y 1 .
The computational domain is [ 0 , 0.25 ] × [ 0 , 1 ] and the ratio of specific heat is γ = 5 / 3 . The reflective boundary conditions are employed for the left and right boundaries. The boundary conditions on the top and bottom boundaries are given as ( ρ , u , v , p ) = ( 1 , 0 , 0 , 2.5 ) and ( ρ , u , v , p ) = ( 2 , 0 , 0 , 1 ) , respectively. Figure 14 gives the density results computed by different schemes with meshes of 240 × 960 until t = 1.95. It is apparent that all schemes are able to maintain the symmetry of the growth of instability. Compared to WENO-Z, the proposed WENO-NZ scheme displays more finer microstructures. Figure 15 further gives the density results computed to t = 1.95 by using denser meshes of 400 × 1600 . Obviously, the WENO-NZ scheme exhibits a greater richness of complex vortex structures and better resolution than the WENO-Z scheme.

5. Conclusions

In this study, we modified the global smoothness indicator of the fifth-order WENO-Z scheme to improve its performance. First, the five-point global stencil was subdivided into four smaller sub-stencils that are usually used to construct the third-order WENO-Z scheme. Based on these lower-order smoothness indicators, a novel global smoothness indicator with fifth-order accuracy was reconstructed, which resulted in the new WENO-NZ scheme. The results from ADR analysis implied that the proposed scheme was a low-dissipation scheme, whatever the value of power q. Accuracy tests showed that the WENO-NZ scheme with either q = 1 or q = 2 could achieve fifth-order accuracy even at critical points. Compared with the WENO-Z scheme, although the WENO-NZ had the same theoretical convergence accuracy, the WENO-Z failed to balance the low-dissipation property and fifth-order convergence. Moreover, a series of one- and two-dimensional benchmark problems demonstrated the WENO-NZ scheme showed lower numerical dissipation and better resolution than WENO-Z. However, for some two-dimensional problems involving small vortex structures, the power q was best to be 2 for the WENO-NZ scheme since the numerical dissipation was too small to suppress oscillation when q = 1.

Author Contributions

Conceptualization, S.H. and M.L.; methodology, S.H. and M.L.; software, S.H.; validation, S.H.; formal analysis, S.H. and M.L.; investigation, S.H.; resources, S.H.; data curation, S.H.; writing—original draft preparation, S.H.; writing—review and editing, S.H. and M.L.; visualization, S.H.; supervision, M.L.; project administration, M.L.; funding acquisition, M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China under grant no. (11971411).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Harten, A. High resolution schemes for hyperbolic conservation laws. J. Comput. Phys. 1983, 49, 357–393. [Google Scholar] [CrossRef]
  2. Toro, E.F.; Billett, S.J. Centred TVD schemes for hyperbolic conservation laws. IMA J. Numer. Anal. 2000, 20, 47–79. [Google Scholar] [CrossRef]
  3. Tang, S.J.; Li, M.J. Construction and application of several new symmetrical flux limiters for hyperbolic conservation law. Comput. Fluids 2020, 93, 104741. [Google Scholar] [CrossRef]
  4. Harten, A.; Osher, S.; Engquist, B.; Chakravarthy, S.R. Some results on uniformly high-order accurate essentially nonoscillatory schemes. Appl. Numer. Math. 1986, 2, 347–377. [Google Scholar] [CrossRef]
  5. Harten, A.; Engquist, B.; Osher, S.; Chakravarthy, S.R. Uniformly high order accurate essentially non-oscillatory schemes, III. J. Comput. Phys. 1987, 71, 231–303. [Google Scholar] [CrossRef]
  6. Shu, C.W.; Osher, S. Efficient implementation of essentially non-oscillatory shock-capturing schemes. J. Comput. Phys. 1988, 77, 439–471. [Google Scholar] [CrossRef]
  7. Liu, X.D.; Osher, S.; Chan, T. Weighted essentially non-oscillatory schemes. J. Comput. Phys. 1994, 115, 200–212. [Google Scholar] [CrossRef]
  8. Jiang, G.S.; Shu, C.W. Efficient implementation of weighted ENO schemes. J. Comput. Phys. 1996, 126, 202–228. [Google Scholar] [CrossRef]
  9. Borges, R.; Carmona, M.; Costa, B.; Don, W.S. An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws. J. Comput. Phys. 2008, 227, 3191–3211. [Google Scholar] [CrossRef]
  10. Tang, S.J.; Li, M.J. Novel functional weights for improving the third-order WENO schemes. Int. J. Numer. Meth. Fluids 2021, 93, 3131–3150. [Google Scholar] [CrossRef]
  11. Henrick, A.K.; Aslam, T.D.; Powers, J.M. Mapped weighted essentially non-oscillatory schemes: Achieving optimal order near critical points. J. Comput. Phys. 2005, 207, 542–567. [Google Scholar] [CrossRef]
  12. Guo, J.; Jung, J.H. Radial basis function ENO and WENO finite difference methods based on the optimization of shape parameters. J. Sci. Comput. 2017, 70, 551–575. [Google Scholar] [CrossRef]
  13. Guo, J.; Jung, J.H. A RBF-WENO finite volume method for hyperbolic conservation laws with the monotone polynomial interpolation method. Appl. Numer. Math. 2017, 112, 27–50. [Google Scholar] [CrossRef]
  14. Abedian, R. A finite difference Hermite RBF-WENO scheme for hyperbolic conservation laws. Int. J. Numer. Meth. Fluids 2022, 94, 583–607. [Google Scholar] [CrossRef]
  15. Abedian, R.; Dehghan, M. The formulation of finite difference RBFWENO schemes for hyperbolic conservation laws: An alternative technique. Adv. Appl. Math. Mech. 2023, 15, 1023–1055. [Google Scholar] [CrossRef]
  16. Zeng, F.J.; Shen, Y.Q.; Liu, S.P. A perturbational weighted essentially non-oscillatory scheme. Comput. Fluids 2018, 172, 196–208. [Google Scholar] [CrossRef]
  17. Castro, M.; Costa, B.; Don, W.S. High order weighted essentially non-oscillatory WENO-Z schemes for hyperbolic conservation laws. J. Comput. Phys. 2011, 230, 1766–1792. [Google Scholar] [CrossRef]
  18. Cheng, X.H.; Feng, J.H.; Zheng, S.P. A fourth order WENO scheme for hyperbolic conservation laws. Adv. Appl. Math. Mech. 2020, 12, 992–1007. [Google Scholar] [CrossRef]
  19. Xu, W.Z.; Wu, W.G. An improved third-order weighted essentially non-oscillatory scheme achieving optimal order near critical points. Comput. Fluids 2018, 162, 113–125. [Google Scholar] [CrossRef]
  20. Wang, Y.H.; Du, Y.L.; Zhao, K.L.; Li, Y. A low-dissipation third-order weighted essentially nonoscillatory scheme with a new reference smoothness indicator. Int. J. Numer. Meth. Fluids 2020, 92, 1212–1234. [Google Scholar] [CrossRef]
  21. Kim, C.H.; Ha, Y.; Yoon, J. Modified non-linear weights for fifth-order weighted essentially non-oscillatory schemes. J. Sci. Comput. 2016, 67, 299–323. [Google Scholar] [CrossRef]
  22. Rathan, S.; Gande, N.R.; Bhise, A.A. Simple smoothness indicator WENO-Z scheme for hyperbolic conservation laws. Appl. Numer. Math. 2020, 157, 255–275. [Google Scholar] [CrossRef]
  23. Liu, S.P.; Shen, Y.Q.; Zeng, F.J.; Ming, Y. A new weighting method for improving the WENO-Z scheme. Int. J. Numer. Meth. Fluids 2018, 87, 271–291. [Google Scholar] [CrossRef]
  24. Acker, F.; Borges, R.; Costa, B. An improved WENO-Z scheme. J. Comput. Phys. 2016, 313, 726–753. [Google Scholar] [CrossRef]
  25. Tang, S.J.; Feng, Y.J.; Li, M.J. Novel weighted essentially non-oscillatory schemes with adaptive weights. Appl. Math. Comput. 2022, 420, 126893. [Google Scholar] [CrossRef]
  26. Don, W.S.; Borges, R. Accuracy of the weighted essentially non-oscillatory conservative finite difference schemes. J. Comput. Phys. 2013, 250, 347–372. [Google Scholar] [CrossRef]
  27. Pirozzoli, S. On the spectral properties of shock-capturing scheme. J. Comput. Phys. 2006, 219, 489–497. [Google Scholar] [CrossRef]
Figure 1. The construction stencils of the new global smoothness indicator.
Figure 1. The construction stencils of the new global smoothness indicator.
Mathematics 11 04449 g001
Figure 2. Numerical solutions of the linear equation with discontinuity computed by different WENO schemes; (a) density distribution; (b) the values of smoothness indicators at t = 0.
Figure 2. Numerical solutions of the linear equation with discontinuity computed by different WENO schemes; (a) density distribution; (b) the values of smoothness indicators at t = 0.
Mathematics 11 04449 g002
Figure 3. The nonlinear weights of different schemes for the linear equation with discontinuity at t = 0; (a) WENO-Z1; (b) WENO-NZ1; (c) WENO-Z2; (d) WENO-NZ2.
Figure 3. The nonlinear weights of different schemes for the linear equation with discontinuity at t = 0; (a) WENO-Z1; (b) WENO-NZ1; (c) WENO-Z2; (d) WENO-NZ2.
Mathematics 11 04449 g003
Figure 4. The dispersion and dissipation properties of WENO-Z and WENO-NZ schemes; (a) dispersion; (b) dissipation.
Figure 4. The dispersion and dissipation properties of WENO-Z and WENO-NZ schemes; (a) dispersion; (b) dissipation.
Mathematics 11 04449 g004
Figure 5. The numerical results for the linear advection problem with two critical points computed by different WENO schemes; (a) L errors; (b) L orders.
Figure 5. The numerical results for the linear advection problem with two critical points computed by different WENO schemes; (a) L errors; (b) L orders.
Mathematics 11 04449 g005
Figure 6. Numerical solutions of linear advection problem with several critical points computed by different WENO schemes.
Figure 6. Numerical solutions of linear advection problem with several critical points computed by different WENO schemes.
Mathematics 11 04449 g006
Figure 7. Numerical solutions of sod problem computed by different schemes; (a) density distribution; (b) a zoom of contact discontinuity.
Figure 7. Numerical solutions of sod problem computed by different schemes; (a) density distribution; (b) a zoom of contact discontinuity.
Mathematics 11 04449 g007
Figure 8. Numerical solutions of lax problem computed by different schemes; (a) density distribution; (b) a zoom of discontinuity.
Figure 8. Numerical solutions of lax problem computed by different schemes; (a) density distribution; (b) a zoom of discontinuity.
Mathematics 11 04449 g008
Figure 9. Numerical solutions of the Shu-Osher problem computed by different schemes; (a) density distribution; (b) a zoom of the high-frequency wave.
Figure 9. Numerical solutions of the Shu-Osher problem computed by different schemes; (a) density distribution; (b) a zoom of the high-frequency wave.
Mathematics 11 04449 g009
Figure 10. Numerical solutions of 2D Riemann problem computed by different schemes; (a) WENO-Z1; (b) WENO-Z2; (c) WENO-NZ1; (d) WENO-NZ2.
Figure 10. Numerical solutions of 2D Riemann problem computed by different schemes; (a) WENO-Z1; (b) WENO-Z2; (c) WENO-NZ1; (d) WENO-NZ2.
Mathematics 11 04449 g010
Figure 11. Numerical solutions of double Mach reflection problem computed by the WENO-Z1 scheme with uniform meshes of 960 × 240 .
Figure 11. Numerical solutions of double Mach reflection problem computed by the WENO-Z1 scheme with uniform meshes of 960 × 240 .
Mathematics 11 04449 g011
Figure 12. A zoom of numerical solutions near the Mach stem for the double Mach reflection problem, which is computed by different schemes with uniform meshes of 960 × 240 ; (a) WENO-Z1; (b) WENO-Z2; (c) WENO-NZ1; (d) WENO-NZ2.
Figure 12. A zoom of numerical solutions near the Mach stem for the double Mach reflection problem, which is computed by different schemes with uniform meshes of 960 × 240 ; (a) WENO-Z1; (b) WENO-Z2; (c) WENO-NZ1; (d) WENO-NZ2.
Mathematics 11 04449 g012
Figure 13. A zoom of numerical solutions near the Mach stem for the double Mach reflection problem, which is computed by different schemes with uniform meshes of 1600 × 400 ; (a) WENO-Z1; (b) WENO-Z2; (c) WENO-NZ1; (d) WENO-NZ2.
Figure 13. A zoom of numerical solutions near the Mach stem for the double Mach reflection problem, which is computed by different schemes with uniform meshes of 1600 × 400 ; (a) WENO-Z1; (b) WENO-Z2; (c) WENO-NZ1; (d) WENO-NZ2.
Mathematics 11 04449 g013
Figure 14. Numerical solutions of Rayleigh–Taylor instability problem computed by different schemes with uniform meshes of 240 × 960 ; (a) WENO-Z1; (b) WENO-Z2; (c) WENO-NZ1; (d) WENO-NZ2.
Figure 14. Numerical solutions of Rayleigh–Taylor instability problem computed by different schemes with uniform meshes of 240 × 960 ; (a) WENO-Z1; (b) WENO-Z2; (c) WENO-NZ1; (d) WENO-NZ2.
Mathematics 11 04449 g014
Figure 15. Numerical solutions of Rayleigh–Taylor instability problem computed by different schemes with uniform meshes of 400 × 1600 ; (a) WENO-Z1; (b) WENO-Z2; (c) WENO-NZ1; (d) WENO-NZ2.
Figure 15. Numerical solutions of Rayleigh–Taylor instability problem computed by different schemes with uniform meshes of 400 × 1600 ; (a) WENO-Z1; (b) WENO-Z2; (c) WENO-NZ1; (d) WENO-NZ2.
Mathematics 11 04449 g015
Table 1. The errors, orders, and CPU times for linear advection problem with smooth initial condition computed by different WENO schemes.
Table 1. The errors, orders, and CPU times for linear advection problem with smooth initial condition computed by different WENO schemes.
MethodN L 1 ErrorOrder L 2 ErrorOrder L ErrorOrderTime
WENO-Z1406.4906 × 10−6---7.3342 × 10−6---1.0792 × 10−5---0.09360
802.0181 × 10−75.012.2523 × 10−75.093.2397 × 10−75.060.31200
1606.3021 × 10−95.007.0126 × 10−95.039.9715 × 10−95.021.23241
3201.9695 × 10−105.002.1894 × 10−105.013.1007 × 10−105.015.44443
6406.1552 × 10−125.006.8396 × 10−125.009.6797 × 10−125.0025.02256
12801.9237 × 10−135.002.1371 × 10−135.003.0235 × 10−135.00117.01635
WENO-Z2406.4581 × 10−6---7.2051× 10−6---1.0281 × 10−5---0.10920
802.0168 × 10−75.002.2462 × 10−75.003.1934 × 10−75.010.24960
1606.3015 × 10−95.007.0101 × 10−95.009.9414 × 10−95.011.21681
3201.9694 × 10−105.002.1893 × 10−105.003.1006 × 10−105.005.55364
6406.1552 × 10−125.006.8396 × 10−125.009.6797 × 10−125.0024.96016
12801.9237 × 10−135.002.1371 × 10−135.003.0235 × 10−135.00120.44837
WENO-NZ1406.4712× 10−6---7.2545 × 10−6---1.0476 × 10−5---0.03120
802.0173 × 10−75.002.2486 × 10−75.013.2110 × 10−75.030.37440
1606.3017 × 10−95.007.0111 × 10−95.0019.9534× 10−95.011.40401
3201.9694 × 10−105.002.1893 × 10−105.003.1006 × 10−105.006.30244
6406.1552 × 10−125.006.8396 × 10−125.009.6797 × 10−125.0029.37499
12801.9237 × 10−135.002.1371 × 10−135.003.0235 × 10−135.00136.98448
WENO-NZ2406.4580× 10−6---7.2028 × 10−6---1.0281 × 10−5---0.15600
802.0168 × 10−75.002.2462 × 10−75.003.1934 × 10−75.010.35880
1606.3015 × 10−95.007.0101 × 10−95.009.9414 × 10−95.011.41961
3201.9694 × 10−105.002.1893 × 10−105.003.1006 × 10−105.006.17764
6406.1552 × 10−125.006.8396 × 10−125.009.6797 × 10−125.0029.46859
12801.9237 × 10−135.002.1371 × 10−135.003.0235 × 10−135.00138.91889
Table 2. The errors, orders, and CPU times for the linear advection problem with two critical points were computed by different WENO schemes.
Table 2. The errors, orders, and CPU times for the linear advection problem with two critical points were computed by different WENO schemes.
MethodN L 1 ErrorOrder L 2 ErrorOrder L ErrorOrderTime
WENO-Z1407.0429 × 10−5---9.5195 × 10−5---2.1744 × 10−4---0.06240
802.4102 × 10−64.873.1785 × 10−64.906.6772 × 10−65.030.28080
1607.8990 × 10−84.931.0142 × 10−74.972.0989 × 10−74.991.12321
3202.5157 × 10−94.973.2185 × 10−94.987.7672 × 10−94.765.60044
6407.8219 × 10−115.001.0239 × 10−104.973.5698 × 10−104.4425.47496
12802.4235 × 10−125.013.2682 × 10−124.971.7250 × 10−114.37118.74796
WENO-Z2406.6140 × 10−5---9.2603 × 10−5---2.2383 × 10−4---0.09360
802.2447 × 10−64.883.0595 × 10−64.926.6812 × 10−65.070.37440
1607.2388 × 10−84.959.6470 × 10−84.992.0987 × 10−74.991.23241
3202.2821 × 10−94.993.0181 × 10−95.006.5525 × 10−95.005.46004
6407.1426 × 10−115.009.4316 × 10−115.002.0464 × 10−105.0026.31737
12802.2327 × 10−125.002.9471 × 10−125.006.3923 × 10−125.00121.25958
WENO-NZ1407.1643 × 10−5---9.5751 × 10−5---2.1375 × 10−4---0.12480
802.2953 × 10−64.963.0766 × 10−64.966.6998 × 10−65.000.48360
1607.2863 × 10−84.989.6755 × 10−84.992.0990 × 10−75.001.27921
3202.3087 × 10−94.983.0290 × 10−95.006.5526 × 10−95.006.22444
6407.2598 × 10−114.999.4765 × 10−115.002.0464 × 10−105.0028.90699
12802.2695 × 10−125.002.9652 × 10−125.006.3923 × 10−125.00137.24968
WENO-NZ2407.3503 × 10−5---9.6607 × 10−5---2.1206 × 10−4---0.07800
802.3293 × 10−64.983.0796 × 10−64.976.7004 × 10−64.980.35880
1607.3012 × 10−85.009.6541 × 10−85.002.0988 × 10−75.001.48201
3202.2852 × 10−95.003.0181 × 10−95.006.5526 × 10−95.006.31804
6407.1435 × 10−115.009.4313 × 10−115.002.0464 × 10−105.0029.28139
12802.2327 × 10−125.002.9471 × 10−125.006.3923 × 10−125.00139.57409
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Han, S.; Li, M. An Improved WENO-Z Scheme for Hyperbolic Conservation Laws with New Global Smoothness Indicator. Mathematics 2023, 11, 4449. https://doi.org/10.3390/math11214449

AMA Style

Han S, Li M. An Improved WENO-Z Scheme for Hyperbolic Conservation Laws with New Global Smoothness Indicator. Mathematics. 2023; 11(21):4449. https://doi.org/10.3390/math11214449

Chicago/Turabian Style

Han, Shuang, and Mingjun Li. 2023. "An Improved WENO-Z Scheme for Hyperbolic Conservation Laws with New Global Smoothness Indicator" Mathematics 11, no. 21: 4449. https://doi.org/10.3390/math11214449

APA Style

Han, S., & Li, M. (2023). An Improved WENO-Z Scheme for Hyperbolic Conservation Laws with New Global Smoothness Indicator. Mathematics, 11(21), 4449. https://doi.org/10.3390/math11214449

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop