MCNP以及使用方法

1、MCNP以及使用方法總 目 錄MC方法解粒子輸運(yùn)問(wèn)題MCNP及使用粒子輸運(yùn)問(wèn)題粒子輸運(yùn)系統(tǒng)輸運(yùn)結(jié)束條件？所求量粒子源輸運(yùn)長(zhǎng)度粒子和物質(zhì)作用的核素、作用類型、作用后粒子能量和方向MC方法解粒子輸運(yùn)問(wèn)題逐一模擬每個(gè)粒子在輸運(yùn)系統(tǒng)中的輸運(yùn)過(guò)程，模擬足夠多的粒子（相對(duì)統(tǒng)計(jì)誤差可接受） （MC原理：大量、隨機(jī)性實(shí)驗(yàn)）計(jì)算粒子對(duì)所求物理量的平均貢獻(xiàn)（ MC原理：平均結(jié)果近似問(wèn)題的解） 一般情況下，當(dāng)模擬足夠多粒子后，認(rèn)為所模擬的粒子的平均行為近似為輸運(yùn)系統(tǒng)中粒子的平均行為；粒子對(duì)所求物理量的平均貢獻(xiàn)即為問(wèn)題的蒙卡解。開(kāi)始數(shù)據(jù)預(yù)處理, 各記錄單元清零取一個(gè)粒子歷史源分布抽樣輸運(yùn)過(guò)程碰撞過(guò)程歷史終止否？統(tǒng)計(jì)處

2、理做完給定歷史數(shù)否？結(jié)果的處理與輸出終止記錄過(guò)程記錄過(guò)程記錄過(guò)程記錄過(guò)程MC輸運(yùn)流程圖 根據(jù)問(wèn)題的幾何條件、物理假定、處理方法，可歸納為以下幾種：粒子從系統(tǒng)逃脫；粒子經(jīng)碰撞被吸收；經(jīng)俄國(guó)輪盤賭后，歷史被終止；粒子能量低于給定能量（閾能）；粒子位置越過(guò)某一界面；粒子飛行時(shí)間超過(guò)給定時(shí)間；粒子權(quán)重小于某個(gè)小量。粒子歷史終止條件總 目 錄MC方法解粒子輸運(yùn)問(wèn)題MCNP及使用MCNP目錄MCNP簡(jiǎn)介MCNP運(yùn)行MCNP輸入文件編寫舉例MCNPA General Monte Carlo Neutron and Photon Transport Code開(kāi)發(fā)： 美國(guó)洛斯阿拉莫斯國(guó)家實(shí)驗(yàn)室（ Los Ala

3、mos National Laboratory ）的蒙特卡羅小組MCNP發(fā)展歷史 MC方法研究（1948年開(kāi)始）MCS （1963年）MCN（1965年）中子，三維幾何，具有獨(dú)立的數(shù)據(jù)庫(kù)MCNG=MCN+MCG （1973）模擬高能光子MCNP=MCNG+MCP （1977年）模擬低能光子（能量達(dá)1KeV)Stan Ulam, John,von Neumann Nicholas MetropolisM.C 方法的命名It was Nicholas Metropolis who noted that Stan had an uncle who would borrow money from re

4、latives because he “just had to go to Monte Carlo” and thus named the mathematical method “Monte Carlo.”MCNP發(fā)展歷史MCNP3 （1983年）MCNP3A （1986年）MCNP3B （1988年）MCNP4 （1990年）MCNP4A（1993年）MCNP4B （1997年）（LAHET+HMCNP）MCNP4C （2000年）MCNPX：質(zhì)子、中子、光子、電子耦合（UNIX 第一個(gè)版本，具有并行功能，可以計(jì)算電子）（記數(shù)分析圖）可以計(jì)算微擾可以更好處理不可分辨共振MCNP5(2003

5、年)(ANSI standard Fortran 77)(ANSI standard Fortran 90)MCNP可處理問(wèn)題neutron onlyphoton onlyelectron onlycombined neutron/photon(n-p)neutron/photon/electronphoton/electronelectron/photonPhotonuclear （mcnp4c不能處理光核反應(yīng)）Keff（裂變系統(tǒng)）MCNP粒子能量限制neutron：10E-11 MeV to 20 MeV(all isotope)up to 150 MeV(some isotope)pho

