模拟计算结果的精度验证是一个系统性工程,通常需要从数值层面和物理层面两个维度进行交叉检验。下面是主要的验证方法:
与实验数据或已知解对比
这是直接的验证手段,核心思路是用"已知答案"来检验模拟的准确性。
与实验数据对比:将仿真结果与物理实验或实测数据进行对比,计算均方根误差(RMSE)、相对误差等量化指标。一般而言,偏差在5%~10%以内且变化趋势一致,通常可视为可信。
与解析解/理论解对比:选择可简化的模型,使其存在已知的数学解析解(如悬臂梁弯曲、一维热传导、管道层流等),将仿真结果与理论解对比,误差在10%以内通常认为合理。
与高精度仿真或基准案例对比:将简化模型的结果与更精细的模型(如3D vs
2D)、行业内公认的标准测试案例、或更高保真度的仿真(如DNS)进行对比。
与经验公式/行业标准对比:将结果与成熟的理论公式、经验图表或行业设计规范进行交叉校验。
网格无关性验证(收敛性分析)
这是数值模拟中基础也关键的一步,目的是排除网格离散化带来的数值误差。
空间网格收敛:逐步加密网格(如从粗到细多套网格),观察关键输出参数(z大应力、流量、温度等)的变化。当继续加密网格后结果变化量小于1%~5%或趋于稳定时,认为达到网格无关性。
时间步长收敛:对于瞬态问题,逐步缩小时间步长,检查结果是否趋于稳定。
迭代收敛性:对于迭代求解算法,检查结果是否随迭代次数增加而收敛,收敛残差是否降至合理水平。
物理合理性检查
即使数值上收敛了,还需要判断结果是否"合理"。
守恒律检查:能量、质量、动量是否守恒(考虑流入流出边界),误差应在合理范围内。
物理趋势检查:应力集中是否发生在预期的几何不连续处?流体是否从高压流向低压?温度梯度方向是否符合热源/冷源位置?对称结构在对称载荷下结果是否对称?
量级合理性:关键结果的数值大小是否在预期的数量级范围内(如金属结构应力是几个MPa还是几百GPa)。
可视化检查:通过云图、动画等手段观察物理量分布,发现异常行为或不连续区域。
敏感性分析与不确定性量化
用于评估模型对输入参数变化的鲁棒性,以及结果的可信范围。
敏感性分析:系统地改变关键输入参数(材料属性、边界条件、摩擦系数、载荷大小等),观察输出结果的变化程度。高度敏感的参数是潜在的主要误差来源,需要更精确的输入数据。
不确定性量化:通过蒙特卡洛模拟、多项式混沌展开等概率方法,为输入参数定义概率分布,输出结果的置信区间(如±X%)而非单一预测值。
模型与参数审查
模型简化合理性:忽略的细节是否对关心的结果影响很小?简化假设是否明确且合理?
参数来源可靠性:材料参数、边界条件等是否来源权威且与实际系统一致?
单元类型与求解器选择:单元阶次是否足够?求解器类型、收敛容差、时间步长等设置是否合理?
重复性与一致性检查
重复性:在相同输入条件下,多次运行仿真,检查结果是否具有可重复性和一致性。
软件间交叉验证:使用不同的仿真软件对同一问题进行计算,对比结果是否一致。