第五百八十四章 如果权威错了呢？（第3/4页）

见此情形。

老郭便立马喊起了陆光达的名字：

“光达！”

陆光达原本正拧着眉毛看着面前的一份文件呢，闻言立马抬起了头。

见到老郭和徐云后他先是微微一怔，旋即眼中便冒出了光，大喜道：

“老郭，韩立同志！哎呀呀，你们可算来了！”

随后在众人的注视下。

老郭跟在超市推购物车似的将徐云推到桌边，朝陆光达面前的稿纸努了努下巴：

“光达，怎么了这是？”

之前由于保密需要，陆光达只派人用简洁的内容通知了老郭理论组这边遇到了一些情况，大致和中子运输有关，但再多老郭就不了解了。

陆光达闻言叹了口气，对老郭解释道：

“老郭，还记得之前你带回来的那套外文期刊吗？”

“里头有一份长友同志从洛斯阿拉莫斯国家实验室中带出来的文件，文件除了一些重要参数，还描述了一个定态次的临界状态模型。”

老郭面色沉重的点了点头：

“当然记得。”

这件事他怎么可能忘的了呢？

毕竟这份文件可是他的至交好友，用生命换来的绝密文件啊……

随后陆光达顿了顿，继续说道：

“我们在翻译好这份文件后立刻进行了全组学习和讨论，在不久前也确实产出了很不错的理论成果。”

“这些理论成果把我们的研究效率硬生生的推了一大截，很多原先停滞的地方也开始出现了松动。”

“只是在后续的中子通量密度计算的时候，我们突然发现了一个情况……”

说到这里。

陆光达抬头看了老郭一眼，从桌上拿起了一份文件：

“根据我们对模型的后续衍生计算，发现有些计算结果存在明显的异常。”

老郭闻言眉头一皱，取过文件看了起来。

过了一会儿。

他的鼻翼间忍不住发出了一道轻咦：

“唔？”

中子输运方程。

这是原子弹研制过程中非常非常重要的一个模块。

上辈子是原子弹的同学应该都知道。

核爆过程中子与核碰撞的概率是一个很复杂的过程——无论是是应用还是计算上都是如此。

原子弹最开始就是搞轰击，然后炸出中子。

中子传播过程中遇到新的核，接着发射新的中子。

这些中子会随机向不同方向运动，再次进行撞击，如此反复……

这么一轮又一轮的过程，必须要在数学上精确到每一轮过程中中子的运动状态。

用术语来描述就是这样的：

初始在堆内某一位置具有某一能量及某一运动方向的中子，稍晚些时候，将运动到堆内的另一位置以另一能量和另一运动方向出现。

这种运动轨迹用数学方程组表示，便是中子输运方程。

但问题是……

链式反应后产生的中子能量分布很广，需要求解多群的玻尔兹曼方程，而且这玩意还没有解析解。

所以呢。

只能离散后再通过多种计算方法求数值解，核武器里面核燃料的形状也比较复杂，所以求解起来更加困难。

后世的计算机算力强，计算这个问题可以直接用蒙卡计算。

但眼下这个时代只能靠手解单群的中子输运方程，这就很麻烦了。

可你不解决这个问题又不行，因为没有具体单解的话，很多应用上的操作是无法进行的。

例如控制棒在哪里插？

高浓缩铀如何达到临界体积？

合适的燃料摆放方式是什么？

没有具体的数值，这些东西是搞不起来的。

因此当初在拿到洛斯阿拉莫斯国家实验室文件的时候，老郭是既悲痛又开心。

悲痛是因为这份文件的获取过程太过坎坷，不止一位同志战友牺牲在了护送途中。

开心则是因为有了这份文件，很多难点应该就可以顺利解决了。

但如今看来……

这件事远远没有那么简单。

例如他手上的这份计算稿纸，这是一轮非常标准的的一般数值的计算过程。

也就是当粒子的平均自由程非常小时。

在扩散条件下通过光学厚胞腔……也就是原子弹应用过程中的一个模块的数值，来求解离散纵坐标。

其中输运方程的形式如下：

ut＋b·▽u=0

这里u=（x，t）。

其中时间变量：t≥0。

空间变量：x=（x1，……，xn）∈Rn。

龙套向量：b=（b1，……，bn）∈Rn，这是一个固定的向量。

接着在边界Γ：Rn×{t=0}上，给定初值，g：Rn→R。

观察上面这个方程，不难发现u沿某个特定方向的导数为0。

这时固定一个任意的点（x，t），并定义z（s）=u（x＋sb，t＋s），s∈R。