Finite Difference WENO Scheme Research Articles

Numerical simulation of 3D seismic wave based on alternative flux finite-difference WENO scheme

Summary High-frequency non-physical oscillations may occur due to shock waves in seismic wavefield and dynamic rupture simulation. In this study, we introduced the alternative flux finite-difference weighted essentially non-oscillatory (AWENO) scheme to address potential shock wave issues in computational seismology effectively. The wavefield of the body-fitted curvilinear domain was accurately computed through conservative grid mapping, ensuring accurate implementation of free surface boundary conditions on irregular surfaces using characteristic boundary conditions and minimizing artificial boundary reflections with exponential decay absorbing layers. Finally, we compared our scheme with the GRTM for flat surfaces and the CGFDM3D-EQR for irregular surfaces to demonstrate its correctness and accuracy, and validated its non-oscillatory characteristics. The aforementioned scheme is anticipated to assume a significant function in simulating more intricate seismic wavefields or dynamic ruptures.

Local Characteristic Decomposition–Free High-Order Finite Difference WENO Schemes for Hyperbolic Systems Endowed with a Coordinate System of Riemann Invariants

Local Characteristic Decomposition–Free High-Order Finite Difference WENO Schemes for Hyperbolic Systems Endowed with a Coordinate System of Riemann Invariants

High Order Structure-Preserving Finite Difference WENO Schemes for MHD Equations with Gravitation in all Sonic Mach Numbers

High Order Structure-Preserving Finite Difference WENO Schemes for MHD Equations with Gravitation in all Sonic Mach Numbers

A new type of increasingly higher order finite difference and finite volume MR-WENO schemes with adaptive linear weights for hyperbolic conservation laws

In this paper, the new fifth-order, seventh-order, and ninth-order finite difference and finite volume multi-resolution weighted essentially non-oscillatory (MR-WENO) schemes with adaptive linear weights are presented for hyperbolic conservation laws on structured meshes. They are termed as high-order finite difference and finite volume ALW-WENO schemes. These ALW-WENO schemes only apply one small stencil and one large stencil in reconstruction processes, which could achieve the desired accuracy in the region of smoothness and non-oscillatory properties in the region of containing strong shocks. The linear weights that sum to one can be automatically adjusted to any positive numbers. This is the first time that arbitrary high-order finite difference and finite volume WENO schemes are designed by using only two unequal-sized central spatial stencils. The structure of these novel WENO schemes is simple, so it is easier for obtaining high-order accuracy and solving multi-dimensional problems in large scale engineering applications. Compared to traditional MR-WENO schemes with same order, the computational efficiency can be further improved. Some benchmark tests indicate that these new ALW-WENO schemes have good robustness and performance.

A sixth-order finite difference WENO scheme for non-linear degenerate parabolic equations: an alternative technique

PurposeThis paper aims to construct a sixth-order weighted essentially nonoscillatory scheme for simulating the nonlinear degenerate parabolic equations in a finite difference framework.Design/methodology/approachTo design this scheme, we approximate the second derivative in these equations in a different way, which of course is still in a conservative form. In this way, unlike the common practice of reconstruction, the approximation of the derivatives of odd order is needed to develop the numerical flux.FindingsThe results obtained by the new scheme produce less error compared to the results of other schemes in the literature that are recently developed for the nonlinear degenerate parabolic equations while requiring less computational times.Originality/valueThis research develops a new weighted essentially nonoscillatory scheme for solving the nonlinear degenerate parabolic equations in multidimensional space. Besides, any selection of the constants (sum equals one is the only requirement for them), named the linear weights, will obtain the desired accuracy.

A New Third-Order Finite Difference WENO Scheme to Improve Convergence Rate at Critical Points

