“两个有效缓发中子份额进行联立,则可以得到中子代时间和燃料需要的密度数……”
“接着一部分中子会被燃料吸收发生辐射俘获,卫东航同志他们的计算结果可以把这部分的相对几率进行修正……”
“假设每次引起裂变反应释放的平均中子数为ν,而且有了一个中子被吸收后裂变的概率,那么我们给这两个做乘积,就会得到燃料每吸收一个中子产生的裂变中子数为η=νσ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]0+λlc0[(1-β)k]1+λlc0[(1-β)k]2+……+λlc0(1-β)k+λlc0=[[(1-β)k][1-βk(1-β)[1-k(1-β)]]+β1-k(1-β)]n0。
然后再引入爆轰方程,就可以得到znd模型了。
非常简单,也非常好理解,有手就行。
而随着znd模型的顺利建立,剩下的便是……
参数的引入与计算了。
当然了。
到了这一步,单纯的人力肯定是不合适的……或者说不稳妥——毕竟事关重大嘛。
因此除了人力计算之外,计算机也是离不开的一个辅助项。
于是很快便有理论组的成员将这些关键数据和公式摘抄下来,前往地面送到了计算机中心进行同步推导——之前在太上项目成立的时候,104机已经被送到了基地。
如今计算机所的几位专家都在地面上待命,随时可以进行着参数的输入与计算。
不过陆光达等人也没有干坐着,而是继续进行了笔算。
多一方计算就多了一个保障,到时候各方把结果一汇总对比就行了。
接着很快。
会议室内便响起了噼里啪啦的算盘声。
“小珠进一,大珠退一……”
“二上三去五……”
“三百五十七乘四点九九,结果应该是……”
不同材质的算盘在不同力度的拨动下发出了不同频率的声音,在徐云耳中组成了一曲轻快动听的乐章。
但很快,徐云便意识到了一个“真相”:
令他感觉动听的并非是算盘的敲击声,而是在打算盘的这些……人。
大于、陆光达、王淦昌、程开甲……这些前辈用心血与忠诚演奏出的曲子,怎么可能不好听呢?
蓦然,徐云又想起了穿越前看到一个新闻:
某个国产手机品牌时隔三年,在海对面的层层封锁下研制出了5g国产芯片。
这个新闻看似商业竞争,背后却是国与国的博弈。