收敛速度 — CAE术语解说
收敛速度
老师,收敛速度是指网格细化时精度提升的方式吗?
是的。当网格的代表尺寸h减半时,误差减少的程度——这就是收敛速度(Convergence Rate)。误差 ∝ h^p 中的指数p称为"收敛阶数",p=1(一阶精度)时,h减半误差也减半;p=2(二阶精度)时,误差减为1/4。FEM的线性一阶单元通常是p=2(能量范数),二阶单元是p=3或更高。验证这个收敛阶数是否实际出现,就是"网格收敛分析",这是验证(Verification)的基础。
定义
实际上怎样确认收敛阶数呢?
逐步细化网格(例如:h=1、0.5、0.25),计算各分辨率下的应力、速度等目标量。绘制图表时,纵轴为误差(与参考解的差异),横轴为h(对数标度),如果成为直线则幂律关系成立。斜率就是p。当参考解不存在时(大多数实际问题),使用Richardson外推法(从两个网格阶段的结果推估真值的方法)。GCI(Grid Convergence Index)法是IAEA和航空航天规范推荐的标准方法。
空间与时间收敛的区别
时间方向也有收敛速度吗?
有的。非定常分析中,除了空间网格尺寸h,时间步长dt也影响收敛。隐式欧拉法(一阶精度)时,dt减半则时间误差减半;Crank-Nicolson法(二阶精度)时误差减为1/4。OpenFOAM的icoFoam有Crank-Nicolson选项,但默认是ddtSchemes: Euler(一阶)。时间精度成为问题的情况是在LES、DNS分析中,当需要精确追踪涡的脱落频率时,此时需要使空间和时间的收敛阶数相匹配,否则会出现"即使细化空间网格也被时间误差主导"的情况。
实际工作中网格收敛确认要做到什么程度?
不需要在所有分析中都进行完整的网格收敛分析。实际判断标准是——①量产设计的筛选分析:可省略网格收敛确认,②重要部件的强度验证、认证试验补充:通常需要用GCI法进行收敛确认(如航空航天的FAA/EASA认证),③新模型类型或边界条件的初次分析:必须确认——这样的感觉。"细化网格后结果在1%以内稳定"是现场的经验性目标,这种状态通常被称为"充分收敛"。
收敛阶数低于理论值时会有什么问题?
网格可能存在某些问题的信号。例如——存在尖锐的角(应力奇点)时无法达到理论高阶精度,单元宽高比极端不均匀,接触或材料不连续处网格不匹配等原因都可能导致问题。还有一种是解本身包含奇点的情况——裂纹尖端的应力在理论上以1/sqrt(r)发散,无论网格多细化,普通单元也无法出现h^2的收敛。对于具有奇异性的问题需要特殊单元(如Singular quarter-point element)。
关联术语
GCI法和收敛阶数的确认方法非常实用易懂!也不知道奇点的影响。
CAE术语的准确理解是团队内沟通的基础。 — Project NovaSolver也考虑支持实务工作者的学习。
CAE的未来,与实务工作者一起思考
Project NovaSolver致力于直面收敛速度等实务课题的本质,构建支撑工程现场的工具。
查看项目最新信息 →相关话题
价值
更详细
错误