蒙特卡洛方法在中子输运中的应用.pdfVIP

《中子输运理论与数值方法》课程作业 ——蒙特卡洛方法 1 目录 2 1. 前言 半个多世纪以来，由于科学技术的发展和 电子计算机的发明，蒙 特卡洛(Monte Carlo)方法作为一种独立的方法被提出来，并首先在核 武器的试验与研制中得到了应用。蒙特卡洛方法是一 计算方法，但 与一般数值计算方法有很大区别。它是以概率统计理论为基础的一 方法。由于蒙特卡洛方法能够比较逼真地描述事物的特点及物理实验 过程，解决一些数值方法难以解决的问题，因而该方法的应用领域日 趋广泛。蒙特卡洛模拟计算是解决中子在介质中输运较为成熟、有效 的方法，对于原子能、辐射防 、剂量学和辐射生物物理学等研究领 3 域实际问题的计算，都可以利用蒙特卡洛方法予以实现。 粒子输运过程可以用玻耳兹曼方程加以描述，然而，以此基础上 发展起来的近似数值方法如扩散近似法、离散坐标方法在处理截面与 能量相关以及散射各向异性介质、复杂几何条件问题时碰到了较大困 难。而蒙特卡洛方法在处理这类问题时得心应手，有很强的解题能力， 并且近似较少，接近于真实情况。 粒子辐射问题计算通常有输运方程法、蒙特卡洛法(MC 法)、实验 测量法以及经验法等几种方法。蒙特卡洛计算法又称随机抽样法或统 计试验法，是基于计算机模拟的思想，抓住物理过程的数量和几何特 征，进行数字模拟试验，该方法是求解辐射输运 问题的一种相当成熟 和有效的方法，而且它对于各种复杂问题，具有良好的通用性，实用 性相当广泛，几乎涉及核科学的各个领域。本文主要介绍蒙特卡洛的 概念、原理和应用及研究现状。 2. 蒙特卡洛方法概述 蒙特卡洛方法又称随机抽样技巧或统计试验方法。半个多世纪以 来，由于科学技术的发展和 电子计算机的发明 ，这种方法作为一 独立的方法被提出来，并首先在核武器的试验与研制中得到了应用。 蒙特卡洛方法是一 计算方法，但与一般数值计算方法有很大区别。 它是以概率统计理论为基础的一种方法。由于蒙特卡洛方法能够比较 逼真地描述事物的特点及物理实验过程，解决一些数值方法难以解决 4 的问题，因而该方法的应用领域日趋广泛。 蒙特卡洛方法的主要组成部分有： (1)概率密度函数(pdf)— 必须给出描述一个物理系统的一组概率 密度函数； (2)随机数产生器—能够产生在区间[0,1]上均匀分布的随机数； (3)抽样规则—如何从在区间[0,1]上均匀分布的随机数出发,随 机抽取服从给定的pdf 的随机变量； (4)模拟结果记录—记录一些感兴趣的量的模拟结果； (5)误差估计—必须确定统计误差（或方差）随模拟次数以及其它 一些量的变化； (6)减少方差的技术—利用该技术可减少模拟过程中计算的次数； (7)并行和矢量化—可以在先进的并行计算机上运行的有效算法 2.1 蒙特卡洛方法的基本思想 可以通俗地说，蒙特卡洛方法是用随机试验的方法计算积分，即 将所要计算的积分看作服从某种分布密度函数 f(r)的随机变量ｇ(r)的 数学期望   g  0 g (r)f (r)dr (0.1) 通过某 试验，得到 N 个观察值 r ，r ，…，r （用概率语