MCNP以及使用方法

1、MCNP以及使用方法總 目 錄MC方法解粒子輸運問題MCNP及使用粒子輸運問題粒子輸運系統(tǒng)輸運結(jié)束條件？所求量粒子源輸運長度粒子和物質(zhì)作用的核素、作用類型、作用后粒子能量和方向MC方法解粒子輸運問題逐一模擬每個粒子在輸運系統(tǒng)中的輸運過程，模擬足夠多的粒子（相對統(tǒng)計誤差可接受） （MC原理：大量、隨機(jī)性實驗）計算粒子對所求物理量的平均貢獻(xiàn)（ MC原理：平均結(jié)果近似問題的解） 一般情況下，當(dāng)模擬足夠多粒子后，認(rèn)為所模擬的粒子的平均行為近似為輸運系統(tǒng)中粒子的平均行為；粒子對所求物理量的平均貢獻(xiàn)即為問題的蒙卡解。開始數(shù)據(jù)預(yù)處理, 各記錄單元清零取一個粒子歷史源分布抽樣輸運過程碰撞過程歷史終止否？統(tǒng)計處

2、理做完給定歷史數(shù)否？結(jié)果的處理與輸出終止記錄過程記錄過程記錄過程記錄過程MC輸運流程圖 根據(jù)問題的幾何條件、物理假定、處理方法，可歸納為以下幾種：粒子從系統(tǒng)逃脫；粒子經(jīng)碰撞被吸收；經(jīng)俄國輪盤賭后，歷史被終止；粒子能量低于給定能量（閾能）；粒子位置越過某一界面；粒子飛行時間超過給定時間；粒子權(quán)重小于某個小量。粒子歷史終止條件總 目 錄MC方法解粒子輸運問題MCNP及使用MCNP目錄MCNP簡介MCNP運行MCNP輸入文件編寫舉例MCNPA General Monte Carlo Neutron and Photon Transport Code開發(fā)： 美國洛斯阿拉莫斯國家實驗室（ Los Ala

3、mos National Laboratory ）的蒙特卡羅小組MCNP發(fā)展歷史 MC方法研究（1948年開始）MCS （1963年）MCN（1965年）中子，三維幾何，具有獨立的數(shù)據(jù)庫MCNG=MCN+MCG （1973）模擬高能光子MCNP=MCNG+MCP （1977年）模擬低能光子（能量達(dá)1KeV)Stan Ulam, John,von Neumann Nicholas MetropolisM.C 方法的命名It was Nicholas Metropolis who noted that Stan had an uncle who would borrow money from re

4、latives because he “just had to go to Monte Carlo” and thus named the mathematical method “Monte Carlo.”MCNP發(fā)展歷史MCNP3 （1983年）MCNP3A （1986年）MCNP3B （1988年）MCNP4 （1990年）MCNP4A（1993年）MCNP4B （1997年）（LAHET+HMCNP）MCNP4C （2000年）MCNPX：質(zhì)子、中子、光子、電子耦合（UNIX 第一個版本，具有并行功能，可以計算電子）（記數(shù)分析圖）可以計算微擾可以更好處理不可分辨共振MCNP5(2003

5、年)(ANSI standard Fortran 77)(ANSI standard Fortran 90)MCNP可處理問題neutron onlyphoton onlyelectron onlycombined neutron/photon(n-p)neutron/photon/electronphoton/electronelectron/photonPhotonuclear （mcnp4c不能處理光核反應(yīng)）Keff（裂變系統(tǒng)）MCNP粒子能量限制neutron：10E-11 MeV to 20 MeV(all isotope)up to 150 MeV(some isotope)pho

6、ton: 1keV to 100GeVelectron: 1keV to 1GeV工作方式Monte Carlo程序的工作方式：用戶代碼型：EGS，GEANT，MARS輸入文件型：MCNP，F(xiàn)LUKAMCNP屬于第二種源碼開放，用戶可以自己修改程序，重新編譯MCNP做題三步驟：已知、求/證明；解；答。MCNP（求解）輸入文件：告訴MC程序已知什么、求什么輸出文件：告訴我們答案是什么輸出文件介紹 程序運行完之后，主要生成以下幾類輸出文件：標(biāo)準(zhǔn)輸出（out*）用于接續(xù)運行(runtp*)統(tǒng)計結(jié)果輸出（mctal）mesh統(tǒng)計結(jié)果輸出（meshtal)輸出文件out*輸入文件重復(fù)區(qū)域、截面、物理設(shè)置