In this work, a new, improved third-order finite difference weighted essentially non-oscillatory scheme is presented for one- and two-dimensional hyperbolic conservation laws and associated problems. The parameter p which is regulate dissipation is introduced in the nonlinear weights in the framework of the conventional WENO-Z scheme, and the higher-order global smoothness indicator is obtained by the idea of Wang [Wang, Y. H., Y. L. Du, K. L. Zhao and L. Yuan. 2020. ‘A Low-dissipation Third-order Weighted Essentially Nonoscillatory Scheme with a New Reference Smoothness Indicator’. International Journal for Numerical Methods in Fluids. 92 (9): 1212–1234.], the sufficient condition of nonlinear weights is proved by using Taylor expansions. Finally, the value range of parameter p is obtained. The proposed scheme is verified to achieve the optimal order near critical points by linear convection equations with different initial values, and the high-resolution characteristic of the present scheme is proved on a variety of one- and two- dimensional standard numerical examples. Numerical results demonstrate that the proposed scheme gives better performance in comparison with the other third-order WENO schemes.

Finite difference alternative unequal-sized weighted essentially non-oscillatory schemes for hyperbolic conservation laws

In this paper, two unequal-sized weighted essentially non-oscillatory (US-WENO) schemes are proposed for solving hyperbolic conservation laws. First, an alternative US-WENO (AUS-WENO) scheme based directly on the values of conserved variables at the grid points is designed. This scheme can inherit all the advantages of the original US-WENO scheme [J. Zhu and J. Qiu, “A new fifth order finite difference WENO scheme for solving hyperbolic conservation laws,” J. Comput. Phys. 318, 110–121 (2016).], such as the arbitrariness of the linear weights. Moreover, this presented AUS-WENO scheme enables any monotone fluxes applicable to this framework, whereas the original US-WENO scheme is only suitable for the more dissipative smooth flux splitting. Therefore, the method in this paper has a smaller L1 and L∞ numerical errors than the original scheme under the same conditions. Second, in order to further improve the computational efficiency of the above AUS-WENO scheme, a hybrid AUS-WENO scheme is proposed by combining a hybrid strategy. This strategy identifies the discontinuous regions directly based on the extreme points of the reconstruction polynomial corresponding to the five-point stencil, which brings the important advantage that it does not depend on the specific problem and does not contain any artificial adjustable parameters. Finally, the performance of the above two AUS-WENO schemes in terms of low dissipation, shock capture capability, discontinuity detection capability, and computational efficiency is verified by some benchmark one- and two-dimensional numerical examples.

A fifth-order finite difference HWENO scheme combined with limiter for hyperbolic conservation laws

In this paper, a simple fifth-order finite difference Hermite WENO (HWENO) scheme combined with limiter is proposed for one- and two-dimensional hyperbolic conservation laws. The fluxes in the governing equation are approximated by the nonlinear HWENO reconstruction which is the combination of a quintic polynomial with two quadratic polynomials, where the linear weights can be artificial positive numbers only if the sum equals one. And other fluxes in the derivative equations are approximated by high-degree polynomials directly. For the purpose of controlling spurious oscillations, an HWENO limiter is applied to modify the derivatives. Instead of using the modified derivatives both in fluxes reconstruction and time discretization as in the modified HWENO scheme (Z. Zhao et al. (2020) [27]), we only apply the modified derivatives in time discretization while remaining the original derivatives in fluxes reconstruction. Comparing with the modified HWENO scheme, the proposed HWENO scheme is simpler, more accurate, efficient, and has better resolution. In addition, the HWENO scheme has a more compact spatial reconstructed stencil and is more efficient than the same order finite difference WENO scheme of Jiang and Shu (1996) [11] and WENO scheme with artificial linear weights of Zhu and Qiu (2016) [29]. Various benchmark numerical examples are presented to show the fifth-order accuracy, efficiency, high resolution, and robustness of the proposed HWENO scheme.

A new finite difference mapped unequal-sized WENO scheme for Hamilton-Jacobi equations

In this paper, a new fifth-order finite difference mapped unequal-sized weighted essentially non-oscillatory (MUS-WENO) scheme is proposed for solving Hamilton-Jacobi equations in multi-dimensions. This new MUS-WENO scheme uses the same widest spatial stencils as that of the classical same order finite difference WENO scheme [19], and could obtain smaller truncation errors and optimal fifth-order accuracy in smooth regions. The MUS-WENO scheme uses a convex combination of a quartic polynomial with two quadratic polynomials. The linear weights can be artificially set as any positive constants on condition that their summation is one. Together with a new mapping function, associated new mapped nonlinear weights are proposed to reduce the difference with the linear weights. So it is a first time that the new finite difference MUS-WENO scheme with a very small ε could reach optimal accuracy order convergence even near extreme points in smooth regions and avoid spurious oscillations near strong discontinuities when solving Hamilton-Jacobi equations. Generally speaking, the main advantages of such new MUS-WENO scheme comparing with the classical WENO scheme [19] are its good convergence, robustness, efficiency and easy extension to multi-dimensions. Numerical experiments are proposed to show the good performance of this new MUS-WENO scheme.

