Abstract

Radionuclide dispersion model, which is of critical importance to the emergency response of severe nuclear accident, is used to estimate the consequences arising from accidental or routine releases and to predict areas of high contamination. It is difficult to evaluate the radioactive consequence accurately and rapidly for the accidental release of radionuclides from marine reactor because of the complex mobility feature in the sea. Based on CFD method, a finite-volume, three-dimensional regional oceanic dispersion model was developed in this paper to simulate the dispersion of radionuclides originating from marine reactor. The simulated dose variation of 137Cs presented good agreement with the monitoring data of marine radioactive pollution caused by Fukushima Dai-ichi nuclear accident, which demonstrated the validity of the method. A severe accident scenario of marine reactor was simulated and analyzed, which indicates that the regional oceanic dispersion model can provide dose assessment for nuclear emergency response.

1. Introduction

The marine reactor, equipped in the offshore mobile platform, is used to provide nuclear propulsion (e.g., the nuclear-powered icebreaker and nuclear merchant ship) or electrical supply (e.g., the floating nuclear power plant, FNPP). Since the first nuclear submarine put into service in 1955, some 700 marine reactors have served at sea on various vessels worldwide and about 200 marine reactors are still in operation [1]. However, risks of accidents caused by extreme events such as storms, tsunamis, flooding, and collisions are unavoidable safety concerns to marine reactors. Differed from the land-based nuclear power plant, the marine reactor always experiences the effects of ship motions or load variation due to sea condition [2], which results in greater accident probability. There have been many emergencies and incidents in the history. From 1959 to 2007, Russian nuclear-powered submarines occurred 165 safety-related events [3]. In 1985, the submarine K431 was being refueled when a spontaneous chain reaction occurred [4], which resulted in violent explosion and radioactive contamination over a significant area in Chazhma Bay [5]. Though no indications of leakage have been observed from the sunken submarines Komsomolets [6, 7] and Kursk [8], the public are acutely sensitive to the potential long-term risks of contamination. The worst accidents, which always caused international disputes, were accompanied by serious ecological and radiological hazards.

Nowadays, 95% of global goods are moved by shipping [1]. Meanwhile, shipping contributes about 7% of the anthropogenic SO2 emissions and 15% of the NOx emissions [9]. Nuclear propulsion is the only proven emissions-free option prone to substitute conventional ship propulsion [10, 11]. In future, the application of marine nuclear propulsion would become more widespread because of the constraints on fossil fuel use in transport [12]. As an economical solution to drive development in remote but resource rich areas such as the Arctic, the first FNPP, Academician Lomonosov, is expected to be delivered to Russia in 2019 [13]. China will prioritize the development of FNPPs in the coming five years to support offshore oil and gas exploitation in the South China Sea [14]. All of them would increase the risk of radionuclides release from marine reactors to ocean, consequently.

Marine dispersion model plays a major role in obtaining emergency decision to minimize the potential hazards to humans and the environment and can serve for further dose assessments. Environmental Modeling for Radiation Safety (EMRAS) project [15], organized by the International Atomic Energy Agency (IAEA) from 2003 to 2007, included a working group on model validation for radionuclides transport in watershed river and estuarine systems. The working group studied five scenarios in aquatic systems [16], but without any marine environment consideration. In the subsequent EMRAS II program (from 2009 to 2011), even no specific working group for the aquatic system was included [17]. The Fukushima Dai-ichi accident, which caused serious marine pollution, has largely reinforced the development of marine nuclear emergency capabilities, and since 2011, a significant number of radionuclide dispersion models have been published. Kawamura et al. [18] applied the particle random-walk model SEA-GEARN to simulate the oceanic spreading of 131I and 137Cs, which became first studies published after the accident. Miyazawa et al. [19] developed a dispersion model of 137Cs based on an Eulerian transport equation. A high-resolution (1 km) regional ocean model was used to simulate the 137Cs concentration from March to May by Tsumune et al. [20]. More recently, a comprehensive intercomparison of marine dispersion models which were applied to the Fukushima Dai-chi nuclear accident were conducted in the frame of MODARIA program by IAEA [21]. At present, the MODARIA II program is being executed, which includes a working group assessing the transportation of radionuclides released in the marine environment [22].

