缘起
今天(其实是去年,2020年,呸!)看到一篇文献,
[1] KHAYRULLINA A, VAN HOOFF T, BLOCKEN B, et al. Validation of steady RANS modelling of isothermal plane turbulent impinging jets at moderate Reynolds numbers [J]. European Journal of Mechanics - B/Fluids, 2019, 75(228-43).
作者研究了一种平面冲击射流(impinging jet),评估了不同的湍流模型在该应用中的性能。文中展示了一种网格无关性验证的方式,叫 GCI index。比较有趣,因为我以前并没有见过。今天我们来模仿一下。
先搜索一下 GCI index,在很多地方都能见到,比如以下这个归属 NASA 的网站和来自 ARFM 的文献[1][2],所以应该是一种比较公认的、靠谱的网格无关性验证方法。
https://www.grc.nasa.gov/www/wind/valid/tutorial/spatconv.htmlwww.grc.nasa.govROACHE P J. QUANTIFICATION OF UNCERTAINTY IN COMPUTATIONAL FLUID DYNAMICS [J]. Annual Review of Fluid Mechanics, 1997, 29(1): 123-60.
动手
先分析一下上面那张图,发现其实 GCI index 是不同网格之间结果的误差。这就和我之前的研究很像了。可以看看下面这篇文章的内容。
https://zhuanlan.zhihu.com/p/93112564zhuanlan.zhihu.com那末,我们首先构造指示 GCI index 阴影图像的函数,老规矩用 MATLAB 快速实现一下,写出如下函数。参数依次为 axes 对象,数据,颜色。
function shaded_error_bar(ax