MCNP使用說(shuō)明--經(jīng)典版

1、MCNP4c3說(shuō)明堆工所21室何濤（翻譯自UWMadison教程）MCNP摘要n簡(jiǎn)介簡(jiǎn)介n輸入文件基礎(chǔ)輸入文件基礎(chǔ)n幾何描述幾何描述n源的描述源的描述n計(jì)數(shù)描述計(jì)數(shù)描述n材料，物理成分和數(shù)據(jù)材料，物理成分和數(shù)據(jù)n高級(jí)計(jì)數(shù)描述高級(jí)計(jì)數(shù)描述n高級(jí)幾何描述高級(jí)幾何描述n各種簡(jiǎn)化各種簡(jiǎn)化n臨界問(wèn)題臨界問(wèn)題MCNP簡(jiǎn)介n對(duì)Unix的熟悉n運(yùn)行MCNPn繪制MCNP幾何圖象nMCNP輸入文件結(jié)構(gòu)學(xué)習(xí)目的MCNP簡(jiǎn)介n懂得怎樣用Unix命令行運(yùn)行MCNPn懂得MNCP文件名的慣例n能夠使用繪圖工具描述幾何形狀運(yùn)行MCNPn% mcnp i= o= optionsnOptions(選項(xiàng)）i 處理輸入文件默認(rèn)值

2、p 繪圖x 處理截面默認(rèn)值r 粒子傳輸默認(rèn)值z(mì) 標(biāo)繪計(jì)數(shù)結(jié)果 標(biāo)繪截面運(yùn)行MCNPn% mcnp i= o= 選項(xiàng)n默認(rèn)文件名inp輸入文件outpASCII輸出文件runtpe二進(jìn)制重啟文件n通過(guò)命令行改變默認(rèn)值% mcnp inp=exl oupt= exlo runexlr% mcnp name=exl 練習(xí)1an運(yùn)行一次% mcnp i=demo1 什么文件被創(chuàng)建？n再運(yùn)行一次% mcnp i=demo1 這時(shí)候什么文件被創(chuàng)建？不要使用默認(rèn)文件名n始終清楚地定義文件名或者i=inName o=outName r=runNamen=baseNamei=inname n=baseOutNa

3、men這將會(huì)防止你覆蓋先前地計(jì)算結(jié)果n這將會(huì)幫助你知道哪一個(gè)結(jié)果是正確的練習(xí)1bn運(yùn)行一次% mcnp n=demo1 什么文件被創(chuàng)建？n再運(yùn)行一次% mcnp i=demo1 ndemo1_ 這時(shí)候什么文件被創(chuàng)建？繪制幾何圖象n計(jì)算機(jī)上的二維幾何圖象顯示n能夠用來(lái)檢查幾何問(wèn)題的很多方面：柵元和表面序號(hào)材料密度材料位置n幾何錯(cuò)誤用紅色虛線顯示n經(jīng)常繪圖檢查幾何結(jié)構(gòu)練習(xí)2an繪制圖象 mcnp i=demo1 n= demo1_ ip概念圖象放大全景顯示改變方向驗(yàn)證材料，柵元，表面，密度等學(xué)習(xí)目標(biāo)：輸入文件基礎(chǔ)n了解MCNP中的物理單位n了解MCNP輸入文件三個(gè)主要部分n了解MCNP輸入文件的格

4、式規(guī)定n了解MCNP輸入文件的簡(jiǎn)寫(xiě)特點(diǎn)MCNP中的物理單位n長(zhǎng)度：cmn能量：MeVn時(shí)間：剎（10-8s）n溫度：MeV(KT)n原子密度：1024原子/cm3n質(zhì)量密度：g/cm3n截面：靶（10-24cm2）MCNP輸入文件n標(biāo)題卡n柵元卡要求空行分隔n曲面卡要求空行分隔n數(shù)據(jù)卡推薦空行作為結(jié)束MCNP輸入文件格式n每行最多80個(gè)字符n不含控制字符，比如：Tabn注釋行：標(biāo)題卡之后的任何位置都可插入第一列是字母“C”，且隨后四個(gè)空格從輸入數(shù)據(jù)之后的$符號(hào)后開(kāi)始n以上三種情況可以單獨(dú)或同時(shí)存在輸入簡(jiǎn)寫(xiě)nnR：表示將它前面的一個(gè)數(shù)據(jù)重復(fù)n遍例如：2 4R = 2 2 2 2 2nnI：表示在

5、與它前后相鄰的兩個(gè)數(shù)之間插入n個(gè)線性插值點(diǎn)。例如：1 5I 7 = 1 2 3 4 5 6 7nxM：表示它前面的數(shù)據(jù)與x之積例如：5 4M = 4 20nnJ：表示從它所在位置跳過(guò)n項(xiàng)不指定的數(shù)據(jù)而使用缺省值。輸入簡(jiǎn)寫(xiě)規(guī)則n如果n(R,I,J)中的n缺省，則假設(shè)n=1。n如果xM中的x缺省，則致命錯(cuò)誤。nnR前面必須放有一個(gè)數(shù)或由R或M產(chǎn)生的數(shù)據(jù)項(xiàng)。nnI前面必須放有一個(gè)數(shù)或由R或M產(chǎn)生的數(shù)據(jù)項(xiàng)，后面還要跟有一個(gè)常數(shù)。nxM前面必須放有一個(gè)數(shù)或由R或M產(chǎn)生的數(shù)據(jù)項(xiàng)。nnJ前面可以放除了I以外的任何內(nèi)容。輸入文件MCNP教程n簡(jiǎn)介簡(jiǎn)介n輸入文件基礎(chǔ)輸入文件基礎(chǔ)n幾何描述幾何描述n源的描述源的描

6、述n計(jì)數(shù)描述計(jì)數(shù)描述n材料，物理成分和數(shù)據(jù)材料，物理成分和數(shù)據(jù)n高級(jí)計(jì)數(shù)描述高級(jí)計(jì)數(shù)描述n高級(jí)幾何描述高級(jí)幾何描述n各種簡(jiǎn)化各種簡(jiǎn)化n臨界問(wèn)題臨界問(wèn)題MCNP幾何n幾何基礎(chǔ)n快速開(kāi)始n曲面n組合曲面nMacrobodiesn柵元特性n例子學(xué)習(xí)目標(biāo)：幾何n懂得四種定義曲面的方法n懂得怎樣由曲面創(chuàng)建柵元n了解Macrobodies的定義細(xì)節(jié)n懂得進(jìn)行曲面變換n懂得何時(shí)使用特殊曲面幾何基礎(chǔ)n“universe”根據(jù)材料和特性被分成不同的區(qū)域n整個(gè)無(wú)窮的universe必須包括在幾何模型之內(nèi)n幾何的基本單位是柵元n所有的柵元都由閉合曲面定義n所有的曲面都能將universe分成兩部分曲面n由方程定義曲

7、面由方程及參數(shù)確定例如：一個(gè)球心在原點(diǎn)半徑為R的球 j so R平行于y軸半徑為R的圓柱j c/y x z R垂直于z軸的平面j pz z柵元中的復(fù)合曲面n柵元中的點(diǎn)和曲面的關(guān)系通過(guò)柵元對(duì)曲面的坐向聯(lián)系起來(lái)：“+”和“-”曲面將universe分為兩個(gè)半?yún)^(qū)n布爾算符將不同的半?yún)^(qū)與創(chuàng)建的柵元聯(lián)系起來(lái)交（Intersection)聯(lián)（Union)余（Complement)坐向n柵元中所有的點(diǎn)都通過(guò)坐向與定義柵元的曲面聯(lián)系起來(lái)。坐向說(shuō)明了柵元中的點(diǎn)在曲面的那一邊n+ 正的坐向?qū)τ陂_(kāi)放曲面（平面)，點(diǎn)在曲面的坐標(biāo)軸正方向；對(duì)于封閉的曲面（球，圓柱等），點(diǎn)在曲面以外。n+ 負(fù)的坐向?qū)τ陂_(kāi)放曲面（平面)，

