Abstract

The use of simulation has spread to all areas of engineering and science, and the use of numerical models based on partial differential equations has thus multiplied. The resolution of these models is generally based on the discretization of the space in which the solutions to the equations under consideration are sought. The finite differences method or the finite elements method are two examples of such a discretization. This discretization simplifies the solving but implies a form of uncertainty on the value of any quantity of interest. To quantify this discretization uncertainty, the grid convergence index (GCI), based on the Richardson extrapolation technique, is now standard in the Verification and Validation (V&V) literature. But alternative approaches were also proposed in the statistical literature, such as Bayesian approaches with Gaussian process models. The objective of this work is to compare on a standard test case from the literature (Timoshenko's beam) the well-established GCI-based approach to the--younger--Bayesian approach for the quantification of discretization uncertainty.