# Source Term Derivation and Radioactive Release Evaluation for JRTR Research Reactor under Severe Accident

Authors: Ned Xoubi
Journal: Science and Technology of Nuclear Installations
 (2020)

---

## Abstract

The source term for the JRTR research reactor is derived under an assumed hypothetical severe accident resulting in generation of the most severe consequences. The reactor core is modeled based on the reactor technical design specifications, and the fission products inventory is calculated by using the SCALE/TRITON depletion sequence to perform burnup and decay analyses via coupling the NEWT 2-D transport lattice code to the ORIGEN-S fuel depletion code. Fifty radioisotopes contributed to the evaluation, resulting in a source term of 3.7 × 1014 Bq. Atmospheric dispersion was evaluated using the Gaussian plume model via the HOTSPOT code. The plume centerline total effective dose (TED) was found to exceed the IAEA limits for occupational exposure of 0.02 Sv; the results showed that the maximum dose is 200 Sv within 200 m from the reactor, under all the weather stability classes, after which it starts to decrease with distance, reaching 0.1 Sv at 1 km from the reactor. The radiation dose plume centerlines continue to the exceed international basic safety standards annual limit of 1 mSv for public exposure, up to 80 km from the reactor.


---

## Body

## 1. Introduction

The utilization of research reactors in a variety of scientific, production, and training applications requires unique design characteristics and operational modes, consequently exerting additional burden on their safety systems and procedures.Although, in comparison to power reactors, research reactors are small and contain smaller inventory of fuel and radionuclides, they have a higher number of reported accidents. The lack of containment building, proximity to densely populated areas, and high burnup cores, coupled with life extension, have raised worldwide interest for in-depth safety evaluation of these systems.In 2008, the international atomic energy agency (IAEA) published a research reactor safety report [1]. It aims to provide safety analysts, regulators, and reactor operation staff and management, with the essential calculation methods and techniques for evaluating the source term and analyzing the radiological consequences of accidents in research reactor.The magnitude, composition, form, and mode of release of radioactive elements released during a reactor accident define the source term. The derivation of the source term is an essential step of the safety analysis [1]. One approach commonly used in research reactors safety analyses is to assume a hypothetical accident that results in generating the most severe consequences [1]; this is the subject of this paper.Recently, several studies for the evaluation of research reactors source term have been conducted. Foudil et al. have estimated the source term and doses for the NUR research reactor under a hypothetical severe accident [2]. Muswema et al. carried out a study to drive the source term and analyze the radiological safety for the TRICOII research reactor in Kinshasa [3]. Charalampos et al. published a study about source term derivation and radiological consequences for the Greek Research Reactor [4]. Birikorang et al. have estimated the source term and ground deposition of radionuclides following a hypothetical release from Ghana Research Reactor [5]. Sadeghi et al. calculated the source term and dose of a postulated accident release from the Tehran Research Reactor [6]. Ullah et al. carried out a study to evaluate the source term for the upgraded LEU Pakistan Research Reactor-1 [7].In this work, the JRTR research and training reactor core inventory is calculated using the SCALE/TRITON [8, 9] depletion sequence to perform burnup and decay analyses via coupling the NEWT 2-D [10] transport lattice code to the ORIGEN-S [11] fuel depletion code.The corresponding source term and the release of radionuclides during postulated severe accident are estimated in accordance with the IAEA [1] recommendation and the US-NRC regulatory guide 1.183 [12]. The atmospheric dispersion of radionuclide concentrations are simulated based on the Gaussian Plume Model (GPM) using the HOTSPOT code [13].


## 2. Reactor Description