6、ton: 1keV to 100GeVelectron: 1keV to 1GeV工作方式Monte Carlo程序的工作方式：用戶代碼型：EGS，GEANT，MARS輸入文件型：MCNP，F(xiàn)LUKAMCNP屬于第二種源碼開(kāi)放，用戶可以自己修改程序，重新編譯MCNP做題三步驟：已知、求/證明；解；答。MCNP（求解）輸入文件：告訴MC程序已知什么、求什么輸出文件：告訴我們答案是什么輸出文件介紹 程序運(yùn)行完之后，主要生成以下幾類輸出文件：標(biāo)準(zhǔn)輸出（out*）用于接續(xù)運(yùn)行(runtp*)統(tǒng)計(jì)結(jié)果輸出（mctal）mesh統(tǒng)計(jì)結(jié)果輸出（meshtal)輸出文件out*輸入文件重復(fù)區(qū)域、截面、物理設(shè)置

7、、動(dòng)態(tài)內(nèi)存信息輸運(yùn)過(guò)程信息輸運(yùn)完成后的結(jié)果綜述tally卡結(jié)果詳細(xì)信息tally卡結(jié)果綜合信息輸運(yùn)粒子數(shù)、時(shí)間信息輸出文件runtp*主要用于接續(xù)運(yùn)行運(yùn)行命令為：mcnp5 c i=*在*文件中第一行為：continue第二行給出運(yùn)行歷史（粒子數(shù)或運(yùn)行時(shí)間）注意：如果用nps，則nps代表總粒子數(shù)如果用ctme，則ctme代表繼續(xù)運(yùn)行時(shí)間MCNP目錄MCNP簡(jiǎn)介MCNP運(yùn)行MCNP輸入文件編寫舉例MCNP運(yùn)行運(yùn)行平臺(tái)Windows平臺(tái)下運(yùn)行Linux平臺(tái)下運(yùn)行常用指令：mcnp5 i=* o=* r=*mcnp5 i=* (輸出文件out*，runtp* 自動(dòng)從e字母開(kāi)始排)mcnp5 i=*

8、 r=runtp* c （接續(xù)運(yùn)行）mcnp5 i=* c (默認(rèn)為runtpe) （接續(xù)運(yùn)行）對(duì)于Mcnp4c版本，mcnp5mcnp，其余不變Windows平臺(tái)下運(yùn)行Linux平臺(tái)下運(yùn)行登錄服務(wù)器 user：flukau password：flukau軟件：SSHSecureShellClientSSH Secure SSH Secure Shell ClientLinux相關(guān)命令網(wǎng)上搜。運(yùn)行SSH Secure 使用方法同ftp添加profile登陸密碼同用戶名：flukau鼠標(biāo)右鍵，添加文件夾，用自己的名字命名；以后上傳、下載文件以及運(yùn)行程序都在自己的文件夾內(nèi)進(jìn)行。按圖示打開(kāi)輸入運(yùn)行命

9、令的窗口；也可以雙擊桌面的SSH Secure Shell Client，但后者要再次登錄。進(jìn)入自己的文件夾手動(dòng)刷新，顯示新生成的文件MCNP目錄MCNP簡(jiǎn)介MCNP運(yùn)行MCNP輸入文件編寫舉例粒子輸運(yùn)問(wèn)題的“已知” 和 “求”已知：輸運(yùn)系統(tǒng)已知：輸運(yùn)方式求：？粒子類型位置方向能量幾何描述介質(zhì)材料輸運(yùn)粒子類型粒子死亡條件運(yùn)算終止條件減方差技巧統(tǒng)計(jì)量統(tǒng)計(jì)方式已知：粒子源輸入文件輸入文件采用卡片結(jié)構(gòu)所謂的“卡片”就是一行輸運(yùn)粒子類型卡幾何相關(guān)卡源定義卡減方差技巧卡統(tǒng)計(jì)卡材料定義卡物理過(guò)程卡問(wèn)題截止卡用戶數(shù)組卡其它卡信息塊 可選空行分隔符可選標(biāo)題卡柵元卡 空行分隔符！面卡 空行分隔符！數(shù)據(jù)卡 空行終

