  A machine-learning approach for data-driven Reynolds-Averaged Navier--Stokes
(RANS) predictions of turbulent flows including estimates of turbulence
modelling uncertainties is developed by combining a Bayesian symbolic
identification methodology for learning customised RANS model corrections for
selected classes of flows and a space-dependent model-aggregation algorithm
that combines the predictions of a set of competing machine-learned RANS models
by means of weighting functions depending on a vector of local flow features.
The customised model corrections are learned by using the SBL-SpaRTA algorithm,
recently proposed by \citet{cherroud2022sparse}, which delivers sparse model
correction terms in analytical form and whose parameters are described by
probability distribution functions. This makes the learned models naturally
interpretable and endowed with a measure of uncertainty. The learned models are
subsequently aggregated by training Random Forests Regressors (RFR), which
associates a model performance score with a set of local flow features. The
scores can be interpreted as the probability that a candidate model performs
better than its competitors, given the flow behavior at a given location.
Predictions of new flows are then formulated as a locally weighted average of
the solutions of a set of machine-learned models. An uncertainty measure is
obtained by propagating through the models their posterior parameter
distributions and by combining competing model predictions to estimate the
inter-model variance. The proposed space-dependent model aggregation procedure
(X-MA) is applied to flows of varying complexity, showing significant
improvement with respect to the baseline model and good generalizability to
unseen flows. The results make X-MA an attractive candidate for the development
of a generalizable data-driven turbulence modelling framework with quantified
uncertainty.