7、、動態(tài)內(nèi)存信息輸運過程信息輸運完成后的結(jié)果綜述tally卡結(jié)果詳細(xì)信息tally卡結(jié)果綜合信息輸運粒子數(shù)、時間信息輸出文件runtp*主要用于接續(xù)運行運行命令為：mcnp5 c i=*在*文件中第一行為：continue第二行給出運行歷史（粒子數(shù)或運行時間）注意：如果用nps，則nps代表總粒子數(shù)如果用ctme，則ctme代表繼續(xù)運行時間MCNP目錄MCNP簡介MCNP運行MCNP輸入文件編寫舉例MCNP運行運行平臺Windows平臺下運行Linux平臺下運行常用指令：mcnp5 i=* o=* r=*mcnp5 i=* (輸出文件out*，runtp* 自動從e字母開始排)mcnp5 i=*

8、 r=runtp* c （接續(xù)運行）mcnp5 i=* c (默認(rèn)為runtpe) （接續(xù)運行）對于Mcnp4c版本，mcnp5mcnp，其余不變Windows平臺下運行Linux平臺下運行登錄服務(wù)器 user：flukau password：flukau軟件：SSHSecureShellClientSSH Secure SSH Secure Shell ClientLinux相關(guān)命令網(wǎng)上搜。運行SSH Secure 使用方法同ftp添加profile登陸密碼同用戶名：flukau鼠標(biāo)右鍵，添加文件夾，用自己的名字命名；以后上傳、下載文件以及運行程序都在自己的文件夾內(nèi)進(jìn)行。按圖示打開輸入運行命

9、令的窗口；也可以雙擊桌面的SSH Secure Shell Client，但后者要再次登錄。進(jìn)入自己的文件夾手動刷新，顯示新生成的文件MCNP目錄MCNP簡介MCNP運行MCNP輸入文件編寫舉例粒子輸運問題的“已知” 和 “求”已知：輸運系統(tǒng)已知：輸運方式求：？粒子類型位置方向能量幾何描述介質(zhì)材料輸運粒子類型粒子死亡條件運算終止條件減方差技巧統(tǒng)計量統(tǒng)計方式已知：粒子源輸入文件輸入文件采用卡片結(jié)構(gòu)所謂的“卡片”就是一行輸運粒子類型卡幾何相關(guān)卡源定義卡減方差技巧卡統(tǒng)計卡材料定義卡物理過程卡問題截止卡用戶數(shù)組卡其它卡信息塊 可選空行分隔符可選標(biāo)題卡柵元卡 空行分隔符！面卡 空行分隔符！數(shù)據(jù)卡 空行終

10、止符 推薦使用其它 可選輸入文件舉例行輸入格式書寫規(guī)則柵元卡和面卡之間必須空一行!面卡和數(shù)據(jù)卡之間必須空一行！每行不超過80列；所有卡都必須從第15列開始；如果15列為空白，則表示它是前一張卡的繼續(xù)卡； 而至少一個空格以后的符號&表示隨后的卡是繼續(xù)卡數(shù)據(jù)項由一個或多個空格分開英文大小寫都可書寫規(guī)則注釋為c（行注釋符）或者$（末尾注釋符）nR代表重復(fù) 1 3R = 1 1 1 1nI代表遞加 1 2I 4 = 1 2 3 4xM代表遞乘 1 3M 3M = 1 3 9nJ代表跳過（默認(rèn)） 比如：DBCN X1 X20 DBCN 15J 8 3J 1(只給第16和第20個參數(shù)賦值，其它默認(rèn)) MC

