WENO Scheme Research Articles

A Moment-Based Hermite WENO Scheme with Unified Stencils for Hyperbolic Conservation Laws

A Moment-Based Hermite WENO Scheme with Unified Stencils for Hyperbolic Conservation Laws

An Improved Fifth-Order WENO Scheme for Solving Hyperbolic Conservation Laws Near Critical Points

An Improved Fifth-Order WENO Scheme for Solving Hyperbolic Conservation Laws Near Critical Points

Numerical study of bifurcation structures in reflected shock-wave/laminar-boundary-layer interaction within an end-wall tube

Abstract The primary aim of this study is to analyze the unsteady characteristics of the interaction between a reflected shock wave and a laminar boundary layer in an end-wall shock tube. Our direct numerical simulations at shock Mach numbers of Ms=1.9, 2.5, and 3.5 using a fifth-order WENO scheme and three-step Runge-Kutta time integration method revealed inhomogeneity and anisotropy in the shock bifurcation. Surprisingly, the upper and lower bifurcated structures maintain a notably asymmetric flow during the forward propagation of the reflected shock bifurcation. The inverse flow in the bifurcation resembles a crooked earthworm structure, exhibiting high-frequency oscillations indicative of instability. However, at higher shock intensities, the earthworm transforms into a stable strip-like configuration, facilitating the entrapment of inverse flow and leading to rapid bifurcation height growth and early convergence. Additionally, isolated islands with high density, temperature, and pressure emerge in the transitional region behind the bifurcated shocks, due to variations in wave propagation speed.

Three-dimensional compact multi-resolution weighted essentially non-oscillatory reconstruction under the Arbitrary Lagrange–Euler framework

A third-order compact multi-resolution weighted essentially non-oscillatory (CMR-WENO) reconstruction method for three-dimensional (3D) hybrid unstructured grids is developed using the Arbitrary Lagrange–Euler framework. The finite volume method is used to discretize the governing equations, and some turbulent and moving boundary problems are simulated. Only one compact center stencil comprising the neighboring cells of each control cell is required to construct the polynomials in the algorithm. As a result, the number of stencils and stencil cells is significantly reduced when compared with the traditional WENO scheme. This simplifies the code and improves the robustness of the algorithm. By ensuring the cell average and first-order derivatives are consistent with that in stencil cells an over-determined system of equations can be used to reconstruct the polynomials. This system can then be solved using the compact least squares method to avoid an ill-conditioned coefficient matrix. Furthermore, a coupled implicit iteration strategy is used to solve for the unknown coefficients, so no extra determination is required for the derivatives of each control cell. The final interpolation function for discontinuities in the flow field is obtained using CMR-WENO to nonlinearly combine polynomials of different orders, which further improves the stability of the algorithm. The CMR-WENO can be implemented on 3D hybrid unstructured grids and can be used to simulate complex problems such as those involving turbulence and moving boundaries. Finally, the algorithm presented here is verified to be third-order accurate and to exhibit good robustness when used on several representative numerical examples.

A New Alternative WENO Scheme Based on Exponential Polynomial Interpolation with an Improved Order of Accuracy

A New Alternative WENO Scheme Based on Exponential Polynomial Interpolation with an Improved Order of Accuracy

An efficient hybrid WENO scheme with minimized dispersion and adaptive dissipation properties for compressible flows

A numerical scheme with good spectral properties is important for the simulation of compressible flows with various of length scales for fine flow scales resolving. The MDAD-HY scheme (Li et al., 2022) using a discontinuity detector and scale sensor achieves the minimized dispersion and adaptive dissipation property. However, the discontinuity detector is devised based on the ratio of the 1st-order and 2nd-order derivatives on two sides of the interface introducing excessive numerical cost. To address this issue, an efficient hybrid WENO scheme with minimized dispersion and adaptive dissipation properties is proposed in this work. Based on the characteristic-decomposition approach, the numerical flux of the present hybrid scheme is achieved by switching between the linear MDAD scheme and the MDAD-WENO scheme according to a new efficient non-dimensional discontinuity detector. The linear flux is reconstructed in a component-wise method to decrease the characteristic-projection operations. To further improve the spectral property of the present scheme, an adaptive parameter controlling the contribution of the optimal linear scheme according to the discontinuity indicator is introduced. Several benchmark test cases involving broadband of length scales and discontinuities are adopted to verify the efficiency and the high-resolution capability of the present scheme.

