Very high-order accurate finite volume scheme on curved boundaries for the two-dimensional steady-state convection-diffusion equation with Dirichlet condition

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.