8、點(diǎn)在曲面的坐標(biāo)軸負(fù)方向；對(duì)于封閉曲面（球，圓柱等），點(diǎn)在曲面以?xún)?nèi)。柵元的復(fù)合曲面：交n同時(shí)滿(mǎn)足兩個(gè)坐向的空間n算符輸入：在兩個(gè)曲面號(hào)中用空格n2 1只表示同時(shí)滿(mǎn)足坐向+2和坐向-1的空間區(qū)域柵元的復(fù)合曲面：聯(lián)n任意滿(mǎn)足兩個(gè)坐向之一的空間n算符輸入：在兩個(gè)曲面號(hào)中用冒號(hào):n2:1表示任意滿(mǎn)足坐向+2和坐向-1之一的空間區(qū)域柵元的復(fù)合曲面：余n表示柵元之外的空間n算符輸入：在曲面號(hào)前用#n#5表示柵元5之外的空間n余以后的區(qū)域可以和 其他區(qū)域進(jìn)行交和并的運(yùn)算n-2 #5代表曲面2之內(nèi)且在曲面5之外的區(qū)域。柵元練習(xí)n截面測(cè)量實(shí)驗(yàn)n中央有一個(gè)圓柱孔(R=2cm)的鈹球(R=25cm)。n孔的中央有一個(gè)

9、氚靶（厚0.5cm，R=2cm）n鈹球嵌在半徑為40cm的混凝土球殼內(nèi)。四類(lèi)MCNP曲面n方程定義曲面n平面，球，圓柱，圓錐，圓環(huán)，任意的二次曲面nMacrobodiesn基于閉合圖元的復(fù)雜曲面n用點(diǎn)定義對(duì)稱(chēng)曲面n平面，線性曲面或二次曲面n由三個(gè)點(diǎn)定義一般平面n一般平面MCNP曲面：方程定義n基本格式（Chapter 3, section III.A, Table 3.1, p 3-12)nj n a listj：曲面號(hào)：19999n：缺省值為0，表示不進(jìn)行任何坐標(biāo)變換 0：用TRn卡對(duì)曲面坐標(biāo)變換 0：曲面j是伴隨曲面n的周期邊界a：方程助記名list：方程描述的數(shù)據(jù)項(xiàng)MCNP曲面：圓錐n圓

10、錐的等式定義了兩個(gè)“葉”。n參數(shù)中額外的條目是用來(lái)區(qū)分“正葉”和“負(fù)葉”的n只有在圓錐平行于軸的時(shí)候才有效MCNP曲面：Macrobodiesn有限的“模塊”構(gòu)成的曲面（Chapter 3, Section III.D, Table 3.1, p 3-12)nBOX任意指向的正交框nRPP直角平行六面體n所有的表面垂直于各自的軸nSPH球n與方程表示的球是一樣的nRCC直圓柱體n軸與底面垂直，但是方向任意nRHP(HEX)直六面棱柱n與RCC相似但底為任意的六邊形MCNP曲面：Macrobodiesn“坐向”與其他封閉曲面相類(lèi)似n+ 正的坐向，點(diǎn)在曲面以外。n+ 負(fù)的坐向，點(diǎn)在曲面以?xún)?nèi)。n能夠

11、與其他類(lèi)型的曲面相復(fù)合n從能夠被分別索引的“面”構(gòu)造Macrobody 的“面”n面是按順序編號(hào)的（見(jiàn)說(shuō)明書(shū)page 3-21）n參考使用Macrobody編號(hào)和“面”編號(hào)nMacrobody RCC的圓柱側(cè)面 j=5n5.1nMacrobody RPP的ymax平面 j=10n10.3MCNP曲面：用點(diǎn)定義對(duì)稱(chēng)曲面n用面上的一到三個(gè)點(diǎn)描述，且面關(guān)于X,Y或Z軸對(duì)稱(chēng)。（見(jiàn)Page3-16）n每一對(duì)坐標(biāo)點(diǎn)定義曲面上的一個(gè)點(diǎn)n第一個(gè)坐標(biāo)：點(diǎn)離軸的距離n第二個(gè)坐標(biāo)：點(diǎn)離軸的半徑MCNP曲面：用點(diǎn)定義對(duì)稱(chēng)曲面n一個(gè)點(diǎn)：定義一個(gè)平面n兩個(gè)點(diǎn)：定義平面或者線性曲面（圓柱、圓錐）n三個(gè)點(diǎn)：定義平面、線性曲面

12、或者二次曲面（球或者一般的二次曲面）n所有的點(diǎn)都在同一葉上n所有的葉必須是可定義的如平面、線性曲面或者二次曲面MCNP曲面：用三點(diǎn)定義平面n任意三點(diǎn)定義一個(gè)平面n所產(chǎn)生的平面方程系數(shù)遵循原則n原點(diǎn)是負(fù)方向n(0,0,)是正方向n(0,0)是正方向n(,0,0)是正方向致命錯(cuò)誤曲面坐標(biāo)變換nTRn坐標(biāo)變換卡（見(jiàn)page3-30）nTRn Ox Oy Oz Bxx Byx Bzx Bxy Byy Bzy Bxz Byz Bzz Mnn：變換號(hào)，與曲面匹配nOx Oy Oz，變換的原點(diǎn)位移矢量nBxx Bzz變換的旋轉(zhuǎn)矩陣（余弦或度 *）nMn1意味著位移矢量是輔助坐標(biāo)系的原點(diǎn)在基本坐標(biāo)系定義的位置

