“两个有效缓发中子份额进行联立，则可以得到中子代时间和燃料需要的密度数……”
“接着一部分中子会被燃料吸收发生辐射俘获，卫东航同志他们的计算结果可以把这部分的相对几率进行修正……”
“假设每次引起裂变反应释放的平均中子数为ν，而且有了一个中子被吸收后裂变的概率，那么我们给这两个做乘积，就会得到燃料每吸收一个中子产生的裂变中子数为η=νσfσγ＋σf=ν11＋α……”
陆光达板书的速度很快，一眨眼就能写下一大行字，但台下大多数人都能跟得上他的思路。
同时与之前一样。
由于没有任何早期经验的缘故，此番即便是陆光达亲自出手，整个推导过程也出现了一些错漏的地方。
不过好在现场的这些‘评审’各个能力超群，基本上在某个环节出现错误后十秒钟内，便会有五六双手同时举了起来。
因此这些错漏虽然偶尔存在，但都很快便被拔除了个干净。
三个小时后。
在时间来到晚上十点半之际。
哗——
陆光达用最后半截粉笔头在黑板上画了个圈，同时看下了台下：
“好了，同志们，这就是我们最后构建出来的znd模型了。”
“其中热中子吸收截面σa=σf＋σγ，有效增殖系数为ηfeppnl，常数源是……”
陆光达很快报出了七八个关键参数，迅速构建出了一个纸面的znd模型。
这个znd模型除了之前的数学计算之外，理论逻辑其实也很简单：
陆光达他们先计算出了一个常数源方程，当k＞1时这个方程没有稳定解，当k=1时上述方程方程稳态解不唯一。
但k＜1时，方程存在见渐近解。
同时在当初陆光达他们计算中子运输方程的时候，理论组曾经得出过一个非常重要的结论：
中子的链式裂变反应装置对吸收截面0.5％的变化响应是非常剧烈的。
在这个基础上。
陆光达他们根据先驱核平衡浓度反推出了一个平衡方程，表达式为dc0dt=0=βνnfσfvn0－λc0。
若截面在t=0时刻发生0.5％的变化，那么在t=0.1s时，瞬发中子的增殖为[（1－β）k]1000。
在每一个增殖间隔l内，裂变产物在衰变时释放λlc个缓发中子，缓发中子源在接下来的第一个增殖间隔内产生（1－β）kλlc，第二个间隔内产生[（1－β）k]2λlc缓发中子。
以此类推。
如果产生1000代瞬发中子增殖间隔内均存在一个缓发中子源，而且假设裂变碎片的浓度保持不变（c=c0），那么0.01s后中子的数目为：
n（1000l）=[（1－β）k]mn0＋λlc0[（1－β）k]m－1＋λlc0[（1－β）k]m－2＋……＋λlc0（1－β）k＋λlc0=[[（1－β）k]m＋[1－βk（1－β）[1－k（1－β）]]＋β1－k（1－β）]n0。
然后再引入爆轰方程，就可以得到znd模型了。
非常简单，也非常好理解，有手就行。
而随着znd模型的顺利建立，剩下的便是……参数的引入与计算了。
当然了。
到了这一步，单纯的人力肯定是不合适的……或者说不稳妥——毕竟事关重大嘛。
因此除了人力计算之外，计算机也是离不开的一个辅助项。
于是很快便有理论组的成员将这些关键数据和公式摘抄下来，前往地面送到了计算机中心进行同步推导——之前在太上项目成立的时候，104机已经被送到了基地。
如今计算机所的几位专家都在地面上待命，随时可以进行着参数的输入与计算。
不过陆光达等人也没有干坐着，而是继续进行了笔算。
多一方计算就多了一个保障，到时候各方把结果一汇总对比就行了。
接着很快。
会议室内便响起了噼里啪啦的算盘声。
“小珠进一，大珠退一……”
“二上三去五……”
“三百五十七乘四点九九，结果应该是……”
不同材质的算盘在不同力度的拨动下发出了不同频率的声音，在徐云耳中组成了一曲轻快动听的乐章。
但很快，徐云便意识到了一个“真相”：
令他感觉动听的并非是算盘的敲击声，而是在打算盘的这些……人。
大于、陆光达、王淦昌、程开甲……这些前辈用心血与忠诚演奏出的曲子，怎么可能不好听呢？
蓦然，徐云又想起了穿越前看到一个新闻：
某个国产手机品牌时隔三年，在海对面的层层封锁下研制出了5g国产芯片。