When radionuclides released from a nonfixed marine reactor, the short-term, near-range prediction of marine radioactivity dispersion is urgent and meaningful for emergency response at the early stage of accident. It helps the rescuers to reduce or even avoid the radiation exposure with an accurate prediction of the near-range radionuclide distribution around the marine reactor under accident. However, the significant uncertainty of the accident region makes it difficult to implement rapid direct monitoring. Though there have been many marine dispersion models, most of them are (1) used for the long-term, large-scale prediction [23] and (2) applied to predict the release from fixed source (e.g., coastal nuclear power plant), not optimized for the mobile platform [21].

Compared with the existing models, computational fluid dynamics (CFD) approaches, which are particularly useful when modeling pollutant dispersion on complex terrain [24], offer a large potential to increase the accuracy of the simulation [25]. The objective of this study was to develop a finite-volume, three-dimensional regional oceanic dispersion model, based on the CFD tools in OpenFOAM [26] package and to simulate the accidental marine dispersion of radionuclides originating from the floating marine reactors accurately and swiftly.

2. Model Description

Apart from the direct discharge of high radioactive liquid wastes to the sea, atmospheric deposition of radionuclides to the sea surface is another pathway from the accident site to the ocean when severe marine reactor accident happens. In this paper, only the direct release is taken into consideration.

The facilities with the marine reactor can be located varying from the offshore waters to the deep ocean. Oceanic dispersion simulations of radionuclides basically consist of a radionuclide dispersion model and an ocean hydrodynamic model [27]. The dispersion model is employed to calculate the movement of radionuclides while the ocean hydrodynamic model is used to model background flow field. In this article, we focus on the radioactive pollutant dispersion at the early stage of accident, namely, a regional small-scale simulation in the short term.

2.1. General Fluid Dynamic Equation

Navier–Stokes equations are applied nowadays to denote collectively the conservation equations of mass, momentum, and energy [28]. CFD approaches, which are capable of modeling small-scale and detailed flow structure, simulate fluid flow by solving the full Navier–Stokes equations [23]. In Cartesian coordinates, the differential form of the generic conservation equation for the transport of a property is expressed as [29]:where is the generic diffusivity for and is the generic source or sink of . Equation (1) describes the nonlinear and coupled governing equations (e.g., momentum, continuity, salinity, and density equations) of fluid flow, which are difficult to solve. In small-scale ocean hydrodynamic model used in this work, the flow field distribution is the main factor affecting the spread of radionuclides and following reasonable assumptions must be made:(1)Seawater is treated as an incompressible, isothermal, and Newtonian fluid [30](2)Radionuclide concentration is considered as a passive scalar, which has no dynamical effect on the seawater motion [31](3)Coriolis force is neglected within small-scale(4)Concentration variations of other solutes are not taken into consideration

The simplified governing equations are sufficient for the description of incompressible, laminar flow; they are still nonlinear and subject to instability [30]. The waters where marine reactors operated are almost always turbulent [32]. Reynolds-averaged approaches to turbulence are introduced to solve the Navier–Stokes equations, which yields the Reynolds-averaged Navier–Stokes (RANS) equations [29].

2.2. Turbulence Model

The standard turbulence model [29] is selected to solve the RANS equations. The averaged continuity and momentum equations can be written as,where is the seawater density, is the mean velocity in the x, y, z direction, is the fluctuating velocity, and means the kinematic viscosity of seawater.

To close the equations, a turbulence model is introduced. Assume that the turbulence effect can be described as an increased viscosity, which leads to the eddy-viscosity model for the Reynolds stress:and the eddy-diffusion for a scalar:

In equation (4), is the turbulent kinetic energy, and is the Kronecker delta. The eddy-viscosity can be expressed as,where stands for the turbulent dissipation rate.

The rate of production of turbulence kinetic energy by the mean flow can be written as,

The standard turbulence model is given by:

