Error quantification in computational fluid dynamics is a subject of increasing interest and research. Although solution verification is conventionally performed using systematic grid refinement, Richardson extrapolation can be restrictive in its applicability, possibly requiring more than three grids to establish monotonic behavior. Ultimately, a single grid error estimation method may prove of greater utility for practical application in computational aerodynamics. In this work, an error transport equation is implemented in a three-dimensional upwind unstructured Navier-Stokes solver. An approach for deriving errors in related quantities of interest from the solution of the error transport equation is presented. Error quantification is demonstrated in two and three dimensions for aerodynamic flows where experimental data are available for comparison. Predicted error bars are found to contain fine grid solutions, test data, and the results of Richardson extrapolation. Furthermore, the error transport equation provides meaningful error quantification for aerodynamic coefficients in cases where Richardson extrapolation cannot be applied.
Nomenclatureverticeŝ n = unit vector normal to boundary surface P = pressure p = order of convergence exponent Q = conserved variable vector R = error transport equation residual vector R = derived result r = refinement ratiô t = unit vector tangential to boundary surface U = velocity magnitude u, v, w = velocity components V = volume x m , y m , z m = coordinates of moment center = angle of attack = ratio of specific heats F = upwind Roe flux vector f = upwind Roe flux in one dimension t = time step x, y = grid spacings " = error transport equation solution vector " = error = viscosity = density = shear stress Subscripts extrap = extrapolated i = grid point n = normal node REF = reference t = turbulent w = wall node x, y, z = components in the x, y, and z directions 1 = freestream Superscripts h = discretized solution of mesh width h L = left state n = time level R = right state