第八章蒙特卡羅方法應用程序介紹-文庫吧資料
2025-07-26 23:06本頁面
【正文】 ) 問題截斷條件 (9) 用戶數(shù)據(jù)數(shù)組 (10)外圍卡 數(shù)據(jù)卡中，標識符必須從前 5列開始填寫。 P 該曲面卡的助記符。每三個數(shù)定義空間一個點， MCNP將把它們轉(zhuǎn)換成所需要的曲面系數(shù)以產(chǎn)生平面： AX+BY+GZ－ D＝ 0 格式： j n P X1 Y1 Z1 X2 Y2 Z2 X3 Y3 Z3 j 曲面號， 1≤ j ≤99999，寫在第 1～ 5 列上。而方程定義的SQ可以自由選取指向。 曲面的指向與方程指定曲面（ SQ除外）是一樣的。 c) 給出三對坐標，則定義的是二次曲面（ PX、 PY、 PZ、SO、 SX、 SY、 SZ、 CX、 CY、 CZ、 KX、 KY、 KZ或 SQ）。 a) 給出一對坐標，則定義一個平面（ PX、 PY或 PZ）。例如在一張 Y卡上可以給出： j Y y1 r1 y2 r2 其中， ( ) 是第 i 點的坐標。 list 1～ 3對點的坐標。 n TRn卡的號，如果沒有坐標變換，則該項缺省。用這些卡描述的曲面必須是分別關(guān)于 X、 Y或 Z軸對稱的，并且如果該曲面是由多葉組成的，則指定的坐標點必須全都在同一個葉上。 參見 MCNP手冊，表 。 a 曲面助記符。 n 對應坐標變換卡 TRn，表示該曲面是在輔助坐標系下描述的，而該輔助坐標系與基本坐標系之間的關(guān)系由 TRn卡給出。 例如： 2 3 1 IMP:N=2 IMP:P=4 3 LIKE 2 BUT TRCL=1 IMP:N=10 3) 曲面描述卡 (1) 由方程定義曲面 格式： j n a list j 曲面號， 1≤ j ≤99999，寫在第 1～ 5 列上。 例如： 10 16 1 2 3 IMP:N=4 IMP:P=8 表示柵元 10由曲面 1的正面、曲面 2的負面和 曲面 3的正面的交集組成，填充質(zhì)量密度為 克 / 厘米 3 的 16號材料。這兒允許的關(guān)鍵詞是：帶有粒子標識符的 IMP、 VOL、 PWT、 EXT、 FCL、WWN、 DXC、 NONU、 PD和 TMP，以及關(guān)于重復結(jié)構(gòu)的 4個柵元參數(shù)卡： U卡、 TRCL卡、LAT卡和 FILL卡。 例如： 3 0 1 2 4 $ 定義柵元 3 3 $ 與下行相同 (1 2 4) 在柵元卡上可定義柵元參數(shù)以代替在輸入文件中數(shù)據(jù)卡部分定義的柵元參數(shù)。 非 運算有兩種形式： (1)n， n是某個柵元號， n表示一個由不在柵元n內(nèi)的點組成的空間區(qū)域。柵元用各相關(guān)曲面的布爾運算表示，布爾算符包括 交 （用空格表示 ）、 并 （用冒號 :表示）和 非 （用 表示）。假定曲面 S 的曲面方程為 f (x,y,z)＝ 0，則對于 f (x,y,z)＞ 0的區(qū)域?qū)τ谇? S 具有正的指向；而對于 f (x,y,z)＜ 0的區(qū)域?qū)τ谇? S 具有負的指向。 list 描述柵元 j和柵元 n之間差別的柵元參數(shù)。 params 任選的柵元參數(shù)說明。 geom 柵元的幾何說明。正值為原子密度，負值為質(zhì)量密度。 m ＝ 0 為真空柵元。 2) 柵元描述卡 格式： j m d geom params 或： j LIKE n BUT list j 柵元號， 1≤ j ≤99999，寫在第 1～ 5 列上。即 MCNP無法查出各柵元之間的重疊和空隙，只有當粒子丟失時，才會發(fā)現(xiàn)幾何錯誤。 對于 MCNP的警告信息，用戶不應忽視，應搞清楚它們的含義。 第一個出現(xiàn)的致命錯誤是真的，而后面的錯誤可能不一定是真的，這取決于前面出現(xiàn)的致命錯誤的情況。 例如：光子截斷卡 CUT:P 3J 表示前 3個參量使用缺省值，只改變第四項參量的值。當省略某張輸入卡時，則該卡上的全部參數(shù)均使用缺省值。 例如：中子重要性卡為 IMP:N 光子重要性卡為 IMP:P (6)缺省值 MCNP的許多輸入?yún)?shù)都有缺省值，因此用戶不需要每次都給出各個輸入?yún)⒘康闹怠＿@些卡片是：IMP、 EXT、 FCL、 WWN、 WWE
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1