第五百八十四章 如果权威错了呢？(6/8)

算这个问题可以直接用蒙卡计算。

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

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

例如控制棒在哪里插？

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

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

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

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

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

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

但如今看来......

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

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

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

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

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

ut+b??u=0

这里u=(x，t)，

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

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

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

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

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

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

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

dz(s)ds=b??u(x+sb，t+s)+ut(x+sb，t+s)?1=0。

从这个表达式不难看出。

对每个点(x，t)，u在穿过(x，t)且方向是(b，1)的直线上是个常数，实际上就是它在t=0时刻的初值。

接着再加上一个扩散方程的增值项，很轻松就可以得到一个指数项是e的正数次的结果。

至少以老郭的数学水平看来，这个推导过程不存在什么明显异常。

但是在看到结果时，他整个人却瞬间愣住了。

只见此时此刻。

最后无穷项级数的求和上，显示的赫然是一个指数项是e的负数次的结果！

看到这里。

老郭勐然抬起头，看向了对面的陆光达。

陆光达则无奈朝他