SIMMER-III Analyses of Local Fuel-Coolant Interactions in a Simulated Molten Fuel Pool: Effect of Coolant Quantity
Studies on local fuel-coolant interactions (FCI) in a molten pool are important for the analyses of severe accidents that could occur for sodium-cooled fast reactors (SFRs). To clarify the mechanisms underlying this interaction, in recent years, several experimental tests, with comparatively larger difference in coolant volumes, were conducted at the Japan Atomic Energy Agency by delivering a given quantity of water into a molten pool formed with a low-melting-point alloy. In this study, to further understand this interaction, interaction characteristics including the pressure buildup as well as mechanical energy release and its conversion efficiency are investigated using the SIMMER-III, an advanced fast reactor safety analysis code. It is found that the SIMMER-III code not only reasonably simulates the transient pressure and temperature variations during local FCIs, but also supports the limited tendency of pressurization and resultant mechanical energy release as observed from experiments when the volume of water delivered into the pool increases. The performed analyses also suggest that the most probable reason leading to such limited tendency should be primarily due to an isolation effect of vapor bubbles generated at the water-melt interface.
The disaster in March 2011 at the Fukushima Daiichi Nuclear Power Plant in Japan has caused many people to realize that severe accidents might occur, even if their probability is extremely low [1–3]. For sodium-cooled fast reactors (SFR), since the core is not generally designed at the most reactive configurations, in the past, the sequences and consequences of hypothetic severe accidents have been extensively investigated [4–7]. It is believed that, by assuming pessimistic conditions (e.g., minimal fuel discharge from the core), the accident might proceed into a transition phase where a large whole-core scale pool containing sufficient fuel capable of exceeding the prompt criticality by fuel compaction can be formed (see Figure 1) [8, 9]. Since, during the pool enlargement, a certain amount of liquid coolant might be entrapped within the pool (e.g., due to the failure of control rod guide tube), local fuel-coolant interaction (FCI) in the fuel pool is regarded as one of the various initiators that could lead to such compactive fluid motions .
Unfortunately, it is instructive to note that, over the past decades, although extensive studies (e.g., CCM , KROTOS , FARO , and TROI ), with an emphasis on ascertaining the mechanisms of steam explosion and debris bed formation in a FCI, have been conducted, most of them were performed in a fuel-injection (FI) mode, namely, by injecting or pouring melt into a coolant pool, while the knowledge and data gained from a coolant-injection (CI) mode has not been accumulated sufficiently, especially in a scenario where a certain amount of liquid coolant is entrapped within a larger melt pool. Park et al.  and Sibamoto et al.  are the few representative authors that independently performed their investigations in a CI mode within various situations (such as energetic or nonenergetic conditions). In their studies, valuable information and data on the jet penetration behavior (e.g., penetration velocity, depth, and cavity characteristics) were collected, while the knowledge on the pressure buildup by local FCIs in the pool is still limited, despite their crucial importance for accident progression evaluation.
To clarify the characteristics of molten fuel pool at various disturbances, in the past years, several series of studies, including specifically designed in-pile and out-of-pile experiments as well as numerical simulations , have been initiated at the Japan Atomic Energy Agency (JAEA). Among those studies, for the local FCIs, as reported in a recent publication , several tests were conducted by delivering a given quantity of water into a molten pool formed with a low-melting-point alloy. From the experimental analyses , it was recognized that, for a given melt and water temperature within the nonfilm boiling range, with the increasing water volume, a limited pressure buildup is achievable for both the melt and cover gas regions. In this study, to further understand the mechanisms underlying this interaction, SIMMER-III, an advanced fast reactor safety analysis code [17–19], which is believed to be a unique simulation code currently in the world for event progression analyses in a whole-core scale , is used for analyses. In Section 2, the physical models and methods of this code are concisely described with an emphasis on the main features related to FCI analysis, while, in Section 3, after a briefing of the performed experiments, the analytical geometry and conditions are determined. Finally, in Section 4, detailed interpretations and discussions on the obtained results (e.g., pressure buildup, mechanical energy release, and conversion efficiency) are further presented.
2. Models and Methods of SIMMER-III
SIMMER-III is a two-dimensional, multivelocity field, multiphase, multicomponent, Eulerian, fluid-dynamics code coupled with a structure (fuel-pin) model and a space- and energy-dependent neutron kinetics model (see Figure 2) . The fluid-dynamics portion, which constitutes about two-thirds of the code, is interfaced with the structure model through heat and mass transfer at structure surfaces, while the neutronics portion provides nuclear heat sources based on the mass and energy distributions calculated by the other code elements [18, 19].
The basic geometric structure of SIMMER-III is a two-dimensional - system (see Figure 3), although optionally an - system can be also used for various fluid-dynamics calculations. The code models five basic SFR core materials: mixed-oxide fuel, stainless steel, sodium, control (B4C), and fission gas. Since a material can exist at different physical states during core disruptive accidents (CDAs), fuel, for example, should be represented by fabricated pin fuel, liquid fuel, a refrozen crust on structure, solid particles, broken fuel pellets (also called chunks), and fuel vapor. As a result, as depicted in Figure 3, a total of 30 density components and 17 energy components are being modeled in the current code .
2.1. Fluid-Dynamics Model
In the SIMMER-III code, the following conservation equations involving fluid mass, momentum, and internal energy are solved: where subscripts , (or ), and denote the components of density, velocity, and energy, respectively; , , , , and represent macroscopic (smeared) density, volume fraction, velocity, specific internal energy, and mass transfer rate, respectively; , , , and are pressure, time, gravity, and Heaviside unit function, respectively; is momentum exchange function between and structure; and are mass transfer rate and interfield momentum exchange function between and , respectively; , , , and are rates of energy interchange due to heat transfer, melting/freezing, nuclear heating rate, and vaporization/condensation, respectively; and, finally, and are virtual mass terms in momentum and internal energy equations, respectively .
The overall fluid-dynamics solution algorithm is based on a time factorization approach originally developed for AFDM (advanced fluid-dynamics model) [19, 20], in which the intracell interfacial area source terms, heat and mass transfer, and momentum exchange functions are determined separately from intercell fluid convection. A semi-implicit procedure is used for the intercell convection. A higher-order differencing scheme is also implemented to improve the resolution of fluid interfaces by minimizing numerical diffusion.
In addition to the constitutive models, an equation-of-state (EOS) model is available to close and complete the fluid-dynamic conservation equations . The analytic EOS model employs flexible thermodynamic functions that can be readily adapted to non-SFR materials based on experimental measurements.
2.2. Flow Regime and Interfacial Area Model
In the SIMMER-III calculation, the binary contact areas must be determined to obtain the mass, momentum, and energy transfer terms for the probable 52 contact interfaces among the 8 fluid energy components (liquid fuel, liquid steel, liquid sodium, fuel particle, steel particle, control particle, fuel chunks, and vapor mixture) and 3 structure surfaces (a fuel pin and left/right can walls) (see Figure 3) . Such binary contact areas are determined based on the convectible interfacial areas and a flow regime which describes the geometry of the multiphase flow.
In SIMMER-III, the multiple flow regimes are modeled for both pool flow, in which the effect of structure is negligible, and channel flow, which is confined by structure . As depicted in Figure 4, bubbly, dispersed, and in-between transition regimes are modeled for the pool flow which is the scenario involved in this study. The upper limit of the bubbly regime and the lower limit of the dispersed regime are defined by user-specified void fractions, and , respectively, with typical values being 0.3 and 0.7. In the SIMMER-III flow-regime modeling, it is generally assumed that a cell consists of two local regions, bubbly and dispersed regions as schematically described in Figure 5. The transition regime shown in Figure 4 is defined by a combination of the two regimes as a kind of interpolated flow regime. This treatment is advantageous because the flow characteristics can be determined continuously over the entire void fraction range, without abrupt change upon flow regime transition. To better represent the transient flow, SIMMER-III also has an interfacial area convection model, which is originally proposed by Ishii : where is the interfacial area of component and is the interfacial area source term (source ).
A total of nine interfacial areas associated to moving energy components are convected. These are the surface areas of rea1 liquids (fuel, steel, and coolant) in the bubbly flow region; the surface areas of real liquids in the dispersed flow region; the surface areas of fuel and steel particles; and the surface area of bubbles in the bubbly flow region . The changes of the interfacial areas due to processes of fragmentation, flashing, turbulence-driven breakup, coalescence, and droplet or bubble production are treated as “source terms” in the interfacial area convection equation . The source term treatment in the interfacial area convection equation better represents the physical phenomena relevant to the change of interfacial areas. The modeling of the source terms is presently based on the AFDM approach except the droplet fragmentation . The present SIMMER-III models only a hydrodynamically induced process as the fragmentation mechanism.
The source term of the hydrodynamic fragmentation is modeled by the following simple linear correlation assuming that the change of interfacial area in time interval is equal to the average change of interfacial area: where the equilibrium interfacial area is related through the equilibrium radius and the volume fraction of droplets:
The equilibrium radius is assumed to be the maximum stable size of droplet governed by the hydrodynamic stability criterion related to a critical Weber number: where is defined as an instantaneous relative velocity difference between continuous and dispersed phases and is the continuous-phase density. A constant value of 12 is simply used for , although strictly this is the value for suddenly induced relative velocity. Since (3) is used to calculate the source terms of the droplet fragmentation in the vapor continuous (dispersed flow) and liquid continuous (bubbly flow) regions, the surface tension of droplets in (5) is calculated by for the dispersed flow region and for the bubbly flow region, where and are the surface tension of the continuous and dispersed liquids, respectively, in contact with a third gas or vapor phase. Equation (6b) is used to evaluate an interfacial tension between two immiscible liquid phases, although this relation is strictly appropriate to nonpolar materials . The dimensional time scale of the fragmentation process is often correlated with a dimensionless time of droplet breakup by Rayleigh-Taylor or Kelvin-Helmholtz instability: where is the density ratio of the continuous dispersed phases. Simpkins and Bales (1972) fitted gas-liquid data for complete droplet fragmentation of the Rayleigh-Taylor instability by the following functional form of the dimensionless time in terms of the instantaneous Bond number :
For liquid-liquid system, Yuen et al. (1994) correlated the dimensionless time based on the experimental data in the mercury/water system by 
Although the coefficient in the above functional form may take different values depending on the definition of the end of the breakup process, these correlations predict that the breakup time in liquid-liquid systems is smaller than that in gas-liquid systems for the same Weber number if the Rayleigh-Taylor instability mechanism dominates the fragmentation process. In SIMMER-III, (8) and (9) are applied to the droplet fragmentation in the dispersed (vapor continuous) and bubbly (liquid continuous) flow regions, respectively. A similar approach using the different breakup times between liquid-liquid and gas-liquid systems was applied to the three-velocity field code IVA3, which solves particle-number-density equations instead of IFA convection equations .
2.3. Heat and Mass Transfer Model
In addition to the interfacial area, the heat transfer coefficients between fluids are required in heat and mass transfer calculations. In general, the definition of the coefficients corresponds to the heat transfer resistance between the bulk and the interface temperatures. For the convective heat transfer processes in droplets, the heat transfer coefficients are simply calculated by using thermal conductivity and heat transfer length assuming steady-state temperature profiles. The convective heat transfer in continuous fluids is modeled based on the Nusselt number empirical correlations.
After the interfacial areas and heat-transfer coefficients are obtained, the conservation equations without convection are solved for intracell heat and mass transfer in two steps . The first step calculates the phase-transition processes occurring at interfaces, described by a nonequilibrium heat-transfer-limited model. This is a nonequilibrium process because the bulk temperature does not generally satisfy the phase-transition condition when the mass transfer occurs at the interface. The second step of mass and energy transfer is through equilibrium process occurring when the bulk temperature satisfies the phase-transition condition.
The basic concept of the nonequilibrium mass transfer model is to determine phase-transition rate by energy balance at the interface between the energy components and as shown below. The heat transfer rates from the interface are where superscript denotes the interface, represents the temperature, and and are the interfacial area and heat transfer coefficient between components and , respectively.
The net energy transfer rate from the interface is defined as
If the net heat flow, , is zero, sensible heat is exchanged without phase transition. If is positive, the energy is lost at the interface, namely, either a liquid component freezes or a vapor component condenses. Then, the mass transfer rate for this case is determined from or where is the specific enthalpy and subscripts , , and are the components of liquid, solid, and vapor phase, respectively, while and represent saturation and solidus points, respectively.
If is negative, on the other hand, the energy is gained at the interface, namely, either a solid component melts or a liquid component vaporizes. Then, the mass transfer rate for this case is determined from or where subscript denotes the liquidus point.
In the above equations (12)~(15), the heat of phase transition (the effective latent heat) is defined as the difference between the enthalpy at the interface and the bulk enthalpy of a component undergoing the phase-transition process. Although more correctly the bulk enthalpy should be replaced by the interfacial one, SIMMER-III does not calculate temperature gradients in liquid and vapors.
3. SIMMER-III Application to Local FCI
To enhance the reliability of SIMMER-III code for reactor safety evaluation, in the past, a comprehensive and systematic assessment program of the code was conducted in two steps [7, 17, 19]: Phase 1, for fundamental assessment of individual code models, and Phase 2, for integral code assessment. Phase 1 assessment applied the code to single and multiphase flow benchmark problems, small-scale experiments with reactor and simulant materials, and physical problems with known solutions, while Phase 2 assessment involved applications of the code to integral and complex multiphase situations relevant to accident analyses, such as boiling and nonboiling pool dynamics [7, 17, 25], fuel relocation and freezing , fuel-coolant interaction [27, 28], core expansion dynamics, and disrupted core neutronics. The above two-step program was successfully conducted in collaboration with the European partners KIT (formerly known as FZK), CEA, and IRSN. It was demonstrated that the SIMMER-III code can reasonably simulate the highly transient multiphase phenomena occurring during CDAs [17, 27], even for the energetic FI-mode steam explosions . Therefore, although further confirmation would be still necessary, to some degree, we do have some confidence for its utilization on the CI-mode FCIs, especially for the nonenergetic pure vapor generation and expansion events involved in this study.
3.1. Overview of Experiments
It is instructive to note that, over the past decades, extensive knowledge and findings on the CDAs have been accumulated, in response to the increase in experimental evidence and upgrading of computer codes . For instance, it becomes gradually clear that the potential of forming a stable vapor film at the fuel-sodium interface during CDAs should not be quite high [29, 30]. Therefore, similar to the experimental study , here, our attention is solely confined on the nonfilm boiling conditions and the primary interest is to ascertain whether, under such a heat-exchange mode, a limited pressurization is maintainable, even if the water mass delivered into the pool is remarkably increased.
Figure 6 shows the representative experimental facility, while the relevant instrumentations are further described in Table 1. It can be seen that most of the apparatuses are contained in an outer vessel which provides protection from FCI events occurring in an interaction vessel located in it.
The interaction vessel is a rigid cylindrical stainless steel vessel of 140 mm in inner diameter. Several thermocouples, dynamic pressure transducers, and strain gauges are installed at different positions (see Figure 6). For all experimental runs, a low-melting-point alloy (Bi 60%-Sn 20%-In 20%), with its density and melting point being around 8500 kg/m3 and 352 K, respectively, is utilized for simulating the fuel material. Before depositing on the bottom of the interaction vessel, the mass of the alloy blocks is well weighed so that the molten pool, formed later by heating the alloy blocks, can keep target heights.
Water is used as simulant material for coolant. The combination of water and the above-mentioned low-melting-point alloy has been elaborately selected after comprehensive evaluation in cost, operation easiness, and feasibility of reliable measurement as well as the capability for simulating the potential local FCIs in accident conditions for SFRs from several aspects including melt density and density of liquid coolant and coolant vapor along with the typical nonfilm boiling heat-exchange modes [16, 29, 30]. Water, with its volume carefully measured, is loaded into a glass flask connected to the bottom end of the drive rod. Currently, different sizes of flask (such as 50 and 100 mL) can be utilized depending on the corresponding water volume required in an experimental run.
An electric heater, wrapped around the external wall of the interaction vessel, is used to heat the alloy and water. To avoid the heat convection within the interaction vessel (thereby ensuring the water temperature not to rise too rapidly), an inner cylinder, with its bottom covered by a thin aluminum foil, is equipped through hanging from the top of the interaction vessel. It should be recognized that the thin aluminum foil, despite being quite effective in preventing the heat convection between the melt and water regions, is susceptible to being pierced through by the drive rod, thereby leading to neglectable interference on the water-delivering process. During the heating, as illustrated in Figure 6, air and argon flows are being delivered, to achieve a more accurate control of water temperature and prevent the oxidation of the alloy, respectively.
After the temperature of water and melt reaches target values, through operating the motor, the flask is transported downwards and ruptured by a crushing cone prepositioned on the bottom of the interaction vessel. Before the initiation of any formal runs, several preliminary tests were performed to confirm that the flask can be thoroughly destroyed . Therefore, the confining role of the flask, if any, is expected to be remarkably reduced during FCIs .
Table 2 lists the detailed condition of performed experimental tests . Based on a film boiling criterion proposed from thermal aspect by Kondo et al. , it can be judged that the initial temperature of melt and water in those cases can guarantee that the nonfilm boiling mode should be attainable, although some uncertainties in hydrodynamic difference between the jet flow and current CI-mode condition might exist. In addition, from past SCARABEE tests , it becomes also known that, during the pool enlargement, coolant might be driven into the molten pool in a scenario where a comparatively larger pressure difference is produced by the usual FCIs (i.e., the FI-mode) at the failure of control rod guide tube. Since, in such a scenario, the coolant subcooling may not decrease too much, therefore, as depicted in Table 2, relatively higher water subcooling (around 40 K) is employed for current tests.
3.2. Determination of Analytical Geometry and Condition
Prior to numerical calculation, it should be mentioned that, during experimental analyses , some experimental limitations and uncertainties have been recognized. Since some of those uncertainties, for example, the exact initial geometry of the water lump formed after the rupture of glass flask, might not be much easier to be modelled, therefore, for the calculation, a comparatively simplified scenario, which is achieved by eliminating the known experimental uncertainties as far as possible, is employed. We believe that, if the calculated results obtained from such a simpler idealized situation are basically consistent with the evidence that was confirmable by a much larger experimental database, it perhaps means that even if some uncertainties exist, the developed facility and method might be overall valid over current range of conditions.
Figure 7 illustrates the analytic geometry adopted. It is evident that, compared to the experimental facility (see Figure 6), a much larger dimension is applied. This treatment is specifically made due to the fact that, with the increasing water volume, the outer edge of water lump (i.e., the potential initial FCI sites) in the experiments might become closer to the pressure measuring points (e.g., the PS2), thereby leading to probably enlarged measurement errors. Further, we notice that although our developed experimental method (such as the water-delivery and subsequent flask-rupture process) provides an effective way to achieve well-controlled water mass as well as water and melt temperatures, it may also induce some pool fluctuations, as a result possibly impairing the accuracy of experimental analyses. To focus on the effect of water volume on local FCIs, a static water lump is assumed in our calculations.
Table 3 lists the conditions of the calculated runs in this work. To reduce the potential uncertainty caused by the initial geometry of water lump during experiments, an identical shape factor (i.e., radius/height = 1) is employed for the water column in all the runs calculated. To achieve such a shape factor more easily, as shown in Figure 7, a region of uniform fine meshes is applied within the interaction zone.
4. Results of Analyses
Since a similar tendency of the pressure and temperature history can be observed for all the runs calculated in Section 4.1, a representative one, the Run #6, in which a volume of around 29.1 cc water is assumed to be entrapped within the melt pool, is randomly selected for illustrating the typical transient behavior predicted by SIMMER-III, while, in Sections 4.2 and 4.3, the interaction characteristics between different runs (e.g., pressure buildup, mechanical energy release, and its conversion efficiency) will be further examined and compared.
4.1. Transient Pressure and Temperature History in a Specific Case
Figure 8 shows the calculated temperature and pressure history for the Run #6, while Figure 9 further shows some representative visual frames. From the change of water temperature, it is evident that, after the water and melt are directly contacted, a short premixing phase is observable. In addition, it can be also recognized that as compared to the melt temperature, the measured water temperature tends to be remarkably lower (<400 K), even after the completion of vapor expansion at around 0.06 s, which perhaps indicates that, for this case, only a limited quantity of liquid water is evaporated during the FCI occurrence. As for the melt part, possibly due to a greater dimension employed (resulting in much larger melt pool formed), different from experiments in which a noticeable decrease of melt temperature (around several tens of Kelvins) can be found , the average temperature of melt pool in the current simulation seems to be nearly constant, while, for the cover gas region, similar to experiments in , a quite rapid temperature increase (to be close to the melt temperature) exists. As pointed out in previous experimental analyses , in addition to the high-temperature steam generated, the rapid temperature increase for the cover gas should be also caused by the splashing of melt. In experimental analyses, this judgment is supportable by the posttest inspection that some solid alloy crusts can be found on the wall of the inner cylinder or/and the cover gas region of the interaction vessel, while, for the simulations, some visual evidence, for example, the ejection of melt (see Figure 9), is further observable.
From the melt pressure history shown in Figure 8, it seems that the pressure history, after the direct contact between melt and water is achieved, can be divided into three stages, namely, an initial mild change stage (or even constant for some specific runs), a rapid increase stage, and a decrease stage. By combining the visual observations from Figure 9, the following three-phase event progression might be identifiable: Phase , premixing of melt and coolant in which the vapor generation and condensation occurs simultaneously, Phase , rapid evaporation leading to short-term pressure buildup, and Phase , vapor expansion. In addition, similar to experimental results , it seems that the history of melt pressure can be also characterized by a two-phase pressure with its peak value and duration of less than 1.0 MPa and tens of milliseconds, respectively. As for the cover gas, possibly due to the significantly enlarged cover gas region (see Figures 6 and 7), a much mild pressure rising is observable as compared to previous experimental analyses .
4.2. Comparative Analyses of Pressure Buildup
Since, in actual reactor accident conditions, the compactive fluid motion is supposed to be triggered by a pressure pulse that pushes the liquid fuel away from the core center toward the core peripheries [4, 7], therefore, similar to previous experimental analyses , for the melt part, the impulse on the melt slug, estimated by integrating the transient pressure over time and the cross-sectional area of the interaction vessel, is used for the following analyses :
As for the cover gas region, the average pressure rising rate, defined as the ratio of the maximum pressure value to the time it needs to reach the maximum value, is used as well.
Figure 10 summarizes the pressure buildup characteristics for both the melt and cover gas regions as obtained from SIMMER-III calculations and experimental measurements. It is evident that, although, possibly due to the remarkable difference encountered between experimental and calculated conditions (such as system dimension, melt mass, and height of melt pool as well as the volume of cover gas region), the absolute values are quite different, both the calculated and experimental results suggest that, with the increasing water volume, a limited pressurization is attainable during the local FCIs. In addition, it can be also found that, despite those differences, a knee point at a water volume of around 20 cc seems to be attainable for both calculations and experiments. Therefore, although further elaborated analyses might be preferable, to some extent, it is imaginable that possibly only a less part of melt (e.g., those close to the water lump) contributes to the rapid heat transfer during local FCIs.
(a) Impulse on melt slug
(b) Cover gas pressure rising rate
Figure 11 shows the calculated steam generating rate from the time Phase (2) (see Figure 8) starts to the time the peak of the melt two-phase pressure is attained. It can be clearly seen that, as water volume increases, the steam generating rate seems to be increasing initially (<20 cc) and then tends to be saturated 20 cc~ (or >20 cc), as a result suggesting the presence of a mechanism that suppresses the water evaporation.
It can be easily understood that a relatively larger volume (i.e., higher heat capacity) of water delivered into the pool (e.g., >20 cc) would probably lead the initial equilibrium temperature in the mixing region of the two fluids to be much decreased (through direct contact), as a result reducing the water evaporating rate and further restraining the followed pressure buildup. However, aside from this aspect, owing to the knowledge obtained from one of our past studies , in which the fuel-to-steel heat transfer characteristics were investigated within a molten fuel/steel mixture system, we believe that, although a stable film boiling contact mode should not appear within current melt and water temperature conditions, there is the possibility that the vapor bubbles generated at the interface to some extent may also hinder the direct contact between water and melt, leading to suppressed water evaporation. In addition, the remaining liquid water, before being driven away from its original site during the latter vapor expansion phase, might play a certain role as well in inhibiting the steam evaporation (e.g., through enhanced cooling or recondensation of the steam close to its surface). However, we must notice that, physically speaking, the cooling effect caused by the remaining liquid water, should be only likely to act as a subordinate role, instead of the dominant one.
Table 4 lists the time period of the premixing phase as abstracted from SIMMER-III calculations for several cases within the saturated-pressure range. It can be easily seen that no evident variation is observable, even if the water volume entrapped within the pool is largely increased. Since, theoretically, such a time period should represent a kind of effective heat-transfer rate, to some extent, it might be deducible that the initial cooling of the surrounding melt by the liquid water released should not play a remarkable role at current melt and water temperature conditions. In other words, although further confirmation would be inevitably necessary, we believe that the most probable reason leading to the limited pressurization in our study should be due to the isolation effect caused by the vapor bubbles generated at the melt-water interface, as might be qualitatively inferable as well from the visual evidence shown in Figure 9.
4.3. Comparative Analyses of Mechanical Energy Release
In general, total mechanical energy released in a FCI contains four parts , namely, the compression work of the cover gas , the kinetic energy and potential energy of the melt slug, and the strain energy of the interaction vessel wall . In our earlier experimental analyses , it has been confirmed that similar to the FI-mode investigations , and are the largest two contributions to the total released mechanical energy, while the fraction of and seems to be negligible (<5%). Therefore, in this study, will be simply approximated to
Since for all the experimental and calculated runs performed, the FCI events are observable to be completed within an extremely short time (less than a scale of second) (see Figure 8), an adiabatic process is assumed for the calculation of : where and are the initial and maximum pressures of the cover gas, respectively; , , and are the volumes of cover gas corresponding to , , and , respectively; and is the ratio of specific heats.
Finally, the mechanical energy conversion efficiency is defined as  where is the total thermal energy used for the interaction. Although, in the past FI-mode investigations, it is more common to calculate from the exothermic (melt) side [12, 13, 32], in our earlier experimental study, regarding the CI-mode FCI , due to the unknown mass of melt involved in the interaction, based on energy conservation, the , which is actually the maximum thermal energy absorbable by water, has been tentatively employed for analyses: where , , , and are the total mass, specific heat, subcooled degree, and vaporization heat of water, respectively, while and are the specific heat of the generated steam and the saturation temperature of water, respectively.
As for the numerical analyses, might be directly estimable from the melt side using the following form: where is the melt specific heat; and are the maximum number of meshes in radial () and axial () directions, respectively; and and are the mass and temperature of melt in mesh (,), respectively, while and denote the initial and final time of local FCIs, respectively.
On the other hand, since the exact mass of steam generated is available from SIMMER-III calculations, might be approximated from the water side as well: where is the mass of steam generated.
Figure 12 shows the calculated mechanical energy release based on (17)~(19). It can be clearly seen that both the experimental and SIMER-III calculated data support that, as water volume increases, a saturated tendency of released mechanical energy is observable. Actually, it should not be surprising because the values are directly calculated based on the pressure data in both the melt and cover gas regions as shown in Figure 10. Another point that can be seen from Figure 12 is that, compared to the pressure characteristics, the difference in absolute values of seems to be less prominent. Again, it should be reasonable because, for a given mass of water, the released mechanical energy might be theoretically constant since, as aforementioned, only a less part of melt contributes to the rapid heat transfer during local FCIs, in spite of a much larger dimension employed in the SIMMER-III calculations as compared to the experimental facility.
Figure 13 shows the calculated total thermal energy based on the above three definitions ((21)~(23)). It is evident that, since only a limited part of water, as analyzed above, is probably evaporated during local FCIs, the conservative estimations based on (21), as employed in our earlier experimental analyses, would lead to rather overestimated values. In addition, the quite less difference between calculations using (22) and (23) seems to confirm our judgment in Section 4.2 that the cooling effect caused by the remaining liquid water should not be much notable.
Figure 14 illustrates the calculated conversion efficiency over varying water quantities. It can be easily seen that, as water volume increases, the calculated seems to be slightly decreasing initially and then tends to be constant. This might be due to the fact that, at much lower water volume condition, the cooling effect caused by the remaining liquid water should be extremely lower or even negligible, as a result making the calculated a little higher. On the other hand, for the calculations with much larger water volumes, the vapor bubbles generated at the melt-water interface, as mentioned above, to some extent may hinder the direct contact between melt and water, thereby resulting in saturated values. Nevertheless, it seems that, for all the runs calculated, the obtained values are varying over a comparatively narrow range (20~27%), which perhaps demonstrates again that the cooling effect induced by the remaining unevaporated water should be rather limited. Further, it is also confirmable that, possibly due to the distinct difference involved in experimental conditions (e.g., material and temperature difference and interaction mechanism (steam explosion or nonsteam explosion) as well as interaction mode (FI or CI)), whatever the water volume is, the calculated values are much higher than the reported results in several FI-mode investigations [12, 13, 32]. To further understand the difference between FI and CI modes, more calculations over a variety of parametric conditions are being planned.
5. Concluding Remarks
Motivated to provide further evidence for understanding the mechanisms of local FCIs in a molten pool, in this study, SIMMER-III, an advanced fast reactor safety analysis code, is used for analyses. Compared to experiments, the calculations were performed in a simpler idealized situation (namely, by lessening the known experimental uncertainties as far as possible). It is found that the SIMMER-III code can reasonably reflect the transient pressure and temperature tendencies during local FCIs as understood from experiments. In addition, from the comparative analyses between different cases, it is confirmable that, similar to experiments, as water volume increases, a limited pressure buildup and resultant mechanical energy release are observable. The good agreement on the evidence observed to some degree suggests that, despite the existence of some uncertainties, the developed experimental facility and method might be overall valid over current range of conditions, thereby providing us with enhanced confidence for future investigations using reactor material.
To estimate the mechanical energy conversion efficiency , several definitions for estimating the total thermal energy are introduced. It is confirmable that the calculated values in our experiments are much larger than the reported results in some FI-mode investigations.
Although a stable film boiling contact mode should not appear within current melt and water temperature conditions, the performed analyses in this work suggest that the most probable reason leading to the limited pressurization and resultant mechanical energy release should be due to an isolation effect of steam bubbles generated at the melt-water interface.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors would like to express their sincere gratitude to M. Sugaya of NESI Inc., for his technical assistance SIMMER-III running and analyses, which played an indispensable role in this study.
W. Maschek, C. Munz, and L. Meyer, “Investigations of sloshing fluid motions in pools related to recriticalities in liquid-metal fast breeder reactor core meltdown accidents,” Nuclear Technology, vol. 98, no. 1, pp. 27–43, 1992.View at: Google Scholar
Y. Tobita, K. Morita, K. Kawada, H. Niwa, and N. Nonaka, “Evaluation of CDA energetics in the prototype LMFBR with latest knowledge and tools,” in Proceedings of the 7th International Conference on Nuclear Engineering (ICONE '99), Tokyo, Japan, April 1999.View at: Google Scholar
H. Yamano, Y.-I. Onoda, Y. Tobita, and I. Sato, “Transient heat transfer characteristics between molten fuel and steel with steel boiling in the CABRI-TPA2 test,” Nuclear Technology, vol. 165, no. 2, pp. 145–165, 2009.View at: Google Scholar
T. Suzuki, K. Kamiyama, H. Yamano et al., “Evaluation of core disruptive accident for sodium-cooled fast reactors to achieve in-vessel retention,” in Proceedings of the 8th Japan-Korea Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS '12), Beppu, Japan, December 2012.View at: Google Scholar
B. W. Spencer, K. Wang, C. A. Blomquist, L. M. Mcumber, and J. P. Schneider, “Fragmentation and quench behavior of corium melt streams in water,” Tech. Rep. NUREG/CR-6133, ANL-93/32, Argonne National Laboratory, Lemont, Ill, USA, 1994.View at: Google Scholar
J. H. Kim, B. T. Min, I. K. Park, and S. W. Hong, “Triggered steam explosions with the corium melts of various compositions in a two-dimensional interaction vessel in the troi facility,” Nuclear Technology, vol. 176, no. 3, pp. 372–386, 2011.View at: Google Scholar
H. S. Park, N. Yamano, Y. Maruyama, K. Moriyama, Y. Yang, and J. Sugimoto, “Study on energetic fuel-coolant interaction in the coolant injection mode of contact,” in Proceedings of the 6th International Conference on Nuclear Engineering (ICONE '98), San Diego, Calif, USA, May 1998.View at: Google Scholar
Y. Sibamoto, Y. Kukita, and H. Nakamura, “Visualization and measurement of subcooled water jet injection into high-temperature melt by using high-frame-rate neutron radiography,” Nuclear Technology, vol. 139, no. 3, pp. 205–220, 2002.View at: Google Scholar
Y. Tobita, S. Kondo, H. Yamano et al., “The development of SIMMER-III, an advanced computer program for LMFR safety analysis, and its application to sodium experiments,” Nuclear Technology, vol. 153, no. 3, pp. 245–255, 2006.View at: Google Scholar
K. Kamiyama, Y. Tobita, T. Suzuki, and K. Matsuba, “Development of an evaluation methodology for the fuel-relocation into the coolant plenum in the core disruptive accident of sodium-cooled fast reactors,” in Proceedings of the 22nd International Conference on Nuclear Engineering (ICONE '22), vol. 3, p. 12, Prague, Czech, July 2014.View at: Publisher Site | Google Scholar
H. Yamano, S. Fujita, Y. Tobita et al., SIMMER-III: A Computer Program for LMFR Core Disruptive Accident Analysis, Version 3. A Model Summary and Program Description, Japan Nuclear Cycle Development Institute, 2003.
W. R. Bohl, D. Wilhelm, F. R. Parker, J. Berthier, L. Goutagny, and H. Ninokata, “AFDM: an advanced fluid dynamics model, vol. I: scope, approach and summary,” Los Alamos National Laboratory LA-11692-MS, 1990.View at: Google Scholar
V. P. Carey, Liquid-Vapor Phase-Change Phenomena: An Introduction to the Thermophysics of Vaporization and Condensation Processes in Heat Transfer Equipment, Hemisphere, Washington, DC, USA, 1992.
H. Jacobs, “Analysis of large-scale melt-water mixing events,” in Proceedings of the CSNI Specialists Meeting on Fuel-Coolant Interactions, Santa Barbara, Calif, USA, January 1993.View at: Google Scholar
G. Kayser, J. Charpenel, and C. Jamond, “Summary of the SCARABEE-N subassembly melting and propagation tests with an application to a hypothetical total instantaneous blockage in a reactor,” Nuclear Science and Engineering, vol. 128, no. 2, pp. 144–185, 1998.View at: Google Scholar