MCNP-5A程序使用說(shuō)明書

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