“两个有效缓发中子份额进行联立，则可以得到中子代时间和燃料需要的密度数……”

“接着一部分中子会被燃料吸收发生辐射俘获，卫东航同志他们的计算结果可以把这部分的相对几率进行修正……”

“假设每次引起裂变反应释放的平均中子数为ν，而且有了一个中子被吸收后裂变的概率，那么我们给这两个做乘积，就会得到燃料每吸收一个中子产生的裂变中子数为η=νσ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后中子的数目为：

n1000l=[1－βk]mn0＋λlc0[1－βk]m－1＋λlc0[1－βk]m－2＋……＋λlc01－βk＋λlc0=[[1－βk]m＋[1－βk1－β[1－k1－β]]＋β1－k1－β]n0。

然后再引入爆轰方程，就可以得到znd模型了。

非常简单，也非常好理解，有手就行。

而随着znd模型的顺利建立，剩下的便是……

参数的引入与计算了。

当然了。

到了这一步，单纯的人力肯定是不合适的……或者说不稳妥——毕竟事关重大嘛。

因此除了人力计算之外，计算机也是离不开的一个辅助项。

于是很快便有理论组的成员将这些关键数据和公式摘抄下来，前往地面送到了计算机中心进行同步推导——之前在太上项目成立的时候，104机已经被送到了基地。

如今计算机所的几位专家都在地面上待命，随时可以进行着参数的输入与计算。

不过陆光达等人也没有干坐着，而是继续进行了笔算。

多一