Extrapolation-based Discretization Error and Uncertainty Estimation in Computational Fluid Dynamics

The solution to partial differential equations generally requires approximations that result in numerical error in the final solution. Of the different types of numerical error in a solution, discretization error is the largest and most difficult error to estimate. In addition, the accuracy of the discretization error estimates relies on the solution (or multiple solutions used in the estimate) being in the asymptotic range. The asymptotic range is used to describe the convergence of a solution, where an asymptotic solution approaches the exact solution at a rate proportional to the change in mesh spacing to an exponent equal to the formal order of accuracy. A non-asymptotic solution can result in unpredictable convergence rates introducing uncertainty in discretization error estimates. To account for the additional uncertainty, various discretization uncertainty estimators have been developed. The goal of this work is to evaluation discretization error and discretization uncertainty estimators based on Richardson extrapolation for computational fluid dynamics problems. In order to evaluate the estimators, the exact solution should be known. A select set of solutions to the 2D Euler equations with known exact solutions are used to evaluate the estimators. Since exact solutions are only available for trivial cases, two applications are also used to evaluate the estimators which are solutions to the Navier-Stokes equations: a laminar flat plate and a turbulent flat plate using the turbulence model. Since the exact solutions to the Navier-Stokes equations for these cases are unknown, numerical benchmarks are created which are solutions on significantly finer meshes than the solutions used to estimate the discretization error and uncertainty. Metrics are developed to evaluate the accuracy of the error and uncertainty estimates and to study the behavior of each estimator when the solutions are in, near, and far from the asymptotic range. Based on the results, general recommendations are made for the implementation of the error and uncertainty estimators. In addition, a new uncertainty estimator is proposed with the goal of combining the favorable attributes of the discretization error and uncertainty estimators evaluated. The new estimator is evaluated using numerical solutions which were not used for development and shows improved accuracy over the evaluated estimators.

[1]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[2]  D. F. Mayers,et al.  The deferred approach to the limit in ordinary differential equations , 1964, Comput. J..

[3]  Hugh W. Coleman,et al.  Comprehensive Approach to Verification and Validation of CFD Simulations—Part 1: Methodology and Procedures , 2001 .

[4]  C. L. Rumsey,et al.  Application of FUN3D and CFL3D to the Third Workshop on CFD Uncertainty Analysis , 2013 .

[5]  M. Nallasamy,et al.  Validation of a Computational Aeroacoustics Code for Nonlinear Flow about Complex Geometries Using Ringleb's Flow , 2005 .

[6]  L. Richardson.  The Approximate Arithmetical Solution by Finite Differences of Physical Problems Involving Differential Equations, with an Application to the Stresses in a Masonry Dam , 1911 .

[7]  Carl Ollivier-Gooch,et al.  Obtaining and Verifying High-Order Unstructured Finite Volume Solutions to the Euler Equations , 2009 .

[8]  L. Eça,et al.  Evaluation of numerical error estimation based on grid refinement studies with the method of the manufactured solutions , 2009 .

[9]  F. Menter.  Improved two-equation k-omega turbulence models for aerodynamic flows , 1992 .

[10]  Christopher J. Roy,et al.  Verification of Euler/Navier–Stokes codes using the method of manufactured solutions , 2004 .

[11]  William J. Rider,et al.  On sub-linear convergence for linearly degenerate waves in capturing schemes , 2008, J. Comput. Phys..

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

[13]  Patrick Roache,et al.  Error Bars for CFD , 2003 .

[14]  Josef Stoer,et al.  Numerische Mathematik 1 , 1989 .

[15]  Christopher J. Roy,et al.  Verification and Validation in Scientific Computing: Planning, management, and implementation issues , 2010 .

[16]  Patrick J. Roache,et al.  Criticisms of the “Correction Factor” Verification Method , 2003 .

[17]  E. Fehlberg,et al.  Low-order classical Runge-Kutta formulas with stepsize control and their application to some heat transfer problems , 1969 .

[18]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[19]  J. Cadafalch,et al.  Verification of Finite Volume Computations on Steady-State Fluid Flow and Heat Transfer , 2002 .

[20]  S. Priebe,et al.  Direct Numerical Simulation of Shockwave and Turbulent Boundary Layer Interactions , 2009 .

[21]  I. Babuska,et al.  A‐posteriori error estimates for the finite element method , 1978 .

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

[23]  Hugh W. Coleman,et al.  Comprehensive Approach to Verification and Validation of CFD Simulations—Part 2: Application for Rans Simulation of a Cargo/Container Ship , 2001 .

[24]  P. Roache.  Perspective: A Method for Uniform Reporting of Grid Refinement Studies , 1994 .

[25]  J. Oden,et al.  A Posteriori Error Estimation in Finite Element Analysis , 2000 .

[26]  Patrick J. Roache,et al.  Verification and Validation in Computational Science and Engineering , 1998 .

[27]  Subrahmanya P. Veluri.  Code Verification and Numerical Accuracy Assessment for Finite Volume CFD Codes , 2010 .

[28]  T. Xing,et al.  Factors of Safety for Richardson Extrapolation , 2010 .

[29]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[30]  Jun Shao,et al.  Discussion: Criticisms of the “Correction Factor” Verification Method 1 , 2004 .

[31]  J. Schetz.  Boundary Layer Analysis , 1992 .

[32]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[33]  Timothy G. Trucano,et al.  Verification and validation. , 2005 .

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

[35]  D. Wilcox.  Turbulence modeling for CFD , 1993 .

[36]  Christopher J. Roy,et al.  Review of code and solution verification procedures for computational simulation , 2005 .

[37]  P. Roache.  Verification of Codes and Calculations , 1998 .

[38]  Brian R. Smith,et al.  Description of a Website Resource for Turbulence Modeling Verification and Validation , 2010 .