11、NP輸入文件中物理量的單位長度厘米能量MeV時間10-8 秒溫度MeV（kT）原子密度1024 個原子 / 厘米3質(zhì)量密度克 / 厘米3截面10-24 厘米2原子量中子質(zhì)量的1.008664967倍阿伏加德羅常數(shù)6.0231023 幾何描述類似于搭積木（面、塊搭成）面的布爾運算(交、并、余) 柵元長方體 六個平面交圓柱體 一個圓柱面和兩個平面交球體 一個球面不規(guī)則體：比如，球內(nèi)挖去一個小圓柱體剩余的部分面卡可定義的面包括：平面、球面、橢球面、雙曲面、拋物面、柱面、錐面、圓環(huán)面四種形式直接用面方程進(jìn)行定義（面編號 關(guān)鍵字 參數(shù)）使用點定義軸對稱面使用三點確定普通平面使用小體積元定義面BOX RP

12、P SPH RCC面方程定義曲面卡助記符類型說 明方 程卡片上的數(shù)據(jù)項PPXPYPZ平面一般垂直X軸垂直Y軸垂直Z軸SOSSXSYSZ球面球心在原點一般球心在X軸球心在Y軸球心在Z軸C/XC/YC/ZCXCYCZ圓柱面平行于X軸平行于Y軸平行于Z軸軸心在X軸軸心在Y軸軸心在Z軸面方程定義曲面卡（續(xù)）助記符類型說 明方 程卡片上的數(shù)據(jù)項K/XK/YK/ZKXKYKZ圓錐面平行于X軸平行于Y軸平行于Z軸軸心在X軸軸心在Y軸軸心在Z軸1只用于單葉錐面SQ橢球面雙曲面拋物面軸平行于X、Y或Z軸GQ圓柱面圓錐面橢球面雙曲面拋物面軸不平行于X、Y或Z軸助記符類型方 程卡片上的數(shù)據(jù)項TXTYTZ橢圓或圓形的

13、圓環(huán)面對稱軸平行于X、Y或Z軸面方程定義曲面卡（續(xù)）曲面的布爾運算運算符：交：空格并：冒號:余：井號#曲面方向： 假定曲面 S 的曲面方程為 f (x,y,z)0 對于f (x,y,z)0的區(qū)域，為+S， “+”號可不寫； 對于f (x,y,z)0的區(qū)域，為-S。小體積元方向： 體積元內(nèi)為-，外為+幾何描述舉例長方體：20*5*10面卡：1 px 02 px 5 3 py 0 4 py 205 pz 0 6 pz 10長方體內(nèi)：1 -2 3 -4 5 -6長方體外：-1:2:-3:4:-5:6xyz圓柱體：半徑5底面圓心(0,25,0)頂面圓心(0,45,0)面卡：1 cy 02 py 25

14、3 py 45 圓柱體內(nèi)：-1 2 -3圓柱體外：1:-2:3柵元卡幾何描述 + 材料定義柵元編號 材料號 材料密度 幾何描述 其它說明其它說明格式：關(guān)鍵詞=*，關(guān)鍵詞IMP、VOL、PWT、EXT、FCL、WWN、DXC、NONU、PD、TMP和關(guān)于重復(fù)結(jié)構(gòu)的U、TRCL、LAT、FILL。對于某些重復(fù)結(jié)構(gòu)可使用簡寫形式： 柵元編號 LIKE 柵元編號n BUT 參數(shù)列表 表示該柵元除了參數(shù)列表中列出的和柵元n不同，其它都同。 注意：柵元不能有相交的區(qū)域，不能有沒定義的區(qū)域。數(shù)據(jù)卡輸運粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計卡物理過程卡（可選）問題截止卡

15、用戶數(shù)組卡（可選）其它卡Mode卡Mode N /P /E，默認(rèn)為N注意： 輸運的粒子類型和粒子源發(fā)射的粒子類型不是一回事。幾何相關(guān)卡VOL （默認(rèn)給出所有柵元的體積）AREA （默認(rèn)給出所有曲面的面積）重復(fù)結(jié)構(gòu)卡UTRCLLATFILLTR數(shù)據(jù)卡輸運粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計卡物理過程卡（可選）問題截止卡用戶數(shù)組卡（可選）其它卡源定義卡SDEF粒子信息：類型、位置、方向、能量格式： SDEF 源變量=說明根據(jù)源變量特點，說明有三種級別的表達(dá)方式：源變量=顯值 （包括單值、默認(rèn)）源變量=Dn （該變量值有個隨機(jī)分布，由SI和SP卡聯(lián)合給出其概

