Assessment of different reconstruction techniques for implementing the NVSF schemes on unstructured meshes

Three new far‐upwind reconstruction techniques, New‐Technique 1, 2, and 3, are proposed in this paper, which localize the normalized variable and space formulation (NVSF) schemes and facilitate the implementation of standard bounded high‐resolution differencing schemes on arbitrary unstructured meshes. By theoretical analysis, it is concluded that the three new techniques overcome two inherent drawbacks of the original technique found in the literature. Eleven classic high‐resolution NVSF schemes developed in the past decades are selected to evaluate performances of the three new techniques relative to the original technique. Under the circumstances of arbitrary unstructured meshes, stretched meshes, and uniform triangular meshes, for each NVSF scheme, the accuracies and convergence properties, when implementing the four aforementioned far‐upwind reconstruction techniques respectively, are assessed by the pure convection of several scalar profiles. The numerical results clearly show that New‐Technique‐2 leads to a better performance in terms of overall accuracy and convergence behavior for the 11 NVSF schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

A new r‐ratio formulation for TVD schemes for vertex‐centered FVM on an unstructured mesh

The r‐ratio is a parameter that measures the local monotonicity, by which a number of high‐resolution and TVD schemes can be formed. A number of r‐ratio formulations for TVD schemes have been presented over the last few decades to solve the transport equation in shallow waters based on the finite volume method (FVM). However, unlike structured meshes, the coordinate directions are not clearly defined on an unstructured mesh; therefore, some r‐ratio formulations have been established by approximating the solute concentration at virtual nodes, which may be estimated from different assumptions. However, some formulations may introduce either oscillation or diffusion behavior within the vertex‐centered (VC) framework. In this paper, a new r‐ratio formulation, applied to an unstructured grid in the VC framework, is proposed and compared with the traditional r‐ratio formulations. Through seven commonly used benchmark tests, it is shown that the newly proposed r‐ratio formulation obtains better results than the traditional ones with less numerical diffusion and spurious oscillation. Moreover, three commonly used TVD schemes—SUPERBEE, MINMOD, and MUSCL—and two high‐order schemes—SOU and QUICK—are implemented and compared using the new r‐ratio formulation. The new r‐ratio formulation is shown to be sufficiently comprehensive to permit the general implementation of a high‐resolution scheme within the VC framework. Finally, the sensitivity test for different grid types demonstrates the good adaptability of this new r‐ratio formulation. Copyright © 2015 John Wiley & Sons, Ltd.     

Efficient high‐resolution relaxation schemes for hyperbolic systems of conservation laws

In this work, we present a total variation diminishing (TVD) scheme in the zero relaxation limit for nonlinear hyperbolic conservation law using flux limiters within the framework of a relaxation system that converts a nonlinear conservation law into a system of linear convection equations with nonlinear source terms. We construct a numerical flux for space discretization of the obtained relaxation system and modify the definition of the smoothness parameter depending on the direction of the flow so that the scheme obeys the physical property of hyperbolicity. The advantages of the proposed scheme are that it can give second‐order accuracy everywhere without introducing oscillations for 1‐D problems (at least with) smooth initial condition. Also, the proposed scheme is more efficient as it works for any non‐zero constant value of the flux limiter ? ? [0, 1], where other TVD schemes fail. The resulting scheme is shown to be TVD in the zero relaxation limit for 1‐D scalar equations. Bound for the limiter function is obtained. Numerical results support the theoretical results. Copyright © 2007 John Wiley & Sons, Ltd.     

Weighted essential non‐oscillatory schemes for tidal bore on unstructured meshes

In this paper, the third‐order weighted essential non‐oscillatory (WENO) schemes are used to simulate the two‐dimensional shallow water equations with the source terms on unstructured meshes. The balance of the flux and the source terms makes the shallow water equations fit to non‐flat bottom questions. The simulation of a tidal bore on an estuary with trumpet shape and Qiantang river is performed; the results show that the schemes can be used to simulate the current flow accurately and catch the stronger discontinuous in water wave, such as dam break and tidal bore effectively. Copyright © 2008 John Wiley & Sons, Ltd.     

A well‐balanced stable generalized Riemann problem scheme for shallow water equations using adaptive moving unstructured triangular meshes

We propose a well‐balanced stable generalized Riemann problem (GRP) scheme for the shallow water equations with irregular bottom topography based on moving, adaptive, unstructured, triangular meshes. In order to stabilize the computations near equilibria, we use the Rankine–Hugoniot condition to remove a singularity from the GRP solver. Moreover, we develop a remapping onto the new mesh (after grid movement) based on equilibrium variables. This, together with the already established techniques, guarantees the well‐balancing. Numerical tests show the accuracy, efficiency, and robustness of the GRP moving mesh method: lake at rest solutions are preserved even when the underlying mesh is moving (e.g., mesh points are moved to regions of steep gradients), and various comparisons with fixed coarse and fine meshes demonstrate high resolution at relatively low cost. Copyright © 2013 John Wiley & Sons, Ltd.     

