常用蒙特卡羅程序介紹ppt課件_第1頁
常用蒙特卡羅程序介紹ppt課件_第2頁
常用蒙特卡羅程序介紹ppt課件_第3頁
常用蒙特卡羅程序介紹ppt課件_第4頁
常用蒙特卡羅程序介紹ppt課件_第5頁
已閱讀5頁,還剩41頁未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、l蒙特卡羅方法應(yīng)用軟件的特點蒙特卡羅方法應(yīng)用軟件的特點l常用的通用蒙特卡羅程序簡介常用的通用蒙特卡羅程序簡介lMCNP程序輸入的描述程序輸入的描述 建立完善的通用蒙特卡羅程序可以避免大量的重復(fù)性工作,并且可以在程序的基礎(chǔ)上,開展對于蒙特卡羅方法技巧的研究以及對于計算結(jié)果的改進(jìn)和修正的研究,而這些研究成果反過來又可以進(jìn)一步完善蒙特卡羅程序。通用蒙特卡羅程序通常具有以下特點:具有靈活的幾何處理能力參數(shù)通用化,使用方便元素和介質(zhì)材料數(shù)據(jù)齊全能量范圍廣,功能強(qiáng),輸出量靈活全面含有簡單可靠又能普遍適用的抽樣技巧具有較強(qiáng)的繪圖功能l MORSE程序程序l 較早開發(fā)的通用蒙特卡羅程序,可較早開發(fā)的通用蒙特卡

2、羅程序,可以解決中子、光子、中子光子的聯(lián)合輸運(yùn)以解決中子、光子、中子光子的聯(lián)合輸運(yùn)問題。采用組合幾何結(jié)構(gòu),使用群截面數(shù)據(jù),問題。采用組合幾何結(jié)構(gòu),使用群截面數(shù)據(jù),程序中包括了幾種重要抽樣技巧,如俄國輪程序中包括了幾種重要抽樣技巧,如俄國輪盤賭和分裂技巧,指數(shù)變換技巧,統(tǒng)計估計盤賭和分裂技巧,指數(shù)變換技巧,統(tǒng)計估計技巧和能量偏移抽樣等。程序提供用戶程序,技巧和能量偏移抽樣等。程序提供用戶程序,用戶可根據(jù)需要編寫源分布以及記錄程序。用戶可根據(jù)需要編寫源分布以及記錄程序。l EGS程序程序l EGS是是Electron-Gamma Shower 的的縮寫,它是一個用蒙特卡羅方法模擬在任意縮寫,它是一

3、個用蒙特卡羅方法模擬在任意幾何中,能量從幾個幾何中,能量從幾個KeV到幾個到幾個TeV的電子的電子-光子簇射過程的通用程序包。由美國光子簇射過程的通用程序包。由美國Stanford Linear Accelerator Center提供。提供。EGS于于1979年第一次公開發(fā)表,提供使用。年第一次公開發(fā)表,提供使用。EGS4是是1986年發(fā)表的年發(fā)表的EGS程序的最新版本。程序的最新版本。l MCNP程序程序l MCNP是美國是美國Los Alamos國家實驗國家實驗室開發(fā)的大型多功能通用蒙特卡羅程序,可室開發(fā)的大型多功能通用蒙特卡羅程序,可以計算中子、光子和電子的聯(lián)合輸運(yùn)問題以以計算中子、光

4、子和電子的聯(lián)合輸運(yùn)問題以及臨界問題,中子能量范圍從及臨界問題,中子能量范圍從10-11MeV至至20MeV,光子和電子的能量范圍從,光子和電子的能量范圍從1KeV至至1000MeV 。程序采用獨特的曲面組合幾何結(jié)。程序采用獨特的曲面組合幾何結(jié)構(gòu),使用點截面數(shù)據(jù),程序通用性較強(qiáng),與構(gòu),使用點截面數(shù)據(jù),程序通用性較強(qiáng),與其它程序相比,其它程序相比,MCNP程序中的減方差技巧程序中的減方差技巧是比較多而全的。是比較多而全的。 MCNP的輸入包括幾個文件,但主要的的輸入包括幾個文件,但主要的一個是由用戶編寫的一個是由用戶編寫的INP文件,該文件包括描文件,該文件包括描述問題所必須的全部輸入信息。文件采

