9.2 验证与确认


文档摘要

9.2 验证与确认 收敛只是过程合格,本节处理结果审计:验证(verification)回答"方程解对了吗"——网格无关性与离散误差带;确认(validation)回答"方程对了吗"——与实验或理论解对表。两问合起来,才构成数字风洞验收单的最后一栏签字。这是全书的方法论顶点:前面七章的每个配置,都会在这一节被问一句"凭什么信"。 验证:网格收敛性研究与 GCI 离散误差是数值解与精确解(同一方程)之差,它的标准审讯方式是网格收敛性研究:三套系统加密的网格各算一遍,观察关键量的收敛行为,再用 Richardson 外推估计精确解并给出网格收敛指数(GCI)。流程固定五步: 第一步,选关键量:验收单签字的量(阻力、压降、换热系数),不是全场云图。

9.2 验证与确认

收敛只是过程合格,本节处理结果审计:验证(verification)回答"方程解对了吗"——网格无关性与离散误差带;确认(validation)回答"方程对了吗"——与实验或理论解对表。两问合起来,才构成数字风洞验收单的最后一栏签字。这是全书的方法论顶点:前面七章的每个配置,都会在这一节被问一句"凭什么信"。

验证:网格收敛性研究与 GCI

离散误差是数值解与精确解(同一方程)之差,它的标准审讯方式是网格收敛性研究:三套系统加密的网格各算一遍,观察关键量的收敛行为,再用 Richardson 外推估计精确解并给出网格收敛指数(GCI)。流程固定五步:

第一步,选关键量:验收单签字的量(阻力、压降、换热系数),不是全场云图。第二步,定加密比:相邻网格的特征尺寸比 r 取一点三左右,网格量约各差一倍以上,三套起步。第三步,配齐各套网格的口径:边界层、加密区按比例整体缩放,只加密一个区域不叫网格研究。第四步,同收敛判据各算一遍。第五步,计算观察收敛阶 p 与 GCI:相邻解的差异 epsilon 加外推,GCI 约等于安全系数(工程取一点二五)乘相对误差除以(r 的 p 次方减一)。工程验收惯例:关键量的 GCI 低于百分之一到五视为网格无关,报告里同时给出误差带。

GCI 的正确用法是当误差带用:"压降为十二点四千帕,GCI 百分之二点一"比"压降为十二点四千帕"多出的不只是数字,是可辩护性——评审问"网格够不够细",答案是一个量化的区间而不是"应该够了吧"。

图 9-2:验证与确认的分工

图 9-2:验证与确认的分工

确认:与真实世界对账

确认的对表对象按可信度排序:专为本任务做的实验(金标准)、公开发表的同构型实验数据、经典理论解(充分发展管流的泊肃叶解、绕圆柱的阻力系数-雷诺数曲线)、行业基准算例库。对表的纪律有三条。其一,对表量选"模型敏感量":用湍流模型就该对分离点与压力分布,而不是只对总阻力。其二,偏差要拆账:数值误差带(GCI 给出)从总偏差里扣除,剩下的才是模型误差——不拆账会把数值问题错诊成模型问题。其三,适用范围声明:确认在什么雷诺数、什么工况段成立,结论不外推出验证过的区间。

完整算例:方管绕流的验证确认全流程

背景。新人培训科目:二维方柱绕流,计算斯特劳哈尔数与阻力系数,与经典实验数据对表。麻雀虽小,两关全有。

操作与结果。验证段:三套网格(粗中细,单元数约按一比二点二比五加密,r 约一点三),阻力系数分别为一点五四、一点四九、一点四七;观察收敛良好,外推估计约一点四六,最细网格 GCI 约百分之三点二——误差带可接受,但中等网格与细网格差距已小于 GCI,生产用中等网格即可,省下的算力投入瞬态采样。确认段:阻力系数一点四七对文献值约一点四七至一点五零(随雷诺数区间),偏差在百分之一内;斯特劳哈尔数零点一三二对文献零点一二至零点一六区间,落位合理;旋涡脱落频率的瞬态统计用了超过一百个脱落周期,统计误差随采样延长收敛。

