CITATION程序及使用簡介

1、CITATION程序及使用簡介目錄 CITATION程序的發(fā)展 程序模型、功能及應(yīng)用 輸入文件格式及說明 計算實例CITATION程序歷史 CITATION 美國橡樹嶺國家實驗室反應(yīng)堆部開發(fā) 最早于1969年發(fā)布 1972年第二版修正后開始被廣泛應(yīng)用 80年代在PC機上實現(xiàn)（原來在IBM360機） 1996年發(fā)布列表輸入版，提升計算能力和效率，對大型反應(yīng)堆實際問題的求解進入實用階段 是堆芯擴散計算的主要標(biāo)準(zhǔn)程序程序模型 求解中子擴散方程 可求解帶有任意能組間散射的多能群問題 應(yīng)對點、一、二、三維空間問題 可處理X-Y-Z，-R-Z，六角形-z 和三角形-z 等幾何形狀 采用顯式的有限差分 直接

2、迭代法，可用雙向行列松弛法加速程序功能 基于宏觀數(shù)據(jù)進行堆芯靜態(tài)分析，臨界計算得到有限增殖系數(shù)和通量等值或進行臨界搜索 燃耗計算和多次循環(huán)燃料管理分析 可計算共軛通量，求解一階微擾問題CITATION輸入格式 文本輸入方式 輸入數(shù)據(jù)按照類型分為不同節(jié) 每節(jié)以三位整數(shù)開頭并作為標(biāo)識符 節(jié)內(nèi)數(shù)據(jù)排布根據(jù)各自規(guī)定進行 不同輸入數(shù)據(jù)節(jié)按照大小排序 以999為結(jié)束符號 整個輸入文件以兩行標(biāo)題開頭 一個輸入文件可以計算多個問題CITATION輸入格式 參數(shù)化輸入 固定格式，每行72列 例：24I3代表24個整數(shù)，每個整數(shù)三位長度 控制類卡片需要輸入的數(shù)據(jù)數(shù)量是固定的 0一般指選用默認(rèn)值，非0即為輸入值 數(shù)

3、據(jù)類卡片數(shù)量根據(jù)實際問題確定001 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0700 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 7 7 0 0 0 0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00輸入說明 CITATION的輸入節(jié)輸入說明 001 總控制信息數(shù)據(jù)節(jié) 包括5個卡片（即下述五行），總計79個參數(shù) 卡片1：001（3

4、） 卡片2：控制選擇（243） 卡片3：編輯選擇（243） 卡片4：迭代次數(shù)和機器運行時間的限制（243） 卡片5：一般限制（612.0）001 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0700 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 7 7 0 0 0 0 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00輸入說明 001節(jié)卡片

5、2 NGC10（28-30）：本征值問題的類型 0為有效增殖因子計算 NGC19（55-57）：宏觀截面選擇 如果0，僅僅使用008 輸入宏觀截面 001節(jié)卡片3 IEDG10（28-30）：按能組給出點通量值 如果0，則給出001 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0700 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 7 7 0 0 0 0 0.00000E+00 0.00000E+00 0.00000E+00 0.

6、00000E+00 0.00000E+00 0.00000E+00輸入說明 003：中子通量問題的描述節(jié) 確定幾何形狀、邊界條件，加速方法，收斂條件等 004：幾何網(wǎng)格描述節(jié) 確定網(wǎng)格數(shù)量，區(qū)域大小等 005：區(qū)域布置描述節(jié) 指定各區(qū)域所用的材料號 008：宏觀截面數(shù)據(jù)節(jié) 給出各材料的群常數(shù)數(shù)據(jù)輸入說明 003節(jié)卡片2 NUAC5（13-15）：幾何選擇 6為二維平板（X、Y） NUAC11（31-33）：左邊界條件 1為反射邊界條件，0為外推邊界條件 卡片3：EPI1（1-12）：通量收斂標(biāo)準(zhǔn) 卡片4：XMIS3（25-36）：堆芯功率，MW003 0 0 0 0 6 0 0 0 0 0 1

7、 1 0 0 0 0 0 0 0 -1 0 0 0 0 1.00000E-04 1.00000E-05 1.00000E-04 1.00000E-04 1.00000E-04 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 0.25000E+00 0.00000E+00輸入說明 004節(jié) 數(shù)據(jù)大小不限，按照一定規(guī)則輸入XYZ各方向網(wǎng)格數(shù)量和區(qū)域大小數(shù)據(jù) X方向分成三段，分別是24cm，32cm和24cm，網(wǎng)格分別是12,16和12個 Y方向條件相同004 12 24.0 16 32.0 12 24.0 12 24.0 16

8、 32.0 12 24.0 輸入說明 005節(jié) 數(shù)據(jù)大小不限 給出XY方向上不同區(qū)域的材料 如有Z方向，則分層給出相應(yīng)數(shù)據(jù)005 3 2 3 2 1 3 3 3 3 輸入說明 008節(jié) 卡片2：能組數(shù)和散射范圍（33 ） 卡片3：截面（26，512.0） 區(qū)域、能群、擴散系數(shù)、吸收截面、中子產(chǎn)生截面、微擾數(shù)據(jù)、功率數(shù)據(jù)等 卡片4：散射截面（612.0） 系統(tǒng)聯(lián)系讀取卡片3和4，直至0數(shù)據(jù) 卡片5：裂變源分布函數(shù)（612.0）輸入說明008 2 1 0 1 1 1.400000000 0.010000000 0.007000000 0.00000E+00 0.00000E+00 0.00000E

9、+00 0.010000000 1 2 0.400000000 0.150000000 0.200000000 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2 1 1.400000000 0.010000000 0.007000000 0.00000E+00 0.00000E+00 0.00000E+00 0.010000000 2 2 0.400000000 0.150000000 0.200000000 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 3 1 1.300000000 0.00

10、8000000 0.003000000 0.00000E+00 0.00000E+00 0.00000E+00 0.010000000 3 2 0.500000000 0.050000000 0.060000000 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0 0 1.00000E+00 0.00000E+00計算實例 TWIGL二維時空動力學(xué)基準(zhǔn)題 工物系反應(yīng)堆物理實驗室次臨界裝置 靜態(tài)IAEA PWR三維基準(zhǔn)題TWIGL基準(zhǔn)題 二維兩群問題 計算模型為1/4堆芯 內(nèi)部為對稱邊界條件，外圍是真空邊界 此動態(tài)基準(zhǔn)題有線性和階躍反應(yīng)性輸入兩種

11、情況，此處僅對初始時刻做臨界計算 輸入文件：次臨界裝置 三維兩群問題 堆芯為燃料棒和有機玻璃慢化劑，四周和下部為石墨反射層 所有邊界均為真空邊界條件 具體尺寸和材料群常數(shù)等見輸入文件：三維IAEA基準(zhǔn)題 該壓水堆堆芯由177個組件構(gòu)成，計算取1/4模型CITATION使用 Windows環(huán)境下命令行程序 指定輸入文件和輸出文件即可運行：.Citation twigl.txt twigl.oTWIGL基準(zhǔn)題運行結(jié)果0.00E+005.00E-041.00E-031.50E-032.00E-032.50E-033.00E-033.50E-034.00E-034.50E-031 3 5 7 9 11 13 15 1