5、用卡述問題所必須的全部輸入信息。文件采用卡片結(jié)構(gòu),每行代表一張卡片,文件由一系列片結(jié)構(gòu),每行代表一張卡片,文件由一系列卡片組成,對于任一特定的問題,只需用到卡片組成,對于任一特定的問題,只需用到INP全部輸入卡片的一小部分。全部輸入卡片的一小部分。MCNP輸入文件中物理量的單位輸入文件中物理量的單位長度厘米能量MeV時間10-8 秒溫度MeV(kT)原子密度1024 個原子 / 厘米3質(zhì)量密度克 / 厘米3截面10-24 厘米2原子量中子質(zhì)量的1.008664967倍阿伏加德羅常數(shù) 6.0231023 l 輸入文件的基本形式輸入文件的基本形式l 信息塊信息塊l 信息塊的卡片放在信息塊的卡片放在

6、INP文件中標(biāo)題卡文件中標(biāo)題卡之前。信息塊給出了之前。信息塊給出了MCNP的一些運(yùn)行信息,的一些運(yùn)行信息,信息塊上各部分的意思和運(yùn)行行信息是一樣信息塊上各部分的意思和運(yùn)行行信息是一樣的,當(dāng)運(yùn)行行信息與信息塊中所指定的信息的,當(dāng)運(yùn)行行信息與信息塊中所指定的信息相矛盾時,則忽略信息塊中相應(yīng)的信息,而相矛盾時,則忽略信息塊中相應(yīng)的信息,而以運(yùn)行行信息為準(zhǔn)。以運(yùn)行行信息為準(zhǔn)。l信息塊是可選的,信息塊的第一張卡信息塊是可選的,信息塊的第一張卡片,必須在第片,必須在第 18 列寫上列寫上“MESSAGE:”,從第一張卡片的第從第一張卡片的第 980 列到后續(xù)卡片的第列到后續(xù)卡片的第 180 列都可填寫運(yùn)行

7、信息。在標(biāo)題卡之前用列都可填寫運(yùn)行信息。在標(biāo)題卡之前用一個空行分隔符結(jié)束信息塊。一個空行分隔符結(jié)束信息塊。l 初始運(yùn)行的輸入文件初始運(yùn)行的輸入文件信息塊空行分隔符選擇項標(biāo)題卡僅一行,占用第 180 列。作為輸出標(biāo)題。柵元卡空行分隔符定義構(gòu)成整個系統(tǒng)的各個基本介質(zhì)單元以及相應(yīng)的物理信息。曲面卡空行分隔符定義組成柵元的曲面信息。數(shù)據(jù)卡空行分隔符其它數(shù)據(jù),包括問題類型、源描述、材料描述、計數(shù)描述,問題截斷條件等。其它選擇項l 接續(xù)運(yùn)行的輸入文件接續(xù)運(yùn)行的輸入文件l接續(xù)運(yùn)行必須在運(yùn)行行信息或信息塊中給接續(xù)運(yùn)行必須在運(yùn)行行信息或信息塊中給出出C項選擇,即項選擇,即Cm,表示從,表示從RUNTPE文件中文

8、件中讀出第讀出第m次轉(zhuǎn)儲的內(nèi)容接著運(yùn)算,如果次轉(zhuǎn)儲的內(nèi)容接著運(yùn)算,如果m未指未指定,則讀最后一次轉(zhuǎn)儲的數(shù)據(jù)。如果不需要定,則讀最后一次轉(zhuǎn)儲的數(shù)據(jù)。如果不需要改變內(nèi)容,則不需要接續(xù)輸入文件,僅需運(yùn)改變內(nèi)容,則不需要接續(xù)輸入文件,僅需運(yùn)行行RUNTPE以及在運(yùn)行行加上以及在運(yùn)行行加上C選擇。選擇。信息塊空行分隔符選擇項CONTINUE寫在第 18 列數(shù)據(jù)卡空行分隔符只允許部分?jǐn)?shù)據(jù)卡。(FQ,DD,NPS,CTME,IDUM,RDUM,PRDMP,LOST,DBCN,PRINT,KCODE,MPLOT,ZA,ZB,和 ZC)其它選擇項l 卡片格式卡片格式l INP輸入文件的每一行稱之為一輸入文件的每