The JRTR is a multipurpose research reactor optimized for research, training, and isotope production. It is light water moderated and cooled open-tank-in-pool type reactor with nominal power of 5 MW upgradable to 10 MW. The reactor core (shown in Figure1) consists of 18 MTR plate type fuel assemblies loaded in a 33 squares grid matrix, the remaining 15 grid positions are loaded with beryllium blocks serving as internal reflector, and the grid central position is designated as a flux trap experimental facility. Heavy water reflectors are placed at the four corners of the core. The fuel assemblies and beryllium blocks are interchangeable and can be shuffled to reconfigure the reactor core.Figure 1
Cross-sectional (x-y) view of the JRTR reactor core configuration.The core is designed to produce the utmost thermal fluxes at the core’s internal and external irradiation facilities. The maximum in-core thermal-flux-to-power ratio is 2.9 × 1013 cm−2s−1 per MW, producing a maximum thermal flux of 1.45 × 1014 cm−2s−1, at full power. The maximum external thermal flux to power ratio is 1.4 × 1013 cm−2s−1 per MW producing a maximum thermal flux of 6.78 × 1013 cm−2s−1 at full power.Forced convection cooling is used to cool the reactor under normal operating conditions. The coolant flows downward at a rate of 170 kg/s, with an inlet temperature of around 37°C and exits the core with an outlet temperature of 44.2°C. The average heat flux is 17.3 W/cm2, with a maximum of 51.9 W/cm2, and the power peaking factor is 2.54 [14–17]. The main design characteristics of the JRTR reactor are listed in Table 1.Table 1
Design characteristics of the JRTR research and training reactor.
Reactor typeTank-in-poolNominal thermal power (MW)5Max thermal flux (n/cm2 s)1.7 × 1014No. of fuel assemblies18Fuel materialU3Si2Uranium enrichment (wt)19.75%No. of control elements (box shape)4Control plate materialHafniumModeratorLight waterInternal reflectorBerylliumExternal reflectorHeavy waterCoolant (flow direction)Light water (downward)Coolant inlet temperature (°C)37.0Water channel thickness (mm)2.35Water channel thickness between two fuel elements (cm)0.390Max. heat flux (W/m2)1,170,000Nominal flow rate (m3/h)612Average core flow velocity (m/s)2.5Core inlet pressure (MPa)0.18


### 2.1. Fuel Assembly

The fuel is MTR plate type, uranium silicide (U3Si2) enriched to 19.75% 235U and in aluminum clad. The core loading is 18 standard fuel assemblies containing approximately 7.0 kg of 235U, as shown in Figure 1. Each assembly has 21 plates with a thickness of 1.27 mm and are placed 2.35 mm a part, allowing cooling water to flow between them as shown in Figure 2. The fuel material is 19.75 wt% enriched uranium silicide (U3Si2) in aluminum matrix sandwiched between two 0.38 mm thick aluminum plates, 680.0 mm long, and 70.7 mm wide.Figure 2
TRITON 2-D model of the JRTR fuel assembly showing the fuel (red), aluminum clad (green), and water coolant (blue).The fuel meat covers an area of 640 × 6.21 mm in each plate with a density of 6.543 g/cm3. Fuel assembly contains 97.2939 g of uranium in each plate, with uranium density of 4.8 g/cm3. The fuel technical specifications and material data are listed in Table 2.Table 2
Technical specifications and material data of the fuel element.
Fuel assembly dimensions (mm)1015 × 76.2 × 8.0No. of plates per fuel element21Fuel materialU3Si2235U enrichment (wt %)19.75Uranium density (g/cm3)4.8Fuel meat density (g/cm3)6.543Fuel plate dimensions (mm)680.0 × 70.7 × 1.27Fuel plate active dimensions (mm)640 × 62.1 × 0.51Cladding material (aluminum alloy)AG3NECladding density (g/cm3)2.700Cladding thickness (mm)0.38Fuel cycle length (EFPD)50Discharge EFPD900


### 2.2. Operating Cycle

The life cycle length of the JRTR is 50 days of full power operation, with an average discharge burnup of 900 effective full power days (EFPD). Fuel reloading consists of adding one new assembly into the F17 position and discharging one assembly from position F09. All of the fuel assemblies are shuffled, such that each assembly is placed in one new position every cycle, covering all 18 positions in the core before being discharged. Thus, at the end of cycle, the core will have assemblies that are burned from one cycle to 18 cycles. Fuel reloading is assumed to take 5 days, during which the reactor is shutdown.


## 2.1. Fuel Assembly