13、n1意味著位移矢量是基本坐標(biāo)系的原點(diǎn)在輔助坐標(biāo)系定義的位置曲面坐標(biāo)變換n有時(shí)候?qū)?biāo)準(zhǔn)曲面進(jìn)行坐標(biāo)變換比直接定義一個(gè)復(fù)雜曲面更加容易n例子：軸平行于(1,1,0)的圓柱n怎樣直接寫(xiě)出等式定義？n替代方法n用等式定義x軸的圓柱n進(jìn)行45度角的變換練習(xí)n繼續(xù)上一章的練習(xí)n使用macrobodies將混凝土球殼替換成立方體n在鈹球內(nèi)使用圓錐形孔代替圓柱孔，在靶處半徑2cm，外表面半徑4cm?；炷燎驓ぬ幦詾榘霃?cm的圓柱孔。n需要用到半徑1cm，與原來(lái)的孔夾角為60度的圓柱孔進(jìn)行檢驗(yàn)n這個(gè)檢驗(yàn)孔是號(hào)角狀的，靠近靶處半徑1cm，靠近鈹球表面處半徑2cm，中間半徑1.2cm。MCNP柵元：柵元基本特性n

14、柵元不僅僅是指幾何形狀，還包括n材料：定義柵元中用來(lái)輸運(yùn)和反應(yīng)的截面n重要性：n基本用途：把無(wú)用的universe和物理模型分離開(kāi)來(lái)n高級(jí)用途：改進(jìn)問(wèn)題的統(tǒng)計(jì)結(jié)果MCNP柵元：材料定義nMn材料卡在輸入文件的數(shù)據(jù)卡部分（見(jiàn)pg 3-108)n提供了材料所含元素或同位素的原子比例或重量百分比nMn zaid1 frac1 zaid2 frac2 zaidn fracnnn材料號(hào)，與柵元卡中條目匹配nzaid根據(jù)原子序數(shù)和原子量定義的同位素IDn通常：zaid=Z*1000+AnA=0，代表天然元素n對(duì)于特殊的截面庫(kù)有可選條目zaid.xsidnfrac元素在材料中的原子比例或重量百分比（若為負(fù)值

15、)nMCNP將會(huì)自動(dòng)歸一化MCNP柵元：材料定義示例nM1 92235 4.5 92238 95.5 8016 13.5n含豐度為4.5%的鈾235的氧化鈾（核燃料）nM25 7000 78 8000 21n空氣的近似nM1 1001 0.5 8016 0.25 6012 0.25n少量的雜質(zhì)一般是不重要的n鋼含有23種不同的元素nM4 26000 88.8 24000 9 74000 2 25000 0.5 14000 0.25 6000 0.1MCNP柵元：重要性n每個(gè)柵元都有重要性n標(biāo)準(zhǔn)的重要性為1n不同的數(shù)值用來(lái)表示重要性的遞減n如果重要性為0，粒子在此柵元中不予考慮n終結(jié)粒子歷史n剩

16、下的universe重要性為0MCNP柵元：重要性的定義nIMP: n, IMP: p, IMP: e, IMP: n,p 等n兩種方法定義重要性n 在柵元描述卡中，曲面描述之后定義n 1 3-8-1 2 # (-3 5) IMP: n=1n 4 0 1 : -4 5IMP: n=0n 作為數(shù)據(jù)卡，每個(gè)柵元對(duì)應(yīng)一個(gè)數(shù)字nIMP: n 1 1 1 0$ For 4 cells整個(gè)歷史的描述n源的定義n這是必須的:沒(méi)有源就沒(méi)有粒子n問(wèn)題截?cái)鄺l件n可選的:沒(méi)有截?cái)鄺l件，問(wèn)題將永遠(yuǎn)計(jì)算下去n計(jì)數(shù)n可選的:沒(méi)有計(jì)數(shù)就不會(huì)知道任何計(jì)算結(jié)果源的描述n通用源卡 SDEF 定義了如下內(nèi)容n粒子在哪里創(chuàng)建n柵元，

17、曲面，(x,y,z)n粒子在何時(shí)創(chuàng)建n粒子的能量和方向n粒子權(quán)重n粒子類(lèi)型n默認(rèn)值n在零點(diǎn)零時(shí)刻創(chuàng)建，E=14MeV，各向同性，權(quán)重為1問(wèn)題截?cái)嗫╪兩種主要的截?cái)鄋粒子數(shù)目NPSn計(jì)算時(shí)間CTME （以分為單位）n個(gè)別粒子的其他截?cái)鄋CUT 對(duì)每個(gè)粒子定義了最大時(shí)間，最小能量等nELPT對(duì)每個(gè)柵元定義了最小能量計(jì)數(shù)n用來(lái)提供物理量的估計(jì)計(jì)數(shù)nFn: a j1 j2 j3 jN Tnn計(jì)數(shù)號(hào)，最后一位定義計(jì)數(shù)類(lèi)型na計(jì)數(shù)粒子標(biāo)志符：n, p, enji計(jì)數(shù)所在的區(qū)域nT可選：所有計(jì)數(shù)區(qū)域合集的計(jì)數(shù)計(jì)數(shù)類(lèi)型n表面計(jì)數(shù)n1穿過(guò)曲面的積分流量n2穿過(guò)曲面的平均通量n柵元計(jì)數(shù)n4一個(gè)柵元的平均通量n6一

18、個(gè)柵元的平均沉積能量n7一個(gè)柵元的平均裂變沉積能量n特殊計(jì)數(shù)n5一個(gè)點(diǎn)或環(huán)探測(cè)器的通量n8一個(gè)脈沖幅度探測(cè)器的通量練習(xí)n以前面的問(wèn)題為基礎(chǔ)建立完整的模型n加入材料定義和柵元重要性n加入默認(rèn)源n計(jì)算100000個(gè)粒子n對(duì)鈹球外表通量計(jì)數(shù)MCNP幾何總結(jié)n四種定義曲面的方法n方程，macrobodies，一般平面，對(duì)稱(chēng)旋轉(zhuǎn)n邏輯算符把曲面與柵元聯(lián)系起來(lái)n邏輯算符：交，并，余n柵元說(shuō)明還需要材料和重要性nmacrobodies是由“面”構(gòu)成的n“面”類(lèi)似其他曲面，但被內(nèi)在描述n曲面可以任意進(jìn)行坐標(biāo)變換MCNP教程n簡(jiǎn)介簡(jiǎn)介n輸入文件基礎(chǔ)輸入文件基礎(chǔ)n幾何描述幾何描述n源的描述源的描述n計(jì)數(shù)描述計(jì)數(shù)描

19、述n材料，物理成分和數(shù)據(jù)材料，物理成分和數(shù)據(jù)n高級(jí)計(jì)數(shù)描述高級(jí)計(jì)數(shù)描述n高級(jí)幾何描述高級(jí)幾何描述n各種簡(jiǎn)化各種簡(jiǎn)化n臨界問(wèn)題臨界問(wèn)題MCNP源的定義n源的粒子特性n獨(dú)立分布n體積源和曲面源n內(nèi)嵌函數(shù)n非獨(dú)立分布n源的檢查學(xué)習(xí)目標(biāo)：源n知道粒子源必要的增殖特性n能夠使用各種獨(dú)立的樣本分布n能夠在笛卡爾坐標(biāo)，柱坐標(biāo)和球坐標(biāo)下定義體積源和球面源n懂得使用內(nèi)嵌的PDF函數(shù)n能夠使用非獨(dú)立樣本分布n能夠檢查使源的產(chǎn)生正確源的描述n通用源卡 SDEF 定義了如下內(nèi)容n粒子在哪里創(chuàng)建n柵元，曲面，(x,y,z)n粒子在何時(shí)創(chuàng)建n粒子的能量和方向n粒子權(quán)重n粒子類(lèi)型n默認(rèn)值n在零點(diǎn)零時(shí)刻創(chuàng)建，E=14MeV，