9、一行稱之為一張卡片都限于使用第張卡片都限于使用第 180 列并構(gòu)成卡片列并構(gòu)成卡片映象。大部分輸入卡片按行填寫;然而,對映象。大部分輸入卡片按行填寫;然而,對數(shù)據(jù)卡允許按列填寫。數(shù)據(jù)卡允許按列填寫。 $ 符號為它所在那行符號為它所在那行數(shù)據(jù)的結(jié)束符,在數(shù)據(jù)的結(jié)束符,在 $ 符號后面的內(nèi)容作為注符號后面的內(nèi)容作為注釋,它可從釋,它可從 $ 符號后面的任一列開始。符號后面的任一列開始。l標(biāo)題卡只占一行,整行都可填入用戶標(biāo)題卡只占一行,整行都可填入用戶需要的信息,也可以是空行。但要注意在其需要的信息,也可以是空行。但要注意在其它地方使用空行是作為結(jié)束符或者分隔符。它地方使用空行是作為結(jié)束符或者分隔符

10、。l 輸入文件中,在標(biāo)題卡之后及最后輸入文件中,在標(biāo)題卡之后及最后的空行結(jié)束卡之前的任何地方都可插入注釋的空行結(jié)束卡之前的任何地方都可插入注釋卡。注釋卡必須是字母卡。注釋卡必須是字母“C寫在寫在 15 列中列中的任意位置,且至少用一個空格隔開后面的的任意位置,且至少用一個空格隔開后面的注釋內(nèi)容。注釋內(nèi)容。l 行輸入格式行輸入格式l 柵元卡、曲面卡和數(shù)據(jù)卡的書寫格柵元卡、曲面卡和數(shù)據(jù)卡的書寫格式是相同的。必須從式是相同的。必須從 15 列開始填寫這些卡列開始填寫這些卡片相應(yīng)的名字片相應(yīng)的名字(或編號或編號)和粒子標(biāo)識符,后面填和粒子標(biāo)識符,后面填寫用空格分隔的數(shù)據(jù)項。假設(shè)寫用空格分隔的數(shù)據(jù)項。假

11、設(shè) 15 列為空,列為空,則表示它是前一張卡片的繼續(xù)卡。如果在一則表示它是前一張卡片的繼續(xù)卡。如果在一行的末尾有一個用空格隔開的符號行的末尾有一個用空格隔開的符號“&”,則,則表示下一行是該行的繼續(xù)卡,數(shù)據(jù)可填寫在表示下一行是該行的繼續(xù)卡,數(shù)據(jù)可填寫在 180 列。一個數(shù)據(jù)項必須在一張卡片上寫完,列。一個數(shù)據(jù)項必須在一張卡片上寫完,不得跨到下一張卡片上。完全空白的一行則不得跨到下一張卡片上。完全空白的一行則為兩組卡片的分隔符。為兩組卡片的分隔符。 對任何給定的帶有粒對任何給定的帶有粒子標(biāo)識符的類型卡只能有一張。需要整數(shù)的子標(biāo)識符的類型卡只能有一張。需要整數(shù)的數(shù)據(jù)項必須填寫整數(shù),其它數(shù)據(jù)

12、可填寫為整數(shù)據(jù)項必須填寫整數(shù),其它數(shù)據(jù)可填寫為整數(shù)或浮點數(shù)以及數(shù)或浮點數(shù)以及MCNP能讀的數(shù)據(jù)。能讀的數(shù)據(jù)。為書寫方便,可以使用四項書寫功能:nR功能,表示將它前面的數(shù)據(jù)重復(fù)n次。例如:2 4R 等同于 2 2 2 2 2nI功能,表示在與其前后相鄰的兩個數(shù)之間,插入n個線性插值點。對于 X nI Y 的結(jié)構(gòu),如果X和Y是整數(shù),且XY剛好是n+1的整倍數(shù),則產(chǎn)生標(biāo)準(zhǔn)的整數(shù)插值,否則產(chǎn)生實數(shù)插值,但Y值直接存儲。例如:1.5 2I 3.0 1.5 2.0 2.5 32.0可能不精確而 1 4I 6 1 2 3 4 5 6都是精確定整數(shù)XM功能,它表示的數(shù)值為前面的數(shù)據(jù)乘上X。例如:1 1 2M