10、止符 推薦使用其它 可選輸入文件舉例行輸入格式書寫規(guī)則柵元卡和面卡之間必須空一行!面卡和數(shù)據(jù)卡之間必須空一行！每行不超過(guò)80列；所有卡都必須從第15列開(kāi)始；如果15列為空白，則表示它是前一張卡的繼續(xù)卡； 而至少一個(gè)空格以后的符號(hào)&表示隨后的卡是繼續(xù)卡數(shù)據(jù)項(xiàng)由一個(gè)或多個(gè)空格分開(kāi)英文大小寫都可書寫規(guī)則注釋為c（行注釋符）或者$（末尾注釋符）nR代表重復(fù) 1 3R = 1 1 1 1nI代表遞加 1 2I 4 = 1 2 3 4xM代表遞乘 1 3M 3M = 1 3 9nJ代表跳過(guò)（默認(rèn)） 比如：DBCN X1 X20 DBCN 15J 8 3J 1(只給第16和第20個(gè)參數(shù)賦值，其它默認(rèn)) MC

11、NP輸入文件中物理量的單位長(zhǎng)度厘米能量MeV時(shí)間10-8 秒溫度MeV（kT）原子密度1024 個(gè)原子 / 厘米3質(zhì)量密度克 / 厘米3截面10-24 厘米2原子量中子質(zhì)量的1.008664967倍阿伏加德羅常數(shù)6.0231023 幾何描述類似于搭積木（面、塊搭成）面的布爾運(yùn)算(交、并、余) 柵元長(zhǎng)方體 六個(gè)平面交圓柱體 一個(gè)圓柱面和兩個(gè)平面交球體 一個(gè)球面不規(guī)則體：比如，球內(nèi)挖去一個(gè)小圓柱體剩余的部分面卡可定義的面包括：平面、球面、橢球面、雙曲面、拋物面、柱面、錐面、圓環(huán)面四種形式直接用面方程進(jìn)行定義（面編號(hào) 關(guān)鍵字 參數(shù)）使用點(diǎn)定義軸對(duì)稱面使用三點(diǎn)確定普通平面使用小體積元定義面BOX RP

12、P SPH RCC面方程定義曲面卡助記符類型說(shuō) 明方 程卡片上的數(shù)據(jù)項(xiàng)PPXPYPZ平面一般垂直X軸垂直Y軸垂直Z軸SOSSXSYSZ球面球心在原點(diǎn)一般球心在X軸球心在Y軸球心在Z軸C/XC/YC/ZCXCYCZ圓柱面平行于X軸平行于Y軸平行于Z軸軸心在X軸軸心在Y軸軸心在Z軸面方程定義曲面卡（續(xù)）助記符類型說(shuō) 明方 程卡片上的數(shù)據(jù)項(xiàng)K/XK/YK/ZKXKYKZ圓錐面平行于X軸平行于Y軸平行于Z軸軸心在X軸軸心在Y軸軸心在Z軸1只用于單葉錐面SQ橢球面雙曲面拋物面軸平行于X、Y或Z軸GQ圓柱面圓錐面橢球面雙曲面拋物面軸不平行于X、Y或Z軸助記符類型方 程卡片上的數(shù)據(jù)項(xiàng)TXTYTZ橢圓或圓形的

13、圓環(huán)面對(duì)稱軸平行于X、Y或Z軸面方程定義曲面卡（續(xù)）曲面的布爾運(yùn)算運(yùn)算符：交：空格并：冒號(hào):余：井號(hào)#曲面方向： 假定曲面 S 的曲面方程為 f (x,y,z)0 對(duì)于f (x,y,z)0的區(qū)域，為+S， “+”號(hào)可不寫； 對(duì)于f (x,y,z)0的區(qū)域，為-S。小體積元方向： 體積元內(nèi)為-，外為+幾何描述舉例長(zhǎng)方體：20*5*10面卡：1 px 02 px 5 3 py 0 4 py 205 pz 0 6 pz 10長(zhǎng)方體內(nèi)：1 -2 3 -4 5 -6長(zhǎng)方體外：-1:2:-3:4:-5:6xyz圓柱體：半徑5底面圓心(0,25,0)頂面圓心(0,45,0)面卡：1 cy 02 py 25