20、各向同性，權(quán)重為1通用源卡基本語(yǔ)法(1)nSDEFvar1=spec1 var2=spec2 (p3-50)標(biāo)量向量(x y z)位置X, Y, ZPOS時(shí)間TME能量ERG方向DIRVEC權(quán)重WGT類(lèi)型PAR幾何約束CEL, SUR通用源卡基本語(yǔ)法(2)n變量說(shuō)明有三種形式n顯值nERG=2n能量為2MeVn分布號(hào)nERG=D1n能量由1號(hào)分布定義n相關(guān)分布nERG FPOS D2n能量是位置的函數(shù)，且從分布2抽取源的粒子類(lèi)型n只允許一種粒子nPAR= (必須清楚，且沒(méi)有分布）n1中子n2質(zhì)子n3電子nMODE卡（問(wèn)題類(lèi)型卡）必須支持粒子類(lèi)型n默認(rèn)值是MODE卡支持的最小數(shù)字獨(dú)立分布n四種概

21、率分布n柱狀圖獨(dú)立分布n四種概率分布n柱狀圖n離散的獨(dú)立分布n四種概率分布n柱狀圖n離散的n分段線性獨(dú)立分布n四種概率分布n柱狀圖n離散的n分段線性n復(fù)合源信息卡nSIntypeI1 I2 Iknn分布號(hào)ntype(類(lèi)型）nH: 柱狀直方分布nL: 離散的源變量值nA: 線性分段值nS: 復(fù)合值nIk依賴(lài)于類(lèi)型的源變量值或分布號(hào)nH: 直方圖分布的箱邊界nL: 離散的源變量值nA: 定義概率分布的點(diǎn)nS: 次分布號(hào)源概率卡nSPnoptionP1 P2 Pknn分布號(hào)noption(說(shuō)明）n空格: H或者L分布的箱概率，A分布的概率密度nD: H或者L分布的箱概率nC: H或者L分布的積累箱概

22、率nV: 對(duì)于柵元分布n概率與柵元體積成正比（乘以Pi）nPi源變量概率直方分布nIk定義直方分布箱的邊界n必須單調(diào)上升nPk定義下列任何一種n每一個(gè)箱的相對(duì)概率n每一個(gè)箱的積累概率nP1必須是0nMCNP將會(huì)對(duì)Pi重新歸一直方分布示例n在(-1,1)之間的均勻分布SI1 -1 1SP1 01n五個(gè)箱且中間的概率最大SI1 H -1 -.6 -.2 .2 .6 1SP1 D 012321n與上相同,使用累計(jì)概率SI1 -1 -.6 -.2 .2 .6 1SP1 C013689離散分布nIk定義離散變量值n不需要單調(diào)上升nPk定義下列任何一種n每一個(gè)值的相對(duì)概率n每一個(gè)值的積累概率nMCNP將會(huì)

23、對(duì)Pk重新歸一離散分布示例n相同的邊界分布SI1 L -1 1SP1 -1 1n五個(gè)值且中間的概率最大SI1 L -1 -.5 0.5.1SP1 D1 2 3 2 1n與上相同,使用累計(jì)概率SI1 L -1 -.5 0.5.1SP1 D1 3 6 8 9n向量分布SI1 L0 0 01 1 12 3 49 8 7SP1 1 2 1 3分段線性分布nIk定義分段線性分布的點(diǎn)n必須單調(diào)上升n分段線性區(qū)間的最小值和最大值nPk定義分段線性區(qū)間每一點(diǎn)上的值nPk和Pk通常為0，但不做要求nMCNP將會(huì)對(duì)Pk重新歸一n在點(diǎn)之間的概率分布是線性?xún)?nèi)插的線性?xún)?nèi)插分布示例n從0到1的傾斜分布SI1 A 0 1S

24、P1 0 1n鋸齒分布SI1 A -1 -.5 0 .5 1SP1 0 1 0 1 0n余弦平方分布近似SI1 A0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1SI1 1 0.90 0.65 0.35 0.10 0 0.10 0.35 0.65 0.90 1定義源的位置n體積源可以在三種坐標(biāo)系統(tǒng)內(nèi)定義n直角坐標(biāo),柱坐標(biāo),球坐標(biāo)n表面源可以定義為n簡(jiǎn)化的體積源,具有曲面特征n位置的取樣以源的類(lèi)型為基礎(chǔ)n源的位置約束在特殊的柵元或曲面內(nèi)直角坐標(biāo)下的體積源n提供X,Y,Z的分布n例子：n在222的箱內(nèi)的均勻源SDEF X=D1 Y=D2 Z=D3SI1-11SP1 0

25、 1SI2-11SP2 0 1SI3-11SP3 0 1直角坐標(biāo)下的體積源n在222的箱內(nèi)特殊的非自然分布SDEF X=D1 Y=D2 Z=D3SI1-101SP1 0 1 2SI2A-11SP2 01SI3A-101SP3 010柱坐標(biāo)下的體積源nSDEF變量nPOS，AXS，RAD，EXTn圓柱的軸沿著AXS方向通過(guò)POSn在一個(gè)以RAD的抽樣值為半徑的圓上均勻抽取粒子的位置n圓位于在離開(kāi)POS的一定距離且垂直于AXS的一個(gè)平面上，這個(gè)距離是EXT的抽樣值柱坐標(biāo)下的體積源n柱坐標(biāo)下源的缺省取樣nEXT 均勻取樣nRAD 冪函數(shù)取樣正比于r1n一般舉例n均勻圓柱表面源n在SIn卡上填寫(xiě)兩端離

26、POS的距離EXTn在SIn卡上填RAD的內(nèi)半徑和外半徑n圓柱體積源n表面的內(nèi)半徑為0柱坐標(biāo)下的體積源舉例n2cm長(zhǎng)的燃料芯塊均勻源，中心在原點(diǎn)，軸為X軸，直徑1cmSDEF POS=0 0 0 AXS=1 0 0 EXT=D1 RAD=D2SI1-11SP1 0 1SI2 00.5SP2 -21 1n上列燃料芯塊外表面的均勻源SDEF POS=0 0 0 AXS=1 0 0 EXT=D1 RAD=D2SI1-11SP1 0 1SI2 0.40.5SP2 -21 1柱坐標(biāo)下的體積源舉例n仍為上列燃料芯塊，中心(2,3,4)，沿(1,1,1)方向SDEF POS=2 3 4 AXS=1 1 1

