Finite Volume Solution Reconstruction Methods For Truncation Error Estimation

The numerical solution to differential equations results in a discrete solution space for the finite volume and finite difference discretization methods. For various reasons, it can be necessary to prolong the solution from a discrete space to a continuous space. The prolongation to a continuous space can be done using various curve-fitting methods which adds an additional level of approximation to the solution. The allowable error of a prolongation operation depends on the specific task required by the user. In this paper we investigate various prolongation methods and identify the minimal requirements specifically for the purpose of truncation error estimation (the difference between the discrete and integral governing equations) for finite volume methods. The reconstruction methods investigated will include k-exact and ENO methods. Truncation error estimation for 1D Burgers’ equation and the k-exact method suggest that the minimum polynomial order is dependent on the highest derivatives in the truncation error expression and, therefore, the discretization scheme. The effect of different reconstruction methods on truncation error estimation is investigated and the minimum polynomial order for accurate truncation error estimation is identified for the Euler equations and is found to be second-order for the weak formulation and third-order for the strong formulation.

[1]  Timothy J. Barth,et al.  Recent developments in high order K-exact reconstruction on unstructured meshes , 1993 .

[2]  Curtis R. Mitchell,et al.  Practical aspects of spatially high-order accurate methods , 1993 .

[3]  Antony Jameson,et al.  Aerodynamic design via control theory , 1988, J. Sci. Comput..

[4]  Chi-Wang Shu Numerical experiments on the accuracy of ENO and modified ENO schemes , 1990 .

[5]  Tom I-P. Shih,et al.  A discrete transport equation for error estimation in CFD , 2002 .

[6]  C. W. Clenshaw,et al.  A method for numerical integration on an automatic computer , 1960 .

[7]  Dominique Pelletier,et al.  VERIFICATION OF ERROR ESTIMATORS FOR THE EULER EQUATIONS , 2000 .

[8]  P. Frederickson,et al.  Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction , 1990 .

[9]  S. Osher,et al.  Uniformly High-Order Accurate Nonoscillatory Schemes. I , 1987 .

[10]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[11]  D. Venditti,et al.  Grid adaptation for functional outputs: application to two-dimensional inviscid flows , 2002 .

[12]  H. Stetter The defect correction principle and discretization methods , 1978 .

[13]  R. Skeel Thirteen ways to estimate global error , 1986 .

[14]  Wanai Li,et al.  High‐order k‐exact WENO finite volume schemes for solving gas dynamic Euler equations on unstructured grids , 2012 .

[15]  C. J. Roy Strategies for Driving Mesh Adaptation in CFD (Invited) , 2009 .

[16]  Curtis R. Mitchell,et al.  K-exact reconstruction for the Navier-Stokes equations on arbitrary grids , 1993 .

[17]  Christopher J. Roy,et al.  SENSEI Computational Fluid Dynamics Code: A Case Study in Modern Fortran Software Development , 2013 .

[18]  V. Pereyra On improving an approximate solution of a functional equation by deferred corrections , 1966 .

[19]  Gecheng Zha,et al.  Improved Seventh-Order WENO Scheme , 2010 .