第八章蒙特卡羅方法應(yīng)用程序介紹(編輯修改稿)
2024-08-16 23:06 本頁面
【文章內(nèi)容簡介】 中打印致命錯誤信息，MCNP不再進行粒子輸運計算，作業(yè)中斷。 第一個出現(xiàn)的致命錯誤是真的，而后面的錯誤可能不一定是真的，這取決于前面出現(xiàn)的致命錯誤的情況。 若在 MCNP運行行上指定 FATAL項，則MCNP忽略致命錯誤，照常運行。 對于 MCNP的警告信息，用戶不應(yīng)忽視，應(yīng)搞清楚它們的含義。 (8)檢查幾何錯誤 MCNP在處理輸入文件的數(shù)據(jù)時，不能檢查一種非常重要的輸入錯誤。即 MCNP無法查出各柵元之間的重疊和空隙，只有當粒子丟失時，才會發(fā)現(xiàn)幾何錯誤。即使如此，可能仍然無法準確判斷錯誤性質(zhì)。 2) 柵元描述卡 格式： j m d geom params 或： j LIKE n BUT list j 柵元號， 1≤ j ≤99999，寫在第 1～ 5 列上。 m 柵元材料號，與材料卡（ Mm）中的序號對應(yīng)。 m ＝ 0 為真空柵元。 d 柵元材料密度。正值為原子密度，負值為質(zhì)量密度。對于真空柵元，該項缺省，不填寫。 geom 柵元的幾何說明。由一系列帶符號的曲面號經(jīng)過布爾運算組成。 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 具有負的指向。正指向的區(qū)域用 +S表示，“ +”號可不寫；負指向的區(qū)域用 S表示。柵元用各相關(guān)曲面的布爾運算表示，布爾算符包括 交 （用空格表示 ）、 并 （用冒號 :表示）和 非 （用 表示）。缺省的運算順序是先 非 ，其次是 交 ，最后是 并 ，使用括號可控制布爾運算的次序。 非 運算有兩種形式： (1)n， n是某個柵元號， n表示一個由不在柵元n內(nèi)的點組成的空間區(qū)域。 (2) ( )，括號內(nèi)是對某一個柵元進行描述的曲面 —— 柵元關(guān)系組，這一形式定義的幾何區(qū)域由不屬于括號內(nèi)描述區(qū)域的點組成的空間。 例如： 3 0 1 2 4 $ 定義柵元 3 3 $ 與下行相同 (1 2 4) 在柵元卡上可定義柵元參數(shù)以代替在輸入文件中數(shù)據(jù)卡部分定義的柵元參數(shù)。格式為：關(guān)鍵詞＝值。這兒允許的關(guān)鍵詞是：帶有粒子標識符的 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 1 2 3 IMP:N=4 IMP:P=8 表示柵元 10由曲面 1的正面、曲面 2的負面和 曲面 3的正面的交集組成，填充質(zhì)量密度為 克 / 厘米 3 的 16號材料。該柵元的中子重要性為 4，光子重要性為 8。 例如： 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 列上。 如果曲面號前有 *號，則該曲面為反射面。 n 對應(yīng)坐標變換卡 TRn，表示該曲面是在輔助坐標系下描述的，而該輔助坐標系與基本坐標系之間的關(guān)系由 TRn卡給出。如果沒有坐標變換，即曲面是在基本坐標系下描述的，則該項缺省。 a 曲面助記符。 list 曲面方程參數(shù)， 1～ 10項，取決于曲面類型。 參見 MCNP手冊，表 。 (2) 用點定義軸對稱曲面 類型為 X、 Y或 Z的曲面卡是用坐標點描述曲面而不是用方程系數(shù)描述。用這些卡描述的曲面必須是分別關(guān)于 X、 Y或 Z軸對稱的，并且如果該曲面是由多葉組成的，則指定的坐標點必須全都在同一個葉上。 格式： j n a list j 曲面號， 1≤ j
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1