“我是今晚的值班员小李。钱主任交代了,二号机十点到十二点归你。”
温卿谢过他,走到二号机前。
机器面板上的指示灯阵列令人眼花缭乱,但她已经熟悉了这些设备的操作。
“你要算什么程序?需要我帮忙输入吗?”小李问。
“我自己来就行。”
温卿从包里取出主程序纸带,“不过可能需要你指导一些操作细节。”
九点五十分,上一组用户结束计算。
小李协助温卿把她的纸带装到纸带读入器上。
“你这个程序不小啊。”
小李看着那卷厚厚的纸带,“至少五千行代码。”
“嗯,简化了很多,但还是复杂。”温卿检查着设备状态。
十点整,机器清空,温卿的程序开始载入。
纸带读入器咔哒咔哒地响着,黄色纸带缓缓移动。
这个过程需要十分钟——五千行代码,一个字符一个字符地读进去。
温卿屏住呼吸。
她最怕的就是读带过程中出错。
一旦某个字符读错,整个程序就可能崩溃,需要重新来过。
还好,这次顺利。
十点十分,程序载入完毕。
温卿开始在控制台输入初始参数。
她要计算的是一个极度简化的非平衡态模型——只考虑一维对称情况,只处理最关键的三个非平衡效应:
速度非平衡、温度非平衡、以及最简单的辐射耦合。
即便如此,模型也已经比传统模型复杂得多。
传统模型只需要解五个方程,她的模型要解十二个,而且很多项是非线性的。
“开始运行。”温卿按下执行键。
机器轰鸣声陡然增大。
面板上的指示灯疯狂闪烁,像一场电子风暴。
喜欢七零疯批女王,踏雪封神请大家收藏:七零疯批女王,踏雪封神更新速度全网最快。
温卿紧盯着输出设备——一台每分钟只能打印十行的慢速打印机。
第一行输出出现在十点十五分:
“初始条件:马赫数=10,高度=60k,开始迭代...”
很好,程序运行正常。
第一个:常规工况
温卿设计的计算方案,是从常规工况逐步过渡到极端工况。
这样既能验证模型在常规情况下的合理性,也能观察非平衡效应如何随条件变化而增强。
十点三十分,第一个工况计算完毕。
温卿拿起打印出的结果,与实验数据对比。
马赫数10、高度60公里——这是相对温和的再入条件。
结果:新模型预测的热流值,与传统模型几乎一致,与实际测量值的误差都在5%以内。
“在常规条件下,非平衡效应确实不显着。”温卿记下结论。
她调整参数,开始计算下一个工况:马赫数12,高度50公里。
机器继续轰鸣。
打印机的速度太慢,温卿等得心焦。
她开始在心里默算一些中间结果——精神力提升后,她的心算能力堪比小型计算机,虽然精度不如机器,但能给出大致趋势。
十一点整,第二个工况结果出炉。
这次出现了微妙差异:
新模型预测的热流值,比传统模型高3%,而实际测量值正好落在两者之间。
“非平衡效应开始显现了。”
温卿在笔记本上标注。
喜欢七零疯批女王,踏雪封神请大家收藏:七零疯批女王,踏雪封神更新速度全网最快。