第172章 数据的复仇:张怀民的建模风暴
清晨五点四十三分,旧车间的灯亮了。
任昭站在控制台前,手里拿着沈知遥发来的电镜报告。纸页边缘有折痕,是他昨晚反复翻看留下的。他把报告放在桌上,旁边是陈老给的那张编号纸条和苏蔓项目组的原始数据拷贝。三支圆珠笔并排别在左口袋,黑色那支刚换过笔芯。
门被推开,张怀民走进来。他背着一台笔记本电脑,屏幕贴着防窥膜。他没说话,走到角落的折叠桌前坐下,打开设备,调出昨晚运行的模型结果。界面显示附着力预测偏差为8%。
“我还是按标准CFD流程算的。”他说,“边界条件都复核过,没有错误。”
任昭走过去,看了一眼屏幕。曲线走势与实测数据明显错位,特别是在进气道肩部区域。
“你用了什么湍流模型?”任昭问。
“SST k-omega。”
“粘性底层处理方式?”
“自动壁面函数。”
任昭点头。这些选择都没问题。正因如此,偏差才更值得警惕。
“你觉得是实验误差?”他问。
“不排除。”张怀民说,“样本表面氧化层会影响测试值。”
“我们已经修正了氧化层影响。”任昭说,“沈知遥的数据是清洗后重测的。”
张怀民没接话,手指在触控板上滑动,放大一组数值。
任昭从本子上撕下一页纸,写下几个参数,递过去:“这是系统推演的结果。你可以拿去对比。”
张怀民接过纸,扫了一眼,眉头皱起:“这应力分布不对。流体不会产生这种局部峰值。”
“不是纯流体。”任昭说,“还有界面作用。”
“什么界面作用?”
“电荷。”
张怀民抬头看他。
“涂层和基材之间,在高速气流下会产生静电积累。”任昭说,“这个力会削弱附着力。”
“没有文献支持这种效应。”张怀民说。
“现在有了。”门口传来声音。
沈知遥走了进来,手里拿着一个密封盒。她把盒子放在桌上,打开,里面是三块金属片,表面涂有纳米涂层。每块片的一端连接着细导线。
“我做了电解液实验。”她说,“把样本浸入模拟大气环境的溶液里,施加不同电压。”
她打开平板,调出图表。横轴是位置坐标,纵轴是界面电势。曲线显示,在曲率突变区,电荷密度显著升高,且与脱落区域完全重合。
“这里。”她指着一处波峰,“电势达到1.7伏,对应第一块脱落点。”
张怀民盯着图表看了十秒,开口:“可能是偶发现象。”
“七次重复实验,六次出现相同分布。”沈知遥说,“概率不支持偶然。”
张怀民合上自己的电脑。
他走到白板前,拿起蓝色马克笔,在左侧画出Navier-Stokes方程的基本形式。然后停下。
“你要加电磁项?”他问。
“不只是电磁。”任昭说,“是双电层理论。”
他在白板右侧写下Gouy-Chapman-Stern模型的简化表达式,标出表面电势与离子浓度的关系。
“把这个耦合进应力张量。”他说,“作为附加项。”
张怀民看着公式,摇头:“计算量会爆炸。而且没有现成求解器支持这种耦合。”
“我们可以改。”任昭说,“先做一维验证。”
沈知遥递给他一组初始参数:“我已经算出了典型工况下的电场分布范围。”
张怀民接过纸,回到电脑前。他新建项目文件夹,命名为“EHD_Adhesion_Model”。然后开始编写新的控制方程模块。
任昭站到他身后,看着代码输入。
“边界条件要重新定义。”他说,“电势不能设为零。”
“我知道。”张怀民敲下一行变量声明,“但这样收敛会很慢。”
“用隐式格式。”任昭说,“时间步长缩小一半。”
两小时后,第一次模拟完成。结果显示附着力预测仍偏高,但偏差从8%降到6.3%。
“方向对了。”沈知遥说。
“网格太粗。”张怀民说,“过渡区需要加密。”
他调整空间离散精度,将关键区域网格细化三倍。
第二次运行开始。等待期间,任昭调出军工复兴系统界面。数学引擎启动,输入当前多场耦合模型。系统判定:此为时代局限性瓶颈,需消耗50军工点获取方向性提示。
他未确认支付。
“不用外力。”他对系统说。
三小时后,第二次结果出来。偏差降至4.1%,曲线形状已接近实测趋势。
“再加一项。”任昭说,“考虑温度梯度引起的离子迁移。”
他在白板上写出热扩散项,标注系数范围。
张怀民沉默地修改代码,加入新变量。
第三次运行中途,程序报错退出。内存溢出。
“变量太多。”张怀民说,“机器撑不住。”
“我去协调超算资源。”任昭说。
他拨通程霄电话,说明需求。二十分钟后,对方回复:已预留四小时优先队列。
“只能跑一次。”任昭说,“必须一次成功。”
沈知遥接手参数优化。她压缩非关键区域网格密度,保留核心过渡区高分辨率。同时设定自适应时间步长策略。
张怀民重写求解器核心循环,采用分块存储技术降低内存占用。
当晚八点十七分,新模型提交至超算中心。
三人守在终端前。
进度条缓慢推进。十二点二十三分,系统提示首次收敛失败。
“初始场设置有问题。”张怀民说,“用实测数据做初值。”
他们导入风洞测试时的温度、压力、速度场数据作为初始条件。
凌晨两点零五分,第二次提交。
进度条走到87%时,再次中断。错误类型:浮点溢出。
“某个节点数值发散。”张怀民调出日志,“找到了,是电势梯度项。”
他限制该变量的上下界,增加稳定性判断逻辑。
第三次提交,时间为凌晨三点四十一分。
这一次,进度条持续前进。
五点零六分,计算完成。
屏幕上弹出结果摘要:附着力预测误差0.3MPa,与实测值偏差小于仪器精度限。
张怀民刷新图像。新模型曲线与实测数据点几乎完全重合,连微小波动都一一对应。
他盯着屏幕,连续看了十一分钟。
任昭没说话,拿起黑色圆珠笔,在记录本上写下:【模型验证通过,可指导工艺重建】。
沈知遥轻声说:“明天就能重启喷涂测试。”
张怀民突然站起来。
他转身面对任昭,双手扶住桌沿,低头。然后弯腰,肩膀下沉,做出一个完整的鞠躬动作。
“任工。”他的声音低,但清楚,“我服了。”
任昭看着他。
“这不是改进。”张怀民直起身,“是革命。”
沈知遥打开新文档,输入标题:《电控界面附着模型操作手册》。
任昭把三支笔重新别好,走到白板前,擦掉旧公式,写下新的建模流程框架。
张怀民坐回电脑前,新建一个共享目录,权限设为全员可读写。他输入第一行文字:【项目名称:第四代战机隐身涂层附着机制联合建模】。
沈知遥将电解液实验数据拖入附件区。
任昭在任务列表中添加一条:【明日九点,召开工艺参数校准会议】。
张怀民打开代码编辑器,开始重构求解器架构,准备将其封装为独立模块。
沈知遥点击保存文档,屏幕右下角时间显示为06:17。
任昭拿起通讯器,拨通实验室值班员号码。
“通知喷涂组。”他说,“今天上午十点,进行首件样件制备。”
通话结束,他看向两人。
“这次。”他说,“让材料听话。”
张怀民的手指停在键盘F5键上方。
沈知遥把激光测距仪放进实验服口袋。
任昭的左手摸了摸眉骨上的旧伤疤。
张怀民按下回车键。