Composition tracking of natural gas–hydrogen mixtures in pipeline flow using high-resolution schemes

A transient pipeline flow model with gas composition tracking is solved for studying the operation of a natural gas pipeline under nonisothermal flow conditions in a hydrogen injection scenario. Two approaches to high-resolution pipeline flow modeling based on the WENO scheme are presented and compared with the implicit finite difference method. The high-resolution models are capable of capturing fast fluid transients and tracking the step changes in the composition of the transported mixture. The implicit method assumes the decoupling of the flow model components in order to enhance calculation efficiency. The validation of the composition tracking results against actual gas transmission pipeline indicates that both models exhibit good prediction performance, with normalized root mean square errors of 0.406% and 1.48%, respectively. Under nonisothermal flow conditions, the prediction response of the reduced model against a high-resolution flow model, with respect to the mass and energy linepack, is at most 3.20%.

Increasingly high-order hybrid multi-resolution WENO schemes in multi-dimensions

In this paper, a new type of increasingly high-order hybrid multi-resolution weighted essentially non-oscillatory (HMR-WENO) schemes is presented in the finite difference framework for solving hyperbolic conservation laws in one, two, and three dimensions. Based on the reconstruction polynomials defined on the one-point, three-point, five-point, seven-point, and nine-point spatial stencils, we reconstruct one zeroth degree reconstruction polynomial, one quadratic reconstruction polynomial, one quartic reconstruction polynomial, one sextic reconstruction polynomial, one octave reconstruction polynomial, and their derivative polynomials together with a new hierarchical bisection method to design a series of new troubled cell indicators which can precisely find all extreme points of associated unequal degree reconstruction polynomials located inside the smallest interval in one dimension. The new troubled cell indicators do not introduce any manual parameters related to different problems. The new hybrid methodology is divided into two parts: if all extreme points of the reconstruction polynomials are nonexistent or outside the smallest interval, the target cell is not a troubled cell and the simple linear upwind schemes are utilized to obtain high-order approximations. Otherwise, the target cell is a troubled cell and the MR-WENO spatial reconstruction procedures with excellent shock-capture ability are adopted. Then a series of HMR-WENO schemes are proposed by using these new troubled cell indicators, which can be easily expanded to arbitrarily high-order accuracies in multi-dimensions. The main benefits of these HMR-WENO schemes are their efficiency, since they could save about 22%-75% CPU time than that of the same order MR-WENO schemes when simulating some benchmark examples in multi-dimensions.

A resolution-enhanced seventh-order weighted essentially non-oscillatory scheme based on non-polynomial reconstructions for solving hyperbolic conservation laws

In high-resolution numerical simulations of flows characterized by both multiscale turbulence and discontinuities, the conflict between spectral characteristics and stability becomes increasingly pronounced as the order of accuracy improves. To address this challenge, we proposed a novel seventh-order weighted essentially non-oscillatory scheme (WENO-K7). This scheme utilizes non-polynomial reconstructions by incorporating kriging interpolation and Gaussian exponential function. Then, a hyper-parameter associated with the Gaussian function is adaptively optimized to achieve higher convergence orders on sub-stencils, reducing numerical errors on global stencils. Additionally, a criterion based on monotone interpolations is devised to automatically identify problematic hyper-parameters, facilitating the transition from non-polynomial to polynomial reconstructions near discontinuities and preserving the essentially non-oscillatory property. Compared to the conventional seventh-order WENO-Z7 scheme, WENO-K7 scheme exhibits smaller computational error and reduced numerical dissipation in smooth regions while maintaining non-oscillatory and high-resolution capabilities around discontinuities. Results from various one- and two-dimensional benchmark cases demonstrate that the proposed WENO-K7 scheme outperforms the widely used WENO-Z7 scheme with only a 12% increase in computational cost. Moreover, the WENO-K7 scheme shares the same sub-stencils as the WENO-Z7 scheme, making it easily applicable to other variants of seventh-order WENO schemes and enhancing their spectral characteristics.

Application of computational algorithms with higher order of accuracy to the modeling of two-dimensional problems on development of hydrodynamic instability

