Research Papers: Techniques and Procedures

Richardson Extrapolation-Based Discretization Uncertainty Estimation for Computational Fluid Dynamics

[+] Author and Article Information
Tyrone S. Phillips

Department of Aerospace and Ocean Engineering,
Virginia Tech,
Blacksburg, VA 24061
e-mail: tphilli6@vt.edu

Christopher J. Roy

Department of Aerospace and Ocean Engineering,
Virginia Tech,
Blacksburg, VA 24061
e-mail: cjroy@vt.edu

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received August 19, 2013; final manuscript received April 1, 2014; published online September 10, 2014. Assoc. Editor: Zvi Rusak.

J. Fluids Eng 136(12), 121401 (Sep 10, 2014) (10 pages) Paper No: FE-13-1505; doi: 10.1115/1.4027353 History: Received August 19, 2013; Revised April 01, 2014

This study investigates the accuracy of various Richardson extrapolation-based discretization error and uncertainty estimators for problems in computational fluid dynamics (CFD). Richardson extrapolation uses two solutions on systematically refined grids to estimate the exact solution to the partial differential equations (PDEs) and is accurate only in the asymptotic range (i.e., when the grids are sufficiently fine). The uncertainty estimators investigated are variations of the grid convergence index and include a globally averaged observed order of accuracy, the factor of safety method, the correction factor method, and least-squares methods. Several 2D and 3D applications to the Euler, Navier–Stokes, and Reynolds-Averaged Navier–Stokes (RANS) with exact solutions and a 2D turbulent flat plate with a numerical benchmark are used to evaluate the uncertainty estimators. Local solution quantities (e.g., density, velocity, and pressure) have much slower grid convergence on coarser meshes than global quantities, resulting in nonasymptotic solutions and inaccurate Richardson extrapolation error estimates; however, an uncertainty estimate may still be required. The uncertainty estimators are applied to local solution quantities to evaluate accuracy for all possible types of convergence rates. Extensions were added where necessary for treatment of cases where the local convergence rate is oscillatory or divergent. The conservativeness and effectivity of the discretization uncertainty estimators are used to assess the relative merits of the different approaches.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Banks, J. W., Aslam, T., and Rider, W. J., 2008, “On Sub-Linear Convergence for Linearly Degenerate Waves in Capturing Schemes,” J. Comput. Phys., 227(14), pp. 6985–7002. [CrossRef]
Roache, P. J., 1994, “Perspective: A Method of Uniform Reporting of Grid Refinement Studies,” ASME J. Fluid Eng., 116(3), pp. 405–413. [CrossRef]
Celik, I. B., Ghia, U., Roache, P. J., Freitas, C. J., Coleman, H. W., and Raad, P. E., 2008, “Procedure for Estimation and Reporting of Uncertainty Due to Discretization in CFD Applications,” ASME J. Fluid Eng., 130(7), pp. 78001–78005. [CrossRef]
Cosner, R. R., Oberkampf, W. L., Rumsey, C. L., Rahaim, C. P., and Shih, T. I.-P., 2006, “AIAA Committee on Standards for Computational Fluid Dynamics: Status and Plans,” Tech. Report No. AIAA-2006-889.
Roache, P. J., 2009, Fundamentals of Verification and Validation, Hermosa, Albuquerque, NM.
Logan, R. W., and Nitta, C. K., 2006, “Comparing 10 Methods for Solution Verification, and Linking to Model Validation,” JACIC, 3, pp. 354–373. [CrossRef]
Cadafalch, J., Perez-Segarra, C. D., Consul, R., and Oliva, A., 2002, “Verification of Finite Volume Computations on Steady-State Fluid Flow and Heat Transfer,” ASME J. Fluid Eng., 124(1), pp. 11–21. [CrossRef]
Xing, T., and Stern, F., 2010, “Factors of Safety for Richardson Extrapolation,” ASME J. Fluid Eng., 132(6), p. 61403. [CrossRef]
Oberkampf, W. L., and Roy, C. J., 2010, Verification and Validation in Scientific Computing, Cambridge University, Cambridge, UK.
Stern, F., Wilson, R. V., Coleman, H. W., and Paterson, E. G., 2001, “Comprehensive Approach to Verification and Validation of CFD Simulations—Part 1: Methodology and Procedures,” ASME J. Fluid Eng., 126(4), pp. 793–802. [CrossRef]
Phillips, T. S., Derlaga, J. M., and Roy, C. J., 2012, “Numerical Benchmark Solutions for Laminar and Turbulent Flows,” Tech. Report No. AIAA-2012-3074.
Roache, P. J., 2003, “Error Bars for CFD,” Tech. Report No. AIAA-2003-408.
Wilson, R., Shao, J., and Stern, F., 2004, “Discussion: Criticisms of the Correction Factor Verification Method,” ASME J. Fluid Eng., 126(4), pp. 704–706. [CrossRef]
Eça, L., and Hoekstra, M., 2002, “An Evaluation of Verification Procedures for CFD Applications,” 24th Symposium on Naval Hydrodynamics, Fukuoka, Japan, July 8–13.
Eça, L., and Hoekstra, M., 2009, “Evaluation of Numerical Error Estimation Based on Grid Refinement Studies With the Method of the Manufactured Solutions,” Comput. Fluids, 38(8), pp. 1580–1591. [CrossRef]
Eça, L., 2010, private communication.
Ainsworth, M., and Oden, J. T., 2000, A Posteriori Error Estimation in Finite Element Analysis, Wiley, New York. [CrossRef]
Phillips, T. S., and Roy, C. J., 2013, “A New Extrapolation-Based Uncertainty Estimator for Computational Fluid Dynamics,” Tech. Report No. AIAA-2013-0260.
Luke, E. A., Tong, X., Wu, J., and Cinnella, P., 2010, “Chem 3.2: A Finite-Rate Viscous Chemistry Solver – The User Guide,” Tetra Research Corporation Tech. Report.
Hoekstra, M., and Eça, L., 1998, “Parnassos: An Efficient Method for Ship Stern Flow Calculation,” Proc. 3rd Osaka Colloquium, Osaka, Japan, pp. 331–357.
Phillips, T. S., and Roy, C. J., 2011, “Residual Methods for Discretization Error Estimation,” Tech. Report No. AIAA-2011-3870.
Ollivier-Gooch, C., Nejat, A., and Michalak, K., 2009, “Obtaining and Verifying Higher-Order Unstructured Finite Volume Solutions to the Euler Equations,” AIAA J., 47(9), pp. 2105–2120. [CrossRef]
Satav, V., Hixon, R., Nallasamy, M., and Sawyer, S., 2005, Validation of a Computational Aeroacoustics Code for Nonlinear Flow About Complex Geometries Using Ringleb's Flow,” Tech. Report No. AIAA-2005-2871.
Eça, L., and Hoekstra, M., “Workshops on CFD Uncertainty Analysis,” http://maretec.inst.utl.pt
Veluri, S. P., 2010, “Code Verification and Numerical Accuracy Assessment for Finite Volume CFD Codes,” Ph.D. thesis, Virginia Tech, Blacksburg, VA.
Veluri, S. P., Roy, C. J., and Luke, E. A., 2012, “Comprehensive Code Verification Techniques for Finite Volume CFD Codes,” Comput. Fluids, 70(30), pp. 59–72. [CrossRef]
Rumsey, C. L., Smith, B. R., and Huang, G. P., 2010, “Description of a Website Resource for Turbulence Modeling, Verification and Validation,” Tech. Report No. AIAA-2010-4742.


Grahic Jump Location
Fig. 2

PARNASSOS Cartesian, stretched, and nonorthogonal grids

Grahic Jump Location
Fig. 1

Comparison of data for two different uncertainty estimators with quadratic fit and bounds for less scatter (top row) and more scatter (bottom row)

Grahic Jump Location
Fig. 3

Samples of the Loci-CHEM computational grids showing the (a) Cartesian tetrahedral grid, (b) curvilinear prismatic grid, (c) highly skewed curvilinear hexahedral grid, and (d) curvilinear hybrid grid

Grahic Jump Location
Fig. 4

Turbulent flat plate setup and the 69 × 49 grid

Grahic Jump Location
Fig. 5

Comparison of the effectivity index for Richardson extrapolation versus various distance measures

Grahic Jump Location
Fig. 6

Error effectivity indices for discretization error estimators

Grahic Jump Location
Fig. 7

Uncertainty effectivity index and conservativeness for discretization uncertainty estimators



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In