High-order Finite Volume Research Articles

An Essentially Non-Oscillatory Spectral Deferred Correction Method for Conservation Laws

We present a computational method based on the Spectral Deferred Corrections (SDC) time integration technique and the Essentially Non-Oscillatory (ENO) finite volume method for the conservation laws (one-dimensional Euler equations). The SDC technique is used to advance the solutions in time with high-order of accuracy. The ENO method is used to define high-order cell edge quantities that are then used to evaluate numerical fluxes. The coupling of the SDC method with a high-order finite volume method (Piece-wise Parabolic Method (PPM)) for solving the conservation laws is first carried out by Layton et al. in [Layton, A. T. and Minion, M. L. [2004] “Conservative multi-implicit spectral deferred correction methods for reacting gas dynamics,” J. Comput. Phys. 194(2), 697–714]. Issues about this approach have been addressed and some improvements have been added to it in [Kadioglu et al. [2012] “A gas dynamics method based on the spectral deferred corrections (SDC) time integration technique and the piecewise parabolic method (PPM),” Am. J. Comput. Math. 1–4, 303–317]. Here, we investigate the implications when the PPM method is replaced with the well-known ENO method. We note that the SDC-PPM method is fourth-order accurate in time and space. Therefore, we kept the order of accuracy of the ENO procedure as fourth-order in order to be able to make a consistent comparison between the two approaches (SDC-ENO versus SDC-PPM methods). We have tested the new SDC-ENO technique by solving several test problems involving moderate to strong shock waves and smooth/complex flow structures. Our numerical results show that we have numerically achieved the formally fourth-order convergence of the new method for smooth problems. Our numerical results also indicate that the newly proposed technique performs very well providing highly resolved shock discontinuities and fairly good contact solutions. More importantly, the discontinuities in the flow test problems are captured with essentially no-oscillations. We have numerically compared the fourth-order SDC-ENO scheme to the fourth-order SDC-PPM method for the same test problems. The results are similar for most of the test problems except in some cases the SDC-PPM method suffers from minor oscillations compared to SDC-ENO scheme being completely oscillation free.

A high-order finite volume method on unstructured grids using RBF reconstruction

This paper proposes a high-order finite volume method based on radial basis function (RBF) reconstruction for the solution of Euler and Navier–Stokes equations on unstructured grids. Unlike traditional polynomial K-exact method, RBF method has stronger adaptability for different reconstruction stencils and more flexibility in choosing interpolating points. We expatiate on the detailed process of flow-field reconstruction by using multiquadric (MQ) basis function for the second-order and third-order schemes on unstructured triangular grids. Subsequently, we validate the accuracy order of RBF method through the numerical test case. Furthermore, the method is used to solve several typical flow fields. Compared with traditional K-exact high-order scheme, RBF method is more accurate and has lower numerical dissipation, which can obtain more elaborate and precise results.

Adjoint-based an adaptive finite volume method for steady Euler equations with non-oscillatory k -exact reconstruction

Abstract In [ An adaptive finite volume solver for steady Euler equations with non-oscillatory k-exact reconstruction, G. H. Hu and N. Y. Yi, J. Comput. Phys., Vol. 312, pp. 235–251 ], a high order finite volume method was proposed for steady Euler equations. The method consists of a Newton iteration method to linearize the governing equations, and a geometrical multigrid method to solve the derived system. The non-oscillatory k -exact reconstruction was developed towards the high quality solution reconstruction, and h -adaptive method was introduced to handle the efficiency issue. To further improve the numerical method, we study the following two issues in detail in this paper. First, we test the new reconstruction patch for the element which locates on the domain boundary, which could improve the convergence of the steady state significantly. Second, the adjoint-based a posteriori error estimation is introduced for the h -adaptive method, which could deliver more efficient adaptive refinement of the mesh. The numerical tests show the effectiveness of the proposed method.

High-order finite volume WENO schemes for Boussinesq modelling of nearshore wave processes

