蒙特卡罗MCNP学习汇总（一）-----MCNP简介及编写第一个程序

辽G.默语

已于 2024-04-16 09:36:58 修改

阅读量3.1w

点赞数 20

分类专栏： MCNP 文章标签：算法

于 2022-04-08 18:35:21 首次发布

本文链接：https://blog.csdn.net/qq_50652615/article/details/124047516

版权

MCNP 专栏收录该内容

5 篇文章 77 订阅

订阅专栏

简介：

什么是MC模拟

一种通过随机抽样解决数学问题的一种数值计算方法。 MC方法解决的主要数学问题 -数值积分问题 -随机物理过程（计算实验）：粒子碰撞、径迹、计数等

MC方法实现粒子的输运：真实模拟一个粒子由产生到消亡的过程不断重复粒子，后去最终的计数（Tally）

1、中子散射、产生光子

2、裂变、产生光子

3、中子俘获

4、中子逃逸

5、光子散射

6、光子逃逸

7、光子俘获

—MC：采用计算机在空间中采用随机方法模拟粒子的传输，估算粒子的数学期望（每个粒子的平均行为）

介绍

A General Monte Carlo N-Particle Transport Code

美国Los Alamos National Laboratory (LANL)开发是当前使用最广泛的粒子输运软件具有很长的历史背景和用户基础采用3-D几何模型具有复杂物理过程模型和核数据驱动涉及粒子：中子、光子和电子内置多种方差减小技术软件中包含燃耗分析、高能物理（基础物理）、重粒子输运等（新版本软件）

应用

临界安全分析（MCNP软件研发的最初目的）辐射屏蔽核安全分析裂变、聚变反应堆设计探测器设计与分析辐射剂量学医学物理核废料存储和处置加速器靶设计石油测井放射医学空间辐射

运行

如何运行MCNP

要运行MCNP，必须：建立一个输入文件（定义）执行命令行使用的软件： MCNP 版本：MCNP5 1.5 OS:WINDOWS

watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBA6L69Ry7pu5jor60=,size_20,color_FFFFFF,t_70,g_se,x_16

注意： MCNP自带的画图需要X11 视窗支持 MCNP有一个独立的可视化程序—VisEd，使用较为广泛，已经逐步替代内置的画图功能

编写第一个程序

格式

Title Card（标题）

One-line problem description(问题描述）

Cell Card （体）

Required blank line as separator（空行，必须、只能一行）

Surface Card（面）

Required blank line as separator（空行，必须、只能一行）

Data Card(粒子)

Recommended blank line as terminator – everything else will be ignored（空行之后的内容都不算）

程序

编写一个大球包裹住小圆柱的程序

This is my first mcnp file$标题
C Cell card
101 1 -3.76 -101 
102 0 101 -102
103 0 102

C Surface card
101 RCC 0 0 0 0 0 7.5 3.75
102 SO 20

C data card
MODE P E
imp:p 1 1 0
imp:e 1 1 0
M1 11000 1 53000 1
SDEF POS=0 0 -1 ERG=0.661
F28:P,E 101
E28 0 0.001 510I 1
F18:P,E 101
E18 0 0.001 510I 1
FT18 GEB 0 0.0567 0
NPS 1000000
PRINT

讲解

$ 注释符

This is my first mcnp file $ 标题
C Cell card $ 体卡（定义自己想构建的几何体，把所有地方都要定义且不能重合）
101 1 -3.76 -101 $ 编号101（自己定义） 1号物质密度3.76 下面定义的编号为101面的里面
102 1 -3.76 101 -102$ 编号102 1号物质密度3.76 101面外面 102面里面
103 0 102$ 编号103 真空 102面外面

C Surface card$面卡（构建自己需要的面，利用面围成体）
101 RCC 0 0 0 0 0 7.5 3.75$编号 101 圆柱面前三位为底面圆心方向向量 Z周方向 7.5高度半径为3.75
102 SO 20$编号102 圆心在原点的圆面半径为20

C data card$体卡（以后说明）
MODE P E
imp:p 1 1 0
imp:e 1 1 0
M1 11000 1 53000 1
SDEF POS=0 0 -1 ERG=0.661
F28:P,E 101
E28 0 0.001 510I 1
F18:P,E 101
E18 0 0.001 510I 1
FT18 GEB 0 0.0567 0
NPS 1000000
PRINT