14、3 py 45 圓柱體內(nèi)：-1 2 -3圓柱體外：1:-2:3柵元卡幾何描述 + 材料定義柵元編號(hào) 材料號(hào) 材料密度 幾何描述 其它說(shuō)明其它說(shuō)明格式：關(guān)鍵詞=*，關(guān)鍵詞IMP、VOL、PWT、EXT、FCL、WWN、DXC、NONU、PD、TMP和關(guān)于重復(fù)結(jié)構(gòu)的U、TRCL、LAT、FILL。對(duì)于某些重復(fù)結(jié)構(gòu)可使用簡(jiǎn)寫形式： 柵元編號(hào) LIKE 柵元編號(hào)n BUT 參數(shù)列表 表示該柵元除了參數(shù)列表中列出的和柵元n不同，其它都同。 注意：柵元不能有相交的區(qū)域，不能有沒(méi)定義的區(qū)域。數(shù)據(jù)卡輸運(yùn)粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計(jì)卡物理過(guò)程卡（可選）問(wèn)題截止卡

15、用戶數(shù)組卡（可選）其它卡Mode卡Mode N /P /E，默認(rèn)為N注意： 輸運(yùn)的粒子類型和粒子源發(fā)射的粒子類型不是一回事。幾何相關(guān)卡VOL （默認(rèn)給出所有柵元的體積）AREA （默認(rèn)給出所有曲面的面積）重復(fù)結(jié)構(gòu)卡UTRCLLATFILLTR數(shù)據(jù)卡輸運(yùn)粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計(jì)卡物理過(guò)程卡（可選）問(wèn)題截止卡用戶數(shù)組卡（可選）其它卡源定義卡SDEF粒子信息：類型、位置、方向、能量格式： SDEF 源變量=說(shuō)明根據(jù)源變量特點(diǎn)，說(shuō)明有三種級(jí)別的表達(dá)方式：源變量=顯值 （包括單值、默認(rèn)）源變量=Dn （該變量值有個(gè)隨機(jī)分布，由SI和SP卡聯(lián)合給出其概

16、率分布）源變量=F 另一個(gè)源變量 Dn（該變量的值依賴于另外一個(gè)變量的分布，由DS卡先說(shuō)明有幾種分布）源變量變量含義缺 省CEL柵元。根據(jù)XXX，YYY，ZZZ和UUU，VVV，WWW確定。SUR曲面。0（柵元源）ERG粒子能量MeV）。14 MeVTME時(shí)間（ 10-8 秒）。0源變量（續(xù)）變量含義缺 省DIR，是VEC和方向UUU，VVV，WWW的夾角余弦。方位角總是在0360之間均勻抽取。體源：在11上均勻抽取 （各向同性源）。面源：在01上的余弦分布 f ()2 中抽取。VECDIR的參考矢量。體源：必須，各向同性除外。面源：垂直于曲面，符號(hào)為NRM。NRM曲面法向的符號(hào)。+1源變量（

17、續(xù)）變量含義缺 省POS抽樣位置的參考點(diǎn)。0,0,0RAD抽樣位置離開(kāi)POS或AXS的徑向距離。0EXT柵元源：沿AXS方向離開(kāi)POS的距離。曲面源：離開(kāi)AXS的角度余弦。0AXS用于EXT和RAD的參考矢量。沒(méi)有方向X抽樣位置的X坐標(biāo)。無(wú)XY抽樣位置的Y坐標(biāo)。無(wú)YZ抽樣位置的Z坐標(biāo)。無(wú)Z源變量（續(xù)）變量含義缺 省CCCCookie-cutter柵元。無(wú)ARA曲面的面積（當(dāng)平面源對(duì)點(diǎn)探測(cè)器要求直接貢獻(xiàn)時(shí)才需要此卡）。無(wú)WGT源粒子權(quán)重。1EFF位置選舍抽樣效率判據(jù)。0.01PAR源發(fā)射的粒子的類型。1中子2光子3電子1：包括中子的問(wèn)題。2：無(wú)中子但有光子問(wèn)題。3：僅有電子的問(wèn)題。SI卡給出源變

