绘图添加网格_自用小玩具：一种网格无关性的验证方法

缘起

今天（其实是去年，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.html​www.grc.nasa.gov

ROACHE 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/93112564​zhuanlan.zhihu.com

那末，我们首先构造指示 GCI index 阴影图像的函数，老规矩用 MATLAB 快速实现一下，写出如下函数。参数依次为 axes 对象，数据，颜色。

function shaded_error_bar(ax