深夜,温卿在计算中心对着空白的代码编辑器发呆。
她忽然想到一个或许可行的替代方案:
用分子动力学模拟。
分子动力学不依赖于宏观本构关系,直接模拟原子间的相互作用。
给定原子间势函数,就能计算材料在任意条件下的响应。
虽然模拟的时空尺度很小,但可以用来研究材料在极端条件下的微观变形机制,从而推断宏观强度行为。
关键还是势函数。
描述钚、铀等锕系元素在极端条件下的原子间相互作用,需要高度精确的多体势,这本身就是前沿课题。
温卿想起记忆中的一个碎片:
关于“机器学习势函数”的概念。
用神经网络拟合第一性原理计算数据,构建既精确又高效的原子间作用模型。
但这个时代,神经网络还只是学术界的玩具,距离工程应用遥不可及。
她只能退而求其次,采用现有的经验势函数,通过调整参数来匹配有限的实验数据。
这无疑是“带着镣铐跳舞”,但总比完全忽略强度效应要好。
炸药与金属的耦合问题,是三个难点中最“工程化”的一个,却也让温卿遭遇了最激烈的争议。
现有的仿真中,炸药爆轰产物与金属壳体的相互作用,被简化为一个“活塞模型”:
爆轰产物提供均匀压力,推动金属向内运动。
这个模型完全忽略了爆轰波阵面的复杂结构、产物气体的非理想性、以及炸药-金属界面的真实物理过程。
温卿从爆炸力学组调来了一组高速摄影照片,那是微缩模型试验中拍摄的。
照片显示,爆轰波撞击金属界面时,会产生极其复杂的反射和透射波系,界面附近出现剧烈的湍流混合。
有些照片甚至捕捉到了金属表面被“侵蚀”的迹象——
高温高压的爆轰产物可能将微量金属熔化甚至气化,形成一层混合等离子体。
“这个混合层会改变能量传递效率。”
温卿在跨组协调会上指出。
“一部分能量消耗在金属的相变上,而不是全部用于做功。
另外,混合层的不稳定性可能导致界面扰动增长,最终影响内爆对称性。”
爆炸力学组组长老秦,一位声音洪亮、作风强势的老专家,直接提出了质疑:
“温工,你考虑的这些效应,在宏观尺度上都是高阶小量。
我们的工程经验表明,只要炸药透镜设计合理,起爆同步性控制得好,界面效应可以忽略。”
“秦老,那是对于现有的大型化设计。”
温卿保持着礼貌但坚定的语气。
“当尺寸缩小到原来的三分之一,特征尺寸减小,表面积与体积比增大,界面效应的影响会被放大。
而且,小型化要求使用能量密度更高的新型炸药,这些炸药的爆轰产物状态更复杂,非理想性更强。”
她调出对比仿真结果:
用传统活塞模型和初步改进的界面模型,分别计算同一个小型化构型。
结果显示,在压缩后期,界面模型预测的压力峰值比活塞模型低15,压力脉动幅度大30。
温卿说。
老秦盯着数据看了很久,最后说:
“你的模型需要验证。我们需要设计专门的界面试验,用实测数据来校核。”
这是一个合理的要求,但也是一个巨大的挑战。
炸药-金属界面过程发生在微秒量级、毫米尺度,诊断极其困难。
需要发展全新的微细观测试技术。
温卿主动请缨,与测试技术组合作,开始了高难度诊断方法的攻关。
这占用了她大量的时间和精力,使她不得不将状态方程和材料强度的研究暂时放缓。
三个月过去了,进展缓慢。
三个方向都遇到了瓶颈:
状态方程缺乏可靠的电子结构数据;
材料强度模型缺乏势函数和验证手段;
界面模型则需要等待新诊断技术的开发。
温卿开始感到一种深切的焦虑。
她每晚只睡三四个小时,灵泉的恢复效果似乎也在高强度消耗下打了折扣。
她的脸色更加苍白,但眼睛却因为长期处于“思维超频”状态而异常明亮。
转机出现在一个凌晨。
温卿在处理一组陈旧的界面测试数据时,发现了异常。
那是十几年前一次失败的试验记录,由于诊断故障,数据残缺不全,一直被封存。
温卿抱着死马当活马医的心态重新分析,却发现了一个被忽略的细节:
在某个特定炸高下,金属的加速历程出现了一个微小的“平台期”——
加速度没有单调增加,而是停顿了零点几微秒。
按照传统模型,这是不可能的。
爆轰压力持续作用,加速度应该持续增加。
温卿盯着那条反常的曲线,大脑飞速运转。
突然,末世记忆中的一个碎片闪现:
那是关于“冲击波与稀疏波相互作用”的图景。
在爆轰波撞击金属后,不仅会产生向金属内部传播的冲击波,还会在爆轰产物中产生向回传播的稀疏波。
这个稀疏波会降低界面压力。
如果炸药与金属的间距合适,稀疏波到达界面的时间,可能正好与金属开始加速的时间耦合,产生暂时的压力抵消效应!
她立刻在草稿纸上推导。
这是一个一维的近似模型,但物理图像清晰。
计算结果显示,在某个临界炸高下,确实会出现一个短暂的压力平台。
“所以界面效应不仅仅是混合和侵蚀,”
温卿兴奋地记录。
“还有波系的动态干涉!这个效应是尺度敏感的——
小型化后炸高相对值变化,这个临界条件更容易被触发!”
这个发现虽然只是冰山一角,却为她打开了一扇窗:
也许不需要完全精确地模拟所有复杂过程,而是抓住最关键的主导物理机制。
对于状态方程,也许不是追求完美的电子结构描述,而是找到影响压缩性的“主导修正项”;
对于材料强度,也许不是构建普适的本构关系,而是识别小型化条件下最可能发生的失效模式。
她调整了研究策略:
从“全面改进”转向“关键突破”,从“追求完美模型”转向“发展实用模型”。
基于这个思路,温卿暂时搁置了第一性原理计算的庞大计划,转而分析现有实验数据中的“异常”规律,尝试用简单的物理模型去解释。
对于状态方程的偏差,她提出了一个唯象的“硬化因子”,这个因子是压力和温度的函数,形式简单,但能很好地拟合已有数据趋势。
对于材料强度,她聚焦于“剪切失稳”这一种最危险的失效模式,建立了简化的判据。
至于界面模型,她将稀疏波干涉效应纳入,发展了一个“波系耦合”的简化界面模型。
虽然粗糙,但物理意义明确,而且能复现那次陈旧试验中的反常现象。