18、量分布區(qū)間的邊界值。格式： SIn選項(xiàng) I1 Ik n 分布號(hào)（n 1999）。選項(xiàng) Ii 的說(shuō)明，允許的值是： 省略或H：直方圖分布的區(qū)間邊界，僅用于標(biāo)量。L：離散的源變量值。A：定義概率分布密度的點(diǎn)。S：分布號(hào)。 I1 Ik 源變量值或分布號(hào)。SI源信息卡SP卡給出概率分布，微分或積分。格式： SPn選項(xiàng) P1 Pk 或： SPn f a b n 分布號(hào)（n 1999）。選項(xiàng) Pi 的說(shuō)明，允許的值是： 省略：對(duì)H或L分布與D相同，對(duì)A分布為概率密度。D： H或L分布的各個(gè)區(qū)間的概率值(不用歸一)。C： H或L分布的各個(gè)區(qū)間的累積概率值。V：僅用于柵元分布。概率與柵元體積成 正比（若有Pi

19、 ，要乘以Pi ）。P1 Pk 源變量概率。 f 內(nèi)部函數(shù)的標(biāo)識(shí)符。a，b 內(nèi)部函數(shù)的輸入?yún)?shù)。SP源概率卡舉例舉例1：位于原點(diǎn)的3MeV點(diǎn)光源，在錐角=120的圓錐內(nèi)均勻發(fā)射粒子，圓錐軸在z軸上。 sdef pos=0 0 0 erg=3 vec=0 0 1 dir=d1 par=2 si1 0.5 1 sp1 0 1舉例2：圓心在原點(diǎn)的球殼源，內(nèi)徑5cm，外徑15cm，各項(xiàng)同性發(fā)射，均勻分布的中子源 sdef pos=0 0 0 rad=d1 par=1 si1 5 15 sp1 -21 2 (球內(nèi)均勻分布的話，半徑符合r2規(guī)律)格式： DSn選項(xiàng) J1 Jk 或： DSnTI1 J1 I

20、k Jk 或： DSnQV1 S1 Vk Sk n 分布號(hào)（n 1999）。選項(xiàng) Ji 的說(shuō)明，允許的值是：省略或 H：連續(xù)分布的源變量值，僅用于標(biāo)量。L：離散的源變量值。S：分布號(hào)。 T 獨(dú)立變量值后面跟相關(guān)變量的值， 這些變量必須是離散的標(biāo)量。 Ii 獨(dú)立變量的值。 Ji 相關(guān)變量的值。DS卡復(fù)合分布SDEF POS=D1 ERG FPOS D2SI1L 5 3.3 6 75 3.3 6SP1.3 .7DS2S 3 4SI3H 2 10 14SP3D 0 1 2SI4-3 a b這是兩個(gè)各向同性點(diǎn)源，其中有的概率在(5,3.3,6 )處，能量服從階梯分布3；而有的概率在(75,3.3,6)

21、處，此時(shí)能量服從Watt 裂變譜分布4。舉例數(shù)據(jù)卡輸運(yùn)粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計(jì)卡物理過(guò)程卡（可選）問(wèn)題截止卡用戶數(shù)組卡（可選）其它卡材料定義卡材料定義需要知道核素構(gòu)成及其所占比例。格式：M材料編號(hào) 核素1 比例1 核素2 比例2 其它說(shuō)明核素表示方式：ZZZAAA.nnX或ZZZAAA AAA=0表示按天然元素取原子量nnX是截面庫(kù)標(biāo)識(shí)，表示從哪個(gè)截面庫(kù)取截面數(shù)據(jù)只有p或者p/e耦合輸運(yùn)，且忽略光核反應(yīng)，可以使用ZZZ000來(lái)簡(jiǎn)單描述；中子截面庫(kù)選擇比較難，最簡(jiǎn)單方法就是只給出ZZZAAA。比例：可歸一，可不歸一；正值表示原子數(shù)比例，負(fù)值表

