Towards higher order discretization error estimation by error transport using unstructured finite-volume methods for unsteady problems

2017 ◽  
Vol 154 ◽  
pp. 245-255 ◽  
Author(s):  
G.K. Yan ◽  
C. Ollivier-Gooch
2020 ◽  
Author(s):  
Lingfa Kong ◽  
Yidao Dong ◽  
Wei Liu ◽  
Huaibao Zhang

Abstract Accuracy of unstructured finite volume discretization is greatly influenced by the gradient reconstruction. For the commonly used k-exact reconstruction method, the cell centroid is always chosen as the reference point to formulate the reconstructed function. But in some practical problems, such as the boundary layer, cells in this area are always set with high aspect ratio to improve the local field resolution, and if geometric centroid is still utilized for the spatial discretization, the severe grid skewness cannot be avoided, which is adverse to the numerical performance of unstructured finite volume solver. In previous work, we explored a novel global-direction stencil and combine it with face-area-weighted centroid on unstructured finite volume methods from differential form to realize the skewness reduction and a better reflection of flow anisotropy. Note, however, that the differential form is hard to achieve higher-order accuracy, and in order to set stage for the method promotion on higher-order numerical simulation, in this research, we demonstrate that it is also feasible to extend this novel method to the unstructured finite volume discretization in integral form. Numerical examples governed by linear convective, Euler and Laplacian equations are utilized to examine the correctness as well as effectiveness of this extension. Compared with traditional vertex-neighbor and face-neighbor stencils based on the geometric centroid, the grid skewness is almost eliminated and computational accuracy as well as convergence rate is greatly improved by the global-direction stencil with face-area-weighted centroid. As a result, on unstructured finite volume discretization from integral form, the method also has a better numerical performance.


Author(s):  
Gary K. Yan ◽  
Carl Ollivier-Gooch

A numerical estimation of discretization error for steady compressible flow solutions is performed using the error transport equation (ETE). There is a deficiency in the literature for obtaining efficient, higher order accurate error estimates for finite volume discretizations using nonsmooth unstructured meshes. We demonstrate that to guarantee sharp, higher order accurate error estimates, one must discretize the ETE to a higher order than the primal problem, a requirement not necessary for uniform meshes. Linearizing the ETE can limit the added cost, rendering the overall computational time competitive, while retaining accuracy in the error estimate. For the Navier–Stokes equations, when the primal solution is corrected using this error estimate, for the same level of solution accuracy the overall computational time is more than two times faster compared to solving the higher order primal problem. In addition, our scheme has robustness advantages, because we solve the primal problem only to lower order.


Sign in / Sign up

Export Citation Format

Share Document