ABSTRACTThis paper presents the formulation and validation of a shock‐capturing hybrid method using the WENO (weighted essential non‐oscillatory)‐HLL (Harten–Lax–Van Leer) scheme for the solution of the two-dimensional extended Boussinesq-type equations. The governing equations are reformulated in a conservative form, using high-order finite volume WENO schemes with HLL Riemann solver to handle the flux terms. The efficient third-order Runge–Kutta time-stepping scheme is used to advance the numerical solution through time. Wave breaking is predicted by applying a wave breaking criterion to automatically degenerate the Boussinesq equations into nonlinear shallow water equations. In order to test whether the present numerical model could accurately reproduce the processes of transformation, breaking, and run-up of nonlinear dispersive waves under different wave conditions in the nearshore region, five test cases have been simulated. The numerical results show that the present shock-capturing model could simulate the nearshore wave processes very well.

High-order upwind finite volume element method for first-order hyperbolic optimal control problems

We present a high-order upwind finite volume element method to solve optimal control problems governed by first-order hyperbolic equations. The method is efficient and easy for implementation. Both the semi-discrete error estimates and the fully discrete error estimates are derived. Optimal order error estimates in the sense of \(L_{2}\)-norm are obtained. Numerical examples are provided to confirm the effectiveness of the method and the theoretical results. doi:10.1017/S1446181116000031

Compact high order finite volume method on unstructured grids I: Basic formulations and one-dimensional schemes

The large reconstruction stencil has been the major bottleneck problem in developing high order finite volume schemes on unstructured grids. This paper presents a compact reconstruction procedure for arbitrarily high order finite volume method on unstructured grids to overcome this shortcoming. In this procedure, a set of constitutive relations are constructed by requiring the reconstruction polynomial and its derivatives on the control volume of interest to conserve their averages on face-neighboring cells. These relations result in an over-determined linear equation system, which, in the sense of least-squares, can be reduced to a block-tridiagonal system in the one-dimensional case. The one-dimensional formulations of the reconstruction are discussed in detail and a Fourier analysis is presented to study the dispersion/dissipation and stability properties. The WBAP limiter based on the secondary reconstruction is used to suppress the non-physical oscillations near discontinuities while achieve high order accuracy in smooth regions of the solution. Numerical results demonstrate the method's high order accuracy, robustness and shock capturing capability. A new reconstruction procedure for high order finite volume scheme is proposed.This procedure can achieve arbitrary order of accuracy using a compact stencil consisting of only face neighboring cells.The spectral properties of the proposed finite volume schemes are analyzed.Some test cases are solved numerically to show the superior performances of the proposed schemes.

Parametric Investigation of Breaking Solitary Wave Over Fringing Reef Based on Shock-Capturing Boussinesq Model

In this paper, a hybrid finite-volume-finite-difference scheme for the solution of the one-dimensional form of the Nwogu's extended Boussinesq-type equations was developed to investigate the breaking of a solitary wave propagating over a fringing reef system. The adoption of a high-order finite volume WENO schemes with HLL Riemann solver for advection terms as well as the third-order Runge—Kutta time-stepping scheme makes the present scheme shock-capturing. From the validation case, it is shown that the present shock-capturing model could simulate the shoreline movement and wave breaking as well as bore propagation very well. In order to comprehensively understand the hydrodynamic characteristics when a solitary wave travels over a fringing reef system, a parametric study has been carried out to analyze the interaction process under various conditions. Effects of different physical parameters on wave breaking axe analyzed and the reflection, transmission and dissipation (RTD) coefficients are calculated based on the integration of kinetic and potential energies in different zones.

An Essentially Nonoscillatory Spectral Deferred Correction Method for Hyperbolic Problems

We present a computational method based on the spectral deferred corrections (SDC) time integration technique and the essentially nonoscillatory (ENO) finite volume method for hyperbolic problems. The SDC technique is used to advance the solutions in time with high-order of accuracy. The ENO method is used to define high-order cell edge quantities that are then used to evaluate numerical fluxes. The coupling of the SDC method with a high-order finite volume method (piece-wise parabolic method (PPM)) is first carried out by Layton et al. [J. Comput. Phys. 194(2) (2004) 697]. Issues about this approach have been addressed and some improvements have been added to it in Kadioglu et al. [J. Comput. Math. 1(4) (2012) 303]. Here, we investigate the implications when the PPM method is replaced with the well-known ENO method. We note that the SDC-PPM method is fourth-order accurate in time and space. Therefore, we kept the order of accuracy of the ENO procedure as fourth-order in order to be able to make a consistent comparison between the two approaches (SDC-ENO versus SDC-PPM). We have tested the new SDC-ENO technique by solving smooth and nonsmooth hyperbolic problems. Our numerical results indicate that the fourth-order of accuracy in both space and time has been achieved for smooth problems. On the other hand, the new method performs very well when it is applied to nonlinear problems that involve discontinuous solutions. In other words, we have obtained highly resolved discontinuous solutions with essentially no-oscillations at or around the discontinuities.

