Quantification of uncertainties in turbulence modeling: A comparison of physics-based and random matrix theoretic approaches

Numerical models based on Reynolds-Averaged Navier-Stokes (RANS) equations are widely used in engineering turbulence modeling. However, the RANS predictions have large model-form uncertainties for many complex flows, e.g., those with non-parallel shear layers or strong mean flow curvature. Quantification of these large uncertainties originating from the modeled Reynolds stresses has attracted attention in the turbulence modeling community. Recently, a physics-based Bayesian framework for quantifying model-form uncertainties has been proposed with successful applications to several flows. Nonetheless, how to specify proper priors without introducing unwarranted, artificial information remains challenging to the current form of the physics-based approach. Another recently proposed method based on random matrix theory provides the prior distributions with maximum entropy, which is an alternative for model-form uncertainty quantification in RANS simulations. This method has better mathematical rigorousness and provides the most non-committal prior distributions without introducing artificial constraints. On the other hand, the physics-based approach has the advantages of being more flexible to incorporate available physical insights. In this work, we compare and discuss the advantages and disadvantages of the two approaches on model-form uncertainty quantification. In addition, we utilize the random matrix theoretic approach to assess and possibly improve the specification of priors used in the physics-based approach. The comparison is conducted through a test case using a canonical flow, the flow past periodic hills. The numerical results show that, to achieve maximum entropy in the prior of Reynolds stresses, the perturbations of shape parameters in Barycentric coordinates are normally distributed. Moreover, the perturbations of the turbulence kinetic energy should conform to log-normal distributions. Finally, the result sheds light on how large the variance of each physical variable should be compared with each other to achieve the approximate maximum entropy prior. The conclusion can be used as a guidance for specifying proper priors in the physics-based, Bayesian uncertainty quantification framework.