27、EXT=D1 RAD=D2SI1-11SP1 0 1SI2 00.5SP2 -21 1球坐標(biāo)下的體積源nSDEF變量nPOS，RADn不允許指定AXS和X,Y,Zn球心在POSnRAD為從球心到粒子均勻取樣位置的距離nRAD默認(rèn)值為正比于r2的冪函數(shù)球坐標(biāo)下的體積源舉例n球心在原點(diǎn)，直徑1cm的均勻源SDEF POS=0 0 0 RAD=D2SI2 00.5SP2 -21 2n上面均勻源的外殼SDEF POS=0 0 0 RAD=D2SI2 0.40.5SP2 -21 2n在三個(gè)不同位置的類(lèi)似均勻源SDEF POS=D1 RAD=D2SI1 L0 0 02 3 4-5 4 3SP1 1 1 1

28、SI2 0.40.5SP2 -21 2曲面源n曲面源可以看作簡(jiǎn)化的體積源n曲面說(shuō)明不是必需的，但是n注意：不要在幾何曲面上創(chuàng)建曲面源而不進(jìn)行說(shuō)明，所以要n指明幾何曲面SUR=n，或者n把曲面源定義在未說(shuō)明的幾何結(jié)構(gòu)之外n對(duì)曲面說(shuō)明后，MCNP將不會(huì)再次檢查是否從該曲面取樣平面曲面源n直角坐標(biāo)下的源分布nX，Y和Z，必須定義平面曲面SUR上的點(diǎn)n應(yīng)該定義源的方向VECn圓形的源分布nPOS定義分布的中心nRAD定義分布半徑nSUR或AXS定義平面平面曲面源舉例n22的箱均勻表面源SDEF X=D1 Y=D2 Z=1 SUR=3.6SI1-11SP1 0 1SI2-11SP2 0 1nR=2的圓柱

29、側(cè)面表面源SDEF SUR=4.3 POS=0 0 0 RAD=D1 SI11$ 默認(rèn)取樣正比于r球表面源n沒(méi)有表面說(shuō)明時(shí)nRAD為定值的簡(jiǎn)化的體積源n不允許指定AXS和EXTn區(qū)域內(nèi)均勻分布n有表面說(shuō)明時(shí)nRAD不需要nAXS和EXT說(shuō)明非均勻分布特點(diǎn)柱表面源n只看作簡(jiǎn)化的體積源n指明RAD的定值n不能用SUR=n否則出現(xiàn)致命錯(cuò)誤n例子：n長(zhǎng)2cm，半徑2cm，中心在原點(diǎn)沿X軸的圓柱測(cè)表面均勻源SDEF POS=0 0 0 AXS=1 0 0 EXT=D1 RAD=2.01SI1-11SP1 0 1點(diǎn)源n兩個(gè)方法說(shuō)明點(diǎn)源的位置nX,Y,ZnPOSn對(duì)于多個(gè)點(diǎn)源,X,Y,Z要求獨(dú)立分布n此時(shí)推

30、薦使用POSnPOS分布：SI1 L0 0 01 1 12 3 49 8 7SP1 1 2 1 3內(nèi)部函數(shù)(Build-in PDFs)n對(duì)于一些源變量，可以使用內(nèi)部函數(shù)nSPn f a bnn 分布號(hào)nf內(nèi)部函數(shù)識(shí)別符na b內(nèi)部函數(shù)參數(shù)n只能用于特定的變量（見(jiàn)p3-60,表3.4）n參見(jiàn)特殊的缺省值和快捷方式（p3-61）源的方向n DIR：cos，是VEC和粒子飛行 方向的夾角余弦nVEC：DIR的參考矢量源的方向默認(rèn)值n對(duì)于體積源（SUR=0)nVEC ：缺省時(shí)產(chǎn)生各向同性分布nDIR：由各向同性均勻抽取n對(duì)于曲面源（SUR不等于0)nVEC ：垂直于曲面，曲面用NRM表征nDIR：在

31、0到1之間由余弦分布抽取的線性樣本源方向舉例n沿X軸的激光束SDEFVEC=1 0 0 DIR=1 n沿X軸方向上有向前的最大輻射的靶SDEFVEC=1 0 0 DIR=D1SI101SP1 -215復(fù)合的獨(dú)立分布nIk定義各級(jí)各個(gè)分布的離散變量值nIk0：使用變量的默認(rèn)值nPk定義各級(jí)各個(gè)分布的概率nMCNP將會(huì)對(duì)Pk重新歸一n各個(gè)分布可以是獨(dú)立的分布，也可以又是另外的復(fù)合分布，最高可以有20級(jí)。復(fù)合分布舉例n輻射能譜通常既有離散能譜又有連續(xù)能譜。復(fù)合分布舉例n輻射能譜通常既有離散能譜又有連續(xù)能譜。SI1SD2D3SP10.60.4SI2A00.511.92.0SP2111.21.50SI2