13、2M 4M 2M 1 1 2 4 16 32nJ功能,表示其后n個數(shù)據(jù)項使用缺省值。例如:DD .1缺省值) 1000 DD J 1000如果nR、nI、及nJ項中缺省n,則假設(shè)n1。l 列輸入格式列輸入格式l列輸入塊的格式:列輸入塊的格式:l Si必須是必須是MCNP卡片名字,它們必須全部卡片名字,它們必須全部是柵元參數(shù)、或者全部是曲面參數(shù)、或者全是柵元參數(shù)、或者全部是曲面參數(shù)、或者全部是其它參數(shù)。部是其它參數(shù)。15 列672 列S1 S2 SmK1D11 D12 D1mK2D21 D22 D2mKnDn1 Dn2 Dnml 粒子標(biāo)識符粒子標(biāo)識符l 幾個輸入卡片都需要粒子標(biāo)識符以幾個輸入卡片

14、都需要粒子標(biāo)識符以區(qū)別中子、光子和電子的輸入數(shù)據(jù)。這些卡區(qū)別中子、光子和電子的輸入數(shù)據(jù)。這些卡片是:片是:IMP、EXT、FCL、WWN、WWE、WWP、WWGE、DXT、DXC、F、F5X、F5Y、F5Z、PHYS、ELPT、ESPLT、CUT和和PERT。粒子標(biāo)識符由上述卡片名字后面的。粒子標(biāo)識符由上述卡片名字后面的冒號、字母冒號、字母N、P或或E組成。組成。l例如:中子重要性卡為例如:中子重要性卡為 IMP:Nl 光子重要性卡為光子重要性卡為 IMP:Pl 缺省值缺省值l MCNP的許多輸入?yún)?shù)都有缺省值,的許多輸入?yún)?shù)都有缺省值,因此用戶不需要每次都給出各個輸入?yún)⒘康囊虼擞脩舨恍枰看?/p>

15、都給出各個輸入?yún)⒘康闹?。?dāng)缺省值符合用戶要求時,便可不在輸值。當(dāng)缺省值符合用戶要求時,便可不在輸入文件中指定。當(dāng)省略某張輸入卡時,則該入文件中指定。當(dāng)省略某張輸入卡時,則該卡上的全部參數(shù)均使用缺省值。如果只想改卡上的全部參數(shù)均使用缺省值。如果只想改變一張卡上的某一個特定參量時,則它前面變一張卡上的某一個特定參量時,則它前面的參量仍需指明,或者用的參量仍需指明,或者用nJ方式跳過前面那方式跳過前面那些使用缺省值的參量。些使用缺省值的參量。l例如:光子截斷卡例如:光子截斷卡 CUT:P 3J -.10l表示前表示前3個參量使用缺省值,只改變第四項個參量使用缺省值,只改變第四項參量的值。參量的值。l

16、 輸入錯誤信息輸入錯誤信息l MCNP對輸入文件出現(xiàn)的錯誤作廣對輸入文件出現(xiàn)的錯誤作廣泛的檢查,如果用戶違反了輸入說明的規(guī)定,泛的檢查,如果用戶違反了輸入說明的規(guī)定,將在終端上以及輸出文件中打印致命錯誤信將在終端上以及輸出文件中打印致命錯誤信息,息,MCNP不再進(jìn)行粒子輸運(yùn)計算,作業(yè)中不再進(jìn)行粒子輸運(yùn)計算,作業(yè)中斷。斷。l 第一個出現(xiàn)的致命錯誤是真的,而第一個出現(xiàn)的致命錯誤是真的,而后面的錯誤可能不一定是真的,這取決于前后面的錯誤可能不一定是真的,這取決于前面出現(xiàn)的致命錯誤的情況。面出現(xiàn)的致命錯誤的情況。l 若在若在MCNP運(yùn)行行上指定運(yùn)行行上指定FATAL項,項,則則MCNP忽略致命錯誤,照

