Abstract
目前國內爐心物理計算普通採用有限差分法。在容許有限的計算機時間和記憶容量範圍內，無法求得準確的結果。採用節點法可解決此困難。由於近代節點法為國外先進技術，無法自國外取得。國內欲使用此法則有賴於根據學術理論，自行撰寫程式。本論文之研究目標為建立二維直角座標節點擴散程式（２DX-Y Nodal Diffusion Code）。為國內首度引進節點法技術，做紮根工作，提供一更精確、更經濟之爐心物理計算工具。本文循序完成了二維直角座標介面流多項式節點擴散法之公式詳細推導及程式撰寫等工作。其中一維中子通率可用二次至四次多項式近似，橫向滲漏可用定值或二次多項式近似。本文並針對不同型式的問題廣泛及有效地測試此程式。結果顯示有效增殖因數Keff及中子通率之計算均相當理想，與參考值如 DIF３D-NDT ，解析解或有限差分法之計算結果比較均十分接近。計算中橫向滲漏採定值或二次近似時，會造成４５度方向中子通率有明顯差異。另外本程式採用粗網再平衡加速法，以加快收斂速度，測試結果證實加速法有其必要，可大輻減少運算時間，增強實用性。由計算結果比較可顯示此多項式節點擴散法較諸有限差分法確能達到快速而準確的預期目標。綜觀本研究成果，參考有限的文獻資料，根據節點擴散法理論，推導詳細之公式，自行撰寫程式，並經測試顯示其可靠性與速度，己達到引進技術之目的。本研究所完成之程式，可做為國內未來進一步發展相關節點法之基石。