Estimation of Intervention Distances for Urgent Protective Actions Using Comparative Approach of MACCS and InterRAS

Abstract

Distances for taking evacuation as a protective measure during early phase of a nuclear accident have been approximated using MELCOR Accident Consequence Code System (MACCS). As a reference data, the source term of Pakistan Research Reactor 1 (PARR-1) and meteorological data of Islamabad, Pakistan, have been considered. Based on comparison with published data and international radiological assessment (InterRAS) code results, it is concluded that MACCS is a rational tool for estimation of urgent protective actions during early phase of nuclear accident by taking into account the variations in meteorological and release concentrations parameters.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Mazzammal Hussain 1 and Salah Ud-Din Khan 2 and Waqar A. Adil Syed 1 and Shahab Ud-Din Khan 3

Academic Editor:Inn Seock Kim

1, Department of Physics, Faculty of Basic and Applied Sciences, International Islamic University, H-10, Islamabad 44000, Pakistan
2, Deanship of Scientific Research, College of Engineering, King Saud University, P.O. Box 800, Riyadh 11421, Saudi Arabia
3, Tokamak Design Division, Institute of Plasma Physics, Chinese Academy of Sciences, 350 Shushanhu Road, Hefei, Anhui 230031, China

Received 19 January 2014; Revised 1 April 2014; Accepted 3 May 2014; 28 May 2014

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

A range of probable accidents are associated with nuclear power plants starting with minor incident to immense disaster. Different researchers had used MACCS code to analyze the nuclear power plants accident consequences. Jongtae Jeong and Wondea Jung studied estimation of early health effects for different combinations of release parameters and meteorological data using MACCS code for Younggwang 3 and 4 nuclear power plants in Korea concluding that, with the same amount of radioactive material released to the atmosphere, a large difference in early health effects from case to case was observed [1]. Haste and coworkers attempted to demonstrate a MELCOR-MACCS capability to simulate the accident scenarios for the whole plant, including the containment response and off-site consequences arising from fission product release from the containment. The results provided a good basis for the NPP analysis foreseen [2]. Ke-Shih performed reevaluation of emergency planning zones (EPZs) for nuclear power plants (NPPs) using MACCS2 code and concluded that the radius identified previously is a reasonable conservative value of EPZs for each of the three operating NPPs in Taiwan [3]. Jeong and Ha studied influence of source term release parameters on health effects for Younggwang 3 and 4 nuclear power plants in Korea using MACCS code and concluded that the research work will be very useful for developing strategies for reducing off-site consequences of accident management if they are combined with influence of weather conditions [4]. Thaning and Baklanov considered a simulated accident at a nuclear power plant that could cause a large release of radioactivity into the atmosphere. The consequence analysis was performed using a 3-dimensional mesoscale model and MACCS code [5].

2. Research Methodology

2.1. Computational Technique: MACCS Code

Sandia National Laboratories developed MELCOR Accident Consequence Code System (MACCS) for the consequence assessment of severe accidents at nuclear power plants. MACCS is organized into three modules, the ATMOS module which performs the atmospheric transport and deposition portion of the calculation, the EARLY module which estimates the consequences of the accident immediately after the accident usually within the first week, and the CHRONC module which estimates the long term consequences of the accident. The flow diagram of MACCS code [6-8] is shown in Figure 1.

Figure 1: Flow diagram of MACCS code.

[figure omitted; refer to PDF]

In MACCS code, Gaussian plume model has been used for the study of atmospheric dispersion of radioactive material and vertical and crosswind distributions. The plume dimensions are defined in vertical and crosswind direction by the standard deviations ( _{σ y} , _{σ z} ) of the normal distribution of material concentration in crosswind and vertical directions. The general form of the Gaussian plume equation is [figure omitted; refer to PDF] where χ ( x , y , z ) is the time integrated air concentration (Bq-s/m³ ) at the downwind location ( x , y , z ) , Q is the source strength (Bq), u - is the mean downwind speed, _{σ y} and _{σ z} are the standard deviations (m) of the normal concentration distribution along crosswind and vertical direction, and h is the release height (m).

Equation (1) is not applicable when plume expands vertically and is bounded by mixing layer or by the ground. To solve this problem, ground and mixing layers are considered as totally reflecting boundaries. This is achieved by adding a mirror image source below the ground and above the inversion layer. By considering this effect in (1), the centerline air concentration χ ( x = 0 , y = 0 , z = H ) and ground concentrations χ ( x = 0 , y = 0 , z = 0 ) , after time of release to the time at which the concentrations become uniform along vertical direction, are given by [figure omitted; refer to PDF] where H = ( h + Δ h ) is the height of the plume centerline (m), h = the initial release height of the plume before plume rise (m), Δ h = the amount of plume rise (m), and L = the height (m) of the inversion layer (mixing height).