This model contains five constant parameters: , , , , and are the most commonly used values.

2.3. Radionuclide Dispersion Model

The most common marine reactor type is the uranium-fueled pressurized water reactor [1]. Various kinds of radionuclides could be emitted from severe marine reactor accidents, among which 137Cs and 131I are representative of volatile, dissolved fission products [33]. Radionuclide behavior in the marine environment is complicated, and it consists of multiple processes, such as dissolution, precipitation, sorption, and suspension. It is unpractical to integrate all the physical and chemical processes in such a short paper. The influence from sorption, precipitation, and suspension is ignored. Therefore, the dispersion of radionuclide contaminants in the sea are predominantly controlled by two factors [34]:(1)Transportation by the ocean current(2)Diffusion caused by the concentration gradient and the turbulence

Turbulent flows are chaotic, diffusive causing rapid mixing. The effective diffusivity in turbulence is at least three orders of magnitude greater than the diffusivity on molecular scale [32]. Consequently, pollutant dispersion, being mainly due to turbulence, is independent of the characteristics of the pollutant itself [34]. In modeling the fluxes of pollutants in coastal waters, the corresponding general solute advection-diffusion equation can be written as,where C stands for the radionuclide concentration at point of (x, y, and z) at each time step; means the turbulent diffusivity of the radionuclide, which also varies in space and time; and represents the sources in the model.

Decay of radionuclides, described by the radioactive decay constant , can be introduced to the advection-diffusion equation by adding the term to the right hand side of equation (9). The relationship between and the half-life () is expressed as the following:

In the following section, the Fukushima nuclear accident data are used to validate the mathematical model mentioned previously. Then, radionuclide dispersion in the open sea was simulated to study spread trend after marine reactor accident.

3. Simulation Configuration

Before the simulation in this case, some configuration in OpenFOAM should be illustrated and explained firstly, which guarantee that the accuracy of simulation. Mesh generation, which impacts substantially on the quality of simulation result, would be the most time-consuming and labor-intensive part in the process of CFD simulation [35]. Consequence assessment in a timely manner is of great importance to nuclear emergency response, which requires fast mesh generation. Automated mesher such as the snappyHexMesh [36], which is part of the OpenFOAM package, has been proved to be sufficiently practical to generate the finite-volume hexahedra meshes. The snappyHexMesh operates an automated and robust meshing process, which is described clearly [35]. The controllable mesh parameters set in the input file can be specified and optimized in advance. When accident happens, snappyHexMesh can generate a high-quality mesh quickly and automatically based on the STL file of the geometry, which can accelerate emergency response. Figure 1 shows a detail of the seabed included in the mesh generated by snappyHexMesh, refined close to the reefs (). In order to obtain solutions independent of mesh, several mesh tests were performed.

The pimpleFoam solver has been secondary developed and applied. The inlet velocity is set as 2 m/s along the x direction. Time step is 0.1 s, and the calculation data are resampled each 0.5 s. In Figure 1, the nuclear merchant ship was firstly located at red dot 1 (20, 230, 10). Since 700ths to 800ths, the ship was drifted directly to the red dot 2 (40, 290, 10) gradually with a constant speed; then, the ship was stuck to the seafloor at dot 2. The release rate of 137Cs is set as a constant (1010 Bq/s) throughout the calculation.

4. Typical Marine Reactor Accident Scenario

When radionuclides released from a marine reactor after the severe nuclear accident, it is urgent to assess the radiation distribution quickly, which gives decision support in emergency rescue to minimize the potential hazards. A hypothetical accident scenario was set to investigate the radionuclides dispersion on the complex seafloor with irregular reefs.

4.1. Accident Scenario Description

The assumed accident is described as following: a nuclear merchant ship loading heavy freight encountered a fierce storm causing navigation system failure. The ship sunk to the seafloor after striking on a rock. The reactor was severely damaged, and radioactivity was detected. The schematic view of the geometric model is showed in Figure 1, and the coordinate origin is represented by a black dot. The seabed is 100 meters under the water. The nuclear merchant ship is noted as a red dot, which is gradually pushed away from point 1 to point 2 by the current. The reefs (denoted from A to G) with a variety of geometries were distributed in the seafloor, which are shown in maroon color. The size of the reefs varies from 50 m (F) to 200 m (G) and the height varies from 15 m (A) to 50 m (G).