17、常運(yùn)行。忽略致命錯誤,照常運(yùn)行。l 對于對于MCNP的警告信息,用戶不應(yīng)的警告信息,用戶不應(yīng)忽視,應(yīng)搞清楚它們的含義。忽視,應(yīng)搞清楚它們的含義。l 檢查幾何錯誤檢查幾何錯誤l MCNP在處理輸入文件的數(shù)據(jù)時,在處理輸入文件的數(shù)據(jù)時,不能檢查一種非常重要的輸入錯誤。即不能檢查一種非常重要的輸入錯誤。即MCNP無法查出各柵元之間的重疊和空隙,無法查出各柵元之間的重疊和空隙,只有當(dāng)粒子丟失時,才會發(fā)現(xiàn)幾何錯誤。即只有當(dāng)粒子丟失時,才會發(fā)現(xiàn)幾何錯誤。即使如此,可能仍然無法準(zhǔn)確判斷錯誤性質(zhì)。使如此,可能仍然無法準(zhǔn)確判斷錯誤性質(zhì)。l 柵元描述卡柵元描述卡格式: jm d geom params或:jLIK

18、E n BUT listj柵元號,1 j 99999,寫在第 15 列上。m柵元材料號,與材料卡(Mm)中的序號對應(yīng)。m 0 為真空柵元。d柵元材料密度。正值為原子密度,負(fù)值為質(zhì)量密度。對于真空柵元,該項缺省,不填寫。geom柵元的幾何說明。由一系列帶符號的曲面號經(jīng)過布爾運(yùn)算組成。params任選的柵元參數(shù)說明。n另一個柵元的名字(編號)。list描述柵元j和柵元n之間差別的柵元參數(shù)。 在柵元的幾何說明中,關(guān)于曲面的指向是一個很重要的概念。假定曲面 S 的曲面方程為 f (x,y,z)0,則對于f (x,y,z)0的區(qū)域?qū)τ谇?S 具有正的指向;而對于f (x,y,z)0的區(qū)域?qū)τ谇?S

19、具有負(fù)的指向。正指向的區(qū)域用+S表示,“+”號可不寫;負(fù)指向的區(qū)域用-S表示。柵元用各相關(guān)曲面的布爾運(yùn)算表示,布爾算符包括交用空格表示 )、并用冒號:表示和非用#表示)。缺省的運(yùn)算順序是先非,其次是交,最后是并 ,使用括號可控制布爾運(yùn)算的次序。非運(yùn)算有兩種形式:#n,n是某個柵元號, #n表示一個由不在柵元n內(nèi)的點組成的空間區(qū)域。 # ( -),括號內(nèi)是對某一個柵元進(jìn)行描述的曲面柵元關(guān)系組,這一形式定義的幾何區(qū)域由不屬于括號內(nèi)描述區(qū)域的點組成的空間。例如:3 0 -1 2 -4$ 定義柵元3 #3$ 與下行相同 #(-1 2 -4) 在柵元卡上可定義柵元參數(shù)以代替在輸入文件中數(shù)據(jù)卡部分定義的柵

20、元參數(shù)。格式為:關(guān)鍵詞值。這兒允許的關(guān)鍵詞是:帶有粒子標(biāo)識符的IMP、VOL、PWT、EXT、FCL、WWN、DXC、NONU、PD和TMP,以及關(guān)于重復(fù)結(jié)構(gòu)的4個柵元參數(shù)卡:U卡、TRCL卡、LAT卡和FILL卡。 在LIKE n BUT格式中,還有兩個關(guān)鍵詞MAT和RHO,分別表示柵元的介質(zhì)號和密度。例如:10 16 -4.2 1 -2 3 IMP:N=4 IMP:P=8表示柵元10由曲面1的正面、曲面2的負(fù)面和 曲面3的正面的交集組成,填充質(zhì)量密度為 4.2 克 / 厘米3 的16號材料。該柵元的中子重要性為4,光子重要性為8。例如:2 3 -3.7 -1 IMP:N=2 IMP:P=4

21、 3 LIKE 2 BUT TRCL=1 IMP:N=10l 曲面描述卡曲面描述卡l 由方程定義曲面由方程定義曲面格式: jn a listj曲面號,1 j 99999,寫在第 15 列上。如果曲面號前有*號,則該曲面為反射面。n對應(yīng)坐標(biāo)變換卡TRn,表示該曲面是在輔助坐標(biāo)系下描述的,而該輔助坐標(biāo)系與基本坐標(biāo)系之間的關(guān)系由TRn卡給出。如果沒有坐標(biāo)變換,即曲面是在基本坐標(biāo)系下描述的,則該項缺省。a曲面助記符。list曲面方程參數(shù),110項,取決于曲面類型。參見MCNP手冊,表3.1。l 用點定義軸對稱曲面用點定義軸對稱曲面l 類型為類型為X、Y或或Z的曲面卡是用坐標(biāo)點描述的曲面卡是用坐標(biāo)點描述