16、率分布）源變量=F 另一個源變量 Dn（該變量的值依賴于另外一個變量的分布，由DS卡先說明有幾種分布）源變量變量含義缺 省CEL柵元。根據(jù)XXX，YYY，ZZZ和UUU，VVV，WWW確定。SUR曲面。0（柵元源）ERG粒子能量MeV）。14 MeVTME時間（ 10-8 秒）。0源變量（續(xù)）變量含義缺 省DIR，是VEC和方向UUU，VVV，WWW的夾角余弦。方位角總是在0360之間均勻抽取。體源：在11上均勻抽取 （各向同性源）。面源：在01上的余弦分布 f ()2 中抽取。VECDIR的參考矢量。體源：必須，各向同性除外。面源：垂直于曲面，符號為NRM。NRM曲面法向的符號。+1源變量（

17、續(xù)）變量含義缺 省POS抽樣位置的參考點。0,0,0RAD抽樣位置離開POS或AXS的徑向距離。0EXT柵元源：沿AXS方向離開POS的距離。曲面源：離開AXS的角度余弦。0AXS用于EXT和RAD的參考矢量。沒有方向X抽樣位置的X坐標(biāo)。無XY抽樣位置的Y坐標(biāo)。無YZ抽樣位置的Z坐標(biāo)。無Z源變量（續(xù)）變量含義缺 省CCCCookie-cutter柵元。無ARA曲面的面積（當(dāng)平面源對點探測器要求直接貢獻(xiàn)時才需要此卡）。無WGT源粒子權(quán)重。1EFF位置選舍抽樣效率判據(jù)。0.01PAR源發(fā)射的粒子的類型。1中子2光子3電子1：包括中子的問題。2：無中子但有光子問題。3：僅有電子的問題。SI卡給出源變

18、量分布區(qū)間的邊界值。格式： SIn選項 I1 Ik n 分布號（n 1999）。選項 Ii 的說明，允許的值是： 省略或H：直方圖分布的區(qū)間邊界，僅用于標(biāo)量。L：離散的源變量值。A：定義概率分布密度的點。S：分布號。 I1 Ik 源變量值或分布號。SI源信息卡SP卡給出概率分布，微分或積分。格式： SPn選項 P1 Pk 或： SPn f a b n 分布號（n 1999）。選項 Pi 的說明，允許的值是： 省略：對H或L分布與D相同，對A分布為概率密度。D： H或L分布的各個區(qū)間的概率值(不用歸一)。C： H或L分布的各個區(qū)間的累積概率值。V：僅用于柵元分布。概率與柵元體積成 正比（若有Pi

19、 ，要乘以Pi ）。P1 Pk 源變量概率。 f 內(nèi)部函數(shù)的標(biāo)識符。a，b 內(nèi)部函數(shù)的輸入?yún)?shù)。SP源概率卡舉例舉例1：位于原點的3MeV點光源，在錐角=120的圓錐內(nèi)均勻發(fā)射粒子，圓錐軸在z軸上。 sdef pos=0 0 0 erg=3 vec=0 0 1 dir=d1 par=2 si1 0.5 1 sp1 0 1舉例2：圓心在原點的球殼源，內(nèi)徑5cm，外徑15cm，各項同性發(fā)射，均勻分布的中子源 sdef pos=0 0 0 rad=d1 par=1 si1 5 15 sp1 -21 2 (球內(nèi)均勻分布的話，半徑符合r2規(guī)律)格式： DSn選項 J1 Jk 或： DSnTI1 J1 I