Cell centered direct Arbitrary-Lagrangian-Eulerian ADER-WENO finite volume schemes for nonlinear hyperelasticity

This paper is concerned with the numerical solution of the unified first order hyperbolic formulation of continuum mechanics proposed by Peshkov & Romenski [Peshkov I, Romenski E. A hyperbolic model for viscous Newtonian flows. Continuum Mechanics and Thermodynamics 2016;28:85–104.], which is based on the theory of nonlinear hyperelasticity of Godunov & Romenski [67][Godunov S, Romenski E. Nonstationary equations of the nonlinear theory of elasticity in Euler coordinates. Journal of Applied Mechanics and Technical Physics 1972;13:868–885.][Godunov S., Romenski E., Elements of Continuum Mechanics and Conservation Laws. Kluwer Academic/ Plenum Publishers; 2003.], further denoted by GPR model. Notably, the governing PDE system is symmetric hyperbolic and fully consistent with the first and the second principle of thermodynamics. The nonlinear system of governing equations of the GPR model is overdetermined, large and includes stiff source terms as well as non-conservative products. In this paper we solve this model for the first time on moving unstructured meshes in multiple space dimensions by employing high order accurate one-step ADER-WENO finite volume schemes in the context of cell-centered direct Arbitrary-Lagrangian-Eulerian (ALE) algorithms.The numerical method is based on a WENO polynomial reconstruction operator on moving unstructured meshes, a fully-discrete one-step ADER scheme that is able to deal with stiff sources [Dumbser M., Enaux C., Toro E., Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws. Journal of Computational Physics . 2008a;227:3971–4001.], a nodal solver with relaxation to determine the mesh motion, and a path-conservative technique of Castro & Parés for the treatment of non-conservative products [Parés C., Numerical methods for nonconservative hyperbolic systems: a theoretical framework. SIAM Journal on Numerical Analysis. 2006;44:300–321.][Castro M, Gallardo J, Parés C. High-order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. applications to shallow-water systems. Mathematics of Computation 2006;75:1103–1134.]. We present numerical results obtained by solving the GPR model with ADER-WENO-ALE schemes in the stiff relaxation limit, showing that fluids (Euler or Navier-Stokes limit), as well as purely elastic or elasto-plastic solids can be simulated in the framework of nonlinear hyperelasticity with the same system of governing PDE. The obtained results are in good agreement when compared to exact or numerical reference solutions available in the literature.

High-order central Hermite WENO schemes: Dimension-by-dimension moment-based reconstructions

In this paper, a class of high-order central finite volume schemes is proposed for solving one- and two-dimensional hyperbolic conservation laws. Formulated on staggered meshes, the methods involve Hermite WENO (HWENO) spatial reconstructions, and Lax–Wendroff type discretizations or the natural continuous extension of Runge–Kutta methods in time. Differently from the central Hermite WENO methods we developed previously in Tao et al. (2015) [34] , the spatial reconstructions, a core ingredient of the methods, are based on the zeroth-order and the first-order moments of the solution, and are implemented through a dimension-by-dimension strategy when the spatial dimension is higher than one. This leads to much simpler implementation of the methods in higher dimension and better cost efficiency. Meanwhile, the proposed methods have the attractive features of the general central Hermite WENO methods such as being compact in reconstruction and requiring neither flux splitting nor numerical fluxes, while being accurate and essentially non-oscillatory. A collection of one- and two-dimensional numerical examples is presented to demonstrate high resolution and robustness of the methods in capturing smooth and non-smooth solutions.