22、示重量份額舉例：光子在水中輸運(yùn) m1 1000 2 8000 1數(shù)據(jù)卡輸運(yùn)粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計(jì)卡物理過(guò)程卡（可選）問(wèn)題截止卡用戶數(shù)組卡（可選）其它卡減方差技巧卡 IMP或WWN卡之一是必須的。助記符卡片類型IMP柵元重要性ESPLT能量分裂和俄國(guó)輪盤賭PWT次級(jí)光子權(quán)重EXT指數(shù)變換VECT方向矢量定義FCL強(qiáng)迫碰撞助記符卡片類型WWE權(quán)重窗的能量或時(shí)間間隔WWN權(quán)重窗的邊界WWP權(quán)重窗的參數(shù)WWG權(quán)重窗生成器WWGE權(quán)重窗生成器的能量或時(shí)間間隔MESH分層重要性網(wǎng)格權(quán)窗生成器PD探測(cè)器貢獻(xiàn)DXCDXTRAN貢獻(xiàn)BBREM韌致輻射偏倚因

23、子格式： IMP:nx1 x2 xi xI n 中子為 N，光子為 P，電子為 E。 N, P 、 P, E 或 N, P , E 也是允許的，如果它們的值相同。xi 柵元 i 的重要性，i 1，2，II 問(wèn)題中的柵元總數(shù)。 缺?。涸谝粋€(gè)MODE N P 問(wèn)題中，若省略了IMP:P卡，則所有柵元的光子重要性都置 1，除非其中子重要性為零，這時(shí)其光子重要性也為零。柵元重要性卡用于輸入各個(gè)柵元的重要性。零重要性可以用來(lái)終止粒子的歷史。IMP卡數(shù)據(jù)卡輸運(yùn)粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計(jì)卡物理過(guò)程卡（可選）問(wèn)題截止卡用戶數(shù)組卡（可選）其它卡統(tǒng)計(jì)卡Fna 統(tǒng)

24、計(jì)類型En能量間隔Tn時(shí)間間隔Cn角度間隔FMn計(jì)數(shù)乘子DEn劑量能量間隔DFn劑量轉(zhuǎn)換系數(shù)EMn能量相乘因子CMn劑量相乘因子FTn統(tǒng)計(jì)卡特殊處理參數(shù)FMesh 網(wǎng)格計(jì)數(shù)Fna統(tǒng)計(jì)類型卡描述Fn單位*Fn單位F1:N/ P/ E面流量/穿透率particles MeVF2:N/ P/ E統(tǒng)計(jì)面通量particles/cm2MeV/cm2F4:N/ P/ E統(tǒng)計(jì)體通量particles/cm2MeV/cm2F5a:N/ P統(tǒng)計(jì)點(diǎn)通量或環(huán)通量particles/cm2MeV/cm2F6:N/ P體能量沉積MeV/g 109 J /gF7:N體裂變能沉積MeV/g 109 J /gF8:P E探測(cè)

25、器脈沖幅度譜pulsesMeV+F8:E電荷沉積charge注意：有些是針對(duì)面，有些是針對(duì)體。舉例舉例1：F1: N 1 3 6 T 分別計(jì)算通過(guò)曲面1、3、6的穿透率以及穿過(guò)三個(gè)曲面的總穿透率。舉例2： F34: P 9 計(jì)算柵元9內(nèi)的體通量。注意：有多個(gè)同類計(jì)數(shù)卡時(shí)，需要在卡號(hào)前加不同數(shù)字進(jìn)行區(qū)分。舉例3：F5: P 計(jì)算點(diǎn)( 0, 20, 15)的點(diǎn)通量。舉例4： *F8: P 10 計(jì)算柵元10內(nèi)的能量沉積。En能量間隔卡分能量段給出Fn的計(jì)數(shù)，否則只給出總計(jì)數(shù)。格式：En E1 Ek注意：若使用此卡，n必須與相應(yīng)的Fn的n一致。舉例：E11 .1 1 20 在四個(gè)能量間隔內(nèi)統(tǒng)計(jì)F11