The fuel is MTR plate type, uranium silicide (U3Si2) enriched to 19.75% 235U and in aluminum clad. The core loading is 18 standard fuel assemblies containing approximately 7.0 kg of 235U, as shown in Figure 1. Each assembly has 21 plates with a thickness of 1.27 mm and are placed 2.35 mm a part, allowing cooling water to flow between them as shown in Figure 2. The fuel material is 19.75 wt% enriched uranium silicide (U3Si2) in aluminum matrix sandwiched between two 0.38 mm thick aluminum plates, 680.0 mm long, and 70.7 mm wide.Figure 2
TRITON 2-D model of the JRTR fuel assembly showing the fuel (red), aluminum clad (green), and water coolant (blue).The fuel meat covers an area of 640 × 6.21 mm in each plate with a density of 6.543 g/cm3. Fuel assembly contains 97.2939 g of uranium in each plate, with uranium density of 4.8 g/cm3. The fuel technical specifications and material data are listed in Table 2.Table 2
Technical specifications and material data of the fuel element.
Fuel assembly dimensions (mm)1015 × 76.2 × 8.0No. of plates per fuel element21Fuel materialU3Si2235U enrichment (wt %)19.75Uranium density (g/cm3)4.8Fuel meat density (g/cm3)6.543Fuel plate dimensions (mm)680.0 × 70.7 × 1.27Fuel plate active dimensions (mm)640 × 62.1 × 0.51Cladding material (aluminum alloy)AG3NECladding density (g/cm3)2.700Cladding thickness (mm)0.38Fuel cycle length (EFPD)50Discharge EFPD900


## 2.2. Operating Cycle

The life cycle length of the JRTR is 50 days of full power operation, with an average discharge burnup of 900 effective full power days (EFPD). Fuel reloading consists of adding one new assembly into the F17 position and discharging one assembly from position F09. All of the fuel assemblies are shuffled, such that each assembly is placed in one new position every cycle, covering all 18 positions in the core before being discharged. Thus, at the end of cycle, the core will have assemblies that are burned from one cycle to 18 cycles. Fuel reloading is assumed to take 5 days, during which the reactor is shutdown.


## 3. Materials and Methods

In this study, the core inventory is calculated using the SCALE/TRITON depletion sequence to perform burnup and decay analyses via coupling the NEWT 2-D transport lattice code to the ORIGEN-S fuel depletion code.The TRITON depletion sequences use the 2-D arbitrary polygonal mesh discrete ordinates transport code NEWT to calculate the flux distribution in the assembly.In the T-DEPL calculation sequence, the cross section processing and neutron transport solution of the NEWT sequence is used to create three-group weighted cross sections based on calculated volume averaged fluxes for each material. These newly weighted cross sections are used to update the ORIGEN-S cross section library via the COUPLE code. The NEWT calculated three-group averaged fluxes for each material and they are used for depletion calculations in the ORIGEN-S code [11]. The three-neutron energy groups are the thermal group below 0.625 eV, the resonance group 0.625 eV to 1 MeV, and the fast group above 1 MeV [11].T-DEPL uses a rigorous iterative procedure to update both fluxes and cross sections as a function of burnup. First, transport calculations are performed to estimate fluxes and prepare weighted cross-sections and other lattice physics parameters. Second, depletion calculations are performed using COUPLE and ORIGEN-S codes to update the nuclide concentrations. At each interval, ORIGEN-S updated nuclide concentrations are input into the assembly lattice transport analysis, and the neutron flux spectrum is recalculated. Cross sections are again updated and used by ORIGEN-S for the following computational step [11].


### 3.1. Assembly Modeling

The TRITON 2-D model was developed for the fuel assembly based on its technical design data provided in Table2. The geometrical dimensions and material compositions of the fuel, clad, and moderator were explicitly modeled, as illustrated in Figure 2.Reactor design operational and power data of 50 effective full power days, refueling cycle, and 5 days shutdown established the burnup step in the model. Burnup calculations were performed using burnup steps of 135.95 MW/MTU, with one additional small step at the beginning of the cycle (BOC) to generate cross sections for effectively unirradiated fuel.The moderator/coolant temperature of 317 K was used as the average water temperature and 325 K as the fuel temperature. The moderator density of water from steam tables at 317 K and 1.4 bars is 0.990 g/cm3.


### 3.2. Core Modeling

Aimed at modeling a realistic worst case scenario, the core is modeled at the end of cycle (EOC) for a fully burned up core. As shown in Figure3, each assembly has been burned for a different number of cycles, and thus it will have a different radionuclide inventory. In each cycle, it is assumed that all assemblies will contribute equally to the reactor power, which results in assembly F09 discharge burnup of 122.36 MWd/kg. The fuel burnup in MWd/kg for each assembly is shown in Figure 3.Figure 3
Fuel burnup (MWD/Kg) of each assembly at the end of cycle of a fully burned up core.


### 3.3. Selected Hypothetical Accident Scenario

The research reactor accident scenario, that is, the accident initiating event and its sequence determine both the magnitude and the nature of radioactive release, and thus the characteristic progressions of different scenarios generally lead to different radioactive release patterns. In this approach, a comprehensive assessment of accident progression for several accident scenarios is performed to calculate several different source terms.Another approach that is commonly used in performing the safety analyses of research reactors, and adopted in this paper, is to assume a hypothetical severe accident that results in a limiting source term. The accident scenario considered in this study assumes that one or more projectiles target the reactor, leading to the destruction the reactor building and pool, resulting in the uncovering of the core and damaging all of the fuel assemblies.Not all the fission products that are released from the damaged fuel assemblies transfer to the atmosphere of the reactor building. The atmospheric release fractions assumed in this study are 1.0 for noble gases, 0.4 for halogens, 0.3 for alkali metals, 0.02 for alkaline earths, 0.05 for tellurium, and 0.0025 for noble metals [1].


### 3.4. Site-Specific Conditions

The JRTR is located in Jordan University of Science and Technology campus, home to more than 20,000 students. The city of Ramtha is approximately 7 km away, with a population of 170,000 people, and approximately two million people reside in Irbid metropolitan located about 15 km from the site.The site has a Mediterranean climate with four seasons. The summer average temperatures range from 27°C to 33°C and drops to an average of 10°C in the winter. The average annual rainfall is 300 mm, spreading over 77 rainy days.The prevailing average hourly wind direction in the JRTR site varies throughout the year. However, for about 11 months of the year, it is most often from the west and it switches to the eastern direction from November 12 to December 20. The average hourly wind speed experiences mild seasonal variation over the year, with a speed of more than 3.2 m/s from November 14 to May 2 and 2.8 m/s for the rest of the year [18]. The frequency of wind speed and direction are shown in Figure 4 [19]. The wind rose diagram for Ramtha shows the annual number of hours the wind blows from the indicated direction.Figure 4
The wind rose for Ramtha showing the annual number of hours the wind blows from each direction [19].The atmospheric stability classes are estimated using a number of meteorological measurements such as weather condition, solar insulation, humidity, and wind speed and are classified fromA (very unstable) to F (moderately stable) [13].


## 3.1. Assembly Modeling

The TRITON 2-D model was developed for the fuel assembly based on its technical design data provided in Table2. The geometrical dimensions and material compositions of the fuel, clad, and moderator were explicitly modeled, as illustrated in Figure 2.Reactor design operational and power data of 50 effective full power days, refueling cycle, and 5 days shutdown established the burnup step in the model. Burnup calculations were performed using burnup steps of 135.95 MW/MTU, with one additional small step at the beginning of the cycle (BOC) to generate cross sections for effectively unirradiated fuel.The moderator/coolant temperature of 317 K was used as the average water temperature and 325 K as the fuel temperature. The moderator density of water from steam tables at 317 K and 1.4 bars is 0.990 g/cm3.


## 3.2. Core Modeling

Aimed at modeling a realistic worst case scenario, the core is modeled at the end of cycle (EOC) for a fully burned up core. As shown in Figure3, each assembly has been burned for a different number of cycles, and thus it will have a different radionuclide inventory. In each cycle, it is assumed that all assemblies will contribute equally to the reactor power, which results in assembly F09 discharge burnup of 122.36 MWd/kg. The fuel burnup in MWd/kg for each assembly is shown in Figure 3.Figure 3
Fuel burnup (MWD/Kg) of each assembly at the end of cycle of a fully burned up core.


## 3.3. Selected Hypothetical Accident Scenario

The research reactor accident scenario, that is, the accident initiating event and its sequence determine both the magnitude and the nature of radioactive release, and thus the characteristic progressions of different scenarios generally lead to different radioactive release patterns. In this approach, a comprehensive assessment of accident progression for several accident scenarios is performed to calculate several different source terms.Another approach that is commonly used in performing the safety analyses of research reactors, and adopted in this paper, is to assume a hypothetical severe accident that results in a limiting source term. The accident scenario considered in this study assumes that one or more projectiles target the reactor, leading to the destruction the reactor building and pool, resulting in the uncovering of the core and damaging all of the fuel assemblies.Not all the fission products that are released from the damaged fuel assemblies transfer to the atmosphere of the reactor building. The atmospheric release fractions assumed in this study are 1.0 for noble gases, 0.4 for halogens, 0.3 for alkali metals, 0.02 for alkaline earths, 0.05 for tellurium, and 0.0025 for noble metals [1].


## 3.4. Site-Specific Conditions

The JRTR is located in Jordan University of Science and Technology campus, home to more than 20,000 students. The city of Ramtha is approximately 7 km away, with a population of 170,000 people, and approximately two million people reside in Irbid metropolitan located about 15 km from the site.The site has a Mediterranean climate with four seasons. The summer average temperatures range from 27°C to 33°C and drops to an average of 10°C in the winter. The average annual rainfall is 300 mm, spreading over 77 rainy days.The prevailing average hourly wind direction in the JRTR site varies throughout the year. However, for about 11 months of the year, it is most often from the west and it switches to the eastern direction from November 12 to December 20. The average hourly wind speed experiences mild seasonal variation over the year, with a speed of more than 3.2 m/s from November 14 to May 2 and 2.8 m/s for the rest of the year [18]. The frequency of wind speed and direction are shown in Figure 4 [19]. The wind rose diagram for Ramtha shows the annual number of hours the wind blows from the indicated direction.Figure 4
The wind rose for Ramtha showing the annual number of hours the wind blows from each direction [19].The atmospheric stability classes are estimated using a number of meteorological measurements such as weather condition, solar insulation, humidity, and wind speed and are classified fromA (very unstable) to F (moderately stable) [13].


## 4. Results and Discussion

### 4.1. Core Inventory

The core fission product inventory calculations have been performed using the SCALE/TRITON depletion sequence, based on reactor maximum full power operation. The activity of 200 fission products has been tracked, consequently accounting for the maximum possible activity in the core. The reactor core radionuclide inventory available for release to the environment is classified into eight groups, based on the NRC regulatory guide 1.183 [12] and is shown in Figure 5 The complete list of the fission products inventory and the their release fractions is included in Supplementary appendix 1 [20].Figure 5
Reactor core radionuclide inventory classified into eight groups.


### 4.2. Source Term

Fifty radioisotopes are selected for the source term evaluation, based on their yield, half-life, and radiotoxicity due to both inhalation and ingestion. The activity released in the form of noble gases and other aerosols is presented in Table3, assuming the core remains exposed and all fission products are released to the atmosphere.Table 3
Source term selected radionuclide inventory and activity release.
RadionuclideGroupHalf-life (t1/2) (ref.)Core inventory (TBq)Release fractionActivity released (TBq)kr83mNoble gas1.83 h7.97E + 021.007.97E + 02kr85Noble gas10.756 y3.49E + 011.003.49E + 01kr85mNoble gas10.776 y1.84E + 031.001.84E + 03kr87Noble gas76.3 m3.71E + 031.003.71E + 03kr88Noble gas2.84 h5.06E + 031.005.06E + 03kr89Noble gas3.15 m6.44E + 031.006.44E + 03xe131mNoble gas11.934 d5.39E + 011.005.39E + 01xe133Noble gas5.2475 d9.91E + 031.009.91E + 03xe135Noble gas9.14 h2.86E + 031.002.86E + 03Xe135mNoble gas15.29 m1.91E + 031.001.91E + 03Xe137Noble gas3.818 m9.53E + 031.009.53E + 03Xe138Noble gas14.08 m9.61E + 031.009.61E + 03i130Halogen12.36 h3.57E + 010.401.43E + 01i131Halogen8.02070 d4.64E + 030.401.86E + 03i132Halogen2.295 h6.93E +030.402.77E + 03i133Halogen20.8 h1.05E +040.404.18E + 03i134Halogen52.5 m1.21E +040.404.84E + 03i135Halogen6.57 h9.83E +030.403.93E + 03Br84Halogen31.80 m1.53E +030.406.13E + 02Cs134Alkali metal2.0648 y2.62E +020.307.87E + 01Cs134mAlkali metal2.903 h1.01E +020.303.03E + 01Cs136Alkali metal13.16 d1.23E +020.303.68E + 01Cs137Alkali metal30.07 y2.88E +020.308.64E + 01Cs138Alkali metal33.41 m1.03E +040.303.09E + 03Cs139Alkali metal9.27 m9.71E +030.302.91E + 03Rb88Alkali metal17.78 m5.11E +030.301.53E + 03Rb89Alkali metal15.15 m6.76E +030.302.03E + 03Ba139Alkaline earths83.06 m9.84E + 030.021.97E + 02Ba140Alkaline earths12.752 d9.40E + 030.021.88E + 02Ba141Alkaline earths18.27 m8.92E + 030.021.78E + 02sr89Alkaline earths50.53 d6.13E + 030.021.23E + 02sr90Alkaline earths28.79 y2.69E + 020.025.39E + 00sr91Alkaline earths9.63 h8.47E + 030.021.69E + 02sr92Alkaline earths2.71 h8.68E + 030.021.74E + 02sb131Tellurium group23.03 m3.99E + 030.052.00E + 02te127Tellurium group9.35 h2.85E + 020.051.42E + 01te127mTellurium group109 d4.01E + 010.052.00E + 00te129Tellurium group69.6 m9.13E + 020.054.57E + 01te129mTellurium group33.6 d1.63E + 020.058.16E + 00te131Tellurium group25.0 m4.08E + 030.052.04E + 02te131mTellurium group30 h7.19E + 020.053.60E + 01te132Tellurium group3.204 d6.85E + 030.053.43E + 02te133mTellurium group55.4 m5.21E + 030.052.61E + 02te134Tellurium group41.8 m1.05E + 040.055.26E + 02Mo99Noble metal65.94 h9.59E + 030.00252.40E + 01rh105Noble metal35.36 h2.16E + 030.00255.39E + 00ru103Noble metal39.26 d4.97E + 030.00251.24E + 01ru105Noble metal4.44 h2.39E + 030.00255.98E + 00tc99mNoble metal6.01 h8.51E + 030.00252.13E + 01Np239Cerium group2.3565 d2.96E + 040.00051.48E + 01Noble gases which can be considered inert and have very low chemical reactivity are the major contributors to the source term, accounting for 62.7% of the released activity, as illustrated in Figure6. The highest release of noble gases is 133Xe with an activity of 9.91 × 1015 Becquerel (Bq); this is followed closely by 138Xe and 137Xe with an activity of 9.61 × 1015 and 9.53 × 1015 Bq, followed by the rest of the xenon radioisotopes. Six krypton isotopes contribute 1.79 × 1016 Bq to the noble gas release, with the largest contribution of 6.44 × 1015 Bq from 89Kr.Figure 6
Source term release activity contribution of each of the eight groups.The second major contributor to the source term is the halogen group consisting of six iodine isotopes in addition to84Br. The group contributes 22% of the source term released activity, mainly from iodine with 1.76 × 1016 Bq. Indeed, the radioisotopes of krypton, xenon, and iodine contribute a total of 6.94 × 1016 Bq accounting for 84% of the source term.Cesium and rubidium which comprise the alkali metal group are the next main contributors with a total activity of 9.80 × 1015 Bq accounting for approximately 12% of the source term.The radioisotopes of the remaining four groups, alkaline earths, tellurium group, noble metals, and cerium group, account for approximately 4% of the released activity, as shown in Figure6.


## 4.1. Core Inventory

The core fission product inventory calculations have been performed using the SCALE/TRITON depletion sequence, based on reactor maximum full power operation. The activity of 200 fission products has been tracked, consequently accounting for the maximum possible activity in the core. The reactor core radionuclide inventory available for release to the environment is classified into eight groups, based on the NRC regulatory guide 1.183 [12] and is shown in Figure 5 The complete list of the fission products inventory and the their release fractions is included in Supplementary appendix 1 [20].Figure 5
Reactor core radionuclide inventory classified into eight groups.


## 4.2. Source Term

Fifty radioisotopes are selected for the source term evaluation, based on their yield, half-life, and radiotoxicity due to both inhalation and ingestion. The activity released in the form of noble gases and other aerosols is presented in Table3, assuming the core remains exposed and all fission products are released to the atmosphere.Table 3
Source term selected radionuclide inventory and activity release.
RadionuclideGroupHalf-life (t1/2) (ref.)Core inventory (TBq)Release fractionActivity released (TBq)kr83mNoble gas1.83 h7.97E + 021.007.97E + 02kr85Noble gas10.756 y3.49E + 011.003.49E + 01kr85mNoble gas10.776 y1.84E + 031.001.84E + 03kr87Noble gas76.3 m3.71E + 031.003.71E + 03kr88Noble gas2.84 h5.06E + 031.005.06E + 03kr89Noble gas3.15 m6.44E + 031.006.44E + 03xe131mNoble gas11.934 d5.39E + 011.005.39E + 01xe133Noble gas5.2475 d9.91E + 031.009.91E + 03xe135Noble gas9.14 h2.86E + 031.002.86E + 03Xe135mNoble gas15.29 m1.91E + 031.001.91E + 03Xe137Noble gas3.818 m9.53E + 031.009.53E + 03Xe138Noble gas14.08 m9.61E + 031.009.61E + 03i130Halogen12.36 h3.57E + 010.401.43E + 01i131Halogen8.02070 d4.64E + 030.401.86E + 03i132Halogen2.295 h6.93E +030.402.77E + 03i133Halogen20.8 h1.05E +040.404.18E + 03i134Halogen52.5 m1.21E +040.404.84E + 03i135Halogen6.57 h9.83E +030.403.93E + 03Br84Halogen31.80 m1.53E +030.406.13E + 02Cs134Alkali metal2.0648 y2.62E +020.307.87E + 01Cs134mAlkali metal2.903 h1.01E +020.303.03E + 01Cs136Alkali metal13.16 d1.23E +020.303.68E + 01Cs137Alkali metal30.07 y2.88E +020.308.64E + 01Cs138Alkali metal33.41 m1.03E +040.303.09E + 03Cs139Alkali metal9.27 m9.71E +030.302.91E + 03Rb88Alkali metal17.78 m5.11E +030.301.53E + 03Rb89Alkali metal15.15 m6.76E +030.302.03E + 03Ba139Alkaline earths83.06 m9.84E + 030.021.97E + 02Ba140Alkaline earths12.752 d9.40E + 030.021.88E + 02Ba141Alkaline earths18.27 m8.92E + 030.021.78E + 02sr89Alkaline earths50.53 d6.13E + 030.021.23E + 02sr90Alkaline earths28.79 y2.69E + 020.025.39E + 00sr91Alkaline earths9.63 h8.47E + 030.021.69E + 02sr92Alkaline earths2.71 h8.68E + 030.021.74E + 02sb131Tellurium group23.03 m3.99E + 030.052.00E + 02te127Tellurium group9.35 h2.85E + 020.051.42E + 01te127mTellurium group109 d4.01E + 010.052.00E + 00te129Tellurium group69.6 m9.13E + 020.054.57E + 01te129mTellurium group33.6 d1.63E + 020.058.16E + 00te131Tellurium group25.0 m4.08E + 030.052.04E + 02te131mTellurium group30 h7.19E + 020.053.60E + 01te132Tellurium group3.204 d6.85E + 030.053.43E + 02te133mTellurium group55.4 m5.21E + 030.052.61E + 02te134Tellurium group41.8 m1.05E + 040.055.26E + 02Mo99Noble metal65.94 h9.59E + 030.00252.40E + 01rh105Noble metal35.36 h2.16E + 030.00255.39E + 00ru103Noble metal39.26 d4.97E + 030.00251.24E + 01ru105Noble metal4.44 h2.39E + 030.00255.98E + 00tc99mNoble metal6.01 h8.51E + 030.00252.13E + 01Np239Cerium group2.3565 d2.96E + 040.00051.48E + 01Noble gases which can be considered inert and have very low chemical reactivity are the major contributors to the source term, accounting for 62.7% of the released activity, as illustrated in Figure6. The highest release of noble gases is 133Xe with an activity of 9.91 × 1015 Becquerel (Bq); this is followed closely by 138Xe and 137Xe with an activity of 9.61 × 1015 and 9.53 × 1015 Bq, followed by the rest of the xenon radioisotopes. Six krypton isotopes contribute 1.79 × 1016 Bq to the noble gas release, with the largest contribution of 6.44 × 1015 Bq from 89Kr.Figure 6
Source term release activity contribution of each of the eight groups.The second major contributor to the source term is the halogen group consisting of six iodine isotopes in addition to84Br. The group contributes 22% of the source term released activity, mainly from iodine with 1.76 × 1016 Bq. Indeed, the radioisotopes of krypton, xenon, and iodine contribute a total of 6.94 × 1016 Bq accounting for 84% of the source term.Cesium and rubidium which comprise the alkali metal group are the next main contributors with a total activity of 9.80 × 1015 Bq accounting for approximately 12% of the source term.The radioisotopes of the remaining four groups, alkaline earths, tellurium group, noble metals, and cerium group, account for approximately 4% of the released activity, as shown in Figure6.