Xtroem-fv: a new code for computational astrophysics based on very high order finite-volume methods – II. Relativistic hydro- and magnetohydrodynamics

xtroem-fv: a new code for computational astrophysics based on very high order finite-volume methods – II. Relativistic hydro- and magnetohydrodynamics

Numerical Discretization of Coupling Conditions by High-Order Schemes

We consider numerical schemes for $$2\times 2$$2×2 hyperbolic conservation laws on graphs. The hyperbolic equations are given on arcs which are one-dimensional in space and are coupled at a single point, the node, by a nonlinear coupling condition. We develop high-order finite volume discretizations for the coupled problem. The reconstruction of the fluxes at the node is obtained using derivatives of the parameterized algebraic conditions imposed at the nodal points in the network. Numerical results illustrate the expected theoretical behavior.

High order finite volume WENO schemes for the Euler equations under gravitational fields

Euler equations with gravitational source terms are used to model many astrophysical and atmospheric phenomena. This system admits hydrostatic balance where the flux produced by the pressure is exactly canceled by the gravitational source term, and two commonly seen equilibria are the isothermal and polytropic hydrostatic solutions. Exact preservation of these equilibria is desirable as many practical problems are small perturbations of such balance. High order finite difference weighted essentially non-oscillatory (WENO) schemes have been proposed in [22], but only for the isothermal equilibrium state. In this paper, we design high order well-balanced finite volume WENO schemes, which can preserve not only the isothermal equilibrium but also the polytropic hydrostatic balance state exactly, and maintain genuine high order accuracy for general solutions. The well-balanced property is obtained by novel source term reformulation and discretization, combined with well-balanced numerical fluxes. Extensive one- and two-dimensional simulations are performed to verify well-balanced property, high order accuracy, as well as good resolution for smooth and discontinuous solutions.

On spatial filtering of flow variables in high-order finite volume methods

A new method of spatial filtering in high-order finite volume methods is presented and assessed. The base of this method is to filter face-averaged variables (fluxes) and then the recovery of cell-averaged ones. Two kinds of filtering method are proposed. The first kind is highly dissipative and appropriate for the numerical regions that need high dissipation, e.g. sponge zones. The second kind, on the other hand, is a precise method and hence is suitable for applying the high-order finite difference filters to the finite volume methods. Applying high-order finite difference filters directly to the high-order finite volume methods without using the proposed method causes stability problems in the large eddy simulation of high Reynolds number flows. The test cases, namely, propagation of the one-dimensional wave packet, advection of a two-dimensional vortical wave and the large eddy simulation of turbulent round jet at Re=104 and 105 are used to examine the accuracy and performance of the proposed methods. The last test case shows the effectiveness of the proposed methods in the large eddy simulation of flows by high-order FVMs. The performance of filtering methods on a curvilinear wavy grid is also investigated.

HIGH-ORDER UPWIND FINITE VOLUME ELEMENT METHOD FOR FIRST-ORDER HYPERBOLIC OPTIMAL CONTROL PROBLEMS

We present a high-order upwind finite volume element method to solve optimal control problems governed by first-order hyperbolic equations. The method is efficient and easy for implementation. Both the semi-discrete error estimates and the fully discrete error estimates are derived. Optimal order error estimates in the sense of $L^{2}$-norm are obtained. Numerical examples are provided to confirm the effectiveness of the method and the theoretical results.

Divergence-Free WENO Reconstruction-Based Finite Volume Scheme for Solving Ideal MHD Equations on Triangular Meshes

AbstractIn this paper, we introduce a high-order accurate constrained transport type finite volume method to solve ideal magnetohydrodynamic equations on two-dimensional triangular meshes. A new divergence-free WENO-based reconstruction method is developed to maintain exactly divergence-free evolution of the numerical magnetic field. In this formulation, the normal component of the magnetic field at each face of a triangle is reconstructed uniquely and with the desired order of accuracy. Additionally, a new weighted flux interpolation approach is also developed to compute the z-component of the electric field at vertices of grid cells. We also present numerical examples to demonstrate the accuracy and robustness of the proposed scheme.

Bound-preserving discontinuous Galerkin methods for relativistic hydrodynamics