This article examines application of computational algorithms with an increased order of accuracy for modeling two-dimensional problems of development of hydrodynamic instabilities. The efficiency of using algorithms to improve the accuracy and reliability of modeling in this area is considered. More specifically, the paper describes a numerical algorithm for solving the problem of development of Richtmayer-Meshkov instability. To construct the algorithm, the authors use the WENO scheme of the fifth order of accuracy Several problems are solved numerically using the developed method. The article models such processes as flows at a time of 4 046 microseconds, a change in the width of the region filled with sulfur hexafluoride, numerical schlieren patterns at a time of 877 microseconds, a change in the width of the region filled with heavy gas. The results are obtained by various methods on grids of different dimensions and compared with experimental data. It is shown that schemes with WENO reconstruction of the 5th order of accuracy demonstrate results closer to full-scale experiments.

Exploring the potential of TENO and WENO schemes for simulating under-resolved turbulent flows in the atmosphere using Euler equations

This paper focuses on the design and analysis of very high-order finite volume methods for the computation of simplified meso- and micro-scale atmospheric flows. In a dry atmosphere, these flows can be represented by the Euler equations with a gravitational source term. Two different approaches are considered here. While one of the approaches is fully conservative for the total energy, the other is formulated in a non-conservative form. The main focus of the paper is to analyze the performance of such models in combination with the traditional WENO reconstruction and the novel TENO reconstruction by examining the spectral properties of these reconstruction methods. The overarching goal is to determine whether the combination of these models and numerical schemes can be used to build an implicit Large Eddy Simulation framework, shedding light on their potential advantages or limitations in representing under-resolved atmospheric processes in the meso- and micro-scales.

An absolutely convergent fixed-point fast sweeping WENO method on triangular meshes for steady state of hyperbolic conservation laws

High order accuracy fast sweeping methods have been developed in the literature to efficiently solve steady state solutions of hyperbolic partial differential equations (PDEs) on rectangular meshes, while they were not available yet on unstructured meshes. In this paper, we extend high order accuracy fast sweeping methods to unstructured triangular meshes by applying fixed-point iterative sweeping techniques to a fifth-order finite volume unstructured WENO scheme, for solving steady state solutions of hyperbolic conservation laws. Similar as other fast sweeping methods, fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in each sweeping order to achieve fast convergence rate to steady state solutions. An advantage of fixed-point fast sweeping methods which distinguishes them from other fast sweeping methods is that they are explicit and do not require inverse operation of nonlinear local systems. This also provides certain convenience in designing high order fast sweeping methods on unstructured meshes. As in the first order fast sweeping methods on triangular meshes, we introduce multiple reference points to determine alternating sweeping directions on unstructured meshes. All the cells on the mesh are ordered according to their centroids' distances to those reference points, and the resulted orderings provide sweeping directions for iterations. To make the residue of the fast sweeping iterations converge to machine zero / round off errors, we follow the approach in our early work of developing the absolutely convergent fixed-point fast sweeping WENO methods on rectangular meshes, and adopt high order WENO scheme with unequal-sized sub-stencils, specifically here a fifth-order finite volume unstructured WENO scheme for spatial discretization. Extensive numerical experiments, including problems with complex domain geometries, are performed to show the accuracy, computational efficiency, and absolute convergence of the presented fifth-order fast sweeping scheme on triangular meshes. Furthermore, the proposed method is compared with the forward Euler time marching method and the popular third order total variation diminishing Rung-Kutta (TVD-RK3) time-marching method for steady state computations. Numerical examples show that the developed fixed-point fast sweeping WENO method is the most efficient scheme among them, and especially it can save up to 70% CPU time costs than TVD-RK3 in order to converge to steady state solutions.

A modified fifth‐order WENO‐Z scheme based on the weights of the reformulated adaptive order WENO scheme

AbstractA modified fifth‐order WENO‐Z scheme is proposed by analogy with the non‐normalized weights of the reformulated fifth‐order adaptive order (AO) WENO scheme. We show that if the original fifth‐order WENO‐AO scheme is rewritten as the form of the conventional WENO combination, the resulting non‐normalized weights can be divided into three parts: a constant one term, a local stencil smoothness measure term and a global stencil smoothness measure term. In order to make use of the latter two terms for constructing a modified WENO‐Z scheme with enhanced performance, we change the form of the third term and introduce an adaptive scaling factor to adjust the contributions from the second and third terms. Numerical examples show that the modified fifth‐order WENO‐Z scheme has the advantage of high resolution in smooth regions and sharp capturing of discontinuities, and it can obtain evidently better results for shocked flows with small‐scale structures compared with the recently developed WENO‐Z+ and WENO‐Z+M schemes.

