Accelerated Bayesian calibration and uncertainty quantification of RANS turbulence model parameters for stratified atmospheric boundary layer flows

TL;DR

This paper introduces a machine learning-accelerated Bayesian inversion method to calibrate and quantify uncertainty in turbulence model parameters for atmospheric boundary layer flows, improving prediction accuracy and uncertainty estimates.

Contribution

It presents a novel Bayesian calibration approach for RANS turbulence models using machine learning acceleration, with comprehensive uncertainty quantification and validation against LES data.

Findings

Uncertainty in key parameters can be reduced by up to 84% through targeted sampling.

Bayesian learned parameters outperform standard deterministic parameters in out-of-sample tests.

Incorporating additional flow quantities further reduces model uncertainty.

Abstract

In operational weather models, the effects of turbulence in the atmospheric boundary layer (ABL) on the resolved flow are modeled using turbulence parameterizations. These parameterizations typically use a predetermined set of model parameters that are tuned to limited data from canonical flows. Using these fixed parameters results in deterministic predictions that neglect uncertainty in the unresolved turbulence processes. In this study, we perform a machine learning-accelerated Bayesian inversion of a single-column model of the ABL. This approach is used to calibrate and quantify uncertainty in model parameters of Reynolds-averaged Navier-Stokes turbulence models. Following verification of the uncertainty quantification methodology, we learn the parameters and their uncertainties in two different turbulence models conditioned on scale-resolving large-eddy simulation data over a range of ABL stabilities. We show how Bayesian inversion of a numerical model improves flow predictions by investigating the underlying mean momentum budgets. Further, we show that uncertainty quantification based on neutral surface layer data recovers the relationships between parameters derived from theoretical modeling, but that learning the parameters based on stable ABL data or data from outside the surface layer can lead to different relationships than neutral surface layer theory. Systematic uncertainty reduction methods reveal that (1) sampling wind speed up to the ABL height can reduce uncertainty in key model parameters by up to 84%, and (2) assimilating fluid flow quantities beyond first-order moment statistics can further reduce uncertainty in ways that wind speed assimilation alone cannot achieve. The parameters learned using Bayesian uncertainty quantification generally yield lower error than standard deterministic parameters in out-of-sample tests and provide uncertainty intervals on predictions.