Higher‐order bounded differencing schemes for compressible and incompressible flows

In recent years, three higher‐order (HO) bounded differencing schemes, namely AVLSMART, CUBISTA and HOAB that were derived by adopting the normalized variable formulation (NVF), have been proposed. In this paper, a comparative study is performed on these schemes to assess their numerical accuracy, computational cost as well as iterative convergence property. All the schemes are formulated on the basis of a new dual‐formulation in order to facilitate their implementations on unstructured meshes. Based on the proposed dual‐formulation, the net effective blending factor (NEBF) of a high‐resolution (HR) scheme can now be measured and its relevance on the accuracy and computational cost of a HR scheme is revealed on three test problems: (1) advection of a scalar step‐profile; (2) 2D transonic flow past a circular arc bump; and (3) 3D lid‐driven incompressible cavity flow. Both density‐based and pressure‐based methods are used for the computations of compressible and incompressible flow, respectively. Computed results show that all the schemes produce solutions which are nearly as accurate as the third‐order QUICK scheme; however, without the unphysical oscillations which are commonly inherited from the HO linear differencing scheme. Generally, it is shown that at higher value of NEBF, a HR scheme can attain better accuracy at the expense of computational cost. Copyright © 2006 John Wiley & Sons, Ltd.     

A semi‐implicit scheme for 3D free surface flows with high‐order velocity reconstruction on unstructured Voronoi meshes

In this paper, we present a computationally efficient semi‐implicit scheme for the simulation of three‐dimensional hydrostatic free surface flow problems on staggered unstructured Voronoi meshes. For each polygonal control volume, the pressure is defined in the cell center, whereas the discrete velocity field is given by the normal velocity component at the cell faces. A piecewise high‐order polynomial vector velocity field is then reconstructed from the scalar normal velocities at the cell faces by using a new high‐order constrained least‐squares reconstruction operator. The reconstructed high‐order piecewise polynomial velocity field is used for trajectory integration in a semi‐Lagrangian approach to discretize the nonlinear convective terms in the governing PDE. For that purpose, a high‐order Taylor method is used as ODE integrator. The resulting semi‐implicit algorithm is extensively validated on a large set of different academic test problems with exact analytical solution and is finally applied to a real‐world engineering problem consisting of a curved channel upstream of two micro‐turbines of a hydroelectric power plant. For this realistic case, some experimental reference data are available from field measurements. Copyright © 2012 John Wiley & Sons, Ltd.     

A gas‐kinetic scheme for the simulation of turbulent flows on unstructured meshes

This paper presents a numerical method for simulating turbulent flows via coupling the Boltzmann BGK equation with Spalart–Allmaras one equation turbulence model. Both the Boltzmann BGK equation and the turbulence model equation are carried out using the finite volume method on unstructured meshes, which is different from previous works on structured grid. The application of the gas‐kinetic scheme is extended to the simulation of turbulent flows with arbitrary geometries. The adaptive mesh refinement technique is also adopted to reduce the computational cost and improve the efficiency of meshes. To organize the unstructured mesh data structure efficiently, a non‐manifold hybrid mesh data structure is extended for polygonal cells. Numerical experiments are performed on incompressible flow over a smooth flat plate and compressible turbulent flows around a NACA 0012 airfoil using unstructured hybrid meshes. These numerical results are found to be in good agreement with experimental data and/or other numerical solutions, demonstrating the applicability of the proposed method to simulate both subsonic and transonic turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.     

High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

On the practical importance of the SSP property for Runge–Kutta time integrators for some common Godunov‐type schemes

We investigate through analysis and computational experiment explicit second and third‐order strong‐stability preserving (SSP) Runge–Kutta time discretization methods in order to gain perspective on the practical necessity of the SSP property. We consider general theoretical SSP limits for these schemes and present a new optimal third‐order low‐storage SSP method that is SSP at a CFL number of 0.838. We compare results of practical preservation of the TVD property using SSP and non‐SSP time integrators to integrate a class of semi‐discrete Godunov‐type spatial discretizations. Our examples involve numerical solutions to Burgers' equation and the Euler equations. We observe that ‘well‐designed’ non‐SSP and non‐optimal SSP schemes with SSP coefficients less than one provide comparable stability when used with time steps below the standard CFL limit. Results using a third‐order non‐TVD CWENO scheme are also presented. We verify that the documented SSP methods with the number of stages greater than the order provide a useful enhanced stability region. We show by analysis and by numerical experiment that the non‐oscillatory third‐order reconstructions used in (Liu and Tadmor Numer. Math. 1998; 79 :397–425, Kurganov and Petrova Numer. Math. 2001; 88 :683–729) are in general only second‐ and first‐order accurate, respectively. Copyright © 2005 John Wiley & Sons, Ltd.