解读。这个科目的教学价值在"流程完整性":误差带与结论同交付、数值与模型误差分账、采样充分性论证。放到工程里,同一套流程只是把方柱换成散热器或翼型,纪律不变——确认不是一次性动作,而是每个新构型开工前的标定环节。

变式。缺乏实验数据的全新产品,确认退而求其次:部件级对表(换热器芯子有数据、整机没有,先对芯子)加跨方法互证(同一子系统用不同模型或网格拓扑复算,结果互相咬合提高信心)。跨方法互证的证据强度弱于实验,报告里要如实降级表述。

验证确认问答与实施表

三套网格之外的工况也要做网格研究吗。不需要全做,但要"锚点工况":在代表性最不利的工况上做完整 GCI,其余工况沿用网格口径并抽查——全做是浪费,不做是裸奔。

GCI 算出来是负的或超预期怎么办。非单调收敛说明网格族之间不成系统:加密不均匀、某套网格质量掉档、或解处于非稳态被稳态硬解。先修网格族的一致性,再谈 GCI。

实验对表差多少能接受。没有统一数字,接受度由误差带与决策敏感度共同决定:结论翻转的临界差值之内必须较真,远离临界差值的偏差按工程精度要求放行——先定"多大算错",再谈"错了多少"。

验证确认实施清单:

动作 输入 输出物
网格收敛性研究 三套以上加密网格 观察收敛阶与 GCI 报告
时间步收敛性 瞬态步长减半复核 统计量差异声明
实验对表 同工况测量数据 偏差拆解(数值与模型)
基准算例复现 公开标准工况 能力证书式记录

关于误差带的使用再补一个评审实战技巧:给出误差带的同时给出"结论对误差带的稳健性"。举例:两方案压降差百分之八、各自 GCI 百分之二——差值大于误差带,排序结论稳健;若差值只有百分之一点五而误差带各百分之二,排序结论就不该进报告,只能写"两方案在数值误差带内无差异"。同一套数据,稳健性声明决定它能不能支撑决策。这也是为什么成熟的仿真报告敢写"无法区分"——敢于承认误差带内的平局,正是数字风洞区别于云图展示的标志。

关于验证资料的沉淀再补一段:网格收敛性研究是典型的"一次性投入、多次复利"的资产。第一次为一个产品线做 GCI 时,你得到的不只是一个误差带,还有"网格口径与误差的对应关系"——同系列几何在相近雷诺数下,这个关系可以复用,后续项目直接按目标误差反推网格规模,省去整套网格研究。把它做成一页"网格规模-误差"速查表放进模板库,产品线里每个新项目开工时先查表再划网格。验证的最高形态不是每个项目都做全流程,而是让历史验证替未来项目付账——这正是第九章纪律与第十一章组织能力的接缝。

还有一个关于对表数据的常见坑要单独提醒:确认用的实验数据自身也有不确定度。风洞的阻塞修正、传感器的标定期、试验工况与仿真工况的微小错位(温度差几度、雷诺数差百分之几),都会以"模型误差"的面目出现在对表偏差里。规范的报告在引用实验数据时同时声明其不确定度,并把仿真误差带与实验误差带合并成总偏差区间——两带不交叠才是真正意义的"对不上"。这个习惯把"我们和实验差了百分之五"这类无主纠纷,变成"仿真带加实验带覆盖了差异"或"确实存在未解释偏差"的清晰判断。对表的对,是对到误差带上,不是对到小数点上。

本节验收点

  • 关键工程量完成三套网格以上的收敛性研究,GCI 误差带入报告。
  • 确认对表选了模型敏感量,数值与模型误差分账拆解。
  • 结论带适用范围声明,未外推出验证区间。
  • 无实验数据任务执行了部件级对表或跨方法互证,证据等级如实标注。

作者与出处
原作者: 灏天文库
来源:灏天文库
整理: 灏天文库整理
由灏天文库平台收录,内容或由平台用户上传,仅供学习交流
发布者: 作者: 灏天文库 转发
评论区 (0)
U