第五百八十四章 如果权威错了呢?(7/10)
来更加困难。
后世的计算机算力强,计算这个问题可以直接用蒙卡计算。
但眼下这个时代只能靠手解单群的中子输运方程,这就很麻烦了。
可你不解决这个问题又不行,因为没有具体单解的话,很多应用上的操作是无法进行的。
例如控制棒在哪里插?
高浓缩铀如何达到临界体积?
合适的燃料摆放方式是什么?
没有具体的数值,这些东西是搞不起来的。
因此当初在拿到洛斯阿拉莫斯国家实验室文件的时候,老郭是既悲痛又开心。
悲痛是因为这份文件的获取过程太过坎坷,不止一位同志战友牺牲在了护送途中。
开心则是因为有了这份文件,很多难点应该就可以顺利解决了。
但如今看来
这件事远远没有那么简单。
例如他手上的这份计算稿纸,这是一轮非常标准的的一般数值的计算过程。
也就是当粒子的平均自由程非常小时。
在扩散条件下通过光学厚胞腔也就是原子弹应用过程中的一个模块的数值,来求解离散纵坐标。
其中输运方程的形式如下:
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)ds=b??u(x+sb,t+s)+ut(x+sb,t+s)?1=0。
从这个表达式不难看出。
对每个点(x,t), u在穿过(x,t)且方向是(b,1)的直线上是个常数,实际上就是它在 t=