High order well-balanced asymptotic preserving finite difference WENO schemes for the shallow water equations in all Froude numbers

In this paper, high order semi-implicit well-balanced and asymptotic preserving finite difference WENO schemes are proposed for the shallow water equations with a non-flat bottom topography. We consider the Froude number ranging from O(1) to 0, which in the zero Froude limit becomes the “lake equations” for balanced flow without gravity waves. We apply a well-balanced finite difference WENO reconstruction, coupled with a stiffly accurate implicit-explicit (IMEX) Runge-Kutta time discretization. The resulting semi-implicit scheme can be shown to be well-balanced, asymptotic preserving (AP) and asymptotically accurate (AA) at the same time. Both one- and two-dimensional numerical results are provided to demonstrate the high order accuracy, AP property and good performance of the proposed methods in capturing small perturbations of steady state solutions.

High Order Semi-implicit WENO Schemes for All-Mach Full Euler System of Gas Dynamics

In this paper, we propose high order semi-implicit schemes for the all Mach full Euler equations of gas dynamics. Material waves are treated explicitly, while acoustic waves are treated implicitly, thus avoiding severe CFL restrictions for low Mach flows. High order accuracy in time is obtained by semi-implicit temporal integrator based on the IMEX Runge-Kutta (IMEX-RK) framework. High order in space is achieved by finite difference WENO schemes with characteristic-wise reconstructions adapted to the semi-implicit IMEX-RK time discretization. The schemes are proven to be asymptotic preserving and asymptotically accurate as the Mach number vanishes. Besides, they can well capture discontinuous solutions in the compressible regime, especially for two dimensional Riemann problems. Numerical tests in one and two space dimensions will illustrate the effectiveness of the proposed schemes.

Numerical study of shock-disturbances interaction in hypersonic inviscid flows with real gas effects using high-order WENO scheme

In the present study, the shock-disturbances interaction in hypersonic inviscid flows with real gas effects is studied by applying a high-order accurate numerical method with the shock capturing technique. To consider real gas effects, the equilibrium air model is utilized here. The strong conservative form of the unsteady compressible Euler equations in the 2D generalized curvilinear coordinates is formulated and the resulting system of equations for the equilibrium air model is discretized by using the fifth-order finite-difference WENO scheme in space and the explicit third-order TVD Runge–Kutta scheme in time to provide a highly accurate and robust equilibrium airflow solver. The solution methodology adopted based on the high-order WENO scheme with the shock capturing technique can reasonably simulate the strong discontinuities in the solution domain without needing any numerical artifact and this feature makes the solution method suitable for the shock-disturbances interaction in hypersonic flows in which real gas effects are to be considered. To study the physical phenomenon of the shock-disturbances interaction problem in high-speed inviscid equilibrium airflows, two test cases are simulated by applying the solution methodology adopted. In the first test case, the interaction between the incoming fast acoustic wave with the bow shock formed around a cylinder with the nose radius of RN=0.0381 m at M∞=8.03 and T∞=182.333 K is simulated. In the second test case, the cylinder radius is taken RN=2.54×10−2 m and the upstream velocity and temperature are u∞=5590m/s and T∞=1833 K, respectively, and all kinds of free stream disturbances, including the entropy wave and the fast and slow acoustic waves, for different wave numbers are considered. The computations are carried out for both the perfect gas and equilibrium air models and the effects of real gas on both the mean and disturbances fields are also studied. The present study introduces the first-known numerical investigation of the shock-disturbances interaction in hypersonic flow over blunt noses with real gas effects in which the numerical results obtained are thoroughly validated by the analytical/theoretical ones and good agreement is found.

Simple high order well-balanced finite difference WENO schemes for the Euler equations under gravitational fields

