实验案例核反应堆屏蔽层设计问题1 实验案例.pdfVIP

实验案例 核反应堆屏蔽层设计问题 1 实验案例 1 1.1 核反应堆屏蔽层设计问题1 1.1.1 问题描述与分析1 1.1.2 模型假设2 1.1.3 中子运动的数学描述2 1.1.4 模拟算法中的问题2 1.1.5 模拟结果2 1.1.6 模拟程序3 1.1.7 思考题5 1 实验案例 1.1 核反应堆屏蔽层设计问题 1.1.1 问题描述与分析 核反应堆屏蔽层是用一定厚度的铅包围反应堆，用以阻挡或减弱反应堆发出 的各种射线。在各种射线中，中子对人体伤害极大，因此，在屏蔽层的设计中， 了解中子穿透屏蔽层的概率，对反应堆的安全运行至关重要。 问题背景 假定屏蔽层是理想的均匀平板。 一个中子进入屏蔽层后运动的物理过程：中子以初速度v0 和方向角α射入屏 蔽层，运动一段距离后与铅核发生碰撞，中子获得新的速度及方向( , 1 )v1 θ 。再游 动一段距离后，与铅核发生第二次碰撞，并获得新的状态 ( , 2 )v2 θ ，如此等等， 经过若干次碰撞后，出现下述情况之一时中子终止运动过程： （1）中子被弹回反应堆； （2 ）中子穿透屏蔽层； （3 ）第n 次碰撞后，中子被屏蔽层吸收。 为使屏蔽层的厚度达到安全设计要求，在计算机上对中子在屏蔽层的运动过 程进行模拟。 1 1.1.2 模型假设 1 假定屏蔽层平行板厚度为D=3d，其中d为两次碰撞之间中子的平均游动距离； 2 假设在第10 次碰撞以后，中子速度下降到为某一很小数值而终止运动（被 引收）.因为每次碰撞后，中子因损失一部分能量而速度下降。 3 假定中子在屏蔽层内相继两次碰撞之间游动的距离服从指数分布； 4 中子经碰撞后的弹射角θ～ U(0, 2π). 1.1.3 中子运动的数学描述弹射角θ —第 i 次碰撞后中子的运动方 i 向与x 轴正向的夹角.x — 第 i 次碰撞后中子所处位置与屏蔽层内壁的距离。 i Ri —中子在第 i 次碰撞前后的游动距离.中子在屏蔽层里随机游动，第 i 次碰 撞以后，按照它的位置坐标 x ，可能有以下三种情况发生： i （1）x 0，中子返回反应堆；（2）x ＞D，中子穿透屏蔽层；（3）0＜x ＜D，若 i i i i10，中子在屏蔽层内继续运动， 若 i=10，中子被屏蔽层吸收。经过第 i 次碰撞，中子在屏蔽层内的位置是 x =x i i +R cosθ ，i=1，2，…，10 ， －1 i i 1.1.4 模拟算法中的问题 由于模拟设计较多中子的模拟，又每个中子碰撞最多达10次，因此模拟100 万各中子，则实际运行最多要处理1000万次碰撞，因此逐个处理则费时，因此 可以利用Matlab的函数功能提高模拟效率。 在Matlab环境下要模拟大量中子时算法设计应注意哪些地方？具体参考后 面代码。关于模拟算法，请自己思考，然后参考下面的实现代码，总结模拟中的 算法。 1.1.5 模拟结果 D 当d 3 时模拟的结果如下表所示。特别注意：d 的取值应该查资料确定， 由于模拟中的d 是随意给定的，所以模拟结果与实际合理的结果差别很大。 表 模拟结果（D d 3 ） 中子数/个 返回/% 穿透/% 吸收/ ％ 1000 83.10