Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Two classes of improved fifth-order weighted compact nonlinear schemes for compressible Euler equations

  • Two enhanced fifth-order weighted compact nonlinear schemes (WCNSs), namely the WCNS scheme with the modified ZN-type weights (WCNS-ZNM) and the perturbed WCNS scheme (WCNS-P), were proposed for compressible Euler equations. The WCNS-ZNM scheme incorporates an enhanced WENO-ZN-type nonlinear interpolation technique to obtain unknown flow values at cell edges. The WCNS-P scheme introduces perturbation terms containing a free parameter into the linear interpolations, effectively reducing numerical errors in smooth regions. A monotone polynomial interpolation approach was designed to enhance numerical stability by filtering out non-smooth regions and realizing the automatic switching between the WCNS-P and WCNS-ZNM schemes in smooth and non-smooth regions. Numerical results exhibit better accuracy and performance in resolving small-scale structures and shock waves, compared to the WCNS scheme with the WENO-Z and WENO-ZN weights.

    Citation: Huanhuan Yang, Yanqun Jiang, Qinghong Tang. Two classes of improved fifth-order weighted compact nonlinear schemes for compressible Euler equations[J]. Electronic Research Archive, 2025, 33(6): 3834-3856. doi: 10.3934/era.2025170

    Related Papers:

    [1] Shuaikang Wang, Yunzhi Jiang, Yongbin Ge . High-order compact difference methods for solving two-dimensional nonlinear wave equations. Electronic Research Archive, 2023, 31(6): 3145-3168. doi: 10.3934/era.2023159
    [2] Jinxiu Zhang, Xuehua Yang, Song Wang . The ADI difference and extrapolation scheme for high-dimensional variable coefficient evolution equations. Electronic Research Archive, 2025, 33(5): 3305-3327. doi: 10.3934/era.2025146
    [3] Jiangshan Wang, Lingxiong Meng, Hongen Jia . Numerical analysis of modular grad-div stability methods for the time-dependent Navier-Stokes/Darcy model. Electronic Research Archive, 2020, 28(3): 1191-1205. doi: 10.3934/era.2020065
    [4] Li Tian, Ziqiang Wang, Junying Cao . A high-order numerical scheme for right Caputo fractional differential equations with uniform accuracy. Electronic Research Archive, 2022, 30(10): 3825-3854. doi: 10.3934/era.2022195
    [5] Huanhuan Li, Lei Kang, Meng Li, Xianbing Luo, Shuwen Xiang . Hamiltonian conserved Crank-Nicolson schemes for a semi-linear wave equation based on the exponential scalar auxiliary variables approach. Electronic Research Archive, 2024, 32(7): 4433-4453. doi: 10.3934/era.2024200
    [6] Zhaoxiang Zhang, Xuehua Yang, Song Wang . The alternating direction implicit difference scheme and extrapolation method for a class of three dimensional hyperbolic equations with constant coefficients. Electronic Research Archive, 2025, 33(5): 3348-3377. doi: 10.3934/era.2025148
    [7] Yunxia Niu, Chaoran Qi, Yao Zhang, Wahidullah Niazi . Numerical analysis and simulation of the compact difference scheme for the pseudo-parabolic Burgers' equation. Electronic Research Archive, 2025, 33(3): 1763-1791. doi: 10.3934/era.2025080
    [8] Changling Xu, Tianliang Hou . Superclose analysis of a two-grid finite element scheme for semilinear parabolic integro-differential equations. Electronic Research Archive, 2020, 28(2): 897-910. doi: 10.3934/era.2020047
    [9] Cheng Wang . Convergence analysis of Fourier pseudo-spectral schemes for three-dimensional incompressible Navier-Stokes equations. Electronic Research Archive, 2021, 29(5): 2915-2944. doi: 10.3934/era.2021019
    [10] Xiu Ye, Shangyou Zhang, Peng Zhu . A weak Galerkin finite element method for nonlinear conservation laws. Electronic Research Archive, 2021, 29(1): 1897-1923. doi: 10.3934/era.2020097
  • Two enhanced fifth-order weighted compact nonlinear schemes (WCNSs), namely the WCNS scheme with the modified ZN-type weights (WCNS-ZNM) and the perturbed WCNS scheme (WCNS-P), were proposed for compressible Euler equations. The WCNS-ZNM scheme incorporates an enhanced WENO-ZN-type nonlinear interpolation technique to obtain unknown flow values at cell edges. The WCNS-P scheme introduces perturbation terms containing a free parameter into the linear interpolations, effectively reducing numerical errors in smooth regions. A monotone polynomial interpolation approach was designed to enhance numerical stability by filtering out non-smooth regions and realizing the automatic switching between the WCNS-P and WCNS-ZNM schemes in smooth and non-smooth regions. Numerical results exhibit better accuracy and performance in resolving small-scale structures and shock waves, compared to the WCNS scheme with the WENO-Z and WENO-ZN weights.



    The compressible Euler equations, a fundamental system in fluid dynamics, have important research and application value in many fields and are expressed as

    Ut+F(U)x+G(U)y=0, (1)

    where

    U=(ρρuρvE),      F(U)=(ρuρu2+Pρuvu(E+P)),      G(U)=(ρvρuvρv2+Pv(E+P)), (2)
    E=Pγ1+12ρ(u2+v2). (3)

    Here, ρ denotes the fluid density, while u and v denote the speed components along the x- and y-axes, respectively. In Eq (3), E represents the total energy per unit fluid volume, with P indicating the fluid pressure, and γ=1.4 representing the specific heat ratio. Due to the hyperbolic property of the compressible Euler equations, their solutions may contain shock waves, discontinuities, and other complex phenomena. Therefore, it is difficult to solve these equations by applying high-order linear methods, because these methods are prone to result in numerical oscillations around jump discontinuities.

    The original weighted essentially non-oscillatory (WENO) scheme [1] is an ideal choice for solving Euler equations among many high-precision numerical schemes, due to its ability to maintain high accuracy in smooth flow regions while effectively suppressing oscillations near discontinuities. This scheme obtains the unknown variables using a nonlinear weighting of all the candidate stencil information. The fifth-order WENO-JS scheme developed by Jiang and Shu [2] introduced novel smoothness indicators for the sub-stencils. Based on that work, a large number of studies have focused on the improvement of the accuracy and computational efficiency of the WENO scheme. For instance, Henrick et al. [3] developed the WENO-M scheme by introducing a mapping weight function to modify JS-type nonlinear weights. Subsequently, Borges et al. [4] established the WENO-Z approach by creating a global smoothness indicator across the five-point stencil and formulating innovative non-oscillatory Z-type weights. By substituting polynomial interpolation with Hermite interpolation, Qiu and Shu [5,6] designed a fifth-order Hermite WENO (HWENO) scheme. Guo and Jung [7,8] presented a new WENO scheme by applying the radial basis function (RBF) reconstruction and designed a monotone polynomial interpolation method to identify non-smooth stencils, where the WENO-JS-type reconstruction was used to maintain non-oscillatory behavior. Building upon the framework of the RBF-WENO scheme established by Guo and Jung [8], researchers have developed several non-polynomial reconstruction approaches, such as the WENO-G [9], TENO-P [10], and HP-WENO [11] schemes. These approaches offer performance optimization through adjustable parameters, showcasing remarkable adaptability. Most recently, Zuo and Zhu [12,13] proposed an innovative high-order hybrid multi-resolution (HMR) WENO method with new troubled-cell indicators. Hu et al. [14] established the WENO-CU6 scheme by implementing an adaptive switching between central and upwind schemes. The WENO-ZQ scheme formulated by Zhu and Qiu [15] applied weighted combinations of polynomial reconstructions at different stencil scales. Building on this, Shen et al. [16] constructed a novel WENO-Z scheme (WENO-ZN) with adaptive weight modification, replacing the Z-type weight constant with a discontinuity-sensitive function. Li et al. [17] later verified that higher-order WENO-ZN variants outperform the WENO-JS and WENO-Z schemes of the same order.

    The weighted compact nonlinear scheme (WCNS), a high-precision scheme that was proposed by Deng and Zhang [18], is widely used for solving Euler equations. Nonomura and Fujii [19] formulated a new WCNS scheme by mixing the numerical flux of grid nodes and midpoints to improve accuracy and stability. Liu et al. [20] demonstrated that the WCNS scheme has the same strong shock-capturing capability as the classical WENO scheme, with low dissipation and high resolution. To improve the resolution of the WCNS scheme, Yan et al. [21] proposed an adaptive Y-type nonlinear weight. The weight is similar to the function in the WENO-ZN scheme and approaches a small value near discontinuities. Wong and Lele [22] proposed a low-dissipative WCNS scheme by using a linear scheme in smooth regions and a hybrid scheme of linear interpolation and WENO-type nonlinear interpolation near discontinuities to improve numerical stability and reduce numerical dissipation while capturing shock waves. To improve computational efficiency, Jiang et al. [23] presented a semi-implicit WCNS formulation with high-order accuracy specifically designed for stiff Navier–Stokes problems, avoiding the limitation of Courant-Friedrichs-Lewy (CFL) conditions required by explicit versions. Jiang et al. [24] extended a fifth-order WCNS scheme with WENO-Z-type weights for solving Hamilton-Jacobi equations that are related to hyperbolic conservation laws. For this kind of equations, Huang et al. [25] formulated a perturbed WCNS scheme (WCNS-P) with a free parameter and boosted the local accuracy by optimizing the parameter in the perturbed linear interpolations constructed on the two-point-wide sub-stencils.

    This paper develops two enhanced fifth-order WCNS schemes, namely the WCNS scheme with the modified ZN-type weights (WCNS-ZNM) and the perturbed WCNS scheme (WCNS-P), for Euler equations. The fifth-order WCNS-ZNM scheme incorporates the improvement on the WENO-ZN-type nonlinear interpolation originally proposed by Li and Shen et al. [17]. The fifth-order WCNS-P scheme is developed by incorporating perturbation terms containing a free parameter into linear interpolations obtained on three-point-wide sub-stencils, extending the WCNS-P approach of third-order accuracy [25]. Optimizing the free parameter in the WCNS-P scheme can minimize numerical errors in the smooth regions, enhancing local accuracy. A local monotone interpolation method is presented to detect non-smooth stencils, enabling automatic switching between the WCNS-P and WCNS-ZNM schemes according to the smoothness of the global stencils, thereby improving numerical stability. This paper is organized as follows: Sections 2 and 3 detail the formulations of WCNS-ZNM and WCNS-P schemes, respectively. Section 4 presents 1D and 2D numerical validations, with conclusions provided in Section 5.

    The fifth-order WCNS-ZNM scheme for 1D Euler equations is derived in this section. Multidimensional extensions can be obtained through the dimension-by-dimension method. A uniform grid (spacing h=(ba)/N) discretizes the domain Ω=[a,b]   (b>a), with cell centers located at xi=a+ih(i=0,1,...,N) and the total number of cells N.

    In classical fifth-order WCNS schemes, flux derivatives Fi at cell nodes xi are calculated by using a sixth-order hybrid difference scheme combining cell-node and cell-edge approximations.

    Fi=6445h(˜Fi+1/2˜Fi1/2)29h(˜Fi+1˜Fi1)1180h(˜Fi+2˜Fi2), (4)

    where ˜Fi±1/2 denote the approximations of the flux F at the cell edges xi±1/2, respectively. To guarantee numerical stability, the characteristic-wise interpolation technique is applied to compute the numerical fluxes at cell edges. The interpolation procedure for ˜Fi+1/2 in the local characteristic space involves these steps:

    1) Estimate the average states at cell edges as Ui+1/2=12(Ui+Ui+1).

    2) Obtain the Jacobian A=(F/U)i+1/2 and its right eigenvectors Ri+1/2=R(Ui+1/2) and left eigenvectors Li+1/2=L(Ui+1/2), i.e., A=Ri+1/2Λi+1/2Li+1/2, where Λi+1/2 is the diagonal matrix and its diagonal elements are λ1=uc,λ2=u,λ3=u+c. Here, u is the fluid speed in the x-direction, and c=γP/ρ is the sound speed.

    3) Apply the Lax–Friedrichs flux splitting, i.e., F=F++F with F±=12(F±αxU) and αx=maxk|λk|. Transform the point values F±i+m to the characteristic space as Q±i+m=Li+12F±i+m(m=2,...,3).

    4) Employ the WENO-ZN interpolation to calculate the left and right states ˜Qi+1/2+ and ˜Qi+1/2 at cell edges in the characteristic space.

    5) Transform the calculated quantities ˜Qi+1/2± from step 4 back to physical space and calculate the cell-edge numerical fluxes as ˜Fi+1/2=Ri+1/2(˜Qi+1/2++˜Qi+1/2).

    The WENO-ZN-type interpolation procedure of ˜Qi+1/2+ is described in detail by a component-by-component manner and by omitting the superscript "+" for clarity. The calculation of ˜Qi+1/2 follows the same procedure.

    To achieve fifth-order accuracy for each component ˜qi+1/2 of ˜Qi+1/2+, we employ a five-point-wide stencilS={xi2,xi1,xi,xi+1,xi+2} composed of three-point-wide sub-stencils Sk={xi+k2,xi+k1,xi+k},k=0,1,2. Using the global stencil S, the linear interpolation of ˜qi+1/2 can be obtained as

    ˜qi+1/2=1128(3qi220qi1+90qi+60qi+15qi+2). (5)

    This equation can alternatively be formulated as a linear weighting of third-order linear interpolations obtained from three sub-stencils Sk, denoted ˜qki+1/2, k=0,1,2.

    ˜qi+1/2=2k=0δk˜qki+1/2, (6)

    where the linear weights are δ0=116,δ1=1016,δ2=516. The expressions of the three third-order interpolations are

    ˜q0i+1/2=18(3qi210qi1+15qi),˜q1i+1/2=18(qi1+6qi+3qi+1),˜q2i+1/2=18(3qi+6qi+1qi+2). (7)

    The high-order linear interpolation (6) may lead to non-physical oscillations around solution discontinuities. To enhance stability and avoid this issue, nonlinear weighting is implemented using three third-order interpolations (7), and the linear weights δk in Eq (6) are replaced with the nonlinear weights ϖk for computing ˜qi+1/2.

    A series of nonlinear weights in the literature (e.g., JS-type [2], Z-type [4], and ZN-type [16,17]) are presented in the WENO versions. The WENO-ZN weights are verified to be superior to the WENO-JS and WCNS-Z weights, in terms of accuracy and resolution, and are expressed as [17]

    ϖk=αk2k=0αk,αk=δk(θ+(τ8βk+ϵ)2),k=0,1,2, (8)

    where ϵ=1012 prevents division by zero and the parameter θ is defined as

    θ=10(β0+β2|β0β2|+ϵ|β0β2|+ϵ)2. (9)

    Here, θ is related to the local smoothness indicators of the sub-stencils Sk, denoted as βk, and is easily affected by numerical disturbances. To resolve this issue, we replace it with a large positive constant. The local smoothness indicators are defined as βk=(hpk)2+(h2pk)2,k=0,1,2, and the specific formulas are

    β0=14(qi24qi1+3qi)2+(qi22qi1+qi)2,β1=14(qi1qi+1)2+(qi12qi+qi+1)2,β2=14(3qi4qi+1+qi+2)2+(qi2qi+1+qi+2)2. (10)

    Define τ8 as the global smoothness indicator for stencil S:

    τ8=(qi24qi1+6qi4qi+1+qi+2)2. (11)

    According to Eqs (7) and (8), the nonlinear weighting of the interpolations of ˜qi+1/2 yields

    ˜qi+1/2=2k=oϖk˜qki+1/2. (12)

    Notice that the WCNS-ZN scheme is composed of a sixth-order central difference scheme (4) for the spatial derivative Fi and a WENO-ZN-type interpolation (12) with the parameter (9) for the numerical fluxes at cell edges. The WCNS-ZNM scheme is designed by replacing the parameter (9) in the WCNS-ZN scheme with a large positive constant.

    Guo and Jung [7,8] proposed the RBF-WENO scheme by selecting multiple quadratic forms and Gaussian RBF with free parameters as the non-polynomial interpolation basis. The local interpolation errors can be controlled by adjusting these parameters. Building on the RBF-WENO concept, Huang et al. [25] established a third-order WCNS-P scheme for Hamilton-Jacobi problems. This work advances the framework to fifth-order accuracy using a five-point stencil S.

    We add perturbation terms with an adjustable parameter ς2 into Eq (7), thereby deriving the following non-polynomial interpolation formula for three-point-wide sub-stencils.

    ˜q0,pi+1/2=(38+a1ς2h2)qi2+(108+b1ς2h2)qi1+(158+c1ς2h2)qi,˜q1,pi+1/2=(18+a2ς2h2)qi1+(68+b2ς2h2)qi+(38+c2ς2h2)qi+1,˜q2,pi+1/2=(38+a3ς2h2)qi+(68+b3ς2h2)qi+1+(18+c3ς2h2)qi+2. (13)

    In Eq (13), the parameter ς2 and the constants al,bl,cl(l=1,2,3) are determined by Taylor series expansion. The flux at cell edges is then evaluated via linear combination of the resulting non-polynomial interpolations, i.e.,

    ˜qi+1/2=2k=0δk˜qk,pi+1/2. (14)

    We expand the right-hand terms of Eq (13) using the Taylor expansion and derive the following expressions:

    ˜q0,pi+1/2=qi+1/2+h2[(a1+b1+c1)ς2qi+1/2]+h3[516q(3)i+1/2+(52a132b112c1)ς2qi+1/2]+h4[(258a1+98b1+18c1)ς2qi+12+45128q(4)i+12]+O(h5), (15)
    ˜q1,pi+1/2=qi+1/2+h2[(a2+b2+c2)ς2qi+1/2]+h3[116q(3)i+1/2+(32a212b2+12c2)ς2qi+1/2]+h4[(98a2+18b2+18c2)ς2qi+1/23128q(4)i+1/2]+O(h5), (16)
    ˜q2,pi+1/2=qi+1/2+h2[(a3+b3+c3)ς2qi+1/2]+h3[116q(3)i+1/2+(12a3+12b3+32c3)ς2qi+1/2]+h4[(18a3+18b3+98c3)ς2qi+1/23128q(4)i+1/2]+O(h5). (17)

    In order to ensure that the non-polynomial interpolations (13) have at least third-order accuracy, all the higher-order terms in Eqs (15)–(17) must be neglected, and the value of the parameter ς2 is obtained as

    ς2=58(5a13b1c1)q(3)i+1/2qi+1/2=124a28b2+8c2q(3)i+1/2qi+1/2=18a3+8b3+24c3q(3)i+1/2qi+1/2. (18)

    Substituting the expression (18) into Eqs (15)–(17) yields

    ˜q0,pi+1/2=qi+1/2+[564(25a1+9b1+c15a13b1c1q(3)i+1/2qi+1/2qi+1/2)+45128q(4)i+1/2]h4+O(h5), (19)
    ˜q1,pi+1/2=qi+1/2+[164(9a2+b2+c23a2b2+c2q(3)i+1/2qi+1/2qi+1/2)3128q(4)i+1/2]h4+O(h5), (20)
    ˜q2,pi+1/2=qi+1/2+[164(a3+b3+9c3a3+b3+3c3q(3)i+1/2qi+1/2qi+1/2)3128q(4)i+1/2]h4+O(h5). (21)

    The derivatives qi+1/2 and q(3)i+1/2 in Eq (18) are numerically approximated by employing the nodal values of stencil S.

    qi+1/2=qi127qi+27qi+1qi+224h, (22)
    q(3)i+1/2=qi1+3qi3qi+1+qi+2h3. (23)

    Substituting Eqs (22) and (23) into expression (18) generates

    ς258(5a13b1c1)24(qi1+3qi3qi+1+qi+2)h2(qi127qi+27qi+1qi+2)124a2+8b28c224(qi1+3qi3qi+1+qi+2)h2(qi127qi+27qi+1qi+2)18a3+8b3+24c324(qi1+3qi3qi+1+qi+2)h2(qi127qi+27qi+1qi+2) (24)

    According to Eqs (19)–(21) and (24), it follows that the constants al,bl,cl(l=1,2,3) satisfy the following conditions:

    a1+b1+c1=0, a2+b2+c2=0, a3+b3+c3=0, (25)
    25a1+9b1+c1=0, 9a2+b2+c2=0, a3+b3+9c3=0, (26)
    5a13b1c10, 3a2+b2c20, a3+b3+3c30, (27)
    58(5a13b1c1)=124a2+8b28c2=18a3+8b3+24c3. (28)

    Solving the above-mentioned linear equations subject to inequality constraints, we have

    a1=5,    b1=15,    c1=10,a2=0,    b2=1,    c2=1,a3=1,    b3=1,    c3=0. (29)

    Based on Eq (29), we plug Eq (13) into the scheme (14) and have

    ˜qi+1/2=1128(3qi220qi1+90qi+60qi+15qi+2)+116ς2h2(5qi215qi1+15qi5qi+1). (30)

    Eq (30) is expanded as a Taylor series, resulting in the truncation error as follows:

    τWCNSP=4732304h5q(5)i+1/2+O(h6). (31)

    Therefore, the WCNS-P scheme composed of the sixth-order central difference scheme (4) for the spatial derivative Fi and the scheme (14) for cell-edge fluxes can attain fifth-order accuracy in smooth regions. If the global stencil S includes discontinuity, the WCNS-P scheme will not maintain the non-oscillating property due to the fact that the estimation of the free parameter ς2 in Eq (24) is carried out based on the global stencil.

    The WCNS-P scheme cannot keep the non-oscillating property in non-smooth regions. Thus, we must find out non-smooth solution regions in which the WCNS-ZNM scheme is employed to enhance stability. Several approaches are used in the literature to identify the discontinuity and then initiate the shock-capturing schemes. For instance, Zuo and Zhu [12,13] developed a series of parameter-independent troubled-cell indicators based on different stencil sizes to determine whether extreme points exist within the interval, thereby accurately assessing the smoothness of the stencil. In this work, we adopt the local monotonicity-preserving interpolation method proposed in [7]. Unlike previous works [12,13], this method does not require the construction of a fourth-degree polynomial for the five-point stencil. Instead, it only necessitates the construction of second-degree polynomials on the three-point sub-stencils. This methodology significantly simplifies the extreme-point detection process while maintaining effectiveness. This subsection will provide a detailed description of this methodology.

    We reconstruct a quadratic polynomial on each sub-stencil. For example, the schematic diagram of the interpolating polynomial on the sub-stencil S0, denoted by F0(x), is shown in Figure 1, where x0c represents the critical point at which F0(x)=0. Notice that the critical point could be generated due to the Gibbs oscillations. Without loss of generality, let xi2=0,xi1=h,xi=2h. The polynomial F0(x) is formulated as

    F0(x)=a0x2+b0x+c0,    x(0,2h). (32)

    The coefficients a0,b0,c0 are expressed as

    a0=qi22qi1+qi2h2, b0=3qi2qi1+qi2h, c0=qi2. (33)

    The critical point x0c is computed as

    x0c=b02a0=3qi24qi1+qiqi22qi1+qih. (34)

    It is obvious that F0(x) is monotone if x0c0 or x0c2h, while the derivative F0(x) changes its sign in the sub-stencil S0, if x0c is located in the interval. If F0(x) is not monotone, it could be oscillatory resulting from the Gibbs phenomenon.

    Figure 1.  The schematic diagram of the interpolating polynomial constructed on S0.

    Let xkc represent the critical points of the quadratic polynomials Fk(x) constructed on Sk, k=0,1,2. If any of the three sub-stencils contains a critical point, the five-point-wide stencil S is assumed to be non-smooth, and then the numerical fluxes obtained based on the information from the stencil S are calculated using the WCNS-ZNM scheme. This discontinuity identification method realizes the adaptive switching between the two schemes in smooth regions and non-smooth regions, i.e.,

    S={WCNSZNMtype stencil,  if x0cI0   or   x1cI1   or   x2cI2,WCNSPtype stencil,      otherwise. (35)

    It follows that the designed WCNS-P scheme is essentially a hybrid scheme. Moreover, the proposed scheme naturally converts to the shock-capturing WCNS-ZNM scheme in non-smooth flow regions.

    High-order shock-capturing schemes with low dissipation and dispersion can accurately simulate the flow field with complicated structures. We apply the approximate dispersion relation (ADR) [26] to analyze the spectral properties of the WCNS-ZN-type (i.e., WCNS-ZN and WCNS-ZNM) and WCNS-P schemes. For the WCNS-ZNM scheme, two different values of the parameter θ are taken into account, i.e., θ=10 and 100. The reduced and modified wave numbers are denoted as ϕ and Φ, respectively. For comparison, the spectral properties of the classical WCNS-Z scheme with the WENO-Z weights [4] are also presented. The spectral characteristics of these schemes are shown in Figure 2. From this figure, our analysis reveals that the WCNS-Z scheme exhibits the largest dissipation, while WCNS-P demonstrates superior dissipative/dispersive properties compared to both linear WCNS and nonlinear WCNS-Z-type (i.e., WCNS-Z, WCNS-ZN, and WCNS-ZNM) schemes across medium wave numbers. The spectral characteristics of the WCNS-ZNM scheme are affected by the choice of the parameter θ, and the scheme with large θ has superior spectral characteristics as compared to the WCNS-ZN scheme. Although the WCNS-ZNM scheme with large θ produces lower numerical dissipation and dispersion, numerical oscillations around discontinuities are easily triggered. In view of spectral properties and numerical stability, parameter θ in WCNS-ZNM is set to be θ=10 for 1D cases and θ=100 for 2D cases.

    Figure 2.  ADR analysis for the WCNS-Z-type and WCNS-P schemes.

    Performance of WCNS-Z-type (including the fifth-order WCNS-Z, WCNS-ZN, and WCNS-ZNM) and WCNS-P schemes are tested by several 1D and 2D numerical examples governed by 1D and 2D Euler equations with different initial and boundary conditions, respectively. Time integration employs third-order TVD Runge-Kutta [2,15], with Δt meeting CFL conditions, e.g., for 1D cases

    ΔtCFLαx/h. (36)

    We set the CFL number to 0.45, with Δth5/3 for accuracy test cases. The numerical errors in the 1D accuracy test are calculated as

    L1=1N+1Ni=0|ρiρexacti|,L=maxi|ρiρexacti|, (37)

    where ρi,ρexacti represent the numerical and exact densities at cell nodes xi, respectively, while the numerical errors in the 2D accuracy test are measured in similar formulas.

    Example 1. Accuracy tests

    Consider the 1D Euler equations with initial data

    ρ(x,0)=1+0.2 sin (πx),u(x,0)=1,P(x,0)=1, (38)

    given by the exact solutions

    ρ(x,t)=1+0.2 sin (π(xut)),u(x,t)=1,P(x,t)=1. (39)

    The computational domain is [0,2] imposed by periodic boundary conditions. The numerical errors and convergence orders of the WCNS-ZN-type and WCNS-P schemes are calculated at time t=1 using uniform grids with different numbers of cells. The results of the WCNS-ZN and WCNS-ZNM schemes are exactly the same. Therefore, only the results of WCNS-ZNM and WCNS-P schemes are presented. Table 1 clearly demonstrates that both schemes maintain fifth-order accuracy as expected, and the numerical errors obtained with the WCNS-P scheme are less than those obtained with the WCNS-ZNM scheme on the same grids.

    Table 1.  L1 and L errors and convergence orders for the 1D case.
    Method N L error Order L1 error Order
    WCNS-ZNM




    20
    40
    80
    160
    320
    7.26E-05
    2.31E-06
    7.26E-08
    2.27E-09
    7.10E-11
    -
    4.97
    4.99
    5.00
    5.00
    4.43E-05
    1.44E-06
    4.57E-08
    1.44E-09
    4.51E-11
    -
    4.94
    4.98
    4.99
    4.99
    WCNS-P




    20
    40
    80
    160
    320
    4.36E-05
    1.50E-06
    4.81E-08
    1.51E-09
    4.73E-11
    -
    4.86
    4.97
    4.99
    5.00
    2.69E-05
    9.37E-07
    3.03E-08
    9.56E-10
    3.00E-11
    -
    4.85
    4.95
    4.98
    4.99

     | Show Table
    DownLoad: CSV

    Consider the 2D Euler equations with initial data

    ρ(x,y,0)=1+0.2 sin (π(x+y)),u(x,y,0)=1,v(x,y,0)=1,P(x,y,0)=1. (40)

    With periodic boundary conditions applied to the domain [0,2]2, the computation runs until t=0.4. Table 2 demonstrates the numerical errors and convergence orders of the WCNS-ZNM and WCNS-P schemes. From this table, it is also observed that the WCNS-ZNM/P maintains fifth-order accuracy, while the numerical errors produced by the WCNS-P scheme are smaller than those produced by the WCNS-ZNM scheme on the same grids.

    Table 2.  L1 and L errors and convergence orders for the 2D case.
    Method M × N L error Order L1 error Order
    WCNS-ZNM




    20 × 20
    40 × 40
    80 × 80
    160 × 160
    320 × 320
    5.89E-05
    1.87E-06
    5.89E-08
    1.84E-09
    5.76E-11
    -
    4.97
    4.99
    5.00
    5.00
    3.75E-05
    1.19E-06
    3.75E-08
    1.17E-09
    3.66E-11
    -
    4.97
    4.99
    5.00
    5.00
    WCNS-P




    20 × 20
    40 × 40
    80 × 80
    160 × 160
    320 × 320
    3.41E-05
    1.18E-06
    3.77E-08
    1.18E-09
    3.71E-11
    -
    4.86
    4.96
    4.99
    5.00
    2.16E-05
    7.49E-07
    2.40E-08
    7.54E-10 2.36E-11
    -
    4.85
    4.96
    4.99
    5.00

     | Show Table
    DownLoad: CSV

    The relationships between L1 errors and CPU time of the fifth-order WCNS-Z-type (i.e., WCNS-Z, WCNS-ZN, and WCNS-ZNM) and WCNS-P schemes are given in Figure 3. Evidently, for both 1D and 2D cases, the fifth-order WCNS-P scheme is computationally more efficient than the three classes of the fifth-order WCNS-Z-type schemes under the same conditions. The computational efficiency of the WCNS-ZN-type schemes is comparable to that of the WCNS-Z scheme.

    Figure 3.  Comparisons of the numerical errors versus CPU time for the WCNS-Z-type and WCNS-P schemes.

    Example 2. 1D Lax problem

    This problem is utilized to evaluate the shock-capturing performance of the presented schemes, with the following initial data

    (ρ,u,P)={(0.445,0.698,3.528),5x<0,(0.5,0,0.571),                     else. (41)

    Here, the domain [5,5] is discretized into 200 cells, with computation terminating at t=1.3. Four fifth-order schemes, WCNS-Z, WCNS-ZN, WCNS-ZNM, and WCNS-P, are selected for the calculation, with Figure 4 presenting numerical results of density distributions. A reference solution is obtained using the standard WCNS-Z scheme with 2000 grid cells. The results demonstrate that the WCNS-P scheme exhibits improvement over the WCNS-Z-type scheme, and the WCNS-ZNM scheme performs better than the WCNS-Z and WCNS-ZN schemes, with less numerical dissipation near the shock waves and stronger shock-capturing capability. Meanwhile, the WCNS-ZN scheme exhibits superior results compared to the WCNS-Z scheme, which agrees with the results in [17].

    Figure 4.  Results of 1D Lax problem at t=1.3 using 200 cells.

    Example 3. 1D Sod problem

    This problem is adopted to evaluate the stability of the presented schemes, with initial data

    (ρ,u,P)={(1.0,0.0,1.0),5x<0,(0.125,0.0,0.1),            else. (42)

    The computational domain is [5,5] discretized by 200 cells. Figure 5 demonstrates the numerical results obtained with four fifth-order schemes (WCNS-Z, WCNS-ZN, WCNS-ZNM, and WCNS-P) at time t=2. We see that the numerical behaviors of the WCNS-P and WCNS-ZN schemes perform better than the other two schemes since the discontinuities are captured more sharply. Slightly better results can be achieved by the WCNS-P scheme in comparison to the WCNS-ZNM scheme.

    Figure 5.  Results of 1D Sod problem at t=2.0 using 200 cells.

    Example 4. Osher-Shu problem

    This problem expounds the interaction of a Mach 3 shock wave with a density wave, and the initial conditions are

    (ρ,u,P)={(27/7,435/9,31/3),5x<4,(1+0.2 sin (5x),0.0,1.0),    else. (43)

    Here, the domain [5,5] is discretized into 400 cells, with computation terminating at t=1.8. Figure 6 illustrates numerical results obtained with four different schemes against a reference solution generated by the fifth-order WCNS-Z scheme using 4000 cells. We see that the WCNS-P and WCNS-ZN schemes exhibit better wave resolution than the WCNS-Z and WCNS-ZN schemes. Comparative analysis reveals the WCNS-P scheme as the least dissipative and the WCNS-Z scheme as the most dissipative among the fifth-order schemes.

    Figure 6.  Results of 1D Osher-Shu problem at t=1.8 using 400 cells.

    Example 5. 2D Riemann problem

    This example represents a classical fluid dynamics problem that involves many complex fluid structures, with initial data described as

    (ρ,u,v,P)={(1.5,0,0,1.5),x0.8,y0.8,(0.5323,1.206,0,0.3),   x<0.8,y0.8,(0.138,1.206,1.206,0.029),x<0.8,y0.8,(0.5323,0,1.206,0.3),x0.8,y0.8. (44)

    The domain [0,1]2 with transmissive boundaries is divided into 400 × 400 cells. Four fifth-order schemes, WCNS-Z, WCNS-ZN, WCNS-ZNM, and WCNS-P, are applied for the computation. Figure 7 plots 35 density contours within the interval of [0.14,1.76] at time t=0.8. Numerical results show that fine-scale structures in the regions of strong velocity gradient are captured by the WCNS-ZNM and WCNS-P schemes more remarkably compared to the WCNS-Z and WCNS-ZN schemes. The WCNS-ZN scheme has a distinct advantage in resolution over the WCNS-Z scheme, which cannot resolve the Kelvin-Helmholtz (KH) instabilities due to the nonlinear shock-capturing mechanism. Compared to the WCNS-ZNM scheme, the WCNS-P scheme can capture much finer small-scale vortex structures incurred by the KH instabilities along the jet's neck and has obvious advantages in dissipation and resolution.

    Figure 7.  Results of 2D Riemann problem at t=0.8 using 400 × 400 cells.

    Example 6. 2D Rayleigh–Taylor instability problem

    The example characterizes interfacial instability between density-stratified fluids under gravitational acceleration, widely adopted for assessing high-order schemes' dissipation and resolution capabilities. The initialization is defined as

    (ρ,u,v,P)={(2.0,0.0,0.025α  cos (8πx),2y+1), 0y<0.5,(1.0,0.0,0.025α  cos (8πx),y+1.5),       else. (45)

    In the above equations, the sound speed is α=γP/ρ with γ=5/3 in this case. The domain [0,0.25]×[0,1] is divided into 240 × 960 cells, with Dirichlet boundary conditions (ρ,u,v,P)=(2.0,0.0,0.0,1.0) at the bottom and (ρ,u,v,P)=(1.0,0.0,0.0,2.5) at the top, and the left and right boundaries are set to reflective boundary conditions. The final time is t=1.95. The fifth-order WCNS-Z, WCNS-ZN, WCNS-ZNM, and WCNS-P schemes are selected for the calculation, and 15 density contours, from 0.85 to 2.25, are shown in Figure 8. We see that the mushroom-cap shape spike is resolved by all of the four schemes. In comparison with the WCNS-Z scheme, the WCNS-ZN scheme has improved resolution but loses the symmetry incurred by the roundoff errors. The WCNS-ZNM and WCNS-P schemes outperform the two other schemes in resolution. Notably, the WCNS-P scheme captures finer curl structures and more small-scale vortices than the WCNS-ZNM scheme, with slightly lower numerical dissipation and better resolution.

    Figure 8.  Results of 2D Rayleigh–Taylor instability problem at t=1.95 using 240 × 960 cells.

    Example 7. 2D Kelvin–Helmholtz instability problem

    This test case is commonly applied to validate the capability of high-order schemes in evolving small shear layers into 2D turbulence, expressed as:

    (ρ,u,v,P)={(2.0, 0.5, 0.01 sin (2πx), 2.5), |y|<0.25,(1.0, 0.5, 0.01 sin (2πx), 2.5),         else. (46)

    We use 400 × 400 cells to discretize the domain [0.5,0.5]2 with periodic boundary conditions, simulating until t=1. Density contour distributions (50 levels, from 0.8 to 2.26) are presented in Figure 9, which shows that all of the four fifth-order schemes (i.e., WCNS-Z, WCNS-ZN, WCNS-ZNM, and WCNS-P) are capable of capturing small-scale vortex structures caused by the KH instabilities, generating 2D turbulence with time. The size of vortices captured by the WCNS-ZNM and WCNS-P schemes is much bigger than that captured by the WCNS-Z and WCNS-ZN schemes. Among the four schemes, the WCNS-P scheme resolves vortices with the biggest size, exhibiting the best dissipation characteristics. Compared to the WCNS-Z scheme with the largest numerical dissipation, the WCNS-ZN scheme has better resolution, though still inferior to the WCNS-ZNM scheme.

    Figure 9.  Results of 2D Kelvin–Helmholtz instability problem at t=1 using 400 × 400 cells.

    Example 8. 2D explosion problem

    The explosion problem evaluates the stability of numerical methods to evolve unstable contacts, with cylindrically symmetric initial conditions

    (ρ,u,v,P)={(1,0,0,1),x2+y2<0.16,(0.125,0,0,0.1),otherwise. (47)

    The domain [0,1.5]2 is divided into 400 × 400 cells. Wall boundary conditions are imposed at the left and bottom boundaries; otherwise, free boundary conditions are imposed. Figure 10 shows the distribution results of 45 density contours in the interval [0.1,0.2] at time t=3.2. The numerical results indicate that the WCNS-ZNM and WCNS-P schemes behave better than the WCNS-Z and WCNS-ZN schemes, since they can capture the curly contact curve and small-scale vortex structures of finer scales. The WCNS-ZNM scheme has better resolution compared to the WCNS-ZN scheme. Among the four fifth-order schemes, the WCNS-Z scheme with excessive numerical dissipation resolves the narrowest mixing layer and the WCNS-P scheme resolves the curliest contact curve.

    Figure 10.  Results of 2D explosion problem at t=3.2 using 400 × 400 cells.

    Example 9. 2D double Mach reflection problem

    The shock wave structure in this example is extremely complex, including multiple reflected shock waves and density discontinuity lines. We use this case to test the performance of the presented schemes in capturing shocks and small-scale vortices. Initially, the pre- and post-shocked values are given by

    (ρ,u,v,P)={(8.0,8.25  cos (π6),8.25 sin (π6),116.5),0<x16+y3,(1.4,0.0,0.0,1.0),                      16+y3<x<1. (48)

    The domain is [0,4]×[0,1] with 960 × 240 cells. Inflow/outflow conditions are set at left/right boundaries. The top boundary uses exact shock motion, while the bottom applies reflective (if x>1/6) or transmissive conditions. We choose four schemes (i.e., WCNS-Z, WCNS-ZN, WCNS-ZNM, and WCNS-P) for calculation. Figure 11 displays the distribution results of 40 density contours in the interval [1.5,21.53] at time t=0.2. We observe that all schemes are capable of providing stable results and capturing the main flow features. The WCNS-ZNM and WCNS-P schemes better resolve the small-scale vortices induced by KH instabilities in the vortex rolling-up region, while they are partly dissipated by using the WCNS-Z scheme. The WCNS-ZN scheme outperforms the WCNS-Z scheme but underperforms compared to the WCNS-ZNM scheme. The WCNS-P scheme exhibits finer vortex detail at slip-line endpoints and in wall-jet zones compared to other fifth-order variants.

    Figure 11.  Results of 2D double Mach reflection problem at t=0.2 using 960 × 240 cells.

    Example 10. 2D forward-facing step problem

    This problem validates the resolution and stability of high-order schemes. The wind tunnel is set in the domain [0,0.6]×[0,1][0.6,3]×[0.2,1]. The left and right boundaries are treated as inflow and outflow conditions, respectively, while all other surfaces are treated as reflective walls. The initial conditions are given by (ρ,u,v,P)=(1.4,3,0,1). The domain is [0,3]×[0,1] with 900 × 300 cells. Density contour distributions (50 levels, from 0.3 to 6.68) are presented in Figure 12. From this figure, we can see that all of the four fifth-order schemes (i.e., WCNS-Z, WCNS-ZN, WCNS-ZNM, and WCNS-P) capture the shocks robustly. Among the four schemes, the WCNS-P scheme captures the most rolling-up structures of the vortex sheet emerging from the triple junction and triggered by the KH instabilities, with the least dissipation and the highest resolution. While WCNS-ZN outperforms WCNS-Z in resolution, it remains inferior to WCNS-ZNM.

    Figure 12.  Results of 2D forward-facing step problem at t=4 using 900 × 300 cells.

    This study developed the enhanced WCNS-ZNM and WCNS-P schemes for the Euler equations. The WCNS-ZNM scheme is composed of the sixth-order hybrid central difference scheme for the flux derivatives and the modified characteristic-wise WENO-ZN-type interpolation technique for numerical fluxes at cell edges. The WCNS-P scheme was constructed by using a linear combination of three at least third-order non-polynomial interpolations with an adjustable parameter. By optimizing this parameter, the local numerical errors can be reduced, and thus, accuracy is improved. To maintain the non-oscillating property of the WCNS-P scheme, we apply the monotone polynomial interpolation approach to detect the non-smooth regions where the WCNS-ZNM scheme is initiated and accordingly realizes the adaptive switching between the two schemes in smooth and non-smooth regions. The ADR analysis indicates that the WCNS-P scheme has the best spectral properties, compared to the fifth-order WCNS-Z-type (i.e., WCNS-Z, WCNS-ZN, and WCNS-ZNM) schemes. Experimental results of accuracy tests show that the WCNS-ZNM and WCNS-P schemes are able to achieve the designed fifth-order accuracy in smooth regions, and the WCNS-P scheme has higher computational efficiency, compared to the three classes of WCNS-Z-type schemes of the same order. 1D and 2D numerical examples verify that the WCNS-ZN scheme shows clear improvement over the WCNS-Z scheme, and the WCNS-ZNM and WCNS-P schemes behave better than the WCNS-ZN scheme, as far as the resolution and shock-capturing capability are concerned. In general, the WCNS-P scheme provides improved numerical results and has better features with respect to high resolution in comparison to the WCNS-ZNM scheme.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The authors declare there is no conflict of interest.



    [1] X. D. Liu, S. Osher, T. Chan, Weighted essentially non-oscillatory schemes, J. Comput. Phys., 115 (1994), 200–212. https://doi.org/10.1006/jcph.1994.1187 doi: 10.1006/jcph.1994.1187
    [2] G. S. Jiang, C. W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys., 126 (1996), 202–228. https://doi.org/10.1006/jcph.1996.0130 doi: 10.1006/jcph.1996.0130
    [3] A. K. Henrick, T. D. Aslam, J. M. Powers, Mapped weighted essentially non-oscillatory schemes: Achieving optimal order near critical points, J. Comput. Phys., 207 (2005), 542–567. https://doi.org/10.1016/j.jcp.2005.01.023 doi: 10.1016/j.jcp.2005.01.023
    [4] R. Borges, M. Carmona, B. Costa, W. S. Don, An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws, J. Comput. Phys., 227 (2008), 3191–3211. https://doi.org/10.1016/j.jcp.2007.11.038 doi: 10.1016/j.jcp.2007.11.038
    [5] J. Qiu, C. W. Shu, Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method: One-dimensional case, J. Comput. Phys., 193 (2004), 115–135. https://doi.org/10.1016/j.jcp.2003.07.026 doi: 10.1016/j.jcp.2003.07.026
    [6] J. Qiu, C. W. Shu, Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method Ⅱ: Two-dimensional case, Comput. Fluids, 34 (2005), 642–663. https://doi.org/10.1016/j.compfluid.2004.05.005 doi: 10.1016/j.compfluid.2004.05.005
    [7] J. Y. Guo, J. H. Jung, Radial basis function ENO and WENO finite difference methods based on the optimization of shape parameters, J. Sci. Comput., 70 (2017), 551–575. https://doi.org/10.1007/s10915-016-0257-y doi: 10.1007/s10915-016-0257-y
    [8] J. Y. Guo, J. H. Jung, A RBF-WENO finite volume method for hyperbolic conservation laws with the monotone polynomial interpolation method, Appl. Numer. Math., 112 (2017), 27–50. https://doi.org/10.1016/j.apnum.2016.10.003 doi: 10.1016/j.apnum.2016.10.003
    [9] B. Jeong, H. Yang, J. Yoon, Development of a WENO scheme based on radial basis function with an improved convergence order, J. Comput. Phys., 468 (2022), 111502. https://doi.org/10.1016/j.jcp.2022.111502 doi: 10.1016/j.jcp.2022.111502
    [10] T. Yang, G. Q. Zhao, Q. J. Zhao, Novel TENO schemes with improved accuracy order based on perturbed polynomial reconstruction, J. Comput. Phys., 488 (2023), 112219. https://doi.org/10.1016/j.jcp.2023.112219 doi: 10.1016/j.jcp.2023.112219
    [11] T. Yang, D. Z. Sun, Q. J. Zhao, G. Q. Zhao, X. Chen, Perturbed polynomial reconstructed seventh-order hybrid WENO scheme with improved accuracy and resolution, Comput. Fluids, 289 (2025), 106548. https://doi.org/10.1016/j.compfluid.2025.106548 doi: 10.1016/j.compfluid.2025.106548
    [12] H. M. Zuo, J. Zhu, Increasingly high-order hybrid multi-resolution WENO schemes in multi-dimensions, J. Comput. Phys., 514 (2024), 113233. https://doi.org/10.1016/j.jcp.2024.1132 33 doi: 10.1016/j.jcp.2024.113233
    [13] H. M. Zuo, J. Zhu, A new type of modified MR-WENO schemes with new troubled cell indicators for solving hyperbolic conservation laws in multi-dimensions, J. Comput. Phys., 508 (2024), 112996. https://doi.org/10.1016/j.jcp.2024.112996 doi: 10.1016/j.jcp.2024.112996
    [14] X. Y. Hu, Q. Wang, N. A. Adams, An adaptive central-upwind weighted essentially non-oscillatory scheme, J. Comput. Phys., 229 (2010), 8952–8965. https://doi.org/10.1016/j.jcp.2010.08.019 doi: 10.1016/j.jcp.2010.08.019
    [15] J. Zhu, J. X. Qiu, A new fifth order finite difference WENO scheme for solving hyperbolic conservation laws, J. Comput. Phys., 318 (2016), 110–121. https://doi.org/10.1016/j.jcp.2016.05.010 doi: 10.1016/j.jcp.2016.05.010
    [16] Y. Shen, K. Zhang, S. Li, J. Peng, A novel method for constructing high accurate and robust WENO-Z type scheme, preprint, arXiv: 2004.07954.
    [17] S. Y. Li, Y. Q. Shen, K. Zhang, M. Yu, High order weighted essentially non-oscillatory WENO-ZN schemes for hyperbolic conservation laws, Comput. Fluids, 244 (2022), 105547. https://doi.org/10.1016/j.compfluid.2022.105547 doi: 10.1016/j.compfluid.2022.105547
    [18] X. G. Deng, H. X. Zhang, Developing high-order weighted compact nonlinear schemes, J. Comput. Phys., 165 (2000), 22–44. https://doi.org/10.1006/jcph.2000.6594 doi: 10.1006/jcph.2000.6594
    [19] T. Nonomura, K. Fujii, Robust explicit formulation of weighted compact nonlinear scheme, Comput. Fluids, 85 (2013), 8–18. https://doi.org/10.1016/j.compfluid.2012.09.001 doi: 10.1016/j.compfluid.2012.09.001
    [20] X. L. Liu, S. H. Zhang, H. X. Zhang, C. W. Shu, A new class of central compact schemes with spectral-like resolution Ⅱ: Hybrid weighted nonlinear schemes, J. Comput. Phys., 284 (2015), 133–154. https://doi.org/10.1016/j.jcp.2014.12.027 doi: 10.1016/j.jcp.2014.12.027
    [21] Z. G. Yan, H. Y. Liu, Y. K. Ma, M. L. Mao, X. G. Deng, Further improvement of weighted compact nonlinear scheme using compact nonlinear interpolation, Comput. Fluids, 156 (2017), 135–145. https://doi.org/10.1016/j.compfluid.2017.06.028 doi: 10.1016/j.compfluid.2017.06.028
    [22] M. L. Wong, S. K. Lele, High-order localized dissipation weighted compact nonlinear scheme for shock- and interface-capturing in compressible flows, J. Comput. Phys., 339 (2017), 179–209. https://doi.org/10.1016/j.jcp.2017.03.008 doi: 10.1016/j.jcp.2017.03.008
    [23] Y. Q. Jiang, S. G. Zhou, X. Zhang, Y. G. Hu, High order all-speed semi-implicit weighted compact nonlinear scheme for the isentropic Navier-Stokes equations, J. Comput. Appl. Math., 411 (2022), 114272. https://doi.org/10.1016/j.cam.2022.114272 doi: 10.1016/j.cam.2022.114272
    [24] Y. Q. Jiang, S. G. Zhou, X. Zhang, Y. G. Hu, High-order weighted compact nonlinear scheme for one- and two-dimensional Hamilton-Jacobi equations, Appl. Numer. Math., 171 (2022), 353–368. https://doi.org/10.1016/j.apnum.2021.09.012 doi: 10.1016/j.apnum.2021.09.012
    [25] X. Q. Huang, Y. Q. Jiang, H. H. Yang, Two classes of third-order weighted compact nonlinear schemes for Hamilton-Jacobi equations, Appl. Math. Comput., 469 (2024), 128554. https://doi.org/10.1016/j.amc.2024.128554 doi: 10.1016/j.amc.2024.128554
    [26] S. Pirozzoli, On the spectral properties of shock-capturing schemes, J. Comput. Phys., 219 (2006), 489–497. https://doi.org/10.1016/j.jcp.2006.07.009 doi: 10.1016/j.jcp.2006.07.009
  • Reader Comments
  • © 2025 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(184) PDF downloads(24) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog