Invariability of neutron flux in internal and external irradiation sites in MNSRs for both HEU and LEU cores

Similar documents
Neutronic Calculations of Ghana Research Reactor-1 LEU Core

The Effect of 99 Mo Production on the Neutronic Safety Parameters of Research Reactors

Consummation of 19.75% UO 2 Fuel Material in the Core of Nigeria Miniature Neutron Source Reactor (MNSR)

Institute of Atomic Energy POLATOM OTWOCK-SWIERK POLAND. Irradiations of HEU targets in MARIA RR for Mo-99 production. G.

Delayed Neutrons Energy Spectrum Flux Profile of Nuclear Materials in Ghana s Miniature Neutron Source Reactor Core

Operational Experience and Programmes for Optimal Utilization of the Nigeria Research Reactor-1

Estimation of Control Rods Worth for WWR-S Research Reactor Using WIMS-D4 and CITATION Codes

1. INTRODUCTION 2. EAEA EXISTING CAPABILITIES AND FACILITIES

Fuel Element Burnup Determination in HEU - LEU Mixed TRIGA Research Reactor Core

Study of Control rod worth in the TMSR

USA HTR NEUTRONIC CHARACTERIZATION OF THE SAFARI-1 MATERIAL TESTING REACTOR

Designing a 2D RZ Venture Model for Neutronic Analysis of the Nigeria Research Reactor-1 (NIRR-1)

RECH-1 RESEARCH REACTOR: PRESENT AND FUTURE APPLICATIONS

Hybrid Low-Power Research Reactor with Separable Core Concept

Effect of WIMSD4 libraries on Bushehr VVER-1000 Core Fuel Burn-up

An experimental testing of coolant flow rate and velocity in the core of Nigeria Research Reactor-1

Study on SiC Components to Improve the Neutron Economy in HTGR

MONTE CALRLO MODELLING OF VOID COEFFICIENT OF REACTIVITY EXPERIMENT

MEASUREMENT AND COMPARISON OF CONTROL ROD WORTH OF BTRR USING INHOUR EQUATION AND PERIOD REACTIVITY CONVERSION TABLE

Chem 481 Lecture Material 4/22/09

Calculation of Control Rod Worth of TRIGA Mark II Reactor Using Evaluated Nuclear Data Library JEFF-3.1.2

VALIDATION OF BATAN'S STANDARD 3-D DIFFUSION CODE, BAT AN-3DIFF, ON THE FIRST CORE OF RSG GAS. Liem Peng Hong"

Available online at ScienceDirect. Energy Procedia 71 (2015 ) 52 61

17 Neutron Life Cycle

Utilization of Egyptian Research Reactor and modes of collaboration

Nuclear Engineering and Design

NEW FUEL IN MARIA RESEARCH REACTOR, PROVIDING BETTER CONDITIONS FOR IRRADIATION IN THE FAST NEUTRON SPECTRUM.

Assessing the Effect of Fuel Burnup on Control Rod Worth for HEU and LEU Cores of Gharr-1

Nuclear Research Reactors

X. Neutron and Power Distribution

Study of boron dilution phenomenon in the core and fuel assemblies of Bushehr VVER-1000 reactor in normal operating conditions

DETAILED CORE DESIGN AND FLOW COOLANT CONDITIONS FOR NEUTRON FLUX MAXIMIZATION IN RESEARCH REACTORS

M.Cagnazzo Atominstitut, Vienna University of Technology Stadionallee 2, 1020 Wien, Austria

Shutdown Margin. Xenon-Free Xenon removes neutrons from the life-cycle. So, xenonfree is the most reactive condition.

SLOWPOKE-2 31 YEARS AND STILL GLOWING

Typical Safety Analysis Application for a Nuclear Research Reactor

ACTIVATION ANALYSIS OF DECOMISSIONING OPERATIONS FOR RESEARCH REACTORS

MA/LLFP Transmutation Experiment Options in the Future Monju Core

Advanced Heavy Water Reactor. Amit Thakur Reactor Physics Design Division Bhabha Atomic Research Centre, INDIA

VERIFICATION OF A REACTOR PHYSICS CALCULATION SCHEME FOR THE CROCUS REACTOR. Paul Scherrer Institut (PSI) CH-5232 Villigen-PSI 2

but mostly as the result of the beta decay of its precursor 135 I (which has a half-life of hours).

JOYO MK-III Performance Test at Low Power and Its Analysis

Validation of Rattlesnake, BISON and RELAP-7 with TREAT

OECD/NEA Transient Benchmark Analysis with PARCS - THERMIX

Utilization of the Irradiation Holes in HANARO

Control of the fission chain reaction

DESIGN OF B 4 C BURNABLE PARTICLES MIXED IN LEU FUEL FOR HTRS

SUB-CHAPTER D.1. SUMMARY DESCRIPTION

THERMAL HYDRAULIC AND SAFETY ANALYSIS OF HEAT TRANSFER AND DISTRIBUTION IN THE GHANA RESEARCH REACTOR-1 (GHARR-1) CORE USING STAR-CCM+ CFD CODE

The investigations in the field of advanced nuclear power systems for energy production and transmutation of RAW in NAS of Belarus

NUCLEAR SCIENCE ACAD BASIC CURRICULUM CHAPTER 5 NEUTRON LIFE CYCLE STUDENT TEXT REV 2. L th. L f U-235 FUEL MODERATOR START CYCLE HERE THERMAL NEUTRON

Belarus activity in ADS field

DOE NNSA B&W Y-12, LLC Argonne National Lab University of Missouri INR Pitesti. IAEA Consultancy Meeting Vienna, August 24-27, 2010

Sensitivity and Uncertainty Analysis Methodologies for Fast Reactor Physics and Design at JAEA

A Method For the Burnup Analysis of Power Reactors in Equilibrium Operation Cycles

Development of 3D Space Time Kinetics Model for Coupled Neutron Kinetics and Thermal hydraulics

EXPERIMENTAL DETERMINATION OF NEUTRONIC PARAMETERS IN THE IPR-R1 TRIGA REACTOR CORE

Lectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 4. Title: Control Rods and Sub-critical Systems

REFLECTOR FEEDBACK COEFFICIENT: NEUTRONIC CONSIDERATIONS IN MTR-TYPE RESEARCH REACTORS ABSTRACT

MIT: Akshay Dave, Lin-wen Hu, Kaichao Sun INL: Ryan Marlow, Paul Murray, Joseph Nielsen. Nuclear Reactor Laboratory

Reactivity Power and Temperature Coefficients Determination of the TRR

Reactor-physical calculations using an MCAM based MCNP model of the Training Reactor of Budapest University of Technology and Economics

VI. Chain Reaction. Two basic requirements must be filled in order to produce power in a reactor:

Critical Experiment Analyses by CHAPLET-3D Code in Two- and Three-Dimensional Core Models

Nuclear Research Facilities in Russia for innovative nuclear development N. Arkhangelskiy, ROSATOM IAEA Consultancy Meeting June 2013

NEUTRONIC CALCULATION SYSTEM FOR CANDU CORE BASED ON TRANSPORT METHODS

Reactivity Effect of Fission and Absorption

NEUTRONIC ANALYSIS STUDIES OF THE SPALLATION TARGET WINDOW FOR A GAS COOLED ADS CONCEPT.

Fundamentals of Nuclear Reactor Physics

Evaluating the Core Damage Frequency of a TRIGA Research Reactor Using Risk Assessment Tool Software

«CALCULATION OF ISOTOPE BURN-UP AND CHANGE IN EFFICIENCY OF ABSORBING ELEMENTS OF WWER-1000 CONTROL AND PROTECTION SYSTEM DURING BURN-UP».

Research Article Calculations for a BWR Lattice with Adjacent Gadolinium Pins Using the Monte Carlo Cell Code Serpent v.1.1.7

THE SLOWPOKE-2 NUCLEAR REACTOR AT THE ROYAL MILITARY COLLEGE OF CANADA: APPLICATIONS FOR THE CANADIAN ARMED FORCES

Reactivity Coefficients

arxiv: v2 [physics.ins-det] 22 Nov 2017

Lectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 5. Title: Reactor Kinetics and Reactor Operation

Analysis of the TRIGA Reactor Benchmarks with TRIPOLI 4.4

Development of education and training programs using ISIS research reactor

Fuel BurnupCalculations and Uncertainties

Induced photonuclear interaction by Rhodotron-TT MeV electron beam

THE NEXT GENERATION WIMS LATTICE CODE : WIMS9

FUEL CYCLE ANALYSIS IN A THORIUM FUELED REACTOR USING BIDIRECTIONAL FUEL MOVEMENT

REACTOR PHYSICS FOR NON-NUCLEAR ENGINEERS

Lesson 14: Reactivity Variations and Control

CALCULATION OF TEMPERATURE REACTIVITY COEFFICIENTS IN KRITZ-2 CRITICAL EXPERIMENTS USING WIMS ABSTRACT

ESTIMATION OF THE DOSE RATE OF NUCLEAR FUEL OF GHANA RESEARCH REACTOR-1 (GHARR-1) A DISSERTATION SUBMITTED TO THE DEPARTMENT OF NUCLEAR

STEK experiment Opportunity for Validation of Fission Products Nuclear Data

Profile SFR-64 BFS-2. RUSSIA

22.06 ENGINEERING OF NUCLEAR SYSTEMS OPEN BOOK FINAL EXAM 3 HOURS

Solving Bateman Equation for Xenon Transient Analysis Using Numerical Methods

Verification of Predicted Powers and Shim Positions for Advanced Test Reactor Cycle 159A-1

Study of the End Flux Peaking for the CANDU Fuel Bundle Types by Transport Methods

Industrial and Commercial Applications of FRM II

REACTOR PHYSICS ASPECTS OF PLUTONIUM RECYCLING IN PWRs

MONTE CARLO SIMULATION OF THE GREEK RESEARCH REACTOR NEUTRON IRRADIATION POSITIONS USING MCNP

Development of sensitivity analysis with the use of TSUNAMI-3D sensitivity coefficients for benchmark uncertainty analysis

Available online at ScienceDirect. Energy Procedia 71 (2015 )

Use of Monte Carlo and Deterministic Codes for Calculation of Plutonium Radial Distribution in a Fuel Cell

CRITICAL AND SUBCRITICAL EXPERIMENTS USING THE TRAINING NUCLEAR REACTOR OF THE BUDAPEST UNIVERSITY OF TECHNOLOGY AND ECONOMICS

Neutronic Issues and Ways to Resolve Them. P.A. Fomichenko National Research Center Kurchatov Institute Yu.P. Sukharev JSC Afrikantov OKBM,

Transcription:

Indian Journal of Pure & Applied Physics Vol. 49, February 2011, pp. 83-90 Invariability of neutron flux in internal and external irradiation sites in MNSRs for both HEU and LEU cores M Albarhoum Department of Nuclear Engineering, Atomic Energy Commission, P. O. Box, 6091, Damascus- Syria Received 24 August 2010; revised 3 January 2011; accepted 11 January 2011 The effect of the control rod position on the values of the neutronic flux in the internal and external irradiation sites for the actual high enriched uranium UAl 4 -Al fuel and U 3 Si 2, U 3 Si, UO 2, and U-9Mo low enriched uranium fuels has been investigated by using a 3-D model for MNSRs. Results show that the values of the neutron flux in these sites depend only on the radial distance of these sites from the control rod axis, but not on the control rod vertical position. The standard deviation and the variance coefficient are all less than 0.0088E12 and 0.812%, respectively except for the UO 2 fuel for which these two parameters are 0.0088 and 1.938%, respectively. These calculated values depend on the fuel type and enrichment for the same power of the core and not on the thickness of the annulus reflector. The invariability of the neutronic flux values with the vertical position of the control rod renders, MNSRs optimum tools for neutron activation analysis. Keywords: Fuel, MNSR, Initial excess reactivity, Control rod, LEU 1 Introduction Miniature Neutron Source Reactors (MNSRs) are research reactors of the tank-in-pool type using Highly Enriched Uranium (HEU) fuel (~90% w/o U 235 ). Their nominal power is about 30 kw and have a fuel cage formed of 10 circles and one central control rod for the regulation and shutdown of the reactor 1,2. Their primary use is the Neutronic Activation Analysis (NAA). This activity uses 10 sites for sample irradiation, five of which are called Internal Irradiation Sites (IIS) (or inner sites) and located inside the Annulus Reflector (AR) and five are called External Irradiation Sites (EIS) (or outer sites), and located outside the annulus reflector as shown in Fig. 1. Maximum values for the thermal neutron fluxes in these sites are: 1.E12 n/cm 2.s, and 0.5E12 n/cm 2.s, respectively at rated power 3. A noticeable physical characteristic of MNSRs is that the increase in temperature and the variation of the Control Rod (CR) position during reactor operation do not affect appreciably the neutron flux distribution in the whole reactor, IIS and EIS included, except in the space surrounding the CR position itself 5. This property (the invariability of the neutron flux in the IIS and EIS) allows to MNSRs to be a stable and reliable tool for NAA. Now since the distribution of the neutron flux in the whole reactor is not so important for the NAA, but only in the IIS and EIS where samples are irradiated, it may be a good question whether this property of MNSRs depend on the fuel type, or on its enrichment (where four low enriched uranium fuels were examines) from one side, and whether this depends on the thickness of reflector or even the presence of the IIS within the AR or not from the other side. This characteristic has been investigated for the aforementioned 4 fuel types namely: U 3 Si 2, U 3 Si, U-9Mo, and UO 2. 2 Methodology The neutronics calculations were performed using a detailed three-dimensional model implemented in the core calculation s code 1-4 CITATION coupled with the lattice cell code WIMSD4 through the BMAC package. The BMAC package was used to generate the Group Constants (GC) through the WIMSD4 code and then it was used to transfer the calculated GC to CITATION code which was used to calculate the neutronic flux and power distributions in the reactor. Four neutronic groups were used for the model. Their upper energy limits were: 10 MeV, 0.821 MeV, 5530 ev, and 0.625 ev, respectively. The GC for the 1 st fuel circle are presented in Table 1 for the 4 fuels. The criticality calculations for one configuration of the reactor were performed and the Initial Excess reactivity (IER) was determined. The position of the CR was then varied vertically along its axis for the whole CR path (23 cm), and the corresponding flux values in the IIS and EIS were found. This was repeated for all LEU fuel types in addition to the

84 INDIAN J PURE & APPL PHYS, VOL 49, FEBRUARY 2011 Fig. 1 Location of inner and outer irradiation sites of PARR-2(Ref.6) Table 1 Calculated group constants for the reflectors of the Syrian MNSR UAl 4 -Al Group 1 Group 2 Group 3 Group 4 Diffusion (cm) 2.419080 0.920384 0.664506 0.202402 Absorption (cm 1 ) 0.835748 10 3 0.533008 10 3 0.918058 10 2 0.102425 U 3 Si 2 Group 1 Group 2 Group 3 Group 4 Diffusion (cm) 2.392110 0.910318 0.633522 0.202627 Absorption (cm 1 ) 0.153834 10 2 0.970878 10 3 0.150914 10 1 0.119892 U 3 Si Group 1 Group 2 Group 3 Group 4 Diffusion (cm) 1.471260 0.603615 1.121810 5.531540 Absorption (cm 1 ) 0.663879 10 2 0.567219 10 6 0.555804 10 4 0.113135 10 2 U-9Mo Group 1 Group 2 Group 3 Group 4 Diffusion (cm) 1.471260 0.603615 1.121810 5.531540 Absorption (cm 1 ) 0.663879 10 2 0.567219 10 6 0.555804 10 4 0.113135 10 2 UO 2 Group 1 Group 2 Group 3 Group 4 Diffusion (cm) 1.471260 0.603615 1.121810 5.531540 Absorption (cm 1 ) 0.663879 10 2 0.567219 10 6 0.555804 10 4 0.113135 10 2 original HEU fuel and for all CR positions. Some statistical parameters (such as the average value of the thermal neutron flux (ANF), the Standard Deviation (STDEV) and the Variance Coefficient (VC)) indicating the variability of the obtained data around their average value were then calculated for each group and compared so as to draw conclusions about the argument. All calculations were made for the nominal power of 30 kw. Calculations were also made for the reflector thickness, which was in turn diminished and new Fuel Rods (FR) were added so that the IIS became outside the AR and in the core of the reactor. The IIS were thus transferred from the AR to the core of the reactor to see the effect of the thickness of the AR on the values of the ANF in the IIS and EIS.

ALBARHOUM: INVARIABILITY OF NEUTRON FLUX IN MNSRs 85 3 Results and Discussion 3.1 UAl 4 -Al Core The thermal (4 th group) neutron fluxes in both the IIS and EIS for the actual UAl 4 -Al HEU fuel are presented in Table 2 and shown in Fig. 2. The average Table 2 Calculated thermal neutron flux values in the IIS and EIS at different CR positions for the UAl 4 -Al fuel CR position Flux (IIS) Flux (EIS) (cm) (n/cm 2.s) (n/cm 2.s) 1.0 1.012500E12 0.479664 E12 2.0 1.01286 E12 0.479684 E12 3.0 1.01337 E12 0.479701 E12 4.0 1.01290 E12 0.479653 E12 5.0 1.01297 E12 0.479543 E12 6.0 1.01324 E12 0.479705 E12 7.0 1.01276 E12 0.479073 E12 8.0 1.01178 E12 0.478805 E12 9.0 1.01068 E12 0.478332 E12 10.0 1.01051 E12 0.478225 E12 11.0 1.01081 E12 0.478055 E12 12.0 1.00924 E12 0.477604 E12 13.0 1.01004 E12 0.477677 E12 14.0 1.00929 E12 0.477431 E12 15.0 1.00868 E12 0.476088 E12 16.0 1.00760 E12 0.476645 E12 17.0 1.00713 E12 0.476385 E12 18.0 1.00686 E12 0.476088 E12 19.0 1.00603 E12 0.475710 E12 20.0 1.00558 E12 0.475578 E12 21.0 1.00527 E12 0.475370 E12 22.0 1.00477 E12 0.475216 E12 23.0 1.00437 E12 0.475096 E12 value of the thermal Neutron Flux (ANF) at the IIS (the sum of all flux values in column 2 of Table 2 divided by their number which is 23) was 1.009532E12 n/cm 2.s and the STDEV of its values was equal to 0.003054E12, and the Variance Coefficient (VC) was 0.301%, while the ANF, STDEV, and the VC at the EIS were equal to 0.477623E12 n/cm 2.s, 0.001675 E12, and 0.350%, respectively. Hence, the flux value in both the IIS and EIS is practically constant since both the STDEV and the VC are very small. The values of the thermal neutron flux for the case of the U 3 Si 2 and U 3 Si LEU fuels were presented in Table 3. 3.2 U 3 Si 2 LEU Core From Table 3, the values of ANF, STDEV and VC for the U 3 Si 2 LEU fuel in the IIS were: 0.943682E12 n/cm 2.s, 0.002640 E12 and 0.280%, respectively. The same parameters for the EIS were: 0.457407 E12 n/cm 2.s, 0.001459E12 and 0.319%, respectively. The case of the U 3 Si 2 LEU fuel was also similar to the UAl 4 -Al HEU fuel. The obtained values for the ANF, STDEV and VC indicate that for this fuel the reactor property was valid as well. 3.3 U 3 Si LEU Core For the case of the U 3 Si LEU fuel, the values of ANF, STDEV and VC in the IIS were: 0.954848E12 n/cm 2.s, 0.002757E12 and 0.288%, respectively (Table 3). The same calculated parameters for the EIS were: 0.462629E12 n/cm 2.s, 0.001515E12 and Fig. 2 Distribution of the thermal neutron flux in the IIS and EIS versus CR position in the core for the UAl 4 -Al HEU fuel

86 INDIAN J PURE & APPL PHYS, VOL 49, FEBRUARY 2011 Table 3 Thermal neutron flux values in the IIS and EIS at different CR positions for the U 3 Si 2 and U 3 Si LEU fuels Fuel type U 3 Si 2 U 3 Si CR position Flux (IIS) Flux (EIS) Flux (IIS) Flux (EIS) (cm) (n/cm 2.s) (n/cm 2.s) (n/cm 2.s) (n/cm 2.s) 1.0 0.946266 E12 0.4592900 E12 0.957892 E12 0.464576 E12 2.0 0.946302 E12 0.459229 E12 0.958373 E12 0.464611 E12 3.0 0.946248 E12 0.459101 E12 0.958318 E12 0.464461 E12 4.0 0.946297 E12 0.459204 E12 0.958335 E12 0.464549 E12 5.0 0.946393 E12 0.459075 E12 0.957829 E12 0.464307 E12 6.0 0.945331 E12 0.458915 E12 0.957055 E12 0.464198 E12 7.0 0.945462 E12 0.458511 E12 0.957295 E12 0.463832 E12 8.0 0.945068 E12 0.458381 E12 0.956271 E12 0.463502 E12 9.0 0.944992 E12 0.458224 E12 0.956792 E12 0.463509 E12 10.0 0.944752 E12 0.458105 E12 0.956371 E12 0.463336 E12 11.0 0.944345 E12 0.457749 E12 0.955862 E12 0.462954 E12 12.0 0.943811 E12 0.457637 E12 0.954880 E12 0.462685 E12 13.0 0.943494 E12 0.457384 E12 0.954696 E12 0.462501 E12 14.0 0.942671 E12 0.457074 E12 0.954427 E12 0.462326 E12 15.0 0.942383 E12 0.456733 E12 0.954033 E12 0.462046 E12 16.0 0.941936 E12 0.456569 E12 0.953381 E12 0.461808 E12 17.0 0.941420 E12 0.456328 E12 0.952613 E12 0.461433 E12 18.0 0.940955 E12 0.456015 E12 0.952404 E12 0.461149 E12 19.0 0.940332 E12 0.455671 E12 0.951817 E12 0.460883 E12 20.0 0.939741 E12 0.455513 E12 0.951308 E12 0.460667 E12 21.0 0.939631 E21 0.455416 E12 0.951060 E12 0.460531 E12 22.0 0.938904 E12 0.455196 E12 0.950060 E12 0.460245 E12 23.0 0.938436 E12 0.455051 E12 0.950435 E12 0.460362 E12 0.3279%, respectively. The STDEV and VC values were so small to consider the CR position not affecting the ANF in both the IIS and EIS. 3.4 UO 2 and U-9Mo LEU Cores The case of the UO 2 LEU and U-9Mo fuels is presented in Table 4. The ANF, STDEV and VC values for the case of the UO 2 LEU fuel in the IIS were: 0.992638E12 n/cm 2.s, 0.002606E12 and 0.262%, respectively as presented in Table 4, while the same parameters for the EIS were: 0.454339E12 n/cm 2.s, 0.008804E12 and 1.93842%, respectively. For this fuel in the EIS, the VC was about 2%, while it was less than unity for all the preceding fuels. However, this value is still small not to make a significant conclusion about this fuel type. From the other side, the ANF, STDEV and VC for the case of the U-9Mo LEU fuel in the IIS were: 0.813240E12 n/cm 2.s, 0.00609476E12 and 0.812%, respectively, while the same parameters for the EIS were: 0.419870E12 n/cm 2.s, 0.001338E12 and 0.3182%, respectively. Although the absolute value of the neutron flux in both the IIS and EIS is about 20% less than that of theual 4 -Al HEU fuel, the values of the STDEV and VC are so small that this fuel can be treated like the other LEU fuel types except for the UO 2 fuel type. The results obtained for the different LEU fuels for the IIS are shown in Fig. 3. It can be concluded that the effect of the CR position on the thermal neutron flux value in the IIS and EIS does not depend on the fuel type. 3.5 Annulus Reflector Thickness effect The effect of the AR thickness on this property has been investigated. This will be done for the actual HEU fuel only. The Internal Radius (IR) of the AR was varied from 11.55 cm (the actual IR) to 19.55 cm and the values of the ANF in the IIS and EIS were calculated. The results are presented in Table 5 and shown in Fig. 4. The values of the statistical parameters, i.e. the ANF, STDEV and VC for the UAl 4 -Al HEU fuel are presented in Table 6. The ANF in the IIS was a bit varying with the internal radius of the AR, but the STDEV and the VC indicated that these flux values were really not affected by the CR position (Table 6). In Table 7, the values of the same

ALBARHOUM: INVARIABILITY OF NEUTRON FLUX IN MNSRs 87 Table 4 Thermal neutron flux values in the IIS and EIS at different control rod positions for the UO 2 and U-9Mo fuels Fuel type UO 2 U-9Mo CR position Flux (IIS) Flux (EIS) Flux (IIS) Flux (EIS) (cm) (n/cm 2.s) (n/cm 2.s) (n/cm 2.s) (n/cm 2.s) 1.0 0.995751 E12 0.457710 E12 0.814485 E12 0.421365 E12 2.0 0.995418 E12 0.457573 E12 0.841869 E12 0.421568 E12 3.0 0.995596 E12 0.457503 E12 0.814282 E12 0.421378 E12 4.0 0.995198 E12 0.457504 E12 0.814454 E12 0.421309 E12 5.0 0.994452 E12 0.457221 E12 0.814263 E12 0.421358 E12 6.0 0.995536 E12 0.457614 E12 0.814380 E12 0.421258 E12 7.0 0.995109 E12 0.457077 E12 0.813100 E12 0.420979 E12 8.0 0.995081 E12 0.457070 E12 0.8142451 E12 0.420865 E12 9.0 0.994665 E12 0.456868 E12 0.813715 E12 0.420723 E12 10.0 0.993838 E12 0.456598 E12 0.813276 E12 0.420468 E12 11.0 0.993790 E12 0.456373 E12 0.812595 E12 0.420224 E12 12.0 0.993068 E12 0.456191 E12 0.812877 E12 0.420066 E12 13.0 0.991887 E12 0.455755 E12 0.812415 E12 0.419949 E12 14.0 0.992309 E12 0.455810 E12 0.811716 E12 0.419625 E12 15.0 0.991993 E12 0.455510 E12 0.811439 E12 0.419452 E12 16.0 0.991629 E12 0.455356 E12 0.811187 E12 0.419162 E12 17.0 0.990900 E12 0.455079 E12 0.810675 E12 0.418950 E12 18.0 0.990278 E12 0.454602 E12 0.810094 E12 0.418669 E12 19.0 0.990060 E12 0.455671 E12 0.809664 E12 0.418394 E12 20.0 0.989007 E12 0.454249 E12 0.809270 E12 0.418137 E12 21.0 0.989105 E21 0.454246 E12 0.808498 E12 0.417857 E12 22.0 0.988034 E12 0.453885 E12 0.808572 E12 0.417825 E12 23.0 0.987976 E12 0.414323 E12 0.807458 E12 0.417442 E12 Fig. 3 Comparison among the actual HEU and the different LEU fuels distribution of the thermal neutron flux in the IIS and EIS versus CR position in the core

88 INDIAN J PURE & APPL PHYS, VOL 49, FEBRUARY 2011 Table 5 Thermal neutron flux values in the IIS at different control rod positions and increasing values for the internal diameter of the annulus reflector X (cm) Internal radius of the annulus reflector(cm) 11.55 12.55 13.55 14.55 15.55 16.55 17.55 18.55 19.55 1. 1.01250 1.12717 1.20273 1.23055 1.22656 1.22306 1.21110 1.22441 1.24302 2. 1.01286 1.12663 1.20243 1.23039 1.22676 1.22361 1.212171 1.22576 1.24423 3. 1.01337 1.12675 1.20271 1.23067 1.22676 1.22433 1.21324 1.22700 1.24530 4. 1.01290 1.12675 1.20267 1.23067 1.22718 1.22445 1.21344 1.22711 1.24541 5. 1.01297 1.12688 1.20274 1.23067 1.22723 1.22404 1.21292 1.22656 1.24486 6. 1.01324 1.12666 1.20236 1.23016 1.22685 1.22297 1.21161 1.22491 1.24340 7. 1.01276 1.12600 1.20175 1.22945 1.22566 1.22242 1.21075 1.22437 1.24290 8. 1.01178 1.12666 1.20236 1.23006 1.22605 1.22234 1.21017 1.22346 1.24217 9. 1.01068 1.12531 1.20094 1.22867 1.22432 1.22039 1.20787 1.22105 1.23980 10. 1.01051 1.12581 1.20146 1.22926 1.22481 1.22075 1.20807 1.22122 1.24001 11. 1.01081 1.12581 1.20148 1.22937 1.22498 1.22068 1.20811 1.22124 1.24009 12. 1.00924 1.12558 1.20127 1.22905 1.22468 1.22057 1.20803 1.22127 1.24006 13. 1.01004 1.12505 1.20070 1.22849 1.22409 1.22016 1.20762 1.22085 1.23970 14. 1.00929 1.12427 1.19984 1.22745 1.22326 1.21929 1.20698 1.22025 1.23898 15. 1.00868 1.12319 1.19854 1.22615 1.22182 1.21811 1.20580 1.21905 1.23779 16. 1.00760 1.12277 1.19820 1.22572 1.22131 1.21737 1.20502 1.21822 1.23699 17. 1.00713 1.12198 1.19722 1.22477 1.22037 1.21657 1.20421 1.21743 1.23610 18. 1.00686 1.12034 1.19536 1.22275 1.21837 1.21459 1.20237 1.21557 1.23410 19. 1.00603 1.11826 1.19306 1.22032 1.21599 1.21228 1.20011 1.21327 1.23169 20. 1.00558 1.11588 1.19044 1.21749 1.21324 1.20947 1.19748 1.21069 1.22897 21. 1.00527 1.11336 1.18763 1.21462 1.21034 1.20666 1.19474 1.20785 1.22606 22. 1.00477 1.11149 1.18554 1.21240 1.21815 1.20458 1.19270 1.20586 1.22391 23. 1.00437 1.11371 1.18815 1.21518 1.21072 1.20707 1.19493 1.20800 1.22624 Fig. 4 Variation of the thermal neutron flux in the IIS versus CR position for different internal radii of the annulus reflector

ALBARHOUM: INVARIABILITY OF NEUTRON FLUX IN MNSRs 89 Table 6 Values of some statistical parameters for the neutron thermal flux as a function of the CR position in the IIS at different control rod positions and increasing values for the internal diameter of the AR Item Internal radius of the annulus reflector(cm) 12.55 13.55 14.55 15.55 16.55 17.55 18.55 19.55 ANF (n/cm 2.s) E-12 1.123 1.198 1.225 1.222 1.218 1.206 1.219 1.238 STDEV 0.005 0.005 0.006 0.005 0.006 0.006 0.006 0.006 VC (%) 0.439 0.458 0.475 0.433 0.500 0.519 0.525 0.525 Table 7 Values of some statistical parameters for the neutron thermal flux as a function of the CR position in the EIS at different CR positions and increasing values for the internal diameter of the AR Item Internal radius of the annulus reflector(cm) 12.55 13.55 14.55 15.55 16.55 17.55 18.55 19.55 ANF (n/cm 2.s) E-12 0.483 0.489 0.482 0.479 0.455 0.431 0.402 0.369 STDEV 0.006 0.006 0.006 0.002 0.002 0.002 0.002 0.002 VC (%) 12.445 12.483 12.556 0.475 4.843 0.485 0.502 0.495 statistical parameters for the EIS are presented. The decrease of the IR of the AR will cause the IER to decrease as well. The addition of new Fuel Rods (FR) will be necessary to recover the lost IER. Now, if the FR were increased to have an acceptable IER it can be shown that 228 FR in the 10 th circle would be required. In this case, the values of the neutron flux in the IIS and EIS will be those as presented in Table 8 with the presence of 3 Dummy Rods (DR) which are made of aluminium and have the same dimensions of the FR. In Table 8, X indicate the CR position. For this configuration of the core the values of the ANF, STDEV and VC were: 0.738987E12 n/cm 2.s, 0.00286E12 and 0.387%, respectively. It can be noticed that the IIS are now in the core of the reactor (and not in the AR) the thermal neutron flux was practically constant with respect to the variation of the CR position. The main neutronic parameters of the reactor are presented in Table 9. Although the fluxes in the IIS and EIS reduced by about 27%, the new core conserved a negative Shut Down Margin (SDM) and an even better Effective Shut Down margin (ESDM). Since this core was acceptable from the neutronic point of view, for the thermal hydraulics it needed some DR to be added to have an acceptable sub-channel shape. This can be satisfied if other 183 FR were added besides to other 369 new DR. This lead to the neutronic characteristics as presented in Table 10. This final configuration satisfied all safety requirements, but had a value of the neutron flux decreased by about 42% with respect to the original reference value (the value of the UAl 4 -Al HEU fuel). The core contained now about 50% more 235 U than the original one. For this configuration, the values of Table 8 Thermal neutron flux values in the IIS and EIS at different control rod positions when 8 cm of Be are substituted by fuel rods X Flux in IIS Flux in EIS (cm) (n/cm2.s) (n/cm2.s) 1. 0.740844 0.378640 2. 0.740378 0.371405 3. 0.740348 0.371376 4. 0.740448 0.371142 5. 0.741474 0.370144 6. 0.741068 0.359970 7. 0.741366 0.369550 8. 0.740296 0.369042 9. 0.740830 0.368723 10. 0.741037 0.368414 11. 0.741088 0.368300 12. 0.740859 0.368172 13. 0.740438 0.368059 14. 0.739792 0.367926 15. 0.739727 0.367626 16. 0.739127 0.367498 17. 0.738059 0.367187 18. 0.736707 0.366804 19. 0.735154 0.366364 20. 0.733460 0.365914 21. 0.732153 0.365611 22. 0.733985 0.365840 23. 0.735483 0.366654 the NF versus CR position are presented in Table 11. Since the ANF, STDEV and VC were: 0.690092E12 n/cm 2.s, 0.002E12 and 0.387%, respectively for the IIS, and 0.381500E12 n/cm 2.s, 0.002E12 and 0.499%, respectively for the EIS, this proves that this property is independent of the fuel type and of the reflector

90 INDIAN J PURE & APPL PHYS, VOL 49, FEBRUARY 2011 Table 9 Main neutronic characteristics of the new and old cores for 228 FR and 3 DR Item Calculated Measured Initial Excess Rea. (IER) with CRO (mk) 3.7449 3.9400 ± 0.02 Neutron Flux in (I.I.S.) (10 12 n/cm 2.s) 0.735483 1.000000 ± 1% Neutron Flux(E.I.S) (10 12 n/cm 2.s) 0.366654 0.500000± 1% Control Rod Worth (mk) 5.6395 7.000000± 0.01 Shut-Down Margin (mk) 1.8946 3.060000± 0.02 Flooding the IIS with water (mk) 0.6368 1.9500 Flooding the EIS with water (mk) 1.1431 1.0000 Effective Shut-Down Margin (mk) 0.1147 0.1100 Mass of 235 U in core (kg) Table 10 Main neutronic characteristics of the new and old cores for 245 FR and 372 DR Item Calculated Measured Initial Excess Rea. (IER) with CRO (mk) 3.7538 3.9400 ± 0.02 Neutron Flux in (I.I.S.) (10 12 n/cm 2.s) 0.687022 1.000000 ± 1% Neutron Flux(E.I.S) (10 12 n/cm 2.s) 0.379765 0.500000± 1% Control Rod Worth (mk) 5.7819 7.000000± 0.01 Shut-Down Margin (mk) 2.0281 3.060000± 0.02 Flooding the IIS with water (mk) 0.4395 1.9500 Flooding the EIS with water (mk) 1.3194 1.0000 Effective Shut-Down Margin (mk) 1.5866 0.1100 Mass of 235 U in core (kg) 1527.16 999.85 Table 11 Thermal neutron flux values in the IIS and EIS at different control rod positions when 8 cm of Be are substituted by fuel rods X Flux in IIS Flux in EIS (cm) (n/cm 2.s) (n/cm 2.s) 1. 0.691699 0.384184 2. 0.691113 0.384730 3. 0.691076 0.384709 4. 0.691262 0.384454 5. 0.692419 0.383343 6. 0.692097 0.383145 7. 0.692530 0.382722 8. 0.691664 0.382210 9. 0.692176 0.381853 10. 0.692323 0.381556 11. 0.692293 0.381371 12. 0.692073 0.381244 13. 0.691646 0.381145 14. 0.690944 0.380968 15. 0.690919 0.380657 16. 0.690381 0.380552 17. 0.689341 0.380234 18. 0.687022 0.379765 19. 0.688062 0.0379839 20. 0.686622 0.379407 21. 0.685029 0.378934 22. 0.683855 0.378676 23. 0.685578 0.378868 thickness and type. The CR position does not affect the values of the neutron flux in the IIS and EIS anyway because the positions of these sites are far from the CR axis although the values of the neutron flux depend on the fuel type for the same power (30 kw). 4 Conclusions The control rod vertical position along its path does not affect the value of the thermal neutron flux in the IIS and EIS because of the considerable distance of these sites from the CR axis, but the value of the neutron flux in these sites is affected by the fuel type and enrichment. The thickness of the reflector does not affect the ANF in both the IIS and EIS as well. The maximum variation of the ANF values occur in the case of the UO 2 LEU fuel. Acknowledgement The author thanks Prof I Othman, Director General of the Atomic Energy Commission of Syria, for his encouragement and continued support. References 1 Albarhoum M, Annals of Nuclear Energy,35 (2008) 1760. 2 Askew J R, Fayers F J & Kemshell P B, J British Nucl Energy Soc, 5 (1966) 564. 3 CIAE Safety Analysis Report (SAR) for the Syrian Miniature Neutron Source Reactor, China (1993). 4 Fowler T B, Vondy D R & Cunningham G W, Nuclear Reactor Core Analysis Code: CITATION ORNL-TM-2496, Rev 2, July (1971). 5 Khattab K, Omar H & Ghazi N, Nucl Engineering and Design, 236 (23) (2006) 2419. 6 Tayyab M, Showket P & Masood, I Annals of Nucl Energy (35) (2008) 1440.