第235章 数据洪流(2/2)
很经典的理论,但……太理想化了。
“实际飞行中,烧蚀可能不均匀。”
她提出自己的看法。
“材料内部有杂质,或者结构不均匀,烧蚀速率就会有差异。烧蚀产物也不是均匀的气体层,可能是湍流、涡旋……”
“你说得对。”
老王叹气,“但我们没有更精确的模型。国外最新的理论根本不公开,我们只能靠自己摸索。”
温卿脑海中闪过一些碎片:
关于非平衡态流体的描述,关于烧蚀表面微观形貌对热流影响的讨论……
“王老师,我能看看我们的数值模拟程序吗?”她问。
老王有些意外:“你看得懂Fortran?”
“培训时学过一点。”温卿说。其实是在末世学的,那时候编程已经是基础技能。
老王犹豫了一下,还是从保险柜里取出一卷纸带:
“这是主程序的源代码。你看可以,但绝不能带出去,也不能抄录。”
纸带很旧了,边缘有些磨损。
温卿小心地展开,在桌上铺开。
代码是用FortranIV写的,有几千行。
她一行行看下去。
程序结构清晰,但物理模型确实简单:
材料烧蚀是均匀的,烧蚀产物是理想气体,等离子体鞘套是稳定层流……
很多假设都是为了简化计算,但可能偏离实际情况太远。
输入数据的工作完成后,温卿开始参与数值模拟。
DJS-130计算机的内存只有32KB,速度慢得令人发指。
喜欢七零疯批女王,踏雪封神请大家收藏:七零疯批女王,踏雪封神更新速度全网最快。
一个完整的再入过程模拟,需要运行十几个小时。
而且经常因为各种原因中断——电压不稳、纸带读错、甚至只是天气潮湿导致接触不良。
温卿很快掌握了操作技巧。
她发现,通过优化算法,可以减少一些不必要的计算量,把运行时间缩短到十小时以内。
“你怎么想到这样改的?”
老王看她修改了输入参数文件,惊讶地问。
“我觉得这部分网格可以粗化,对结果影响不大,但计算量能减一半。”
温卿解释,“还有这个迭代收敛标准,可以放宽一点,也能节省时间。”
老王试运行了一次,结果和之前几乎一样,但时间节省了40%。
“可以啊。”
他拍拍温卿的肩膀,“以后这活儿交给你了。”
但更核心的问题依然没解决:模拟结果和实验数据的偏差。
温卿把多次试验的数据汇总,绘制成图表。
曲线图上,理论模拟线平滑上升,达到峰值后缓慢下降;实际测量线却像锯齿一样,剧烈波动。
偏差最大的时刻,出现在再入后第45秒到55秒之间。
这个时间段,弹头经历最剧烈的气动加热。
“这个现象,我们叫它‘热流脉动’。”
孙研究员在组内会议上说。
“目前国际上也报道过类似现象,但机理不明。美丽国人认为是边界层转换引起的,苏国人则归因于烧蚀产物的不稳定性。”
他转向温卿:“小温,你这段时间一直在处理数据,有什么想法?”
全组人的目光集中在温卿身上。
她是最年轻的,资历最浅,但也是培训成绩最好、进步最快的。
温卿站起身,走到白板前。
喜欢七零疯批女王,踏雪封神请大家收藏:七零疯批女王,踏雪封神更新速度全网最快。