No abstract

To cite this version:Stéphane Clain, Raphaël Loubère, Gaspar Machado. a posteriori stabilized sixth-order finite volume scheme for one-dimensional steady-state hyperbolic equations. 2016. a posteriori stabilized sixth-order finite volume scheme for one-dimensional steady-state hyperbolic equationsWe propose a new family of finite volume high-accurate numerical schemes devoted to solve one-dimensional steady-state hyperbolic systems. High-accuracy (up to the sixth-order presently) is achieved thanks to polynomial reconstructions while stability is provided with an a posteriori MOOD method which control the cell polynomial degree for eliminating non-physical oscillations in the vicinity of discontinuities. Such a procedure demands the determination of a chain detector to discriminate between troubled and valid cells, a cascade of polynomial degrees to be successively tested when oscillations are detected, and a parachute scheme corresponding to the last, viscous, and robust scheme of the cascade. Experimented on linear, Burgers', and Euler equations, we demonstrate that the schemes manage to retrieve smooth solutions with optimal order of accuracy but also irregular solutions without spurious oscillations.

Accuracy may be dramatically reduced when the boundary domain is curved and numerical schemes require a specific treatment of the boundary condition to preserve the optimal order. In the finite volume context, Ollivier-Gooch and Van Altena (2002) has proposed a technique to overcome such limitation and restore the very high-order accuracy which consists in specific restrictions considered in the least-squares minimization associated to the polynomial reconstruction. The method suffers from several drawbacks, particularly, the use of curved elements that requires sophisticated meshing algorithms. We propose a new method where the physical domain and the computational domain are distinct and we introduce the Reconstruction for Off-site Data (ROD) where polynomial reconstructions are carried out on the mesh using data localized outside of the computational domain, namely the Dirichlet condition situated on the physical domain. A series of numerical tests assess the accuracy, convergence rates, robustness, and efficiency of the new method and show that the boundary condition is fully integrated in the scheme with a very high-order accuracy and the optimal convergence rate is achieved.

Summary Obtaining very high‐order accurate solutions in curved domains is a challenging task as the accuracy of discretization methods may dramatically reduce without an appropriate treatment of boundary conditions. The classical techniques to preserve the nominal convergence order of accuracy, proposed in the context of finite element and finite volume methods, rely on curved mesh elements, which fit curved boundaries. Such techniques often demand sophisticated meshing algorithms, cumbersome quadrature rules for integration, and complex nonlinear transformations to map the curved mesh elements onto the reference polygonal ones. In this regard, the reconstruction for off‐site data method, proposed in the work of Costa et al, provides very high‐order accurate polynomial reconstructions on arbitrary smooth curved boundaries, enabling integration of the governing equations on polygonal mesh elements, and therefore, avoiding the use of complex integration quadrature rules or nonlinear transformations. The method was introduced for Dirichlet boundary conditions and the present article proposes an extension for general boundary conditions, which represents an important advance for real context applications. A generic framework to compute polynomial reconstructions is also developed based on the least‐squares method, which handles general constraints and further improves the algorithm. The proposed methods are applied to solve the convection‐diffusion equation with a finite volume discretization in unstructured meshes. A comprehensive numerical benchmark test suite is provided to verify and assess the accuracy, convergence orders, robustness, and efficiency, which proves that boundary conditions on arbitrary smooth curved boundaries are properly fulfilled and the nominal very high‐order convergence orders are effectively achieved.

No abstract