In MACCS, only first five terms are considered and the rest of the terms are neglected. When a uniform vertical distribution is attained, the following equation is used to calculate centerline air concentration: [figure omitted; refer to PDF] For the estimation of radiation doses during early phase of nuclear emergency, the following pathways are taken into consideration:

(i) external and internal exposure from cloud shine,

(ii) exposures from ground shine,

(iii): internal exposure from resuspension inhalation and skin doses from deposited material onto the skin,

(iv) acute and lifetime doses from early exposures.

2.2. Meteorological Conditions

For the execution of MACCS code, the hourly-based data for the year 2010 of Islamabad collected from the ground station (with an approximate height of 10 meters) has been used. As the near surface wind speed increases with altitude, the same data if used for the stake height (61 m) will overestimate the plume rise. This could produce significant underestimation of the radiation doses. To incorporate the effect of wind speed (m/s) with altitude, the following theoretical formula was used to estimate the wind speed at higher altitudes: [figure omitted; refer to PDF] where u = wind speed at height z (m/s), _{u o} = wind speed at surface (m/s), and p = parameter (dimensionless) that varies with stability class and surface roughness.

The wind rose for the reported year is presented in Figure 2. The percentage value of wind direction remained 3.7% in north, 8.7% in northeast, 8.4% in east, 13.8% in southeast, 10.8% in south, 38.3% in southwest, 6.4% in west, and 10% in northwest, respectively. Data presented in Figure 2 represents that dominant wind direction was southwest (38.3%) during year 2010. The dominant atmospheric stability class as presented in Figure 3 was stability class-F.

Figure 2: Wind rose for sixteen directions.

[figure omitted; refer to PDF]

Figure 3: Dominant stability class.

[figure omitted; refer to PDF]

3. Initial Conditions and Assumptions

3.1. Source Term and Release Fraction

The source term for the fission product in the reactor core was taken from the international published data [6] for a postulated accidental airborne release from Pakistan Research Reactor 1 (PARR-1), Islamabad, on upgraded power of 10 MW. The fractions of releases were based on USNRC document NUREG-1150. The release fractions of 1, 0.4, 0.3, 0.05, and 0.02 for noble gases, halogens, alkali metals, the tellurium group, and the Ba-Sr group, respectively, have been considered.

3.2. Meteorological Options

MACCS code provides different options for the selection of the meteorological data. The following four meteorological options (MO) have been used to estimate the radiation doses:

(i) MO-1: user specified day and hour,

(ii) MO-2: weather bin sampling,

(iii): MO-3: constant weather boundary data,

(iv) MO-4: special case (maximum wind speed).

3.3. Release Scenarios

Scenario 1. In this scenario, the releases and meteorological data have been considered at height of 61 meters (stake releases). The total effective dose equivalent (TEDE) for whole body modeled using four different meteorological options is presented in Figure 4.

Figure 4: Scenario 1 mean TEDE (7 days).

[figure omitted; refer to PDF]

Scenario 2. In this scenario, the releases and meteorological data have been considered at ground level (10 m). The TEDE for whole body modeled by using four different meteorological options is presented in Figure 5.

Figure 5: Scenario 2 mean TEDE (7 days).

[figure omitted; refer to PDF]

3.4. Comparison of Doses for Constant Weather Condition

The comparison of radiation doses of the scenarios has been made with the literature data (using the meteorological conditions of the site [6]) for 61-meter height using constant weather conditions. The comparison is presented in Figure 6.

Figure 6: Comparison of mean TEDE (7 days) for constant weather.

[figure omitted; refer to PDF]

4. Comparison of MACCS and InterRAS Codes

Comparison of MACCS code with InterRAS code [9, 10] for releases at the height of 61 meters using constant weather conditions has been made. The InterRAS code estimate radiation doses maximum up to 48 hours. The doses estimated using InterRAS code are comparable to doses estimated using MACCS code. The comparison of MACCS and InterRAS code results is presented in Figure 7.

Figure 7: Comparison of MACCS and InterRAS codes.

[figure omitted; refer to PDF]

4.1. Mean TEDE Trends (One Week)

Mean TEDE for different release heights and meteorological conditions has been analyzed and presented in Figures 4-6. Through analysis, it was found that radiation doses reduce exponentially over the distance.

It was found that the protective action "evacuation" was required at different distances in different accident situations. The radiation doses fall below the intervention level for evacuation (50 mSv) at 10 kilometers for Scenario 1, 10.75 kilometers for Scenario 2, and 2.75 kilometers for constant weather condition using site specific data.