The compressible Euler equations coupled with the gravitational source terms admit a hydrostatic equilibrium state where the gradients of the flux terms can be exactly balanced by those in the source terms. This property of exact preservation of the equilibrium is highly desirable at the discrete level when they are numerically solved. In this study, we design the simple high order well-balanced finite difference weighted essentially non-oscillatory (WENO-R/I) schemes, which base on the WENO reconstruction and interpolation procedures respectively, for this system with an ideal gas equation of state and a prescribed gravitational field. The main idea in achieving the well-balanced property is to rewrite the source terms in a special way and discretize them by using the nonlinear WENO differential operators with the homogenization condition to those for the flux terms. The proposed well-balanced schemes can be proved mathematically to preserve the hydrostatic isothermal and polytropic equilibria states exactly and at the same time maintain genuine high order accuracy. Moreover, the resulting schemes are high efficiency of computation and can be implemented straightforwardly into the existing finite difference WENO code with the Lax-Friedrichs numerical flux for solving the compressible Euler equations. Last but not least, the proposed well-balanced framework also works for a class of moving equilibrium state. Extensive one- and two-dimensional numerical examples are carried out to investigate the performance of the proposed schemes in term of high order accuracy, well-balanced property, shock capturing essentially non-oscillatory and resolving the small perturbation on the coarse mesh resolution.

A Modified Fifth Order Finite Difference Hermite WENO Scheme for Hyperbolic Conservation Laws

In this paper, we develop a modified fifth order accuracy finite difference Hermite WENO (HWENO) scheme for solving hyperbolic conservation laws. The main idea is that we first modify the derivatives of the solution by Hermite WENO interpolations, then we discretize the original and derivative equations in the spatial directions by the same approximation polynomials. Comparing with the original finite difference HWENO scheme of Liu and Qiu (J Sci Comput 63:548–572, 2015), one of the advantages is that the modified HWENO scheme is more robust than the original one since we do not need to use the additional positivity-preserving flux limiter methodology, and larger CFL number can be applied. Another advantage is that higher order numerical accuracy than the original scheme can be achieved for two-dimensional problems under the condition of using the same approximation stencil and information. Furthermore, the modified scheme preserves the nice property of compactness shared by HWENO schemes, i.e., only immediate neighbor information is needed in the reconstruction, and it has smaller numerical errors and higher resolution than the classical fifth order finite difference WENO scheme of Jiang and Shu (J Comput Phys 126:202–228, 1996). Various benchmark numerical tests of both one-dimensional and two-dimensional problems are presented to illustrate the numerical accuracy, high resolution and robustness of the proposed novel HWENO scheme.

Kinetic theory based multi-level adaptive finite difference WENO schemes for compressible Euler equations

In this paper we proposed the kinetic framework based fifth-order adaptive finite difference WENO schemes abbreviated as WENO-AO-K schemes to solve the compressible Euler equations, which are quasi-linear hyperbolic equations that can admit discontinuous solutions like shock and contact waves. The formulation of the proposed schemes is based on the kinetic theory where one can recover the Euler equations by applying a suitable moment method strategy to the Boltzmann equation. The kinetic flux vector splitting strategy is used in WENO-AO framework, which produces the computationally expensive error and exponential functions. Thus, to reduce the computational cost, a physically more relevant peculiar velocity based splitting strategy is used, which is more efficient than the kinetic flux vector splitting. High order of accuracy in time is achieved using the third-order total variation diminishing Runge–Kutta (TVD-RK) scheme. Several one- and two-dimensional test cases are solved for the compressible Euler equations using the proposed fifth-order WENO-AO-K schemes and the results are compared with conventional WENO-AO scheme. Proposed schemes capture the complex flow features in a smooth region accurately, and discontinuity is also well resolved. Error analysis of the proposed schemes shows optimal convergence rates in various norms.

Effect of Flux Evaluation Methods on the Resolution and Robustness of the Two-Step Finite-Difference WENO Scheme

Effect of Flux Evaluation Methods on the Resolution and Robustness of the Two-Step Finite-Difference WENO Scheme

On the conservation of finite difference WENO schemes in non-rectangular domains using the inverse Lax-Wendroff boundary treatments