## 5. Atmospheric Dispersion and Dose Calculation

Atmospheric dispersion of 50 selected radionuclides has been evaluated using the Gaussian plume model (GPM) via the HOTSPOT 3.1.0 code [13]. The radiation dose was calculated using FGR 13 library [21] and the ICRP-recommended absorption types [22]. Site-specific meteorological conditions were applied to Pasquill stability classes A through D in this study.The plume calculations are performed for the 50 radionuclides released into the atmosphere from the reactor at an effective height of 5 m, a receptor height of 1.5 m, a wind speed of 1 m/s, and a sampling time of 10 minutes. Two deposition velocities are used: 0.3 cm/sec for respirable particles and 8 cm/sec for nonrespirable particles.The results of the total effective dose (TED) as a function of distance from the reactor are illustrated in Figure7, which shows that the radiation dose may reach people up to 80 km from the reactor, under neutrally stable (D) weather conditions. Under very unstable (A) weather conditions, the dose will be received by people within 10 km from the reactor.Figure 7
Total effective dose (Sv) as a function of downwind distance (km) from the reactor, for different weather stability classes (A–F).The maximum dose is about 200 Sv for most weather stability classes (A–D), and it extends from 15  m to 100 m from the reactor. For people within 300 m of the reactor, the TED dose is above 1 Sv, which is higher than the IAEA limits for occupational exposure of 0.02 Sv per year for radiation workers [23]. Radiation doses exceeding international basic safety standards annual limit of 1 mSv for public exposure [23] extend up to 80 km as shown in Figure 7, depending on the distance from the reactor and the weather stability class.The results of plume centerline ground deposition are depicted in Figure8, which shows that the maximum ground deposition is between 9 and ×108 kBq/m2 and falls within 300 meters from the reactor depending on the weather stability class. Contaminant spreads to 100 km depositing more than 1000 kBq/m2 for all stability classes.Figure 8
Plume centerline ground deposition (kBq/m2) as a function of distance (km) from the JRTR reactor for different weather stability classes (A–F).


## 6. Conclusions

In order to drive the source term, the JRTR reactor was modeled in the SCALE/TRITON code based on its design technical specifications. Burnup and decay analyses for a fully burned-up core using NEWT transport lattice and ORIGEN-S fuel depletion code were performed to calculate the core radionuclides inventory available for release to the environment.A hypothetical accident scenario assuming that all of the fuel assemblies are damaged and the reactor building is destroyed is selected. The fraction of fission products released to the atmosphere of the reactor building is calculated in accordance with the IAEA recommendation, resulting in the release 191 radioisotopes. Fifty radioisotopes are selected based on their yield, half-life, and radiotoxicity, to evaluate the source term, resulting in an activity of 3.7 × 1014 Bq.Although the reactor core radionuclide inventory available for release to the environment comprises of 191 radioisotopes, it was found that krypton, xenon, and iodine account for 84% of the source term. Noble gases, halogens, and alkali metal groups consisting of 27 radioisotopes are the vast contributors to the source term, accounting for more than 96% of the released activity.The Gaussian Plume Model (GPM) via the HOTSPOT code is used to simulate the atmospheric dispersion of radionuclide concentrations. The results show the contaminants spread to 100 km depositing more than 1000 kBq/m2.Communities as far as 80 km could be exposed to radiation doses higher than the international basic safety standard annual limit of 1 mSv for public exposure. Radiation workers in the vicinity of the reactor could be exposed to 50 times the recommended IAEA limits for occupational exposure of 0.02 Sv per year.



---

*Source: 1868965-2020-07-20.xml*