32、L0.60.891.32SP22110相關(guān)源分布n一個(gè)源變量可能與另一個(gè)源變量值的選取有關(guān)nDSnoption J1 J2 Jknn分布號(hào)noption(說(shuō)明）n省略或H: 直方分布nL: 離散值nS: 復(fù)合值nJk依耐于類(lèi)型的源變量值或分布號(hào)nH: 直方圖分布的箱邊界nL: 離散的源變量值nS: 次分布號(hào)相關(guān)源分布nDSn T I1 J1 I2 J2 Ik Jknn分布號(hào)nIk獨(dú)立變量的離散值nJk相關(guān)變量的離散值n如果找不到匹配的獨(dú)立變量的抽樣值，那么相關(guān)變量取它的缺省值n當(dāng)相關(guān)變量取缺省值時(shí)，通常使用T選擇相關(guān)源分布nDSn Q V1 S1 V2 S2 Vk Sknn分布號(hào)nVk獨(dú)立變量

33、的箱邊界n單調(diào)上升nVk的最大值必須不小于獨(dú)立變量可能的最大值nSk相關(guān)變量的分布號(hào)n如果獨(dú)立變量的樣本抽取值在Vk-1和Vk之間時(shí),由Sk分布抽取相關(guān)變量的值n如果Sk=0,使用變量的缺省值n獨(dú)立變量使用內(nèi)部函數(shù)時(shí)必須采取Q選擇相關(guān)源分布舉例n能量不同的兩個(gè)點(diǎn)源SDEFPOS=D1ERG=FPOS=D2SI1L 5 0 0-5 0 0SP1 1 1DS2L 5 10n基于能量的不同角分布SDEFERG=D1DIR=FERG=D2SI1 H026SP1 0 0.2 1DS2S 3 4柵元舍棄n怎樣從不規(guī)則形狀的柵元體中均勻抽樣？n直角坐標(biāo)下體積源：CEL=3Cookie-cutter舍棄n怎樣

34、從不規(guī)則形狀的表面區(qū)域均勻抽樣？n定義一個(gè)特殊的柵元以適當(dāng)?shù)男螤钆c曲面相交nCCC=n關(guān)于源的舍棄n在取樣以后(使用X,Y,Z或POS,RAD,EXT確定位置）n如果(x,y,z)在柵元之外，舍棄并重新取樣n源變量EFF給定了取樣效率的下限n如果粒子接受率小于EFF,問(wèn)題放棄n默認(rèn)值：0.01nCookie-cutter柵元不能取為問(wèn)題中描述的真的幾何形狀源的檢查nPRINTX1X2X3數(shù)據(jù)卡n如果PRINT卡（輸出打印卡）不給出n則產(chǎn)生最小的基本輸出n如果Xi不給出n則產(chǎn)生最大輸出n如果Xi大于0n則產(chǎn)生最小基本輸出和Xi列表n如果Xi小于0n則產(chǎn)生最大輸出和Xi列表n見(jiàn)說(shuō)明書(shū)3-135的表

35、格n表110給出了前50種源粒子MCNP源的總結(jié)nMCNP給出了大量的靈活的定義源的方法n不同方式的獨(dú)立分布n復(fù)合分布n相關(guān)分布n內(nèi)部函數(shù)n柵元舍棄和Cookie-cutter舍棄n源的定義很復(fù)雜，必須經(jīng)常檢查MCNP教程n簡(jiǎn)介簡(jiǎn)介n輸入文件基礎(chǔ)輸入文件基礎(chǔ)n幾何描述幾何描述n源的描述源的描述n計(jì)數(shù)描述計(jì)數(shù)描述n材料，物理成分和數(shù)據(jù)材料，物理成分和數(shù)據(jù)n高級(jí)計(jì)數(shù)描述高級(jí)計(jì)數(shù)描述n高級(jí)幾何描述高級(jí)幾何描述n各種簡(jiǎn)化各種簡(jiǎn)化n臨界問(wèn)題臨界問(wèn)題MCNP計(jì)數(shù)n基本計(jì)數(shù)n計(jì)數(shù)箱n計(jì)數(shù)分段n計(jì)數(shù)乘子n計(jì)數(shù)箱乘子/歸一化n計(jì)數(shù)標(biāo)記n計(jì)數(shù)輸出優(yōu)化學(xué)習(xí)目標(biāo):MCNP計(jì)數(shù)n知道怎樣選擇正確的計(jì)數(shù)類(lèi)型n能定義計(jì)數(shù)箱

36、n能在幾何上對(duì)計(jì)數(shù)分段n知道使用計(jì)數(shù)乘子和找到乘子因子n知道使用計(jì)數(shù)標(biāo)記n能夠優(yōu)化計(jì)數(shù)輸出n能夠在運(yùn)行時(shí)輸出計(jì)數(shù)Fn計(jì)數(shù)卡n計(jì)數(shù)標(biāo)簽給出了物理量的估值nFn: a j1 (j2 j3) jN Tnn計(jì)數(shù)號(hào),其最后一位定義類(lèi)型na計(jì)數(shù)粒子類(lèi)型n,p,enjN 計(jì)數(shù)區(qū)域nT可選：所有計(jì)數(shù)區(qū)域的總數(shù)或平均數(shù)n括號(hào)內(nèi)的項(xiàng)給出這些區(qū)域的總數(shù)或平均計(jì)數(shù)計(jì)數(shù)類(lèi)型n表面計(jì)數(shù)n1穿過(guò)曲面的積分流量n2穿過(guò)曲面的平均通量n柵元計(jì)數(shù)n4一個(gè)柵元的平均通量n6一個(gè)柵元的平均沉積能量n7一個(gè)柵元的平均裂變沉積能量n特殊計(jì)數(shù)n5一個(gè)點(diǎn)或環(huán)探測(cè)器的通量n8一個(gè)脈沖幅度探測(cè)器的通量計(jì)數(shù)單位計(jì)數(shù)類(lèi)型 Fn單位*Fn單位F1：曲

37、面流量粒子MeVF2：曲面通量粒子/cm2MeV/cm2F4：柵元通量粒子/cm2MeV/cm2F5：探測(cè)器通量粒子/cm2MeV/cm2F6：能量沉積MeV/gJerk/gF7：裂變能量沉積MeV/g Jerk/gF8：脈沖幅度pulseMeV計(jì)數(shù)和Macrobodiesn不能在整個(gè)macrobody曲面上計(jì)數(shù)n必須使用單獨(dú)的“面”或者面的集合n如果曲面5是macrobody的RCC(直圓柱體）nF1: n 5將是錯(cuò)誤的nF1: n (5.1 5.2 5.3）才是正確的n在某些情況下，消除的曲面（eliminated surface)不能用來(lái)計(jì)數(shù)（或者描述源）n如果曲面2.3因?yàn)榕c另一個(gè)曲面

38、相同而被消除，那么nSUR=2.3在SDEF卡中錯(cuò)誤nF1: n 2.3在計(jì)數(shù)中錯(cuò)誤定義計(jì)數(shù)箱n在缺省情況下，計(jì)數(shù)是對(duì)所有能量、時(shí)間、角度、位置的平均積分n用戶(hù)可以定義下列每個(gè)量的計(jì)數(shù)箱n能量箱能譜n余弦/角度箱角度分布n時(shí)間箱瞬時(shí)特性n分段箱空間分布能量、時(shí)間和余弦箱En E1 E2 EkCn C1 C2 CkTn T1 T2 Tknn計(jì)數(shù)號(hào)nn0時(shí)對(duì)于所有計(jì)數(shù)給出箱的缺省值n余弦箱只在類(lèi)型1時(shí)有效nEi 第i個(gè)能量箱的上限nCi第i個(gè)余弦箱的上限（ C1-1） nTi第i個(gè)時(shí)間箱的上限nMCNP將會(huì)自動(dòng)給出時(shí)間箱和能量箱計(jì)數(shù)總和n在余弦卡末尾填寫(xiě)T也可得到計(jì)數(shù)總和計(jì)數(shù)分段n通過(guò)幾何劃分對(duì)計(jì)

39、數(shù)進(jìn)行空間分段nFSnS1 S2 Sk Tnn計(jì)數(shù)號(hào)nSi標(biāo)記的幾何曲面nT可選的總的計(jì)數(shù)箱n應(yīng)用于計(jì)數(shù)n的每一個(gè)柵元或曲面n可以要求SDn卡定義分段體積或面積理解計(jì)數(shù)分段nFS（分段計(jì)數(shù)）卡上的k個(gè)曲面將會(huì)創(chuàng)建k1個(gè)分段體積或曲面1.相對(duì)于曲面S1的坐向與數(shù)符相同的部分2.相對(duì)于曲面S2的坐向與數(shù)符相同的部分但在先前分段已記錄過(guò)的部分除外k.相對(duì)于曲面Sk的坐向與數(shù)符相同的部分，但在先前分段已記錄過(guò)的部分除外k+1所有其余部分k+2 如果在FSn卡上有字符T，給出整個(gè)柵元或曲面的計(jì)數(shù)計(jì)數(shù)分段的面積/體積SDn (D11 D12 D1m) (D21 D22 D2m) (Dk1 Dk2 Dkm)

40、 nn計(jì)數(shù)號(hào)nk計(jì)數(shù)卡上的柵元或曲面數(shù)目nm計(jì)數(shù)卡上的分段數(shù)目nDij第I個(gè)柵元或曲面上的第j個(gè)分段的面積，體積或質(zhì)量n在MCNP不能自動(dòng)計(jì)算體積或面積時(shí)必須使用此卡計(jì)數(shù)乘子nFMn(bin set1) (bin set2) Tnn計(jì)數(shù)號(hào)n(bin set i):箱組i,(乘子組1)(乘子組2)衰減器組)nT：可選項(xiàng)，給出所有箱的總計(jì)數(shù)n用來(lái)計(jì)算以下列形式的通量函數(shù)表示的量C(E)R(E)dEn在這里R(E)是一個(gè)算子，表示加或乘上一個(gè)響應(yīng)函數(shù)（例如截面）計(jì)數(shù)乘子組n乘子組nC m(reactor list 1) (reactor list 2) nC增殖常數(shù)n對(duì)于計(jì)數(shù)類(lèi)型4且C-1，則使用C

