《计算物理(本)》[第11章].pptVIP

第十一章 蒙特卡罗方法 §11.1 蒙特卡罗方法概述 蒙特卡罗(Monte Carlo)方法：利用随机数统计地去计算和模拟给定的问题。 说明： 1、也称统计模拟法、随机抽样法或统计试验法。 2、Monte Carlo方法的命名：世界上著名的赌城摩洛哥的Monte Carlo。 3、方法：先构造一个与物理问题等价的随机过程，当完成大量的随机试验后，结果由多次事件的平均值给出。 4、随机数的抽样：现都用计算机程序来产生，一般不用物理方法抽样。计算机产生的是伪随机数。 5、求解确定性物理问题：如微分方程、定积分、线性方程等。同其它数值计算方法相比是速度慢。 6、求解复杂物理问题：如果物理学问题的严格算法不知道，或非常复杂，则蒙特卡罗方法有意想不到的成功。特别在分子运动学、输运现象、布朗运动、放射性衰变等问题中，由于本身有一定的统计规律性，这种方法很奏效。 7、计算机模拟：蒙特卡罗方法广泛应用于“计算机模拟”，在计算机上模拟真实过程。 §11.2 蒙特卡罗方法的原理 用蒙特卡罗方法寻找某个未知量x时，利用计算机产生的随机变量ξ，得到的期望值E(ξ) x= E(ξ) 为了估算ξ的平均值，构造随机变量ξ的若干独立的实验数序列{ξi|i=1，2，3，……n}，得 § 11.3 伪（赝）随机变量的抽样 实际上，大多数伪随机变量不满足[0，1]均匀分布，而是具有符合分布密度函数为f(x)的分布。 抽取符合f(x) 随机数的步骤如下： 一、离散型分布随机变量的直接抽样法 例1.γ光子与物质相互作用的抽样问题。 物理过程：γ光子与物质作用有康普顿效应、光电效应和电子对效应三种类型 。其中光电效应和电子对效应为光子吸收过程。 设三种过程的碰撞截面分别σs 、σe和σp，则总截面σT = σs +σe+σp 。 根据给定的[0，1]之间均匀分布的随机数γ ，求应产生那种效应？ 二、连续型分布随机变量的直接抽样法 一个γ对应一个ξ。如果积分式能求出解析形式的反函数，则得到变换公式 ξF=F -1(γ) 11.12 例1.均匀分布：在区间（a，b）内找出符合均匀分布的随机变量 。 例2．指数分布：指数分布的一般形式为 f(x)=λe-λx x ≥0其中λ0 ，求符合此分布的随机变量ξ § 11.4 Monte Carlo方法求圆周率π 求π的物理模型如图所示，边长为1的正方形内切一个圆。于是正方形和圆的面积分别为 S正=1，S圆= π /4 §11.5 蒙特卡罗方法求解确定性问题 (定积分) 蒙特卡罗方法求解确定性问题的基本思想是：对于给定的确定性问题，设计一个概率模型，然后采用一定的抽样方法按照所设计的概率模型抽样，最后把由这个模型产生的一个数学特征作为该确定问题的近似解。 一．一维定积分计算 2 ．随机投点法 二．多重定积分计算 §11.6 蒙特卡罗方法对链式反应的模拟 （核临界安全模拟计算） 实际链式反应过程是一个随机过程，可以利用蒙特卡罗方法来研究链式反应。 核爆炸：一般235U核的半衰期很长，约7.07*109 年。因此在某一时刻内，只有非常少量的235U发生裂变，释放的能量很小。但是当裂变一旦引起链式反应后，将释放出大量能量，导致核爆炸。 临界条件定量分析 残存因子：N次核裂变，Nin个中子被铀核吸收 用计算机模拟具有一定大小及形状的体积内发生的大量随机裂变，然后计算放出中子被吸收引起的裂变数Nin，求出相应的f 值。 为简单起见，考虑铀的几何形状为长方形a*a*b，如图所示。 蒙特卡罗法研究这种形状的核材料临界问题： （1）核裂变的随机位置 （X0，Y0，Z0） 随机点的坐标范围为 -0.5a X0 0.5a ； -0.5a Y0 0.5a ； -0.5b Z0 0.5b 立体角：dΩ=sinθdθdφ=-dφdcosθ 按立体角内均匀分布是指: φ角在0～2π之间均匀分布；θ在０～π之间，以cosθ的值在-1～1之间均匀分布。 (3)中子平均自由程 平均自由程：中子与铀核碰撞所走的平均路程。 在平均自由程内与核碰撞几率相同：例如平均自由程为1cm，中子在核块内飞越的距离为0.3cm，则中子有30%几率的与铀核发生碰撞被吸收。 中子飞越距离d可以用[0,1]之间的随机数表示。 以下给出计算残存因子的流程图的说明： M=V ：铀核