Due to high performance of modern computers, Uncertainty Quantification is becoming an important part of engineering design. Every non intrusive Uncertainty Quantification method requires a considerable number of evaluations of the model, meaning that the design process is more expensive in terms of computational resources/time. In Computational Fluid Dynamics, the usual practice is to reduce the computational time by reducing the number of nodes of the used mesh. Each coarsening of the mesh leads to the increase of the error measured as the difference between the real solution and the solution provided by the computational model. In this work, an approach for quantification of the global error around the stochastic domain, in a mesh reduction process, is described and results obtained for a test case are detailed. The method is based on a comparison of the high accurate mesh against coarse mesh with lower accuracy, but less expensive in terms of computational time. The global error is defined as a volume difference between surrogate models created in the stochastic domain. The stochastic domain is given by pre-specified input variables with appropriate boundaries. Surrogate models are used and a non intrusive polynomial chaos model is created with response samples from high and low accuracy mesh. For the chosen test case, the input variables, related to the stochastic space, were the free stream pressure and free stream Mach number. A hypersonic flow solver developed at the von Karman Institute, Cosmic, was used to compute properties of a flow around the reentry spacecraft. A computational expensive mesh was used as a reference mesh. Due to computational resources, it was impossible to use expansive mesh for Monte Carlo simulation or high order Polynomial Chaos. Therefore, the global error estimation approach was applied to find an accurate and relatively inexpensive mesh for Uncertainty Quantification in hypersonic simulation. Multiple meshes with different coarsening were tested, based on expert knowledge of the problem. The global error estimation method allowed for finding a final mesh, with an error on the mean value 0.48% and on the standard deviation 5.89%, which was 4 times faster than the reference mesh