41、=anm在材料卡上標(biāo)識(shí)的材料號(hào)nreactor list i：反應(yīng)表i，由空格（表示乘）和或冒號(hào)（表示加）分開(kāi)的反應(yīng)號(hào)n乘在加之前運(yùn)算，且反應(yīng)表內(nèi)不允許括號(hào)n(R1 R2:R3) = (R1*R2)+R3n(R1 R2:R1 R3) = (R1*R2)+(R1*R3) = R1*(R2+R3)n計(jì)數(shù)乘子是C.R一些計(jì)數(shù)乘子反應(yīng)號(hào)中子光子-1 總的截面-1 非相干散射截面-2 吸收截面-2 相干散射截面-3 彈性散射截面-3 光電截面-4 平均加熱數(shù)-4 對(duì)生成截面-5 產(chǎn)生截面-5 總截面-6 總裂變截面-6 光子加熱數(shù)-7 裂變-8 裂變Q其它反應(yīng)號(hào)見(jiàn)附錄G16 (n,2n)的ENDF反應(yīng)號(hào)

42、17 (n,3n)18 (n,f)計(jì)數(shù)衰減器組n衰減器組nC -1 m1 px1 m2 px2 nC增殖常數(shù)nmi材料號(hào) npxi密度乘以衰減材料厚度 n正值表示原子密度n負(fù)值表示質(zhì)量密度n用計(jì)數(shù)乘子Ce-.px模擬粒子在衰減器中行為n只有衰減材料較薄，不考慮散射效應(yīng)時(shí)這樣考慮特殊的計(jì)數(shù)乘子n特殊乘子組nc -knc增殖常數(shù)nk特殊計(jì)數(shù)乘子n11/權(quán)重n給出徑跡數(shù)n21/速度n給出中子總數(shù)乘子箱EmnME1 ME2 MEkCMnMC1 MC2 MCkTMnMT1 MT2 MTknn 計(jì)數(shù)號(hào)nn0時(shí)對(duì)于所有計(jì)數(shù)給出箱的缺省值n余弦箱只在類(lèi)型1時(shí)有效nMEi：適用于能量箱i的乘子nMCi：適用于余

43、弦箱i的乘子(c1-1)nMTi：適用于時(shí)間箱i的乘子箱的歸一化n乘子箱的一個(gè)主要用途是通過(guò)箱的寬度進(jìn)行歸一化n標(biāo)準(zhǔn)箱對(duì)于每一個(gè)箱提供了全部的直方分布n通過(guò)箱的寬度進(jìn)行劃分可以產(chǎn)生下列結(jié)果：n每單位能量能量箱n每單位立體角余弦箱n每單位時(shí)間時(shí)間箱計(jì)數(shù)標(biāo)記nCFnC1 C2 CknSFnS1 S2 Sknn計(jì)數(shù)號(hào)n當(dāng)n=0時(shí)對(duì)于所有計(jì)數(shù)使用默認(rèn)標(biāo)記nCi柵元標(biāo)記編號(hào)nSi曲面標(biāo)記編號(hào)n當(dāng)任何一個(gè)粒子離開(kāi)柵元Ci或者穿過(guò)曲面Si時(shí)，把粒子軌跡打上標(biāo)記n可以標(biāo)記的光子是由帶標(biāo)記的中子產(chǎn)生的光子修改計(jì)數(shù)輸出（打印層次卡）nFQna1 a2 aknn計(jì)數(shù)號(hào)n計(jì)數(shù)箱類(lèi)型編碼（按照默認(rèn)順序)nF柵元，曲面或

44、探測(cè)器nD直接或者帶標(biāo)記nU用戶(hù)nS分片斷nM乘子nC余弦nE能量nT時(shí)間n最后兩項(xiàng)所對(duì)應(yīng)的計(jì)數(shù)箱將做成一個(gè)表，分別對(duì)應(yīng)表的豎向和橫向。修改計(jì)數(shù)輸出n如果指定了計(jì)數(shù)箱類(lèi)型編碼的子集，子集將會(huì)被放置到卡的末尾，而在它前面按照缺省次序安排未指明的字母nFCn計(jì)數(shù)注釋卡n可在輸出文件中打印注釋信息n在計(jì)數(shù)被修改以后特別有用交互式計(jì)數(shù)顯示nMCNP可以暫停以檢查計(jì)數(shù)的中間狀態(tài)（見(jiàn)附錄B）nTallynn將計(jì)數(shù)n設(shè)為當(dāng)前計(jì)數(shù)以測(cè)繪nFreeqn設(shè)置變量qF D U S M C E T為X軸nFixedqnn測(cè)繪變量q F D U S M C E T的計(jì)數(shù)箱nnLINLIN或者LOGLIN：切換Y軸的對(duì)數(shù)

45、坐標(biāo)關(guān)于計(jì)數(shù)的更多主題n環(huán)形探測(cè)器n探測(cè)器響應(yīng)函數(shù)n微擾n計(jì)數(shù)漲落繪制n計(jì)數(shù)特殊處理MCNP計(jì)數(shù)總結(jié)nMCNP提供了一些物理量的計(jì)數(shù)n能量箱，時(shí)間箱和角度箱便于進(jìn)行更細(xì)致的分析n計(jì)數(shù)分段允許存在簡(jiǎn)單的空間相關(guān)性n計(jì)數(shù)乘子可以計(jì)算更多的物理量n計(jì)數(shù)可以在計(jì)算過(guò)程中檢測(cè)，而且其結(jié)果可以以不同順序輸出MCNP教程n簡(jiǎn)介簡(jiǎn)介n輸入文件基礎(chǔ)輸入文件基礎(chǔ)n幾何描述幾何描述n源的描述源的描述n計(jì)數(shù)描述計(jì)數(shù)描述n材料，物理成分和數(shù)據(jù)材料，物理成分和數(shù)據(jù)n高級(jí)計(jì)數(shù)描述高級(jí)計(jì)數(shù)描述n高級(jí)幾何描述高級(jí)幾何描述n各種簡(jiǎn)化各種簡(jiǎn)化n臨界問(wèn)題臨界問(wèn)題MCNP中的各種約簡(jiǎn)n各種約簡(jiǎn)的語(yǔ)法n截?cái)鄋數(shù)目控制n調(diào)整取樣n部分確定

