论文标题
评估合奏方法,用于量化具有小合奏尺寸的稳态CFD应用中的不确定性
Evaluation of ensemble methods for quantifying uncertainties in steady-state CFD applications with small ensemble sizes
论文作者
论文摘要
贝叶斯不确定性量化(UQ)对行业和学术界感兴趣,因为它提供了一个框架,可通过合并可用数据来量化和降低计算模型中的不确定性。例如,对于计算成本非常高的系统,计算流体动力学(CFD)问题,传统的,精确的贝叶斯方法(例如马尔可夫链蒙特卡洛)是棘手的。为此,基于合奏的贝叶斯方法已用于CFD应用程序。但是,到目前为止,它们对UQ的适用性尚未得到充分的分析和理解。在这里,我们评估了三种广泛使用的基于迭代集合的数据同化方法的性能,即集合卡尔曼滤波器,集合随机最大似然方法以及与UQ问题多个数据同化的集合Kalman滤波器。我们从优化的角度介绍了三种集合方法的推导。此外,使用标量案例来证明三种不同方法的性能,重点是小合奏尺寸的影响。最后,我们评估了三种集合方法,用于量化涉及湍流平均流量的稳态CFD问题中的不确定性。具体而言,Reynolds平均Navier-Stokes(RANS)方程被视为正向模型,并且通过合并观察数据来量化传播速度中的不确定性。结果表明,整体方法无法准确捕获真正的后验分布,但是即使使用非常有限的集合尺寸,它们也可以很好地估算不确定性。根据比较的总体性能和效率,在此处评估的三种集合方法中,合奏随机最大似然方法被确定为近似贝叶斯UQ方法的最佳选择。
Bayesian uncertainty quantification (UQ) is of interest to industry and academia as it provides a framework for quantifying and reducing the uncertainty in computational models by incorporating available data. For systems with very high computational costs, for instance, the computational fluid dynamics (CFD) problem, the conventional, exact Bayesian approach such as Markov chain Monte Carlo is intractable. To this end, the ensemble-based Bayesian methods have been used for CFD applications. However, their applicability for UQ has not been fully analyzed and understood thus far. Here, we evaluate the performance of three widely used iterative ensemble-based data assimilation methods, namely ensemble Kalman filter, ensemble randomized maximum likelihood method, and ensemble Kalman filter with multiple data assimilation for UQ problems. We present the derivations of the three ensemble methods from an optimization viewpoint. Further, a scalar case is used to demonstrate the performance of the three different approaches with emphasis on the effects of small ensemble sizes. Finally, we assess the three ensemble methods for quantifying uncertainties in steady-state CFD problems involving turbulent mean flows. Specifically, the Reynolds averaged Navier--Stokes (RANS) equation is considered the forward model, and the uncertainties in the propagated velocity are quantified and reduced by incorporating observation data. The results show that the ensemble methods cannot accurately capture the true posterior distribution, but they can provide a good estimation of the uncertainties even when very limited ensemble sizes are used. Based on the overall performance and efficiency from the comparison, the ensemble randomized maximum likelihood method is identified as the best choice of approximate Bayesian UQ approach~among the three ensemble methods evaluated here.