20、k Jk 或： DSnQV1 S1 Vk Sk n 分布號（n 1999）。選項 Ji 的說明，允許的值是：省略或 H：連續(xù)分布的源變量值，僅用于標(biāo)量。L：離散的源變量值。S：分布號。 T 獨立變量值后面跟相關(guān)變量的值， 這些變量必須是離散的標(biāo)量。 Ii 獨立變量的值。 Ji 相關(guān)變量的值。DS卡復(fù)合分布SDEF POS=D1 ERG FPOS D2SI1L 5 3.3 6 75 3.3 6SP1.3 .7DS2S 3 4SI3H 2 10 14SP3D 0 1 2SI4-3 a b這是兩個各向同性點源，其中有的概率在(5,3.3,6 )處，能量服從階梯分布3；而有的概率在(75,3.3,6)

21、處，此時能量服從Watt 裂變譜分布4。舉例數(shù)據(jù)卡輸運粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計卡物理過程卡（可選）問題截止卡用戶數(shù)組卡（可選）其它卡材料定義卡材料定義需要知道核素構(gòu)成及其所占比例。格式：M材料編號 核素1 比例1 核素2 比例2 其它說明核素表示方式：ZZZAAA.nnX或ZZZAAA AAA=0表示按天然元素取原子量nnX是截面庫標(biāo)識，表示從哪個截面庫取截面數(shù)據(jù)只有p或者p/e耦合輸運，且忽略光核反應(yīng)，可以使用ZZZ000來簡單描述；中子截面庫選擇比較難，最簡單方法就是只給出ZZZAAA。比例：可歸一，可不歸一；正值表示原子數(shù)比例，負(fù)值表

22、示重量份額舉例：光子在水中輸運 m1 1000 2 8000 1數(shù)據(jù)卡輸運粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計卡物理過程卡（可選）問題截止卡用戶數(shù)組卡（可選）其它卡減方差技巧卡 IMP或WWN卡之一是必須的。助記符卡片類型IMP柵元重要性ESPLT能量分裂和俄國輪盤賭PWT次級光子權(quán)重EXT指數(shù)變換VECT方向矢量定義FCL強迫碰撞助記符卡片類型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韌致輻射偏倚因

23、子格式： IMP:nx1 x2 xi xI n 中子為 N，光子為 P，電子為 E。 N, P 、 P, E 或 N, P , E 也是允許的，如果它們的值相同。xi 柵元 i 的重要性，i 1，2，II 問題中的柵元總數(shù)。 缺?。涸谝粋€MODE N P 問題中，若省略了IMP:P卡，則所有柵元的光子重要性都置 1，除非其中子重要性為零，這時其光子重要性也為零。柵元重要性卡用于輸入各個柵元的重要性。零重要性可以用來終止粒子的歷史。IMP卡數(shù)據(jù)卡輸運粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計卡物理過程卡（可選）問題截止卡用戶數(shù)組卡（可選）其它卡統(tǒng)計卡Fna 統(tǒng)

24、計類型En能量間隔Tn時間間隔Cn角度間隔FMn計數(shù)乘子DEn劑量能量間隔DFn劑量轉(zhuǎn)換系數(shù)EMn能量相乘因子CMn劑量相乘因子FTn統(tǒng)計卡特殊處理參數(shù)FMesh 網(wǎng)格計數(shù)Fna統(tǒng)計類型卡描述Fn單位*Fn單位F1:N/ P/ E面流量/穿透率particles MeVF2:N/ P/ E統(tǒng)計面通量particles/cm2MeV/cm2F4:N/ P/ E統(tǒng)計體通量particles/cm2MeV/cm2F5a:N/ P統(tǒng)計點通量或環(huán)通量particles/cm2MeV/cm2F6:N/ P體能量沉積MeV/g 109 J /gF7:N體裂變能沉積MeV/g 109 J /gF8:P E探測

25、器脈沖幅度譜pulsesMeV+F8:E電荷沉積charge注意：有些是針對面，有些是針對體。舉例舉例1：F1: N 1 3 6 T 分別計算通過曲面1、3、6的穿透率以及穿過三個曲面的總穿透率。舉例2： F34: P 9 計算柵元9內(nèi)的體通量。注意：有多個同類計數(shù)卡時，需要在卡號前加不同數(shù)字進(jìn)行區(qū)分。舉例3：F5: P 計算點( 0, 20, 15)的點通量。舉例4： *F8: P 10 計算柵元10內(nèi)的能量沉積。En能量間隔卡分能量段給出Fn的計數(shù)，否則只給出總計數(shù)。格式：En E1 Ek注意：若使用此卡，n必須與相應(yīng)的Fn的n一致。舉例：E11 .1 1 20 在四個能量間隔內(nèi)統(tǒng)計F11