4.2. Radionuclide Dispersion

The transportation of radionuclide mainly controlled by the flow of seawater. Figure 2 shows the velocity field around the reefs, which corresponds to a plane of 10 m above the seafloor. The flow field is turbulent and inhomogeneous, and wake flows are generated behind the blunt reefs. A reserve flow region behind the reef G is presented clearly in Figure 2, where the flow is moving toward the reef G.

The distribution of pollutant concentration in vertical plane () is shown in Figure 3, which demonstrates the moment of 300ths after the accidental release. The process of dispersion from the bottom of the sea to the surface is clear and intuitive.

The flow, which determines the distribution of radionuclides around the reefs, is affected by the size and shape of them. In Figure 3, it is obvious to see that reef A and G with sharp back edges caused a dip of the radionuclide concentration. To have a better understand of the effect of reefs and the process of dispersion, the horizontal distribution at four different moments (100 s, 300 s, 500 s, and 700 s) of radioactive pollutant concentration is shown in Figure 4. 137Cs disperses continuously in temporal and spatial.

The mobility is a significant feature of the marine reactor, which increases the difficulty in assessing the radiation distribution. Since 700ths referring to the initial time of the accident, the nuclear merchant ship drifted to dot 2 within 100 s by the ocean current and after that the ship stuck in the seabed. The changing processes of radionuclide distribution during () and after () the motion of the ship are shown in Figure 5.

The concentration at the downstream of the marine reactor changes with the motion of the source. High concentration of 137Cs was diluted gradually by the flush of currents. The process of dilution (from to in Figure 5) is inhomogeneous in temporal and spatial, for the turbulence and vortices in the flow field (see Figure 2).

5. Discussion and Conclusion

A high-resolution, small-scale oceanic dispersion model was developed to simulate the initial spreading of radionuclides discharged from the marine reactor, which is nonfixed and floating in the sea. In this study, given that the close similarity to the marine reactor accident scenario, the direct discharge of highly contaminated water from the 1F NPP reactor was used to validate the reliability and accuracy of the oceanic dispersion model. The results showed that it performs satisfactorily in simulating the distribution of 137Cs in the complex coastal terrain, and the calculation results were in good agreement with the monitoring data. Compared with the observed data, the deviations of the simulated results are mainly from the uncertainties in source estimation, the complexity in boundary conditions, and the omission of absorption by bottom sediments.

A marine reactor accident scenario, accompanied by the moving of radioactive source, was simulated to analyze the dispersion of radionuclides in the irregular ocean floor terrain. High-resolution flow field was generated, which displayed vortices around the reefs. The radionuclide concentration from both fixed source and nonfixed source is calculated. In this paper, the vertical dispersion from the seabed to the surface and the horizontal distribution of radioactive pollutant released from a fixed source were presented. In addition, the dispersion of radionuclide from nonfixed source was analyzed. It is suggested that the regional radionuclide dispersion model can be used as an important tool to study and predict the oceanic radionuclide dispersion and to provide dose assessment for emergency response. All the simulated results shown an obvious influence of ocean floor terrain to the distribution of the pollutant, which provides a scientific basis for a rapid development of rescue routes under the marine nuclear accident. Especially, the nonfixed source introduced an extrasignificant uncertainty to the monitoring of the radionuclide concentration; the simulation results supply a strong evidence to the decision of nuclear emergency response.

In this study, the processes of sorption, precipitation, and suspension between radionuclides and the environment were not taken into consideration. In addition, the process of current movements caused by winds was ignored. More studies are needed to improve the capability of the regional oceanic dispersion model.

Data Availability

No data are available on request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Project of Anhui Nuclear Emergency Command Platform Based on Pengze Nuclear Power Plant, the Chinese Academy of Sciences’ Defense Science and Technology Innovation Fund (no. CXJJ-17-M175).