第1293节（1 / 2）

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

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

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