MCNP-5A程序使用說明書

MCNP-5A2010.7 1.4卡片書寫格式7...1.5粒子標(biāo)識符9... 4.2.1VOL柵元體積卡 23 4.3.2ESPLT能量分裂與輪盤賭卡 284.3.3TSPLT時間分裂或輪盤賭卡 304.3.4PWT光子權(quán)重卡 314.3.5EXT指數(shù)變換卡 31 4.3.7FCL強迫碰撞卡 33 4.3.9WWP權(quán)窗參數(shù)卡 34 4.3.12PDN探測器貢獻(xiàn)卡 364.3.13DXCDXTRAN貢獻(xiàn)卡 3.7. 4.4.1SDEF通用源卡 37 4.4.7SSW寫曲面源卡 48 4.4.11用戶提供的子程序SOURCE及SRCDX 53 4.5.3EN記數(shù)的能量卡6.3..4.5.4TN記數(shù)的時間卡6.3..4.5.5CN記數(shù)的余弦卡6.3..4.5.6FQN打印層次卡6.4..4.5.7FMN記數(shù)乘子卡6.4..4.5.8DEN及DFN劑量轉(zhuǎn)換卡6.7. 4.5.12CFN柵元標(biāo)記卡(記數(shù)類型1,2,4,6,7) 69 4.5.14FSN分段記數(shù)卡(記數(shù)類型1,2,4,6,7) 704.5.15SDN分段除數(shù)卡(記數(shù)類型1,2,4,6,7) 714.5.16FUNTALLYX輸入卡 714.5.17子程序TALLYX用戶提供的子程序 7.2 4.5.19DDN探測器診斷卡 744.5.20DXTDXTRAN卡 1FTN記數(shù)特殊處理卡 76 4.6.1MN材料卡8.0..4.6.2MPNN光核子核素選擇3DRXS離散反應(yīng)截面卡4TOTNU總裂變卡8.2..4.6.5NONU裂變關(guān)閉卡8.3..4.6.6AWTAB原子量卡8.3..4.6.7XSN截面文件卡8.4..4.6.8VOID取消材料卡8.4..4.6.9PIKMT光子產(chǎn)生偏倚卡10MGOPT多群特征卡8.5..4.7能量及熱處理方式的指定1PHYS能量物理截斷卡2TMP自由氣體模型熱溫度卡3THTME自由氣體模型熱溫度時間卡MTMS(Α,Β)材料卡問題截斷卡9.0.. 4.8.3NOTRN中性粒子對點探測器的直接貢獻(xiàn) NPS歷史截斷卡 92 4.9用戶數(shù)據(jù)數(shù)組 92 4.10.1PRDMP打印及轉(zhuǎn)儲周期卡 934.10.2LOST丟失粒子控制卡 944.10.3RAND隨機數(shù)發(fā)生器卡 944.10.4DBCN調(diào)試參量卡 954.10.5FILES文件生成卡 974.10.6PRINT輸出打印卡 974.10.7TALNP記數(shù)輸出控制卡 994.10.8MPLOT運行中圖示記數(shù) 994.10.9PTRAC粒子軌跡輸出 100 附錄 110 到全部卡片中的一小部分。全部輸入數(shù)據(jù)卡片類型的見2-**和附錄A。在這本手冊中“卡”NP。MCNP的所有功能都應(yīng)在弄懂的基礎(chǔ)上謹(jǐn)慎使用。尤其是對于探測器記數(shù)及降低方差技MCNP們建議讀者最好閱讀第二章的相應(yīng)部分。輸入文件(INP)可以有兩種形式：初始運行和接續(xù)運行。無論哪種形式，INP文件中都1.1信息塊僅僅是給出MCNP一個執(zhí)行信息的方法。其次，它也是一種避免重復(fù)地輸入常用字符的簡便大寫的、小寫的或者大小寫混合的，在信息塊與標(biāo)題卡之間用一個空行作為塊間的分隔符。**-**也討論的常規(guī)的執(zhí)行行信息一樣。任何更改文件名、程序模塊化的選擇、執(zhí)行行信息的鍵盤輸入都比信息塊沖突消息的優(yōu)先級高。在信息塊中INP=文件名的輸入是1.2初始運行信息塊卡選擇項(空行分隔)柵元描述卡(空行分隔)(空行分隔)數(shù)據(jù)卡(空行結(jié)束)推薦其它卡選擇項符的字母可以是大寫的、小寫的或者大小寫混合的。加資料，比如原文本信息或者問題的另一種描述，結(jié)束符就可以阻止讀取這些附加行。1.3接續(xù)運行接續(xù)運行文件是用于繼續(xù)計算一個早先被終止計算的問題，，例如一個程序先運行了2打印計算結(jié)果。對于接續(xù)運行，必須在信息塊或鍵入的命令行中給出選擇運行條件c或者cn，以此指明當(dāng)前是一個接續(xù)運行方式。接續(xù)運行必須在執(zhí)行行信息中給出Cm，表示從的結(jié)果。在PRDMP卡(page)的介紹中討論了如何選擇轉(zhuǎn)儲時間。下面我們將談到，用戶可從這些轉(zhuǎn)儲中的任何一次重新運行。據(jù)覆蓋掉原有的轉(zhuǎn)儲容，這樣就避免了RUNTPE文件大小的迅速增長，RUNTPE文件的增大也該文件中的continue卡是必須有的，且從第一列寫起，或者在選擇的信息塊和它的空信息塊卡(空行分隔)CNTIMUE卡...(空行結(jié)束)個子集。接續(xù)運行數(shù)據(jù)卡可以是：FQ，DD，NPS，CTME，IDUM，RDUM，PRDMP，還有一個便捷的特性是，如果上述的文件都沒有改變，而且計算環(huán)境允許執(zhí)行行信息，DD.11000我們將一個程序已經(jīng)運行了一分鐘，如果想要模擬更多的粒子，在執(zhí)行行中執(zhí)行c或cn。叫做RUNTPE的重啟文件。E粒子數(shù)NPS，該數(shù)應(yīng)包括以前已完成的粒子數(shù)在。在一次接續(xù)運行中的CTME表示還要運行一方法是十分便利的。1.4卡片書寫格式所有的輸入都限制于80列以，不區(qū)分大小寫。大部分輸入卡片按行填寫，然而數(shù)據(jù)卡容?？瞻卓ㄗ鞣纸绶蚪K止符使用，數(shù)據(jù)輸入文件由一個或多個空白卡分成幾部分。n將作為源概率分布n的表頭文字印出。一個獨立的數(shù)據(jù)卡必須在一卡片上寫完。1.行輸入格式號和粒子標(biāo)識符(與后面的數(shù)據(jù)輸入之間用至少一個空格隔開)，而且可以寫在1-5列的任至少一個空格結(jié)束的一行表示下一行是繼續(xù)卡。繼續(xù)卡上的數(shù)據(jù)也是寫在1-80列?？瞻卓ㄗ鳛檩斎胛募膸讉€部分之間的分隔符。一個獨立的數(shù)據(jù)項必須在一卡片上寫完。對于一個MCNP進(jìn)行適當(dāng)?shù)霓D(zhuǎn)換。實際上，任何非整數(shù)數(shù)據(jù)類型，使用四種書寫功能：I23456。2m4m2m2m便相當(dāng)于112481664128256。省值。例如，下面兩卡片的作用是一樣的。DDJ1000BCNJ13MI4=133.5413M3M=13912R2I2.5=1l1l.52.02.5IR2M=112IRR=11112I43M=12341212I42I10=123468103J4R不合法14I3M不合法14IJ不合法2.列輸入格式列輸入格式對于柵元參數(shù)和源分布是特別有用的。按行排列的柵元重要性或體積可讀性的數(shù)據(jù)逐個放在每一行上。源的選擇(除了缺省的情況)在下一行輸入清楚。不再需要&繼續(xù)符，如果出現(xiàn)也會被忽略。“?！碧柺欠旁诳ㄆ诌@一行的l－5列的某個列上。卡片名字必須全部都是柵元參數(shù)卡片名字，或全部都是曲面參數(shù)卡片名字，或全部都是其它卡片名字。在一個輸入文件，如果片名字，或全部都是曲面參數(shù)卡片名字，或全部都是其它卡片名字。在一個輸入文件，如果一個卡片名字在一個#卡上出現(xiàn)，則在該輸入文件，這個卡片決不能再按行格式輸入。如果在#卡上輸入的數(shù)據(jù)比卡片名字長，第一個數(shù)據(jù)輸入柵元或者曲面的編號。如果輸入了任一據(jù)塊等。如果使用了很多柵元數(shù)據(jù)選擇項，則需要額外的列數(shù)據(jù)塊。RMIJ中使用恰當(dāng)。＃SS…S12mK1DD…D11121mK2DD…D21222m···KD…Dij…Dii1im···KDD…Dnn1n2nm3.除了隨K變動以外，每一列j(S-D)表示一個行輸入卡，這里的1小于n。ij1jj5.D-D必須是S卡的有效項。但D~D當(dāng)中的某些J項或者J項后面跟有一些空1jnjj1+1,jnjiji,j-1i,j-17.S不能再在該輸入文件的其它地方出現(xiàn)。i8.K是選擇項，它們是柵元的名字(程序編號)，如果有任一K不是空格，則全部KjiiKjii否則，忽略K并且K不應(yīng)出現(xiàn)。iiiij1.5粒子標(biāo)識符性卡，用imp:p。為多種粒子指定同樣的數(shù)據(jù)，可以用一卡片代替多。Fn中子點探測器進(jìn)行能量記數(shù)。在熱記數(shù)中的例子中，可能出現(xiàn)雙粒子標(biāo)識符。p1.6缺省值MCNP的很多輸入?yún)⒘渴怯腥笔≈档?，它們列在第四章中。因此用戶并非每次都必須明確給出各個輸入?yún)⒘康闹?。?dāng)缺省值符合用戶要求時，便可不在輸入文件中指定。當(dāng)省略某輸入卡時，則該卡上的全部參量均使用缺省值。但是當(dāng)只想改變一卡片上的某一特定缺省參量時，那么在它前面的參量仍需指明，或者用nJ方式跳過前面那些使用缺省值的參量也可。例如在光子截止卡片使用CUT:P3J-.10可以簡便的跳過前三個使用缺省值改變第四個參數(shù)的數(shù)值。MCNP誤作了廣泛地檢查，當(dāng)存在違反了基本規(guī)定的錯誤時，將在終端和OUTP文件中給出致命性(FATAL)錯誤的信息，并且MCNP將不再做粒子輸運計也可能是假的，這要根據(jù)前面出現(xiàn)的致命性錯誤的情況而定。若在MCNP執(zhí)行命令行上指定MCNP還給出許多警告性的而不是致命的錯誤信息，程序?qū)⒗^續(xù)執(zhí)行。警告是試圖告知忽視，在進(jìn)行重要計算之前，應(yīng)當(dāng)搞清楚它們的意義，再作決斷。除了致命錯誤和警告信息以外；除了致命錯誤和警告信息以外；MCNP在遇到任何嚴(yán)重?fù)p壞之前(如被零除)立即出現(xiàn)“BADTROUBLE”信息，并中斷運行。1.8檢查幾何錯誤P這個幾何系統(tǒng)，從中發(fā)現(xiàn)錯誤，以便糾正。另一種檢查幾何錯誤的方法是建立和運行一個簡短的問題，對這個問題，用外源的粒子軌跡注滿這個系統(tǒng)。即用一個球?qū)⒄麄€幾何系統(tǒng)包起來；該球面為球面型源，進(jìn)行一次短時ID掉任何的FM卡。(2)增加一個柵元和一個大的球面，這個球面包圍整個幾何系統(tǒng)，并且這個系統(tǒng)的外好使非零重要性相等。如果這個系統(tǒng)是一維或二維無限，使用一個或一個以上平面代替一個引起某些粒子丟失。新運行這個歷史以便在OUTP文件上產(chǎn)生某些特殊的輸出。在重新運行期間打印事件記錄。子因。原點放在丟失粒子的坐標(biāo)點上，然后將繪圖平面的水平軸選為丟失粒子的方向余弦(U，V，W)。這樣錯誤的原因很可能就在圖上出現(xiàn)破折號，或者畫圖和所想要的圖之間有些第二章柵元描述卡或者jlikenbutlistm：柵元的材料號，它是材料卡(Mm)中相應(yīng)材料的序號。真空柵元，m=0。d：柵元材料的密度。填入正值時，表示是原子密度(單位為1024個原子/cm3)；填入負(fù)gcm。對于真空柵元，該項不填，直接列出下一項容。geom：柵元的幾何描述。它列出界定該柵元的所有曲面號(帶有數(shù)符，表示坐向)，及它必須作為界定該概元的一個曲面列出(同樣帶有坐向數(shù)符)。n：另一個柵元的序號。list：一組關(guān)鍵詞=數(shù)值，它說明定義的n柵元和m柵元之間的不同屬性。t的數(shù)字，認(rèn)為它是柵元號。這是對一個柵元幾何描述的簡寫。#3$與下行相同#(-12-4)的FILL的輸入可能很長又互相關(guān)聯(lián)，如果放在柵元卡將會簡化它們的概念。在柵元卡的第二和第三項。TMP卡和WWN卡有兩種輸入形式。一種是關(guān)鍵詞=數(shù)值，可以使用TMP1=數(shù)值TMP2=數(shù)IMPN4IMP:P=8EXT:N=-.4X注意以下幾點：2.在對一個柵元進(jìn)行幾何描述時不要使用不常用的復(fù)雜操作符去添加不必要的曲面。額外的曲面會讓程序運行起來變慢而且可能破壞體積和面積計算的必要條件。4.用外源的粒子軌跡注滿這個系統(tǒng)來檢查幾何錯誤。5.如果添加或移動?xùn)旁乃械臇旁獏?shù)卡。如果柵元參數(shù)卡使用列輸入格式則2.1柵元描述的簡寫格式3LIKE2BUTTRCL=1IMP:N=10這表明柵元3除了位置不同(TRCL=1)和不同的中子重要性之外其他方面都和柵元2第三章曲面描述卡3.1用方程定義曲面格式:jnalistj:曲面編號，l≤j≤99999。對于反射曲面和白邊界，只需在曲面描述卡的曲面序號a:曲面助記名(見表2-1)。list：按照規(guī)定次序給出的描述該曲面的數(shù)據(jù)項(1-10項)。xyz的計算值是正的，則稱該點對于這個曲面具有正的坐向。代替以防止這個曲面的坐向相反。DXTRAN2.周期性平面不允許曲面變換。但性的。4.周期性平面只能界定其他周期性平面或者上下界面來界定。5.周期性平面的一邊一定有一個零重要性柵元。6.所有的周期性平面都有一個旋轉(zhuǎn)向量法線，指向幾何體的上面和下面。7.不能使用下一次事件估計例如探測器和DXTRAN球。注：方程1.(yy)2+(zz)2t(xx)=0方程2.(xx)2+(zz)2t(yy)=0方程3.(xx)2+(yy)2t(zz)=0方程6.(xx)2/B2+((yy)2+(zz)2A)2/C21=0方程7.(yy)2/B2+((xx)2+(zz)2A)2/C21=0方程8.(zz)2/B2+((xx)2+(yy)2A)2/C21=0jGQ5.750-.8660-12-23.46439j7cx1*TR761-1.73203060 (r,s)柱坐標(biāo)系統(tǒng)圍繞S軸旋轉(zhuǎn)的橢圓(圖2-1b)。其原點位于原坐標(biāo)系的點(x,y,z)：s2+(ra)2=1bc2〈|r=(xx)2+(zz)2當(dāng)a<c時，出現(xiàn)退化情況，如果0<a<c，得到的是外曲面(圖2-1c)；-c<a<0時，得到的是曲面(圖2-1d)。橢圓環(huán)面的坐標(biāo)變換是局限于輔助坐標(biāo)系的每個軸都要平行于基本坐標(biāo)系的一個軸。3.2用幾何定義曲面坐標(biāo)點必須全都在同一個葉上。jYy1rly2r2這里(y,r=x2+y2)是點i的柱坐標(biāo)。iiii給出的坐標(biāo)對數(shù)不同，描述的曲面類型也不同：(1)給出一對坐標(biāo)，則定義的是平面，即PX、PY或PZ。(2)給出兩對坐標(biāo)，定義的是線性曲面，如PX，PY，PZ，CX，CY，CZ，KX，KY或KZ。當(dāng)用兩個點定義一個錐面時，只生成一個單葉錐面。用這個規(guī)定所確定的這些曲面(SQ除外)的坐向與由方程指定曲面得到的坐向是等同的。而對SQ來講，是這樣定義坐向的，即離對稱軸足夠遠(yuǎn)的點有正的坐向。這與方程定義jSQ-.0833331100068.52083-26.500這描述的是一個以(0，4，0)為中心，以1cm為半徑的球。例5：101-23$cell11Y-32212Y2333423Y213142義的是關(guān)于y軸對稱的環(huán)形柵元，該柵元的截面如圖2-2所示。第一個曲面通過點 (-3,2)及(2,1)，第二個曲面通過點(2,3)、(3,3)及(4,2)。由點(2,1)、(3,1)和(4,2)定義最后一個曲面．這些幾何點坐標(biāo)均指的是(y,r)。3SQ1-1.51000-0.62502.503.3由三個點定義一般平面PxyzxyzxyzAX＋BY＋CZ-D=0面方程系數(shù)遵循如下原則：(1)應(yīng)使坐標(biāo)原點關(guān)于該平面有負(fù)坐向；(2)當(dāng)該平面通過坐標(biāo)原點時(D=0)，則應(yīng)使得點(0,0,∞)對該平面有正坐向；(3)若上兩項都無法做到(即D=C=0)，則應(yīng)使點(0,∞,0)對該平面有正坐向；(4)若上三項都無法做到(即D=C=B=0)，則應(yīng)使點(∞,0,0)對該平面有正坐向；倘若第四項也失敗，那么三個點給在一條直線上，則將印出警告信息。3.4用宏體定義割、其他柵元的定義、SDEF源等等。但不能在SSR/SSW卡、賦值的正負(fù)性不變，并且要加上適當(dāng)?shù)淖⑨?。參看后面的例子。下面列出了可使用的幾何體以及它們的完整描述：BOX任意方向的長方體RPP長方體SPH球RCC直圓柱REC直橢圓柱ELL橢圓體ARB任意多面體BOXBOX:：任意方向為直角的盒子(所有角都是90°)BOXVxVyVzA1xA1yA1zA2xA2yA2zA3xA3yA3z其中VxVyVz=角的坐標(biāo)(x，y，z)A1xA1yA1z=第一個面的向量A2xA2yA2z=第二個面的向量A3xA3yA3z=第三個面的向量一個以原點為中心，棱長為2cm，各面平行于主坐標(biāo)軸的立方體。RPP：長方體RPPXminXmaxYminYmaxZminZmaxSPHVxVyVzR其中VxVyVz=角的坐標(biāo)(x，y，z)，R=半徑RCC：直圓柱RCCVxVyVzHxHyHzR其中VxVyVz=底面中心(x，y，z)HxHyHz=圓柱軸向量R=半徑RHPv1v2v3h1h2h3r1r2r3s1s2s3t1t2t3vvv的底面中心坐標(biāo)h2，h3=00hr1r2r3=從中心軸指向第一個面中心的向量s1s2s3=從中心軸指向第二個面的中心的向量t1t2t3=從中心軸指向第三個面的中心的向量REC：直橢圓柱RECVxVyVzHxHyHzV1xV1yV1zV2xV2yV2z其中VxVyVz=圓柱底面中心的坐標(biāo)(x，y，z)HxHyHz=圓柱軸向的高度向量V1xV1yV1z=橢圓長軸向量(H1向量的法向)V2xV2yV2z=橢圓短軸向量(H、V1向量都成直角)這是一個以y軸為對稱軸，高度為10cm的直橢圓柱。底面中心在(0，5，0)，長軸長TRC直截錐TRCVxVyVzHxHyHzR1R2其中VxVyVz=直截錐的底面中心坐標(biāo)(x，y，z)HxHyHz=截錐的軸向高度向量R1=截錐下底面半徑R2=截錐上底面半徑心(5，0，0)，半徑2cmELL：橢圓體ELLV1xV1yV1zV2xV2yV2zRmV1xV1yV1z=第一個焦點坐標(biāo)V2xV2yV2z=第二個焦點坐標(biāo)Rm=長軸長度V1xV1yV1z=橢圓中心坐標(biāo)V2xV2yV2z=長軸向量(模=長軸長度)Rm=短軸長度ELL000003-2xWED：楔狀體WEDVxVyVzV1xV1yV1zV2xV2yV2zV3xV3yV3zV1xV1yV1z=三角形底面第一個邊的向量V2xV2yV2z=三角形底面第二個邊的向量V3xV3yV3z=高度向量WEDARB任意多面體ARBaxayazbxbybzcxcycz...hxbyhzN1N2N3N4N5N6必須輸入ARB的八個頂點坐標(biāo)(每組三個數(shù)據(jù))，盡管一些可能用不上，甚至一組也用不上。這八組數(shù)據(jù)后面還有六項輸入N，遵循以下規(guī)定：每一項都是四個阿拉伯整數(shù)來定義ARB的一個面，這四個數(shù)字分別代表這個面上角的頂點。例如，輸入1278就定義了這個平面由第1、2、7、8角的坐標(biāo)來界定。但是我們知道三個點就可以定義一個平面，所以這個是由MCNP算出的。如果距離的絕對值大于1.e-6，系統(tǒng)就會給出錯誤信息，并且將這個距離和應(yīng)該在這個平面上的點的坐標(biāo)(x，y，z)打印在OUTP文件中。如果第四個數(shù)字是0，B例子：ARB-5-10-5-5-1055-10-55-1050120000000000123412501350245034500這是一個五個面的多面體，頂點坐標(biāo)分別為(-5,-10,-5),(-5,-10,6),(5,-10,-5),幾何繪圖器中用MBODY=OFF就可以看到表面編號的圖解。xSPH:當(dāng)作常規(guī)曲面，沒有多面體的面zREC:1.橢圓柱曲面TRC:1.錐體表面ELL:當(dāng)作常規(guī)曲面，沒有多面體的面WED:1.包含上下直角斜邊的傾斜平面V下面的輸入文件描述了五個柵元，圖解了一個由多種宏體、柵元/曲面組合的描述。曲面號用斜體寫在它們定義的平面旁邊。注意柵元或曲面的序號不一定從1開始也不用連。4.1e電子輸運問題類型(MODE)卡9)能量及熱中子的處理；11)用戶數(shù)據(jù)組；10)截斷條件；12)外圍卡。：第四章數(shù)據(jù)卡1)問題類型卡；3)坐標(biāo)變換；2)柵元和曲面參數(shù)；4)源的描述；5)多群特征描述；6)光子產(chǎn)生偏倚描述；7)記數(shù)指定；8)材料描述；4.2.1VOL柵元體積卡x=xino=xx…x12inoxx…x12i序給出致命性錯誤信息?？墒褂胣j功能跳過不想填入體積的柵元。柵元的體積按cm3單位的)的旋轉(zhuǎn)柵元，以及多面體，MCNP都將自動算出它們的體積。在計算體積的同時，也可4.2.2AREA曲面面積卡格式：areaxx…x…x12inx=曲面i的面積；i缺?。篗CNP會試圖計算所有曲面的面積。當(dāng)一個曲面的AREA卡沒有輸入數(shù)據(jù)時，就如果MCNP能夠計算出一個曲面的任何一側(cè)的所有柵元的體積，則該曲面的面積也將被計出，該面積又在記數(shù)中要被運用時，程序給出致命性錯誤信息。AREA卡將涉及幾何描述中的所有曲面，并非只涉及一個特定的記數(shù)界面。如果需要一4.2.3重復(fù)結(jié)構(gòu)卡1)UUNIVERSE卡e填充的任何柵元里的全部空間?？紤]被填充柵元和用來填充柵元的Universe之間關(guān)系的一種方式是被填充的柵元是一個往里面窺視第二級的窗口，就象墻上的一個窗戶提供戶外景象一樣。在第二級的柵元可以基本坐標(biāo)系下可有自己的原點，與上一級無關(guān)。然而，如果被填充柵元和用來填充柵元的這二者的坐標(biāo)系轉(zhuǎn)換關(guān)系。Universe。換句話說，一個Universe中的那些柵元不是正好地裝進(jìn)被填充的柵元。最好解決的方法是或者縮小被Universe尺寸，但是這種方法運行慢一些。101-23-45-6fill＝l207-89-10u=1fill=2lat＝130-11u=-24011u=250-l:2:-3:4:-5:6lpx0.0012px49.99993py-9.99994py9.99995pz-56pz57px08pxl09py010py1011s5504就會使問題運行的快一些。該減號指出忽略在較高級柵元的跟蹤。一個問題中的每個柵元或者是問題的實際領(lǐng)域的一部分或者是某個Universe的一部分，但是一個問題中的曲面限制較少?？捎脝蝹€曲面描述在多個Universe中的一些柵元。來填充某個級上的柵元的某個Universe的那些柵元，而不放在較高級，但在一個柵格的情況下除外。子也可以用宏體來定義：10－20fill＝120－30u＝1fill＝2lat＝130－11u＝－24011U＝2502020rpp050－1010－5530rpp010010－5511s55042)TRCL柵元變換卡TRCL卡使得只需一次描述界定在形狀和尺寸是等同的，只是在幾何位置上不同的幾個Unverse，那么，一個單獨的變換就能完全定義所有這些被填充柵元的部，因為該Universe中的柵元將繼承它們填充那些柵元的變換。換。填寫變換的規(guī)定同TR卡。如果使用*trcl符號，轉(zhuǎn)換矩陣項是角度(以度為單位)而個曲面的名字等于原有曲面的名字加上1000倍這個柵元的名字。這個公式給出產(chǎn)生曲面的次完整地描述。3)LAT柵格卡lat=l意味著這個柵格是一個六面體；lat=2意味著柵格是一個六棱柱。LAT卡的非零項意味著相應(yīng)柵元是一個(0,0,0)柵格元素。一個柵格柵元的柵元描述有二個主要用途。一是這個描述是標(biāo)準(zhǔn)的MCNP柵元描述，二是這個柵元曲面描述的順序確定位于這個柵元每個曲面外邊的柵格元素的順序。當(dāng)確定柵格以后，就要決定(0,0,0)柵格元素，并且決定這個柵格的三個柵格指標(biāo)在規(guī)定在柵元卡上填寫界定(0,0,0)元素的曲面。對一個六面體柵格柵元，在第一個列的曲面的外面是(1,0,0)元素，在第二個列的曲面的外面是(-1,0,0)元素；然后按照那個順?biāo)?。這個方法提供柵格排列的順序，因此當(dāng)你指定元素(7,9,3)時，程序就會知道這個元素是哪一個。對六棱柱柵格柵元，在第一個列的曲面外面是(l,0,0)元素，在第二個列的曲面外面是(-1,0,0)元素；然后依次是(0,l,0)，(0,-1,0),(-l,l,0),(l,-1,0)， (0,0,1),(0,0,-l),如圖2-3(b)所示。描述(0,0,0)元素的最后兩個曲面必須是棱柱(a)六面體柵元排列順序(b)六棱柱柵元排列順序好填滿空間，這意味著柵格的對邊必須是相等并平行的。六面體柵格柵元可以在一個或兩個柵格是左旋或右旋都可以。一個柵格必須在它的Universe里。一個問題的實際領(lǐng)域本身可以是一個柵格。如果一個粒子離開實際領(lǐng)域的最后一個柵元(限制柵格的圍)，這個粒子死4)FILL卡充柵元的坐標(biāo)系和在輔助坐標(biāo)系下用來填充柵元的Universe之間的變換。如果沒有指定變換，則這個Universe繼承被填充柵元的變換。如果旋轉(zhuǎn)角是以度為單位而不是余弦，則使度為單位的角度而又要按余弦來填寫一些角度，那么所有FILL及*FILL數(shù)據(jù)都必須放在輸ILL這個柵陣的每一個柵元都是用同一個Universe填充。如果是一個數(shù)組，填充由這個數(shù)組所覆蓋柵陣的那部分；柵陣的其余部分不存在?？捎貌煌腢niverse去填充這個柵陣的各種數(shù)、負(fù)數(shù)或者零表示，并用：號隔開。每個柵格元素的指標(biāo)是由它相對于(0,0,0)元素的該柵陣柵元卡上前兩個曲面定義第一個柵格指標(biāo)的方向。指標(biāo)的數(shù)值圍取決于這個柵格的 元素的位置的是不同的。在柵格柵元卡上第三和第四列的曲面定義第二個柵格指標(biāo)的方向。那些柵陣柵元存在。每一個數(shù)組元素的值是填充相應(yīng)柵格的Universe號。在這個數(shù)組里可陣元素不用任何Universe去填充。如果零Universe(問題的實際領(lǐng)域)是一個柵格，這個變換號或者變換本身。fill1:20:1442040433040(0,1,1)和(l,2,1)都是用Universe4填充；元素(2,l,0)用Universe2填充；元素(l,l,1)和(2,l,l)用Universe3填充。4.2.4TRN坐標(biāo)變換卡格式：trnOOOBBBBBBBBBM123123456789O，O，O：變換的位移矢量；123B~B：變換的旋轉(zhuǎn)矩陣；19M：M=l意味著位移矢量是輔助坐標(biāo)系的原點在基本坐標(biāo)系定義的位置。M=-l意味著位移矢量是基本坐標(biāo)系的原點在輔助坐標(biāo)系定義的位置。缺?。簍rn0001000100011使用：可選擇卡MCNP所有類型的曲面都有某種對稱性，當(dāng)對稱軸平行于某些坐標(biāo)軸時，即使較為復(fù)雜參數(shù)難度就大多了。在這種情況下，用戶可以在一個較為方便的輔助坐標(biāo)系來描述該曲面，是容易的，但這時在該曲面描述卡上要有變換號n，在相應(yīng)的TRn卡指定這個變換。環(huán)面 (Torus)不能做坐標(biāo)變換。在MCNP的輸入文件中，基本上是在輔助坐標(biāo)系下給出較為簡單的曲面描述卡。當(dāng)輔助一個問題中最多可有10個TRn卡。B矩陣說明兩個坐標(biāo)系軸的方向之間的關(guān)系。B是i在基本坐標(biāo)系(XYZ)軸與輔助坐標(biāo)系(x’y’z’)軸之間的夾角余弦(若在TRn前加了元素BBBBBBBBB123456789B的意義與M無關(guān)。用戶不必給出所有B值，可用kJ跳過不指定的項，B可按以下任iiii)B全部給出(9個值)；iii)給出矩陣中列或行方向上三個矢量中的兩個(6個值)。程序?qū)⒂伤鼈兊牟娉水a(chǎn)生iii)給出矩陣中列與行方向上各一個矢量，其中的公共分量必須小于1(共5個值)。用歐拉角方法填完該矩陣；iv)給出矩陣中列或行方向上的一個矢量(3個值)，程序?qū)茨撤N隨機方法生成其他一個矢量由三個元素組成，這些元素可以按矩陣的橫向或豎向排列，但不能走對角錢。0)的矢量將保持原樣。如果正交性大于0.001弧度時，則給出一個警告信息。轉(zhuǎn)面時，使用方式iv)是合適的。曲面不是旋轉(zhuǎn)面(一般情況)時，方式ii)與iii)大系，則使用方式i)是合適的。文件中這個問題的坐標(biāo)系統(tǒng)之間的關(guān)系。tr47.91.30-10001-1004.3降低方差卡下面是定義降低方差卡的參數(shù)的卡片。助記名卡片容助記名卡片容IMP柵元重要性卡ESPLT能量分裂與輪盤賭卡TSPLTEXT時間分裂或輪盤賭卡PWT指數(shù)變換參量VECT光子權(quán)重卡FCL強迫碰撞參量WWE粒子權(quán)窗能量WWNWWN粒子權(quán)窗邊界WWP粒子權(quán)窗參量WWG粒子權(quán)窗產(chǎn)生器WWGE粒子權(quán)窗產(chǎn)生能量MESH網(wǎng)格權(quán)窗產(chǎn)生器的網(wǎng)格層次重要性PDnDXCDXTRAN(球)柵元貢獻(xiàn)BBREM軔致輻射偏倚柵元或曲面，需要提供適當(dāng)?shù)娜笔≈诞?dāng)作空格間隔(可使用nR重復(fù)或者nJ跳躍)。當(dāng)IMP這些卡既需要蒙特卡羅方法的知識，也需要一定的降低方差技巧。4.3.1IMP柵元重要性卡n=對于中子是N；對于光子是P；對于電子是E。如果不同粒子的重要性相同，也允許I=問題中的柵元數(shù)柵元的重要性有兩個用途：a)當(dāng)重要性為0時，用它來終止粒子的歷史；b)用于幾何分裂與輪盤賭技巧，促迸粒子軌跡朝向幾何上較為重要的區(qū)域運動。在中子-光子耦合輸運問題中，用戶可能希望對光子的分裂與中子有所不同，所以可能需要用到兩組重要性(中子和光子)。如果在運行方式為np的問題中，未給出imp:p卡，1。imp:nl22m0120r應(yīng)記住，凡進(jìn)入“0”重要性柵元的軌跡都將被殺死。對探測器進(jìn)行記數(shù)或向DXTRAN。4.3.2ESPLT能量分裂與輪盤賭卡格式：esplt:nRE…RE112020使用：可選擇卡；不能在多組計算中使用。重要性作用設(shè)成相同，并附加一個R=1,代碼中的分裂與輪盤賭都要改變。0權(quán)窗E：粒子實行分裂的能量限(Mev)。i對于能量下降過程中經(jīng)過E的粒子：iR=(R>1)粒子要分裂成的軌跡數(shù)；(0<R<1)俄羅斯輪盤賭的概率。iii對于能量上升過程中(例如裂變)經(jīng)過E的粒子(R<0參看“其他”)：ii1/R=(1/R>1)粒子要分裂成的軌跡數(shù)；(0<1/R<1)輪盤賭的概率。iii其他：如果第一輸入項是負(fù)的，也就是R<0，能量上升時不能做分裂或輪盤賭游戲。i缺?。菏÷赃@卡的粒子沒有能量重要性。T引起的粒子權(quán)重的改變會造成中止能量和權(quán)窗能量的能量修正?？ㄉ系妮斎霐?shù)據(jù)包括最多20對的能量-重要性比率參數(shù)RE。如果能量E不是單調(diào)的，iiiRERiiiiiii粒子會分裂。若0<1/R<1，R代表輪盤賭游戲中的存活概率。如果粒子能量下降到E以iii個能量分界，用Ri得到的值決定輸出。如果能量下降到E以下，一定會發(fā)生特定的分裂或賭。如果能量上升到E以上，通常ii可能并不希望做相反的分裂或輪盤賭游戲(尤其在裂變占優(yōu)勢的問題中)。為此設(shè)置了一種選擇。當(dāng)分裂參數(shù)R給成負(fù)值時，只對能量下降時做分裂或輪盤賭游戲，而在能量上升時，i權(quán)窗如果指定了權(quán)窗，能量的分裂或賭就由權(quán)窗獨自完成。這時的R是權(quán)窗的能量修正。i如果能量下降到E以下，權(quán)窗除以R劃分。相反的，如果能量上升到E以上，權(quán)窗乘以R。iiii如果經(jīng)過了多個能量分界，權(quán)窗的修正由R得到的值決定。i例如卡片：ESPLT:N2.12.01.25.001WWP:N535000JJ過分界E會發(fā)生分裂或賭，同沒有權(quán)窗時一樣。缺省值對很多問題都適用。i例如卡片：ESPLT:N2.12.01.25.001WWP:N535000J14.3.3TSPLT時間分裂或輪盤賭卡nN子為P；電子為E缺省：如果省略了這個卡片表示這些粒子不發(fā)生時間分裂。使用：可選擇卡；不能在多組計算中使用。相反，TSPLT中的Ri表示上升能量的分裂或輪盤賭。)對于小于最小時間T的能量，絕對重i0產(chǎn)生x粒子的粒子在晚期都是重要的。因此一個MODENP問題，例如指定一個TSPLT：PR=(R>1)粒子要分裂成的軌跡數(shù)；(0<R<1)俄羅斯輪盤賭的概率。iiiT：粒子經(jīng)歷分裂或賭的時間，而且T必須是單調(diào)增或單調(diào)減的。ii缺?。菏÷赃@卡的粒子沒有時間重要性。的分裂與輪盤賭可以使用IMP卡。在MODEPEF8加上高度記數(shù)和多組問題不能使用卡上的輸入數(shù)據(jù)包括最多20對的時間-重要性比率參數(shù)RT。如果時間T不是單調(diào)iii的，會發(fā)出警告信息。R可以不是整數(shù)也可以是0和1之間的數(shù)。R>1，粒子會發(fā)生ii分裂。0<R<1，R代表輪盤賭游戲中的存活概率。T的單位是剎(10-12s)，如果粒子時iii間大于T，抽樣進(jìn)行指定的分裂或賭。iN權(quán)窗如果指定了權(quán)窗，能量的分裂或賭就由權(quán)窗獨自完成。這時的R是權(quán)窗的時間修正。i如果粒子時間超過T，權(quán)窗除以R劃分。如果經(jīng)過了多個能量分界，權(quán)窗的修正由R得到iiiNWWP:N535000JJiiii過分界T會發(fā)生分裂或賭，同沒有權(quán)窗時一樣。缺省值對很多問題都適用。iNWWP:N535000J1一分為二；當(dāng)時間超過10000剎時，將以0.2的殘存概率進(jìn)行輪盤賭。另外，這個例子指4.3.4PWT光子權(quán)重卡格式：pwtww…w…w12iIW:中子在第i個柵元碰撞所產(chǎn)生的光子權(quán)重的值；iI的柵元總數(shù)。缺?。喝绻笔〈丝?，則所有柵元的W均被置為-l。ie對每個柵元i,可以指定一個最小的光子權(quán)W(＞0)，中子誘發(fā)的光子將以高于w*I/Iiiis的權(quán)出生，其中I和I分別表示中子從第s個柵元出發(fā)在第i個柵元碰撞的兩個柵元的重si對于每一個具有負(fù)的w(<0)項的柵元i，僅當(dāng)中子誘發(fā)的光子權(quán)重大于-w*w*I/I時，iisis則產(chǎn)生光子，其中w是導(dǎo)致當(dāng)前光子產(chǎn)生事件中的源中子的初始權(quán)，I和I分別表示中子ssi當(dāng)PWT卡上指定的w為0，那么在光子產(chǎn)生是可能的情況下，則導(dǎo)致中子在i柵元中的i每次碰撞都考慮產(chǎn)生光子。想要關(guān)掉i塊中的光子產(chǎn)生過程時，應(yīng)置w=-1.0e6。iPWT卡通過誘發(fā)光子總的權(quán)重與在PWT卡上指定的有關(guān)的權(quán)重閾值的比較來控制中子在碰撞時誘發(fā)光子的數(shù)目。權(quán)重閾值與柳元重要性有關(guān)，并且如果w＜0，也與源中子權(quán)重有i對使用光子?xùn)旁匾?imp:p)，而不使用權(quán)窗(WWNn:P)的問題，PWT卡上的數(shù)據(jù)項最好為缺省項，即W=-1，或者把每個柵元的w是設(shè)置為平均源的權(quán)。ii4.3.5EXT指數(shù)變換卡格式：ext:pAA…A…A12iIA:第i個柵元的數(shù)據(jù)項，每個A的格式均為A=QV,其中Q為拉伸的長度，iimV定義拉伸的方向；m缺省：此即非偏移輸運的情況，則所有柵元的A=0。i在強迫碰撞的柵元中或者在沒有權(quán)窗控制的棚元中不使用指數(shù)變換。僅當(dāng)粒子通量有一個指數(shù)一個指數(shù)分布時，如強吸收問題，使用指數(shù)變換合適。指數(shù)變換方法是通過調(diào)整總截面在所感興趣的方向上拉伸碰撞之間的徑跡長度?？偨孛?=(1P)ttt∑：真的總截面；tQP0，不使用指數(shù)變換；Qp0＜P＜l，固定的拉伸參數(shù)；ataat根據(jù)EXT卡上的每個A項的V的三個選擇來定義拉伸的方向。im(1)A項沒有V部分，也就是對一個給定的柵元僅填寫拉伸度A=Q。這樣只在粒子的imi上作隱式俘獲，在這種情況下人A=Q=S，而且抽取的距離是散射距離而不是抽取碰撞的距離。i (2)指定拉伸方向為VmVm是從碰撞點到點(X，Y，Z)的方向，其中(X，Y，Z)mmmmmm是在VECT卡上指定的。方向μ是粒子方向和拉伸方向V之間的夾角余弦，A的符號給出拉mm伸的方向是朝向(X，Y，Z)還是遠(yuǎn)離(X，Y，Z)。mmmmmm (3)指定拉伸方向為V=X，Y或Z，那么方向μ分別為粒子方向與X-，Y-或Z-軸之間m的夾角余弦，A的符號給出拉伸的方向是朝向X，Y或Z軸，還是遠(yuǎn)離這三個軸。i例子：ext:n00.7VS-SV-.6V0.5VSZ-.4X2299vectV000V111929AiS-SV2-.6V9-.4XQSSSVV2.-V2-V9V9Z拉申參數(shù)P=.7P=∑/∑atP=∑/∑atP=.6P=.5P=∑/∑atP=.4方向指向點(1,1,1)粒子方向遠(yuǎn)離點(1,1,1)遠(yuǎn)離原點指向原點沿著Z軸4.3.6VECT向量輸入卡格式：vectvxyz…vxyz…mmmmnnnnm，n：標(biāo)記向量v，v任意號，這個號是任意的；mnx，y；z：為對應(yīng)向量v的定義坐標(biāo)。mmmmVECT卡的數(shù)據(jù)項是四個為一組，這些組定義指數(shù)變換中的若干個向量或者用戶修補的4.3.7FCL強迫碰撞卡格式：fcl:nxx…x…x12iIlxlii缺省值為x=0，無強迫碰撞；I為問題中的柵元數(shù)。i應(yīng)用了強迫碰撞功能的柵元不再進(jìn)行權(quán)窗游戲。如果x=0，就關(guān)掉了對這一柵元的強迫碰撞功能。用適當(dāng)?shù)臋?quán)重控制把進(jìn)入柵元i的全i部粒子分為碰撞部分和不碰撞部分。如果|xi|<1以|xi|的存活概率對碰撞部分施行俄羅斯輪盤賭以防止碰撞的歷史數(shù)變得太大。如果一些強迫碰撞的柵元彼此鄰近，則建議x選為小i當(dāng)x<0時，強迫碰撞過程僅適用于進(jìn)入柵元i的粒子。強迫碰撞之后，不考慮權(quán)重截i當(dāng)x>0時，強迫碰撞過程適用于進(jìn)入柵元i的粒子以及權(quán)截斷或權(quán)窗游戲幸存的碰撞i粒子。粒子將繼續(xù)被分為非碰撞和碰撞兩部分(概率為|xi|)，直到這些粒子被權(quán)截斷或權(quán)一般情況下令x=l或x=-1。當(dāng)一些強迫碰撞柵元彼此鄰近或強迫碰撞粒子產(chǎn)生的數(shù)ii比所要的粒子數(shù)高些，應(yīng)置x為小數(shù)。i當(dāng)權(quán)窗邊界把一個進(jìn)入柵元i的典型權(quán)重分類時，選x>0。當(dāng)權(quán)窗界把強迫碰撞粒子的i典型權(quán)重分類時，選x>0。i當(dāng)應(yīng)用重要性時，選x>0，因為x<0關(guān)掉權(quán)截斷游戲。ii在使用指數(shù)變換的柵元中，不可做強迫碰撞。4.3.8權(quán)窗卡再充分細(xì)分的需要。權(quán)窗是對指定的空間和能量重要性函數(shù)提供重要性(IMP:n卡)和能量分裂(ESPLT：n卡)的另一種方法。權(quán)窗的優(yōu)點是：(1)權(quán)窗提供一個空間和能量的重要性函數(shù)；(2)權(quán)窗控制粒子的權(quán)重；(3)與其他一些降低方差技巧不矛盾；(4)在曲面跨越、碰撞或者這兩者都可以使用權(quán)窗；(5)能夠控制分裂或者輪盤賭的程度；(6)在指定的空間或者能量區(qū)域上關(guān)掉權(quán)窗；(7)用權(quán)窗產(chǎn)生器自動產(chǎn)生權(quán)窗。權(quán)窗的缺點是：權(quán)窗不如重要性簡單；當(dāng)改變源的權(quán)重時，權(quán)窗必須重新歸一。在重復(fù)結(jié)構(gòu)中，柵元重要性和權(quán)窗還有一點不同。對于柵元重要性(IMP卡)，早univercity中的柵元重要性是被填充柵元的重要性的乘積，各種行為也以這個重要性為依復(fù)結(jié)構(gòu)中推薦使用網(wǎng)格權(quán)窗。格式：wwe:nEE…E…E,j≦9912ijE為能量間隔，E<E<E，E＝0；ii-1i0格式：wwni:nww…w…wi1i2ijiJw為柵元j的權(quán)窗下限和由WWE卡定義的能量間隔E<E<E，E=0。如果沒有WWE卡，iji-1i0使用：當(dāng)使用重要性卡時，不要求此卡。除此之外，要求權(quán)窗(WWN和WWP)卡。如果W＜0，任何進(jìn)入柵元j的粒子被殺死，因為負(fù)輸入相當(dāng)于0重要性。如果負(fù)項用ij于一個能量區(qū)間，那么它們也將用于同一個柵元的其它能量區(qū)間。如果W＞0，進(jìn)入柵元j或在柵元碰撞的粒子將根據(jù)WWP卡的選擇進(jìn)行分裂或輪盤賭。ij如果W＝0，關(guān)閉柵元j的能量箱i的權(quán)窗游戲同時打開權(quán)截斷游戲。應(yīng)在CUT卡上ij把權(quán)截斷指定為這個問題所要求的最低的所允許的權(quán)重。否則，進(jìn)入W＝0柵元的大部分粒ij子將被用權(quán)截斷殺死。當(dāng)用單個重要性函數(shù)或者用一組權(quán)窗邊界不能表示各種能量及空間區(qū)顧及到其它一些降低方差技巧引起的權(quán)變化，對每一個柵元每一能量間隔必須指定權(quán)的下例1：定義中子問題的三個能量間隔和四個柵元的權(quán)窗邊界。wwe:nwwn1:nwwn2:nwwn3:n這個例子沒有WWEE1w11w21w31w11EE23www121314www222324www323334ww1213只用WWN:P卡對一個有三個柵元的問題定義了與能量無關(guān)的光4.3.9WWP權(quán)窗參數(shù)卡pnWUPNWSURVNMXSPLNMWHEMSWITCHNMTIMEETSPLTii)WSURVN：對于經(jīng)輪盤賭存活下來的粒子，它的權(quán)將成為WSURVN與權(quán)窗下限值之積v)SWITCHN：當(dāng)它的值大于0時；權(quán)窗的下限值將被指定為SWITCHN除以粒子所在柵ⅶ)WNORM>0權(quán)窗低分界的乘數(shù)。SWITCHN=0,MTIME=0,ETSPLT=0使用：沒有重要性卡需要使用權(quán)窗的時候使用。利用現(xiàn)成的柵元重要性來定義權(quán)窗的下限值當(dāng)采用與能量無關(guān)的權(quán)窗處理時，則可以利用重要性卡上的數(shù)據(jù)給出權(quán)窗參量，此時應(yīng)將WWP卡上的SWITCHN置一個正的常數(shù)C。在這種處理方式下，柵元的權(quán)窗下限值4.3.10權(quán)窗產(chǎn)生器卡度相關(guān)而有為抽樣所說明的適當(dāng)幾何。這樣平均重要性不能代表這個柵元。不適當(dāng)幾何說明也會引起鄰近柵元之間的重要性有大的區(qū)別。恰好產(chǎn)生器能提供為抽樣使用的幾何說明是否適當(dāng)?shù)男畔?。A)WWG權(quán)窗生器卡格式：wwgIIWJJJJItcgEI：問題的記數(shù)號(Fn卡的n)。由TFn卡定義特殊記數(shù)箱，對這個特殊記數(shù)箱權(quán)窗產(chǎn)tI：參照柵元，典型的參照柵元是源柵元；c>0表示使用以I為參照柵元，使用柵元權(quán)窗生成器cW：產(chǎn)生柵元I的權(quán)窗下限值。如果W=0，它將是平均源權(quán)重的一半；gcgE缺?。喝笔〈丝ǎ划a(chǎn)生權(quán)窗值。WWG卡引起記數(shù)I的最佳重要性函數(shù)的產(chǎn)生。對于柵元權(quán)窗生成器，重要性函數(shù)寫在tB)WWGE權(quán)窗產(chǎn)生能重卡格式：wwge：nEE…E…E；j≤1512ijE：產(chǎn)生權(quán)窗組的能量上限，E＞E。ii+1i缺省：如果沒有此卡，并且又使用權(quán)窗，則相應(yīng)于要運行問題的能量/時間圍產(chǎn)生一個Ji生成一個能量/時間和一個時間/能量權(quán)窗。E4.3.11MESH3-47-3-494.3.12PDN探測器貢獻(xiàn)卡格式：PdnPP…P…P12iIn：記數(shù)號；P：柵元i對探測器貢獻(xiàn)的概率。iiPDn卡可減少對指定探測器相對不太重要的柵元對探測器記數(shù)的貢獻(xiàn)，節(jié)約了計算時間。此卡上對每一柵元給出一個概率值P，它表示對于柵元i中的任一碰撞點或源點，將以i概率P判定它實行對探測器的記數(shù)。為保證無偏性，當(dāng)判定實行記數(shù)時，原有的記錄將乘i上修正因子1/P。用戶通過對離測點較遠(yuǎn)(以平均自由程計)的柵元置小于1的P值，便ii可提高運行速度。也可以將某些柵元的P值置為0，有選擇地遏制了這些柵元對探測器的貢i獻(xiàn)。獻(xiàn)缺省此卡，所有柵元的P都取缺省值1。但對所有的探測器記數(shù)，使用PD0卡建立一組i4.3.13DXCDXTRAN貢獻(xiàn)卡P1P2…P…PiIP：柵元i對貢獻(xiàn)的概率，其缺省值為1；i缺省：m=0，P=1i4.3.14BBREM軔致輻射偏倚BBREMbbb...b49m1m2…mnb=任意正數(shù)1b2...b49=軔致輻射能譜的偏倚參數(shù)m…m=調(diào)用偏倚的材料列表1n4.4源的描述卡對每一個MCNP問題都有四個源中的一個。四個源包括：①通用源(SDEF卡)；②曲面源(SSR卡)；③臨界源(KCODE卡)；④用戶提供的源。除了臨界源之外，其它源都能使用名名EFDSn讀曲面源KCODEKSRC臨界計算的源點的描述應(yīng)包括如下信息：i)能量(ERG)ii開始時間(TME))iv)幾何位置(xxx，yyy，zzz)v)粒子類型(IPT)vi)粒子統(tǒng)計權(quán)重(WGT)vii)粒子開始柵元(ICL)viii)粒子開始曲面(JSU)，如果開始點不在曲面上該值為0。4.4.1SDEF通用源卡缺?。耗J(rèn)為零時刻在點(0，0，0)的14Mev點源，權(quán)重為1。使用：對于通用源問題必須使用，對于臨界源問題為選擇性卡。等號可寫可不寫。源變量不必和MCNP其他地方設(shè)置的源變量完全相同。很多都是中間或者一個缺省分布(例如，各向同性角分布)時，第一級存在。當(dāng)用一個概率分布給出一個源變量時，第二級出現(xiàn)。這個級要求SI，SP卡或者其一。當(dāng)一個源變量與另一個源變量有況下，使用顯值或者使用分布代替有關(guān)的缺省值。省值1CELLxxxyyyzzz能23ERG能量(Mev)4TME時間(剎)05DIR (0，1)的余弦分布6VEC體源：必須，各向同性除外7NRM8POS位置抽樣的參考點0，0，09AD到抽樣位置的徑向距離0EXT0AXS有方向XYZCCARAWGT粒子權(quán)重1EFF位置采樣中排除有效判據(jù)PAR源發(fā)射的粒子類型=光子(MODEP或PE)3=電子(MODEE)RCC位置的坐標(biāo)(x，y，z)。用這些變量可以指定三種不同類型的體分布及三種不同類型的面分布能夠近似比較復(fù)雜的一些分布。三個體分布是長方體、球及圓柱。對每一個體分布SUR的值為0(缺省值)。一個體分棄這個點，再重新抽取另一個點。如果使用這個技巧(柵元舍棄)，必須確保抽樣區(qū)域包含舍棄代替柵元或者和柵元舍棄聯(lián)合。元的面要垂直于坐標(biāo)軸。義為從該球的中心到粒子抽樣位置的距離。然后在半徑為RAD的球面上均勻抽取粒子的位面上一個圓對稱源。展到柵元部的離。這樣，既不會對結(jié)果有明顯的影響，也防止了粒子丟失。位置。但在這種情況下，必須確保這個點在這個曲面上，因為MCNP不作檢查。如果沒有指狀可以是回轉(zhuǎn)橢圓體、球或者平面。柵元舍棄抽樣在這里還適用。但是使用Cookie-Cutter舍棄能夠做柵元舍棄要做的任何事情。圓柱體分布必須指定為退化體分布。義的橢球軸必須平行于一個坐標(biāo)軸。目前還沒有提供有關(guān)橢球曲面的簡單的不均勻抽樣或者偏倚抽樣?？梢允褂靡恍〤ookie－Cutter柵元的一個分布以產(chǎn)生粗糙的不均勻的位置分布。如果SUR的值是球的曲面名字，則在這個曲面上抽取粒子的位置。如果沒有指定矢量方向和球的中心到粒子位置點的矢量之間的夾角余弦。方位角在0o~360o之間的均勻抽取。比較多的粒子從球離記數(shù)區(qū)域的最近的一邊出發(fā)，通常來講，EXT用指數(shù)分布函數(shù)(-31)樣，在上述情況下a＝1是缺省值。Cookie-Cutter柵元太復(fù)雜。對一個曲面源Cookie-Cutter柵元要與源的曲面相交。用一些平面及柱面界定的均勻截面的無限長柵元是合適作一個平面源的Cookie-Cutter柵元。注意：柵元舍棄或者CCC舍棄(Cookie-Cutter)與位置偏倚抽樣結(jié)合幾乎是不正確的游戲。如果使用這種組合，必須確保這種組合是一個合理的游戲，因為MCNP沒有能力檢查種錯誤。源變量SUR、VEC、NRM及DIR是用來確定源粒子的初始飛行方向。根據(jù)參考方向VEC間均勻抽取。如果對一個體分布(SUR=0)沒有指定VEC和DIR，則按缺省值產(chǎn)生各向同性該曲面具有由NRM所確定的符號。如果在一個曲面源的分布上沒有指定DIR，則按缺省值使用余弦分布P(DIR)=2*DIR，0＜DIR＜l。使用DIR的偏倚分布可使較多的粒子沿著朝向幾何記數(shù)區(qū)域的方向出發(fā)。一般來講指數(shù)分布最合適。中