Solution of the Neutron Transport Equation in Two Dimensions via the Finite Difference Method


GÖKER A. İ., ERGİN E.

Bilecik Şeyh Edebali Üniversitesi Fen Bilimleri Dergisi, vol.5, no.1, pp.8-17, 2018 (Peer-Reviewed Journal) identifier

Abstract

It is of paramount importance to determine the distribution of neutrons which trigger the fission reaction in nuclear reactor cores in order to prevent the nuclear accidents. The distribution of neutrons is determined via the neutron transport equation. In this work, the neutron transport equation is solved in a 2 dimensional grid via the finite difference method by using Gaussian and Lorentzian shaped external neutron sources for various absorption cross sections. The profile of the obtained neutron flux values within the reactor core is determined. When the absorption cross section is increased, it was observed that the neutron flux decreases within the entire core for both the Gaussian and Lorentzian shaped external neutron sources. Due to the narrower profile of the Gaussian shaped external neutron source, neutron flux values decrease within the entire reactor in this case relative to the Lorentzian shaped external neutron source. The required conditions for the reactor to stay critical in light of the obtained results are discussed.
Nükleer reaktör çekirdeklerinde fisyon reaksiyonunu tetikleyen nötronların dağılımının tespit edilebilmesinükleer kazaları engellemek için kritik önem taşımaktadır. Nötronların dağılımı ise nötron taşınma denklemiyoluyla tayin edilmektedir. Bu çalışmada, nötron taşınma denklemi orijinde yer alan Gauss ve Lorentz biçimli dışnötron kaynaklarını kullanılarak 2 boyutlu bir çerçevede sonlu farklar yöntemiyle değişik emilme kesit alanları içinçözülmüştür. Elde edilen nötron akı değerlerinin reaktör çekirdeği içindeki dağılımının profili tayin edilmiştir. Emilmekesit alanı artırıldığında hem Gauss hem Lorentz şeklindeki dış nötron kaynaklarında nötron akısının çekirdeğintamamında azaldığı gözlenmiştir. Gauss şeklindeki dış nötron kaynağının daha dar profilinden dolayı bu durumda nötronakı değerleri Lorentz şeklindeki dış nötron kaynağına göre çekirdeğin tamamında azalmaktadır. Elde edilen sonuçlarışığında, reaktörün kritik seviyede kalması için gerekli şartları tartışılmıştır.