22、曲面而不是用方程系數(shù)描述。用這些卡描述的曲面必曲面而不是用方程系數(shù)描述。用這些卡描述的曲面必須是分別關(guān)于須是分別關(guān)于X、Y或或Z軸對稱的,并且如果該曲面是軸對稱的,并且如果該曲面是由多葉組成的,則指定的坐標(biāo)點必須全都在同一個葉由多葉組成的,則指定的坐標(biāo)點必須全都在同一個葉上。上。格式: jn a listj曲面號,1 j 99999,寫在第 15 列上。nTRn卡的號,如果沒有坐標(biāo)變換,則該項缺省。a字母X、Y或Z。list13對點的坐標(biāo)。 每一對坐標(biāo)點定義這個曲面上的一個點。例如在一張Y卡上可以給出:j Y y1 r1 y2 r2其中,( ) 是第 i 點的坐標(biāo)。給出的坐標(biāo)點對數(shù)的不同,描述

23、的曲面類型也不同。給出一對坐標(biāo),則定義一個平面PX、PY或PZ)。給出二對坐標(biāo),則定義的是線性曲面PX、PY、PZ、CX、CY、CZ、KX、KY或KZ)。給出三對坐標(biāo),則定義的是二次曲面PX、PY、PZ、SO、SX、SY、SZ、CX、CY、CZ、KX、KY、KZ或SQ)。當(dāng)用兩點定義一個錐面時,只生成一個單葉錐面。曲面的指向與方程指定曲面SQ除外是一樣的。對于SQ,遠(yuǎn)離對稱軸的點具有正指向。而方程定義的SQ可以自由選取指向。22,iiiizxryl 由三個點定義一般平面由三個點定義一般平面lMCNP對用戶指定的對用戶指定的 P 型曲面,型曲面, 將檢查所給的數(shù)將檢查所給的數(shù)據(jù)個數(shù),若是據(jù)個數(shù),

24、若是 4項,則作一般斜置平面方程的系數(shù)理項,則作一般斜置平面方程的系數(shù)理解,若多于解,若多于 4 項時,便作為三維空間點的坐標(biāo)值理解。項時,便作為三維空間點的坐標(biāo)值理解。每三個數(shù)定義空間一個點,每三個數(shù)定義空間一個點,MCNP將把它們轉(zhuǎn)換成所將把它們轉(zhuǎn)換成所需要的曲面系數(shù)以產(chǎn)生平面:需要的曲面系數(shù)以產(chǎn)生平面:lAX+BY+GZD0格式: jn P X1 Y1 Z1 X2 Y2 Z2 X3 Y3 Z3j曲面號,1 j 99999,寫在第 15 列上。nTRn卡的號,如果沒有坐標(biāo)變換,則該項缺省。P該曲面卡的助記符。 (Xi, Yi, Zi)定義該平面的點坐標(biāo)。l 數(shù)據(jù)卡數(shù)據(jù)卡l在信息卡、柵元描述

25、卡和曲面描述卡在信息卡、柵元描述卡和曲面描述卡之后輸入的是數(shù)據(jù)卡,數(shù)據(jù)卡可分為之后輸入的是數(shù)據(jù)卡,數(shù)據(jù)卡可分為10類:類:l 問題類型問題類型l 幾何卡幾何卡l 減方差減方差l 源描述源描述l 計數(shù)描述計數(shù)描述l 材料及截面描述材料及截面描述l 能量及熱處理能量及熱處理l 問題截斷條件問題截斷條件l 用戶數(shù)據(jù)數(shù)組用戶數(shù)據(jù)數(shù)組l 外圍卡外圍卡l數(shù)據(jù)卡中,標(biāo)識符必須從前數(shù)據(jù)卡中,標(biāo)識符必須從前5列開始列開始填寫。填寫。l 問題類型問題類型MODE卡卡l如果不給出如果不給出MODE卡,則缺省形式是卡,則缺省形式是MODE N,即缺省值是中子輸運(yùn)問題。,即缺省值是中子輸運(yùn)問題。格式: MODEx1 x