4.2. Mean Distance for Intervention Levels

The intervention level for evacuation (50 mSv) was achieved at different distances for different emergency scenarios and for different meteorological options. Mean distances (Km) for taking intervention (evacuation) for different meteorological options (MO) are presented in Table 1.

Table 1: Evacuation intervention distance for different MET options.

Scenarios	Intervention (evacuation) distances (Km)
MO-1	MO-2	MO-3	MO-4
Scenario 1	2.00	7.75	10.00	0.25
Scenario 2	2.75	6.75	10.75	4.75
Literature [6]	--	--	2.75	--

A comparison of MACCS and InterRAS code output results (mean effective dose equivalent to whole body for two days) for constant weather conditions has been presented in Figure 7. From the trend analysis, it is found that the maximum of 5.10 × ^{10 3} mSv has been observed at mean distance of 0.25 km. At the nearby distances, the doses estimated using InterRAS code are comparable to doses estimated using MACCS code.

5. Summary and Conclusion

The most frequent wind direction southwest, that is, 38.3%, and the prevailing stability class F, that is, 58.16%, were recorded during year 2010. Intervention distance for evacuation remained in the range 2.0 to 2.75 kilometers for user specified MET data; 6.75 to 7.75 kilometers for weather bin sampling method; 2.75 to 10.75 kilometers for constant weather data; and 0.25 to 4.75 kilometers for maximum wind speed. With site specific meteorological conditions, intervention distance range is very low, that is, 2.0 to 2.75 kilometers, which is very close to the already published value, that is, 1 to 2 kilometers, using "InterRAS" code.

The seven-day TEDE has been modeled and compared with the reference levels [11] for evacuation, that is, 50 mSv for not more than seven days. It is found that the protective action "evacuation" was required at different distances in different accident scenarios.

MACCS code during early phase of nuclear emergency reasonably estimates the radiation doses by taking into account the variation in meteorological data, release duration, and concentrations. MACCS code could be an efficient tool to be used during early phase of nuclear emergencies if integrated with pre- and postprocessor for handling and presentation of data.

Acknowledgment

The authors would like to extend their sincere appreciation to the deanship of scientific research at King Saud University for its funding of this research through research group Project no. RGP VPP-255.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] J. Jeong, W. Jung, "The estimation of early health effects for different combinations of release parameters and meteorological data," Journal of Korean Nuclear Society , vol. 33, no. 6, pp. 557-565, 2001.

[2] T. Haste, J. Birchley, E. Cazzoli, J. Vitazkova, "MELCOR/MACCS simulation of the TMI-2 severe accident and initial recovery phases, off-site fission product release and consequences," Nuclear Engineering and Design , vol. 236, no. 10, pp. 1099-1112, 2006.

[3] J. Wu, Y. Yang, I. Chen, H. Chen, K. Chuang, "Reevaluation of the emergency planning zone for nuclear power plants in Taiwan using MACCS2 code," Applied Radiation and Isotopes , vol. 64, no. 4, pp. 448-454, 2006.

[4] J. Jeong, J. Ha, "The influence of source term release parameters on health effects," Journal of The Korean Nuclear Society , vol. 31, no. 3, pp. 294-302, 1999.

[5] L. Thaning, A. Baklanov, "Simulation of the atmospheric transport and deposition on a local/meso- and regional scale after hypothetical accidents at the Kola nuclear power plant," Science of the Total Environment , vol. 202, no. 1-3, pp. 199-210, 1997.

[6] "MELCOR Accident Consequence Code System (MACCS)," User's Guide, NUREG/CR-4691-Vol. 1

[7] "MELCOR Accident Consequence Code System (MACCS)," Model Description, NUREG/CR-4691-Vol. 2

[8] "MELCOR Accident Consequence Code System (MACCS)," Programmer's Reference Manual, NUREG/CR-4691-Vol. 3

[9] C. Lin, L. Peng, "Study of interRAS application in early radiation consequences assessment of attacked nuclear power plant," in Proceedings of the 11th China Symposium on Computer Application in Modern Science and Technology, pp. 105-108, Beijing, China, 2003.

[10] M. Hussain, S. U. Khan, W. A. Syed, "Estimation of emergency planning zones for nuclear research reactor using plume dispersion code," in Proceedings of the 20th International Conference on Nuclear Engineering Collocated with the ASME, Power Conference (ICONE20-POWER '12), pp. 479-483, Anaheim, Calif, USA, July-August 2012.

[11] "Regulations on Management of a Nuclear or Radiological Emergency-PAK/914," 2008