26、計(jì)數(shù)：截止能量0.1 MeV；1 MeV；120 MeV；整個(gè)能量范圍，即截止能量20 MeV關(guān)于F8不論 粒子類型取何值，F(xiàn)8都將記錄光子和電子，也就是說(shuō)F8:p、 F8:e、 F8:p e結(jié)果一樣。注意：選擇能量間隔的時(shí)候，應(yīng)該把零區(qū)間和(1E-5)區(qū)間包括在里面。 舉例：E80 1E-5 1E-3 1E-1 零區(qū)間用來(lái)處理由于撞出電子而產(chǎn)生負(fù)計(jì)數(shù)的情況；(1E-5) 區(qū)間內(nèi)將記錄穿過(guò)柵元但沒(méi)有沉積能量的粒子。計(jì)算出通量后，根據(jù)DF卡給出的通量-劑量轉(zhuǎn)化系數(shù)，求得劑量。格式： DEn A E1 Ek DFn B F1 Fk n 計(jì)數(shù)號(hào)。 Ei 第 i 個(gè)能量點(diǎn)（MeV）。 Fi 對(duì)應(yīng)第 i

27、 個(gè)能量點(diǎn)的劑量函數(shù)值。 A LOG或LIN，能量表按對(duì)數(shù)或線性插值。 B LOG或LIN，劑量函數(shù)表的插值方法。缺省：若沒(méi)有填寫A或B，則按對(duì)數(shù)插值處理。劑量計(jì)算DE、DF卡格式： FTn ID1 P1,1 P1,2 P1,3 ID2 P2,1 P2,2 P2,3 n 計(jì)數(shù)號(hào)。 IDi 第 i 項(xiàng)特殊處理的標(biāo)識(shí)符。 FRV 用于計(jì)數(shù)類型1角度間隔的參考方向。GEB 高斯能量展寬。TMC 時(shí)間卷積。INC 指定碰撞次數(shù)。ICD 指定要記錄對(duì)探測(cè)器作貢獻(xiàn)的柵元。ELC 電流計(jì)數(shù)。 Pi,j 第 i 項(xiàng)特殊處理的參數(shù)。FTn 計(jì)數(shù)特殊處理卡數(shù)據(jù)卡輸運(yùn)粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料

28、定義卡減方差技巧卡（可選）統(tǒng)計(jì)卡物理過(guò)程卡（可選）問(wèn)題截止卡用戶數(shù)組卡（可選）其它卡物理過(guò)程卡助記符卡片類型PHYS能量物理截?cái)郥MP自由氣體模型熱溫度THTME熱時(shí)間卡MTmS(,) 材料卡用于控制MCNP的能量以及其它物理狀況。數(shù)據(jù)卡輸運(yùn)粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計(jì)卡物理過(guò)程卡（可選）問(wèn)題截止卡用戶數(shù)組卡（可選）其它卡問(wèn)題截止卡在初始運(yùn)行或接續(xù)運(yùn)行的輸入文件中均可使用，用于終止粒子的歷史或中斷計(jì)算。助記符卡片類型CUT粒子截?cái)郋LPT柵元能量截?cái)郚PS歷史數(shù)截?cái)郈TME計(jì)算時(shí)間截?cái)郈UT粒子截?cái)嗫ǜ袷剑篊UT:n T E WC1 WC2 S

29、WTMWC1、WC2、SWTM涉及到減方差技巧，先不講。T 指截?cái)鄷r(shí)間（10-8 秒），缺省值：中子-無(wú)限大，光子和電子都是等于中子的截?cái)鄷r(shí)間。E 指截?cái)嗄芰浚∕eV），缺省值：中子，光子和電子都是。舉例：CUT:p j 0.01 3j 截?cái)嗄芰繛椋溆嗳笔≈?。格式：NPS N N 粒子的歷史數(shù)。 該卡指定要計(jì)算的粒子歷史數(shù)，當(dāng)計(jì)算的粒子歷史數(shù)達(dá)到指定值時(shí)，MCNP將中止計(jì)算。NPS歷史數(shù)截?cái)嗫ǜ袷剑篊TME x x 計(jì)算時(shí)間（分鐘）。 該卡指定該問(wèn)題在計(jì)算機(jī)上花費(fèi)的時(shí)間限制，超過(guò)該時(shí)間，MCNP將中止計(jì)算。 所有中止計(jì)算的條件只要達(dá)到一個(gè)就中止計(jì)算。CTME計(jì)算時(shí)間截?cái)嗫∕CNP目錄MCNP簡(jiǎn)介MCNP運(yùn)行