网格收敛·Richardson外推·GCI详解
🙋
什么是网格收敛?我听说CAE中网格越细计算结果越会改变,到底要细到什么程度呢?
🎓
网格收敛就是指:网格再细化计算结果基本不再改变的状态。比如汽车车盖变形计算,粗网格得到10mm,细网格得到12mm,再细化到12.1mm……这样的情况下,解就收敛到12mm附近了。这个模拟器中"维数"参数改为1D、2D、3D时,你会发现需要的网格数量增长方式完全不同。
🙋
原来如此!但是每次都用超细的网格计算太耗时了。能不能根据粗一点的结果来预测答案呢?
🎓
这正是Richardson外推法的作用!用3个不同粗度的网格计算,根据这3个结果推测"理论上无限细网格的答案"。在这个工具里输入粗、中、细三种网格下得到的关键数值(比如最大应力),就能计算出推估的"真实解"。实际工作中,这个推估值有时能让你省略掉最细网格的计算。
🙋
明白了!那这个推估有多可信呢?有没有数字标准来判断"已经收敛"?
🎓
有的,那就是GCI(网格收敛指数)。它表示"推估值和现在得到的解之间,如果再细化网格会改善百分之几"。工具里会根据GCI值用颜色表示:2%以下绿色、5%以下橙色、5%以上红色。实务中常见的做法是让GCI_fine降到3%以下来证明网格收敛。试着调整参数看GCI颜色的变化吧。
物理模型与主要公式
网格尺寸$h$与数值解$\phi$的关系,由理论精度阶数$p$按如下方式建模。这是Richardson外推的基础:
$$\phi = \phi_{exact}+ g h^p + \cdots$$
其中$\phi_{exact}$是真实解(连续体解),$g$是常数,$h$是代表网格尺寸。$p$是离散化格式的理论精度阶数(如中心差分为2阶),但实际问题复杂,所以用3个解计算的"表观精度阶数"代替。
从三个不同粗度的网格(细:1,中:2,粗:3)计算表观精度阶数$p$、外推解$\phi_{ext}^{21}$和GCI。
$$ p = \frac{1}{\ln(r_{21})}\left| \ln \left| \frac{\epsilon_{32}}{\epsilon_{21}}\right| + q(p) \right| $$
$$ \phi_{ext}^{21}= \phi_1 + \frac{\phi_1 - \phi_2}{r_{21}^p - 1}$$
$$ GCI_{fine}^{21}= F_s \frac{|\epsilon_{21}|}{r_{21}^p - 1}$$
$r_{21}=h_2/h_1$为细化比,$\epsilon_{21}=\phi_2-\phi_1$为解差,$F_s$为安全系数(通常1.25)。$q(p)$是当$p$未知时的修正项。工具基于这些公式,从输入的3个解自动计算$p$、外推解和GCI,评估收敛性。
常见问题
基本做法是用3段以上系统细化的网格,保持细化比$r$常数(如r=2)。例如粗网格单元数N、中为2N、细为4N。比值不均匀会导致表观精度阶数$p$计算不稳定,外推结果失真。
不一定。外推解只有在解处于渐近区(网格足够细的区域)时才可信。要确认GCI较小、表观精度阶数$p$接近理论值。渐近区外外推解可能偏离真值很远。
一般来说GCI < 5%即解不依赖网格。但因应用不同而异。关键是粗、中、细各网格的GCI应呈现逐级递减趋势。若GCI在某个阶段突然反弹,说明网格尚未进入渐近区。
可能原因:(1)网格还太粗,未进渐近区;(2)网格质量不佳(歪斜、宽高比差);(3)边界条件或物理模型有间断点;(4)数值格式未正确实现。首先应尝试进一步细化网格观察。
实际应用
结构强度分析(汽车、飞机):评估部件最大应力和变形时,必须消除网格依赖性,得到可信结果。轻量化设计中安全系数较小,GCI定量验证往往是设计决策的基础。
流体分析(CFD):翼型升力、阻力系数和热交换器换热、压降等预测,广泛用此法提升精度。复杂湍流模型时理论阶数$p$未知,靠计算的表观阶数判断收敛。
电磁场分析:天线辐射方向图、共振频率、器件寄生参数等网格敏感量,Richardson外推在保持计算成本低的同时实现高精度估计。
学术研究·论文发表:数值计算结果可信性的重要验证手段,已成为可重复研究的标配。建议文章中同时报告网格尺寸和GCI值。
常见误区与注意事项
第一点,"网格越细越接近真值"不总是对。几何奇异点(尖角、接触点)附近,精细网格也可能应力发散。此时应从"求最大应力值"改为"研究应力分布",重新定义CAE目标。
第二点,3个网格必须"系统地"细化。别随意选1mm、0.7mm、0.5mm,而要保持细化比$r$恒定(推荐$r \ge 1.3$),如1mm、0.77mm(1/1.3)、0.59mm(1/1.3²)。否则表观阶数$p$计算失稳,外推无信度。
第三点,不要被GCI迷蛊。GCI_fine < 2%不代表所有物理量都收敛。常见情况是位移收敛但应力未收敛。必须对关键输出量逐一验证收敛性。