An Energy Stable Alternative WENO Scheme with a New Smoothness Indicator for Hyperbolic Conservation Laws

An Energy Stable Alternative WENO Scheme with a New Smoothness Indicator for Hyperbolic Conservation Laws

Efficient Alternative Finite Difference WENO Schemes for Hyperbolic Systems with Non-conservative Products

Efficient Alternative Finite Difference WENO Schemes for Hyperbolic Systems with Non-conservative Products

Efficient WENO Schemes for Nonuniform Grids

A set of arbitrarily high-order WENO schemes for reconstructions on nonuniform grids is presented. These non-linear interpolation methods use simple smoothness indicators with a linear cost with respect to the order, making them easy to implement and computationally efficient. The theoretical analysis to verify the accuracy and the essentially non-oscillatory properties are presented together with some numerical experiments involving algebraic problems in order to validate them. Also, these general schemes are applied for the solution of conservation laws and hyperbolic systems in the context of finite volume methods.

Distinct-positivity-preserving methods for fifth-order finite volume ALW-MR-WENO schemes with a bigger sufficient CFL number for extreme problems

In this paper, the new distinct-positivity-preserving (DPP) methods are proposed for fifth-order finite volume multi-resolution WENO schemes with adaptive linear weights (which are termed as the ALW-MR-WENO schemes) [J. Comput. Phys., 493 (2023), 112471] for solving some extreme problems on structured meshes. Associated WENO spatial reconstruction procedures only require two hierarchical unequal-sized central spatial stencils for achieving fifth-order accuracy in smooth regions and keeping the essentially non-oscillatory property in non-smooth regions. One redefines five cell averages after performing such spatial reconstructions and then designs one quartic polynomial and two cubic reconstruction polynomials based on them in one dimension. After that, a new detective process is used to examine the positivity of the point values of such cubic reconstruction polynomials at some checking points inside the target cell. If the negativity happens, a new over-determined compression limiter is carried out to ensure that the compressed polynomials can achieve the positivity in the whole target cell instead of only at some discrete points inside it. The quartic polynomial could use its four quadrature point values by using a four-point Gauss-Lobatto quadrature formula when performing numerical integrations. Since it is easy to rewrite the point values of one quartic polynomial with the point values of two cubic polynomials timing with associated linear weights, both cubic polynomials could use a four-point Gauss-Lobatto quadrature formula or a three-point Simpson's quadrature formula without losing any algebra precisions. By doing this novelty switching methodology between two quadrature formulas, it is the first time to theoretically prove and increase a sufficient CFL number from 1/12 to 1/7.2 for the fifth-order WENO schemes in one dimension. This methodology can be easily extended to multi-dimensions. Finally, the novelty DPP methods for fifth-order ALW-MR-WENO schemes with a larger practical CFL number of 0.6, instead of the sufficient CFL number of 1/7.2, are also simple and robust enough for some extreme problems without timely halving the CFL number to avoid the appearance of negative density and negative pressure.

Numerical simulation of 3-D 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 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.

Application of the CABARET and WENO Schemes for Solving the Nonlinear Transport Equation in the Problem of Simulating the Propagation of a Sonic Boom Wave in the Atmosphere

Application of the CABARET and WENO Schemes for Solving the Nonlinear Transport Equation in the Problem of Simulating the Propagation of a Sonic Boom Wave in the Atmosphere

Fourth-Order Conservative Non-splitting Semi-Lagrangian Hermite WENO Schemes for Kinetic and Fluid Simulations

AbstractWe present fourth-order conservative non-splitting semi-Lagrangian (SL) Hermite essentially non-oscillatory (HWENO) schemes for linear transport equations with applications for nonlinear problems including the Vlasov–Poisson system, the guiding center Vlasov model, and the incompressible Euler equations in the vorticity-stream function formulation. The proposed SL HWENO schemes combine a weak formulation of the characteristic Galerkin method with two newly constructed HWENO reconstruction methods. The new HWENO reconstructions are meticulously designed to strike a delicate balance between curbing numerical oscillation and introducing excessive dissipation. Mass conservation naturally holds due to the weak formulation of the semi-Lagrangian discontinuous Galerkin method and the design of the HWENO reconstructions. We apply a positivity-preserving limiter to maintain the positivity of numerical solutions when needed. Abundant benchmark tests are performed to verify the effectiveness of the proposed SL HWENO schemes.