26、i xi N,中子輸運(yùn)。P,光子輸運(yùn)。E,電子輸運(yùn)。l 幾何卡幾何卡l幾何卡有以下幾類:幾何卡有以下幾類:助記符卡片類型VOL柵元體積AREA曲面面積UUniverseTRCL柵元變換LAT格子FILL填充卡TR坐標(biāo)變換坐標(biāo)變換卡格式:TRn O1,O2,O3,B1,B2,B3,B4,B5,B6, B7,B8,B9,Mn變換號, 1 n 999 ,*TRn表示 Bi是角度而非角度的余弦。O1,O2,O3坐標(biāo)變換向量的位移。 B1至B9坐標(biāo)變換的坐標(biāo)旋轉(zhuǎn)矩陣。元素B1, B2,B3,B4, B5,B6,B7,B8,B9軸x,x y,x z,x x,y y,y z,y x,z y,z z,z M1

27、,表示位移是輔助坐標(biāo)系原點相對于基本坐標(biāo)系的位移。1,表示位移是基本坐標(biāo)系原點相對于輔助坐標(biāo)系的位移。l 減方差減方差lMCNP運(yùn)用以下卡片來減小方差:運(yùn)用以下卡片來減小方差:助記符卡片類型IMP柵元重要性ESPLT能量分裂和俄國輪盤賭PWT次級光子權(quán)重EXT指數(shù)變換VECT方向矢量定義FCL強(qiáng)迫碰撞助記符卡片類型WWE權(quán)重窗的能量或時間間隔WWN權(quán)重窗的邊界WWP權(quán)重窗的參數(shù)WWG權(quán)重窗生成器WWGE權(quán)重窗生成器的能量或時間間隔MESH重疊重要性網(wǎng)格權(quán)重生成器PD探測器貢獻(xiàn)DXCDXTRAN貢獻(xiàn)BBREM韌致輻射偏倚因子l 源定義源定義助記符卡片類型SDEF通用源SIn源的信息SPn源的概率

28、SBn源的偏倚DSn相關(guān)的源SCn源的注釋SSW寫曲面源SSR讀曲面源KCODE臨界源KSRC臨界計算的源起始點ACODE特征值源通用源卡:格式l 計數(shù)描述計數(shù)描述l下列卡片用來記錄計算結(jié)果:下列卡片用來記錄計算結(jié)果:助記符卡片類型Fna計數(shù)類型FCn計數(shù)注釋En計數(shù)能量間隔Tn計數(shù)時間間隔Cn計數(shù)方向余弦間隔FQn計數(shù)打印層次FMn計數(shù)乘子DEn/DFn劑量能量/劑量函數(shù)EMn計數(shù)能量乘子助記符卡片類型TMn計數(shù)時間乘子CMn計數(shù)余弦乘子CFn計數(shù)柵元標(biāo)志SFn計數(shù)曲面標(biāo)志FSn計數(shù)片段劃分SDn計數(shù)片段的體積/面積FUn子程序TALLYX輸入TFn計數(shù)漲落打印DD探測器和DXTRAN診斷D

29、XTDXTRAN球參數(shù)FTn計數(shù)特殊處理計數(shù)類型卡Fna格式:助記符類型說明Fn單位*Fn單位F1: (N、P、E)面流粒子MeVF2: (N、P、E)面通量粒子/cm2MeV/cm2F4: (N、P、E)體通量粒子/cm2MeV/cm2F5a: (N、P)點或環(huán)探測器通量粒子/cm2MeV/cm2F6: (N、P、N,P)平均沉積能量MeV/克109 J/克F7: N平均裂變沉積能量MeV/克109 J/克F8: (P、E、P,E)脈沖MeV+F8: E沉積電荷電荷無l 材料描述材料描述l這組卡片用于指定在柵元中所使用的材料這組卡片用于指定在柵元中所使用的材料成分和使用那些截面數(shù)據(jù)。成分和使用那些截面數(shù)據(jù)。助記符卡片類型Mm材料成分DRXS離散反應(yīng)

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論