We discuss the issue of conservation of the total mass for finite difference WENO schemes solving hyperbolic conservation laws on a Cartesian mesh using the inverse Lax-Wendroff boundary treatments in arbitrary physical domains whose boundaries do not coincide with grid lines. The numerical fluxes near the boundary are suitably modified so that strict conservation of the total mass is achieved and the high order accuracy and non-oscillatory performance are not compromised. The key point is a suitable definition of the total mass, which is consistent with the high order accuracy finite difference framework over an arbitrary domain with a boundary not necessarily coinciding with grid lines. Extensive numerical examples are provided to demonstrate that our modified method is strictly conservative, and is high order accurate and has as good performance as the original high order WENO schemes with the Lax-Wendroff boundary treatments, for both smooth problems and problems with discontinuities, in both one- and two-dimensional problems involving both scalar equations and systems.

Positivity-Preserving Finite Difference WENO Scheme for Ten-Moment Equations with Source Term

We develop a positivity-preserving finite difference WENO scheme for the Ten-Moment equations with body forces acting as a source in the momentum and energy equations. A positive forward Euler scheme under a CFL condition is first constructed which is combined with an operator splitting approach together with an integrating factor, strong stability preserving Runge–Kutta scheme. The positivity of the forward Euler scheme is obtained under a CFL condition by using a scaling type limiter, while the solution of the source operator is performed exactly and is positive without any restriction on the time step. The proposed method can be used with any WENO reconstruction scheme and we demonstrate it with fifth order accurate WENO-JS, WENO-Z and WENO-AO schemes. An adaptive CFL strategy is developed which can be more efficient than the use of reduced CFL for positivity preservation. Numerical results show that high order accuracy and positivity preservation are achieved on a range of test problems.

Improvement of third-order finite difference WENO scheme at critical points

ABSTRACTIn the paper, an improved third-order finite difference weighted essentially non-oscillatory scheme is presented for achieving the optimal order near critical points. A new global smoothness indicator is obtained with the way of Taylor expansion for the local smoothness indicator. In the framework of the conventional WENO-Z scheme, we present an improved third-order finite difference WENO scheme. The proposed scheme is verified to achieve third-order accuracy by several benchmarks and the behaviour of the present scheme is proved on a variety of one- and two- dimensional standard numerical examples. Numerical results demonstrate that the proposed scheme gives better performance in comparison with the other third-order WENO schemes.

An efficient class of WENO schemes with adaptive order for unstructured meshes

Recent advances in finite-difference WENO schemes for hyperbolic conservation laws have resulted in WENO schemes with adaptive order of accuracy. For instance, a WENO-AO(5,3) scheme can provide up to fifth order of accuracy when the smoothness of the solution in the fifth order stencil warrants it, and yet, it can adaptively drop down to third order of accuracy when the higher order is not warranted by the solution on the mesh. Having an analogous capability for finite-volume WENO schemes for hyperbolic conservation laws, especially on unstructured meshes, can be very valuable. The present paper documents the design of finite volume WENO-AO(4,3) and WENO-AO(5,3) schemes for unstructured meshes. As with WENO-AO for structured meshes, the key advance lies in realizing that there is a favorable basis set, which is very easily constructed, and in which the computation is dramatically simplified. As with finite-difference WENO, we realize that one can make a non-linear hybridization between a large, centered, very high accuracy stencil and a lower order central WENO scheme that is, nevertheless, very stable and capable of capturing physically meaningful extrema. This yields a class of adaptive order WENO schemes that work well on unstructured meshes. On both the large and small stencils we have been able to make the stencil evaluation step very efficient owing to the choice of a favorable Taylor series basis. By extending the Parallel Axis Theorem, we show that there is a significant simplification in the finite volume reconstruction. Instead of solving a constrained least squares problem, our method only requires the solution of a smaller least squares problem on each stencil. This also simplifies the matrix assembly and solution for each stencil. The evaluation of smoothness indicators is also simplified. Accuracy tests show that the method meets its design accuracy. Several stringent test problems are presented to demonstrate that the method works very robustly and very well. The test problems are chosen to show that our method can be applied to many different meshes that are used to map geometric complexity or solution complexity.