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

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

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

        随后在众人的注视。

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

        “光达，怎么了这是？”

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

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

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

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

        老郭面沉重的了：

        “当然记得。”

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

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

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

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

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

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

        说到这里。

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

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

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

        过了一会儿。

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

        “唔？”

        中输运方程。

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

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

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

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

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

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

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

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

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

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

        但问题是....

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

        所以呢。

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

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

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

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

        例如控制棒在哪里？

        缩铀如何达到临界积？

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

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

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

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

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

        但如今看来......

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

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

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

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

        其中输运方程的形式如：

        ut  b??u=0

        这里u=(x，t)，

        其中时间变量：t≥0.，

        空间变量：x=(x1，...，xn)∈rn。

        龙套向量：b=(b1，...，bn)∈rn，这是一个固定的向量.。

        接着在边界y:rnx{t=0}上，给定初值，g:rn→r。

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

        这时固定一个任意的(x，t)，并定义z(s)=u(x  sb，t  s)，s∈r。

        利用一开始的方程就可以得到一个表达式：

        dz(s)