In this paper, we develop a discontinuous Galerkin (DG) method to solve the ideal special relativistic hydrodynamics (RHD) and design a bound-preserving (BP) limiter for this scheme by extending the idea in X. Zhang and C.-W. Shu, (2010) [56]. For RHD, the density and pressure are positive and the velocity is bounded by the speed of light. One difficulty in numerically solving the RHD in its conservative form is that the failure of preserving these physical bounds will result in ill-posedness of the problem and blowup of the code, especially in extreme relativistic cases. The standard way in dealing with this difficulty is to add extra numerical dissipation, while in doing so there is no guarantee of maintaining the high order of accuracy. Our BP limiter has the following features. It can theoretically guarantee to preserve the physical bounds for the numerical solution and maintain its designed high order accuracy. The limiter is local to the cell and hence is very easy to implement. Moreover, it renders L1-stability to the numerical scheme. Numerical experiments are performed to demonstrate the good performance of this bound-preserving DG scheme. Even though we only discuss the BP limiter for DG schemes, it can be applied to high order finite volume schemes, such as weighted essentially non-oscillatory (WENO) finite volume schemes as well.

Compact high order finite volume method on unstructured grids II: Extension to two-dimensional Euler equations

Abstract In this paper, the compact least-squares finite volume method on unstructured grids proposed in our previous paper is extended to multi-dimensional systems, namely the two-dimensional Euler equations. The key element of this scheme is the compact least-squares reconstruction in which a set of constitutive relations are constructed by requiring the reconstruction polynomial and its spatial derivatives on the control volume of interest to conserve their averages on the face-neighboring cells. These relations result in an over-determined linear equation system. A large sparse system of linear equations is resulted by using the least-squares technique. An efficient solution strategy is of crucial importance for the application of the proposed scheme in multi-dimensional problems since both direct and iterative solvers for this system are computationally very expensive. In the present paper, it is found that in the cases of steady flow simulation and unsteady flow simulation using dual time stepping technique, the present reconstruction method can be coupled with temporal discretization scheme to achieve high computational efficiency. The WBAP limiter and a problem-independent shock detector are used in the simulation of flow with discontinuities. Numerical results demonstrate the high order accuracy, high computational efficiency and capability of handling both complex physics and geometries of the proposed schemes.

A fifth-order finite volume weighted compact scheme for solving one-dimensional Burgers’ equation

In the present paper, a high-order finite volume compact scheme is proposed to solve one dimensional Burgers’ equation. The nonlinear advective terms are computed by the fifth-order finite volume weighted upwind compact scheme, in which the nonlinear weighted essentially non-oscillatory weights are coupled with lower order compact stencils. The diffusive terms are discretized by using the finite volume six-order Padé scheme. The strong stability preserving third-order Runge–Kutta time discretizations is used in this work. Numerical results are compared with the exact and some existing numerical solutions to demonstrate the essentially non-oscillatory and high resolution of the proposed method. The numerical results are shown to be more accurate than some numerical results given in the literature.

A high-order local time stepping finite volume solver for one-dimensional blood flow simulations: application to the ADAN model.

In recent years, the complexity of vessel networks for one-dimensional blood flow models has significantly increased, because of enhanced anatomical detail or automatic peripheral vasculature generation, for example. This fact, along with the application of these models in uncertainty quantification and parameter estimation poses the need for extremely efficient numerical solvers. The aim of this work is to present a finite volume solver for one-dimensional blood flow simulations in networks of elastic and viscoelastic vessels, featuring high-order space-time accuracy and local time stepping (LTS). The solver is built on (i) a high-order finite volume type numerical scheme, (ii) a high-order treatment of the numerical solution at internal vertexes of the network, often called junctions, and (iii) an accurate LTS strategy. The accuracy of the proposed methodology is verified by empirical convergence tests. Then, the resulting LTS scheme is applied to arterial networks of increasing complexity and spatial scale heterogeneity, with a number of one-dimensional segments ranging from a few tens up to several thousands and vessel lengths ranging from less than a millimeter up to tens of centimeters, in order to evaluate its computational cost efficiency. The proposed methodology can be extended to any other hyperbolic system for which network applications are relevant. Copyright © 2016 John Wiley & Sons, Ltd.