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

实验案例 核反应堆屏蔽层设计问题 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）中子被弹回反应堆； （2）中子穿透屏蔽层； （3）第n次碰撞后，中子被屏蔽层吸收。 为使屏蔽层的厚度达到安全设计要求，在计算机上对中子在屏蔽层的运动过程进行模拟。 模型假设 1 假定屏蔽层平行板厚度为D=3d，其中d为两次碰撞之间中子的平均游动距离； 2 假设在第10 次碰撞以后，中子速度下降到为某一很小数值而终止运动（被引收）.因为每次碰撞后，中子因损失一部分能量而速度下降。 3 假定中子在屏蔽层内相继两次碰撞之间游动的距离服从指数分布； 4 中子经碰撞后的弹射角θ～ U(0, 2π). 中子运动的数学描述 弹射角θi —第i 次碰撞后中子的运动方向与x 轴正向的夹角. xi — 第i 次碰撞后中子所处位置与屏蔽层内壁的距离。 Ri —中子在第 i 次碰撞前后的游动距离. 中子在屏蔽层里随机游动，第 i 次碰撞以后，按照它的位置坐标 xi，可能有以下三种情况发生： （1）xi0，中子返回反应堆； （2）xi＞D，中子穿透屏蔽层； （3）0＜xi＜D，若i10，中子在屏蔽层内继续运动， 若i=10，中子被屏蔽层吸收。 经过第i 次碰撞，中子在屏蔽层内的位置是 xi=xi－1+Ricosθi ，i=1，2，…，10 ， 模拟算法中的问题 由于模拟设计较多中子的模拟，又每个中子碰撞最多达10次，因此模拟100万各中子，则实际运行最多要处理1000万次碰撞，因此逐个处理则费时，因此可以利用Matlab的函数功能提高模拟效率。 在Matlab环境下要模拟大量中子时算法设计应注意哪些地方？具体参考后面代码。关于模拟算法，请自己思考，然后参考下面的实现代码，总结模拟中的算法。 模拟结果 当时模拟的结果如下表所示。特别注意：d的取值应该查资料确定，由于模拟中的d是随意给定的，所以模拟结果与实际合理的结果差别很大。 表 模拟结果（） 中子数/个 返回/% 穿透/% 吸收/％ 1000 83.10 10.70 6.20 5000 80.90 13.02 6.08 10000 81.36 12.08 6.56 图 模拟结果图示（，中子数为10000） 模拟程序 function sim_zhongzi2 %模拟核反应堆屏蔽层中中子的运动 %为了提高程序运行速度，这里没有逐个动态显示中子的运动， %所以一开始就产生所有中子的初始状态；当然，也可以很容易修改为动态显示程序。 % n=input(中子个数：);%中子个数 xx=[];%存储第10次碰撞后中子的横坐标 yy=[];%存储第10次碰撞后中子的纵坐标 N=10;%每个中子最多碰撞次数 d = 2; D = 3*d; H =10*10; c=zeros(1,3);%c(1)返回反应堆数量，c(2)穿透,c(3)吸收 %产生中子 x=zeros(1,n); y=H*rand(1,n)/10; %中子运动 for j=1:N,%考虑至多10次碰撞 if isempty(x), break; end R = -d*log(rand(1,length(x))); seta = 2*pi*rand(1,length(x)); x = x + (R.*cos(seta)); y = y + (R.*sin(seta)); %c(2)穿透 t=find(xD); c(2) = c(2) + length(t); if length(t)0, xx=[xx,x(t)]; yy=[yy,y(t)];