核辐射防护软件：MCNP二次开发_（9）.材料定义与赋值.docx

PAGE1

PAGE1

材料定义与赋值

在MCNP（MonteCarloN-ParticleTransportCode）中，材料定义是模拟核辐射传输过程的基础。MCNP使用材料卡片（M-card）来定义材料的组成和属性。本节将详细介绍如何在MCNP中定义材料，并赋予相应的属性，以便进行准确的核辐射传输模拟。

材料卡片的基本格式

材料卡片的基本格式如下：

MnZA1W1ZA2W2...

Mn：材料卡片的标识，其中n是材料的编号。

ZA1：第一个核素的原子序数和质量数的组合，格式为Z*1000+A，例如1001表示氢-1。

W1：第一个核素的重量分数。

ZA2：第二个核素的原子序数和质量数的组合。

W2：第二个核素的重量分数。

...：更多的核素及其重量分数。

材料定义的示例

单一核素材料

定义一个由纯氢-1组成的材料：

M110011.0

多核素材料

定义一个由氢-1和氧-16组成的水材料：

M210010.11189880160.888102

密度赋值

在定义材料时，通常还需要指定材料的密度。MCNP使用TRCL卡片来定义材料的密度，格式如下：

TRClD

TRCl：密度卡片的标识，其中l是材料的编号。

D：材料的密度，单位为g/cm3。

例如，定义水材料的密度为1.0g/cm3：

TRC21.0

温度赋值

材料的温度也是影响核辐射传输的重要因素。MCNP使用TMP卡片来定义材料的温度，格式如下：

TMPnT

TMPn：温度卡片的标识，其中n是材料的编号。

T：材料的温度，单位为K（开尔文）。

例如，定义水材料的温度为300K：

TMP2300.0

材料的混合

在实际应用中，材料可能由多种成分混合而成。MCNP支持多种混合方法，包括按重量分数混合和按原子分数混合。

按重量分数混合

按重量分数混合是最常见的方法。例如，定义一个由50%的氢-1和50%的氧-16组成的材料：

M310010.580160.5

按原子分数混合

按原子分数混合时，需要使用-1来标识权重为原子分数。例如，定义一个由2个氢原子和1个氧原子组成的水材料：

M410012.0-180161.0-1

材料库的使用

MCNP提供了多种材料库，可以在输入文件中直接引用。材料库的使用方法如下：

MnLIBZA

Mn：材料卡片的标识，其中n是材料的编号。

LIB：材料库的标识，例如PE表示塑料，H2O表示水。

ZA：材料库中的核素标识。

例如，使用MCNP的水材料库定义水材料：

M5H2O

材料的微观截面数据

为了进行准确的模拟，MCNP需要核素的微观截面数据。这些数据通常包含在XS库文件中。在输入文件中指定XS库文件的方法如下：

*LIBRARY

ала

*LIBRARY：指定XS库的卡片。

ала：XS库文件的名称，例如ncolib。

例如，指定使用ncolib库文件：

*LIBRARY

ncolib

材料的多组分赋值

在某些复杂的应用中，材料可能由多个组分组成。MCNP支持多组分材料的定义，方法是在材料卡片中列出所有组分及其相应的重量分数。

例如，定义一个由水、塑料和空气组成的多组分材料：

M6H2O0.7PE0.2AIR0.1

材料的动态赋值

在某些情况下，材料的属性可能随时间变化。MCNP支持动态赋值，方法是使用时间依赖的密度或温度卡片。格式如下：

TRCl(t)D(t)

TMPn(t)T(t)

TRCl(t)：时间依赖的密度卡片。

D(t)：时间依赖的密度函数。

TMPn(t)：时间依赖的温度卡片。

T(t)：时间依赖的温度函数。

例如，定义一个材料的密度随时间线性变化：

TRC7(t)1.0+0.01*t

材料的微观截面数据生成

在某些情况下，可能需要生成特定材料的微观截面数据。MCNP提供了工具和方法来生成这些数据。具体步骤如下：

准备核素数据：确保有所有需要的核素数据。

生成截面数据：使用MCNP的工具生成截面数据文件。

引用截面数据：在输入文件中引用生成的截面数据文件。

生成截面数据的示例

假设我们需要生成一个由氢-1和氧-16组成的水材料的微观截面数据。步骤如下：

准备核素数据：确保ncolib库文件中有氢-1和氧-16的数据。

生成截面数据：使用MCNP的make_xsec工具生成截面数据文件。

make_xsec-mH2O-omy_xsec.nc

-mH2O：指定材料为水。

-omy_xsec.nc：指定输出文件为my_xsec.nc。

引用截面数据：在MCNP输入文件中引用生成的截面数据文件。

*LIBRARY

my_xsec.nc

材料的物理属性

除了基本的核素组成和密度外，材