46、性n常見(jiàn)問(wèn)題能量和時(shí)間截?cái)鄋CUT：n T E WC1 WC2 SWTMnn粒子類(lèi)型：N為中子，P為光子，E為電子nT最高截?cái)鄷r(shí)間，以剎為單位nE最低截?cái)嗄芰浚訫eV為單位n當(dāng)粒子的壽命超過(guò)T時(shí)就被殺死n當(dāng)粒子的能量小于E時(shí)就被殺死n缺省值n中子：T為無(wú)限大值，E0.0MeVn光子：T為中子壽命，E0.001MeVn電子：T為中子壽命，E0.001MeV權(quán)重截?cái)鄋CUT：n T E WC1 WC2 SWTMnn粒子類(lèi)型：N為中子，P為光子，E為電子nWC1 WC2權(quán)重截?cái)鄋SWTM 源粒子最小權(quán)重n當(dāng)粒子權(quán)重WGT1時(shí)實(shí)行分離n0Ni1時(shí)按概率實(shí)行輪盤(pán)賭n最多五組n當(dāng)粒子的能量上升到Ei 之

47、上時(shí)，執(zhí)行與上相反的“過(guò)程”權(quán)重窗nWWE: n E1 E2 Ej(j100)nWWNi: n Wi1 Wi2 Wij(j0 為SWITCHN除以粒子所在柵元的重要性n通常使用SWITCHN0nMTIMEWWE卡上的權(quán)窗類(lèi)型 0n0能量n1時(shí)間權(quán)窗產(chǎn)生器nMCNP自動(dòng)產(chǎn)生優(yōu)化后的權(quán)窗n權(quán)窗的產(chǎn)生是為了提供相對(duì)于特定柵元（源柵元）的特定計(jì)數(shù)n產(chǎn)生的權(quán)窗可以基于柵元，也可以基于幾何形狀重疊的網(wǎng)格n結(jié)果寫(xiě)入文件中以給WWP卡作參考SWITCHN=-1權(quán)窗產(chǎn)生器nWWG It Ic Wg J J J J IEnIt：對(duì)于優(yōu)化問(wèn)題的計(jì)數(shù)號(hào)nIc：權(quán)窗的參照柵元n0 問(wèn)題柵元號(hào)n0網(wǎng)格產(chǎn)生器nWg：參照柵

48、元Ic的權(quán)窗下限值n0取值為平均源權(quán)的一半nJ：不使用nIE：切換能量（0）和時(shí)間（1）網(wǎng)格產(chǎn)生器nMESHmesh variable=specificationnGEOMnxyz，rec直角坐標(biāo)nrzt，cyl柱坐標(biāo)nREF 參考點(diǎn)的x，y，z坐標(biāo)nORIGINMCNP幾何中多重網(wǎng)格的x，y，z坐標(biāo)nAXS 園柱網(wǎng)格的軸方向向量nVEC 定義的向量，沿AXS方向， 平面角0n網(wǎng)格的邊界應(yīng)該在問(wèn)題幾何之外n詳見(jiàn)3-44網(wǎng)格產(chǎn)生器nMESH mesh variable=specificationnIMESH, JMESH, KMESH：在(x,y,z)或者(r,)方向上的網(wǎng)格點(diǎn)的近似位置nIIN

49、TS, JINTS, KINTS：在(x,y,z)或者(r,)方向上的網(wǎng)格點(diǎn)之間的精確距離n網(wǎng)格的邊界應(yīng)該在問(wèn)題幾何之外n詳見(jiàn)3-44指數(shù)變換nEXT: n A1 A2 Ai AInVECTVa Xa Ya Za Vb Xb Yb Zbnn粒子類(lèi)型nAi第I個(gè)柵元的數(shù)據(jù)項(xiàng)，格式為QVCnQ0：不進(jìn)行變換nQp：0p1，固定的拉伸參數(shù)nQS：pa/t， a是俘獲截面nI問(wèn)題的柵元數(shù)na，b 標(biāo)記矢量Va,Vb的序號(hào)nXa Ya Za定義矢量Va的三個(gè)一組的坐標(biāo)強(qiáng)迫碰撞nFCL: n x1 x2 xi xInn粒子類(lèi)型nxi第i個(gè)柵元的強(qiáng)迫碰撞控制參數(shù)； | xi|1n如果| xi|1，對(duì)碰撞粒子

50、用幸存概率| xi|做輪盤(pán)賭 nI問(wèn)題的柵元總數(shù)n當(dāng)粒子進(jìn)入強(qiáng)迫碰撞的柵元時(shí)，在表面上將不做權(quán)窗處理nxi0，強(qiáng)迫碰撞重復(fù)進(jìn)行直至粒子被殺死源偏倚nSBn option B1 B2 BknSBn f a bnn分布號(hào)noption與SPn卡相同n空格: H或者L分布的箱概率，A分布的概率密度nD: H或者L分布的箱概率nC: H或者L分布的積累箱概率nV: 對(duì)于柵元分布n概率與柵元體積成正比（乘以Pi）nBi源變量概率nf內(nèi)部函數(shù)識(shí)別符na，b 內(nèi)部函數(shù)輸入?yún)?shù)基本的源偏倚n能量偏倚nSDEF ERGD1SI1 0 1 2 3 4 5SP1012121SB101111 10n方向偏倚nSDEF VEC1 0 0 DIRD2SI1 1 0 1 SP1011SB1015偏倚內(nèi)部函數(shù)n沒(méi)有偏倚，通常不給出SI卡nSI卡可以用來(lái)提供偏差nMCNP近似函數(shù)列出了高達(dá)300個(gè)等概率箱nSDEFERGD1SP15aSI10.0050.120SB1 C 0 0.5 1點(diǎn)和環(huán)形探測(cè)器nFn5：p x y z R0nFn5a：p a0 r R0nn5以5結(jié)束的計(jì)數(shù)號(hào)np粒子類(lèi)型nR0探測(cè)點(diǎn)處鄰域球半徑，正號(hào)表示以cm為單位，負(fù)號(hào)表示以平均自由程為單位n點(diǎn)探測(cè)器nX Y Z點(diǎn)探測(cè)器的位置n環(huán)探測(cè)器na定