26、計數(shù)：截止能量0.1 MeV；1 MeV；120 MeV；整個能量范圍，即截止能量20 MeV關(guān)于F8不論 粒子類型取何值，F(xiàn)8都將記錄光子和電子，也就是說F8:p、 F8:e、 F8:p e結(jié)果一樣。注意：選擇能量間隔的時候，應(yīng)該把零區(qū)間和(1E-5)區(qū)間包括在里面。 舉例：E80 1E-5 1E-3 1E-1 零區(qū)間用來處理由于撞出電子而產(chǎn)生負(fù)計數(shù)的情況；(1E-5) 區(qū)間內(nèi)將記錄穿過柵元但沒有沉積能量的粒子。計算出通量后，根據(jù)DF卡給出的通量-劑量轉(zhuǎn)化系數(shù)，求得劑量。格式： DEn A E1 Ek DFn B F1 Fk n 計數(shù)號。 Ei 第 i 個能量點（MeV）。 Fi 對應(yīng)第 i

27、 個能量點的劑量函數(shù)值。 A LOG或LIN，能量表按對數(shù)或線性插值。 B LOG或LIN，劑量函數(shù)表的插值方法。缺省：若沒有填寫A或B，則按對數(shù)插值處理。劑量計算DE、DF卡格式： FTn ID1 P1,1 P1,2 P1,3 ID2 P2,1 P2,2 P2,3 n 計數(shù)號。 IDi 第 i 項特殊處理的標(biāo)識符。 FRV 用于計數(shù)類型1角度間隔的參考方向。GEB 高斯能量展寬。TMC 時間卷積。INC 指定碰撞次數(shù)。ICD 指定要記錄對探測器作貢獻(xiàn)的柵元。ELC 電流計數(shù)。 Pi,j 第 i 項特殊處理的參數(shù)。FTn 計數(shù)特殊處理卡數(shù)據(jù)卡輸運粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料

28、定義卡減方差技巧卡（可選）統(tǒng)計卡物理過程卡（可選）問題截止卡用戶數(shù)組卡（可選）其它卡物理過程卡助記符卡片類型PHYS能量物理截斷TMP自由氣體模型熱溫度THTME熱時間卡MTmS(,) 材料卡用于控制MCNP的能量以及其它物理狀況。數(shù)據(jù)卡輸運粒子類型mode卡幾何相關(guān)卡（可選）源定義卡材料定義卡減方差技巧卡（可選）統(tǒng)計卡物理過程卡（可選）問題截止卡用戶數(shù)組卡（可選）其它卡問題截止卡在初始運行或接續(xù)運行的輸入文件中均可使用，用于終止粒子的歷史或中斷計算。助記符卡片類型CUT粒子截斷ELPT柵元能量截斷NPS歷史數(shù)截斷CTME計算時間截斷CUT粒子截斷卡格式：CUT:n T E WC1 WC2 S

29、WTMWC1、WC2、SWTM涉及到減方差技巧，先不講。T 指截斷時間（10-8 秒），缺省值：中子-無限大，光子和電子都是等于中子的截斷時間。E 指截斷能量（MeV），缺省值：中子，光子和電子都是。舉例：CUT:p j 0.01 3j 截斷能量為，其余缺省值。格式：NPS N N 粒子的歷史數(shù)。 該卡指定要計算的粒子歷史數(shù)，當(dāng)計算的粒子歷史數(shù)達(dá)到指定值時，MCNP將中止計算。NPS歷史數(shù)截斷卡格式：CTME x x 計算時間（分鐘）。 該卡指定該問題在計算機(jī)上花費的時間限制，超過該時間，MCNP將中止計算。 所有中止計算的條件只要達(dá)到一個就中止計算。CTME計算時間截斷卡MCNP目錄MCNP簡介MCNP運行