Numerical simulation of a Hypothetical Core Disruptive Accident in a small-scale model of a nuclear reactor

In the case of a Hypothetical Core Disruptive Accident (HCDA) in a Liquid Metal Fast Breeder Reactor, it is assumed that the core of the nuclear reactor has melted partially and that the chemical interaction between molten fuel and liquid sodium has created a high-pressure gas bubble in the core. The violent expansion of this bubble loads and deforms the reactor vessel and the internal structures, thus endangering the safety of the nuclear plant.The MARA 10 experimental test simulates a HCDA in a 1/30-scale mock-up schematising a reactor block. In the mock-up, the liquid sodium cooling the reactor core is replaced by water and the argon blanket laying below the reactor roof is simulated by an air blanket. The explosion is triggered by an explosive charge.This paper presents a numerical simulation of the test with the EUROPLEXUS code and an analysis of the computed results. In particular, the evolution of the fluid flows and the deformations of the internal and external structures are analysed in detail. Finally, the current computed results are compared with the experimental ones and with previous numerical results computed with the SIRIUS and CASTEM-PLEXUS codes.     

Numerical interpretation of the 8 experiment simulating a hypothetical core disruptive accident

In the case of a hypothetical core disruptive accident (HCDA) in a liquid metal fast breeder reactor (LMFBR), it is assumed that the core of the nuclear reactor has melted partially and that the chemical interaction between the molten fuel and the liquid sodium has created a high-pressure gas bubble in the core. The violent expansion of this bubble loads and deforms the reactor vessel, thus endangering the safety of the nuclear plant. The experimental test 8 simulates the explosive phenomenon in a mock-up included in a flexible vessel with a flexible roof. This paper presents a numerical simulation of the test and a comparison of the computed results with the experimental results and previous numerical ones.     

Computation of a core disruptive accident in the MARS mock-up

A hypothetical core disruptive accident in a liquid metal fast breeder reactor (LMFBR) results from the interaction between molten fuel and liquid sodium, which creates a high-pressure bubble of gas in the core. The violent expansion of this bubble loads and deforms the vessel and the internal structures. The MARS experimental test simulates a HCDA in a small-scale mock-up containing all the significant internal components of a fast breeder reactor. The mock-up is filled with water, topped by an argon blanket, and the explosion is generated by an explosive charge.This paper presents a numerical simulation of the test with the EUROPLEXUS code. The top closure is represented by massive structures and the main internal structures are described by shells. The current numerical results are described and compared with the experimental ones, and previous computations with the CASTEM-PLEXUS code.     

A fundamental approach to specify thermal and pressure loadings on containment buildings of sodium cooled fast reactors during a core disruptive accident

Reactor Containment Building (RCB) is the ultimate barrier to the environment against activity release in any nuclear power plant. It has to be designed to withstand both positive and negative pressures that are credible. Core Disruptive Accident (CDA) is an important event that specifies the design basis for RCB in sodium cooled fast reactors. In this paper, a fundamental approach towards quantification of thermal and pressure loadings on RCB during a CDA, has been described. Mathematical models have been derived from fundamental conservation principles towards determination of sodium release during a CDA, subsequent sodium fire inside RCB, building up of positive and negative pressures inside RCB, potential of in-vessel sodium fire due to failed seals and temperature evolution in RCB walls during extended period of containment isolation. Various heating sources for RCB air and RCB wall and their potential have been identified. Scaling laws for conducting CDA experiments in small-scale water models by chemical explosives and the rule for extrapolation of water leak to quantify sodium leak in reactor are proposed. Validation of the proposed models and experimental simulation rules has been demonstrated by applying them to Indian prototype fast breeder reactor. Finally, it is demonstrated that in-vessel sodium fire potential is very weak and no special containment cooling system is essential.     

Safety Strategy of JSFR Eliminating Severe Recriticality Events and Establishing In-Vessel Retention in the Core Disruptive Accident

In the Japan Sodium Cooled Fast Reactor (JSFR) design, elimination of severe power burst events in the Core Disruptive Accident (CDA) is intended as an effective measure to ensure retention of the core materials within the reactor vessel. The design strategy is to control the potential of excessive void reactivity insertion in the initiating phase by selecting appropriate design parameters such as maximum void reactivity on one hand, and to exclude core-wide molten-fuel-pool formation, which has been the main issue of CDA, by introducing an inner duct on the other hand. The effectiveness of these measures is evaluated based on existing experimental data and computer simulation with validated analytical tools. It is judged that the present JSFR design can exclude severe power burst events. Phenomenological consideration of general characteristics and preliminary evaluations for the long-term material relocation and cooling phases gave the perspective that in-vessel retention would be attained with appropriate design measures.     

Liquid entrainment by an expanding core disruptive accident bubble—a Kelvin/Helmholtz phenomenon

The final stage of a postulated energetic core disruptive accident (CDA) in a liquid metal fast breeder reactor is believed to involve the expansion of a high-pressure core-material bubble against the overlying pool of sodium. Some of the sodium will be entrained by the CDA bubble which may influence the mechanical energy available for damage to the reactor vessel. The following considerations of liquid surface instability indicate that the Kelvin–Helmholtz (K–H) mechanism is primarily responsible for liquid entrainment by the expanding CDA bubble. First, an instability analysis is presented which shows that the K–H mechanism is faster than the Taylor acceleration mechanism of entrainment at the high fluid velocities expected within the interior of the expanding CDA bubble. Secondly, a new model of liquid entrainment by the CDA bubble is introduced which is based on spherical-core-vortex motion and entrainment via the K–H instability along the bubble surface. The model is in agreement with new experimental results presented here on the reduction of nitrogen-gas-simulant CDA bubble work potential. Finally, a one-dimensional air-over-water parallel flow experiment was undertaken which demonstrates that the K–H instability results in sufficiently rapid and fine liquid atomization to account for observed CDA gas-bubble work reductions. An important byproduct of the theoretical and experimental work is that the liquid entrainment rate is well described by the Ricou–Spalding entrainment law.     

Two-dimensional numerical simulation of single bubble rising behavior in liquid metal using moving particle semi-implicit method

Gas-lift pump in liquid metal cooling fast reactor (LMFR) is an innovative conceptual design to enhance the natural circulation ability of reactor core. The two phase flow characteristics of gas–liquid metal make significant improvement of the natural circulation capacity and reactor safety. It is important to study bubble flow in liquid metal. In present study, the rising behaviors of a single nitrogen bubble in 5 kinds of common stagnant liquid metals (lead bismuth alloy (LBE), liquid kalium (K), sodium (Na), potassium sodium alloy (Na–K) and lithium lead alloy (Li–Pb)) and in flowing lead bismuth alloy have been numerically simulated using two-dimensional moving particle semi-implicit (MPS) method. The whole bubble rising process in liquid was captured. The bubble shape, rising velocity and aspect ratio during rising process of single nitrogen bubble were studied. The computational results show that, in the stagnant liquid metals, the bubble rising shape can be described by the Grace's diagram, the terminal velocity is not beyond 0.3 m/s, the terminal aspect ratio is between 0.5 and 0.6. In the flowing lead bismuth alloy, as the liquid velocity increases, both the bubble aspect ratio and terminal velocity increase as well. This work is the fundamental research of two phase flow and will be important to the study of the natural circulation capability of Accelerator Driven System (ADS) by using gas-lift pump.     

液态金属内单个气泡上升行为的MPS法数值模拟

液态金属冷却核反应堆采用气泡泵的概念设计来提升堆芯自然循环能力。液态金属内气液两相流动特征将直接影响核反应系统一回路的自然循环能力及堆芯安全。本研究通过采用移动粒子半隐式（MPS）方法,对液态金属中单个上升气泡的气泡动力学行为进行数值模拟。分析了铅铋合金中3种初始直径不同的单个氮气泡在上升过程中的气泡形状和速度的变化趋势;对比了初始直径相同的单个氮气泡在液钾、液钠、铅铋合金、钾钠合金和锂铅合金5种液态金属中的上升行为;同时将模拟得到的气泡形状与Grace经验关系图进行了对比,验证了MPS方法数值模拟结果的正确性。     

Thermal transient and the temperature profile in a HELICA mock-up simulated by a new finite element homogenous model

This paper deals with a numerical approach for simulating the thermal and mechanical behaviour of pebble beds used as breeder and neutron multiplier in breeding blanket of nuclear fusion reactor. The model of the pebble beds is based on the results of a theoretical and experimental research activity performed by the Authors on ceramic pebble beds (lithium ortosilicate and lithium metatitanate). The results of this activity permitted to determine the effective thermal conductivity of the beds, versus the temperature and the axial pressure and to implement a homogenous model of pebble bed in a FEM code.This paper illustrates an application of the implemented model, considering pebble beds under several cycles of heating and cooling. The examined geometry corresponds to the HELICA mock-up tested by ENEA in the research centre Brasimone. The experimental tests performed on HELICA have been used as a benchmark problem in order to assess the different approaches for simulating pebble beds. In this paper, the simulations performed with two-dimensional models are illustrated. Moreover the numerical results are compared with the experimental ones. Finally, a discussion on results obtained by other authors involved in the benchmark is reported.     

液态锂铅合金的氚释放行为

为了完成聚变堆液态锂铅包层鼓泡提氚系统的工程设计和建造,以金属与氢的作用理论为基础,建立了氚从液态锂铅中的动力学释放行为的数学模型.计算和分析了温度、饱和器氚分压、氦流量对解吸器顶部气相中的氚分压的影响以及氚在液态锂铅中的传质系数、解吸率和吸附率.结果表明:在633～723 K的解吸温度范围内,氚从液态锂铅到气相的整个释放过程虽然包含了氚在熔融合金气泡中的扩散与对流、氚通过与气-液界面相连合金层的扩散、在界而发生的氚原子重组多相反应、氚通过气相边界层的扩散和气相中氚的扩散与对流5个子过程,但起决定作用的是氚在合金内的扩散和气.液界面的多相反应重组,其他子过程意义不大.     

Feasibility study of large MOX fueled FBR core and applicability of various coolants and fuels aimed at the Self-Consistent Nuclear Energy System

The potential of a large MOX fueled fast breeder reactor (FBR) is evaluated with regard to its ability to transmute radioactive nuclides and its safety when incorporated in the self-consistent nuclear energy system (SCNES). The FBR's annual production amounts of selected long -lived fission products (LLFPs), Se-79, Tc-99, Pd-107, I-129, Cs-135 and Sm-151, can be transmuted by using a radial blanket region and part of a lower axial blanket region without any significant impact on the reactor's nuclear and safety characteristics. The other LLFPs are confined in the system. The hazard index level of the LLFPs per one ton of spent fuel from the system after 1000 years is as small as that of a typical uranium ore. To realize self-controllability (passive safety), the proposed FBR core concept employs gas expansion modules and a sodium plenum above the core. To realize self-terminability, even if MOX fuel melting should cause a core compaction, re-criticality of the core can be avoided by a fuel dilution and relocation module. The results show the MOX fueled FBR core has potential applicability to the SCNES. The fundamental applicability of various coolants and fuels is evaluated based on neutron balance toward the final goal of the ideal SCNES. The results show that gas coolant has a potential for increasing the transmutation efficiency of LLFPs. And an improved SCNES with several conventional FBRs and a FP transmutation reactor is also studied.     

Comparison of different LMFBR primary containment codes applied to a benchmark problem

The CONT benchmark calculation exercise is a project sponsored by the Containment Loading and Response Group, a subgroup of the Safety Working Group of the Fast Reactor Coordinating Committee — CEC. A full-size typical pool-type LMFBR undergoing a postulated Core Disruptive Accident (CDA) has been defined by Belgonucleaire-Brussels under a study contract financed by the CEC and has been submitted to seven containment code calculations. The results of these calculations are presented and discussed in this paper.     

Conceptual Design Study of 180 MWt Small-Sized Reduced-Moderation Water Reactor Core

Conceptual design of a Small-sized Reduced-Moderation Water Reactor (S-RMWR) core, which has the thermal output of 180 MW, the conversion ratio of 1.0 and the void reactivity coefficient of negative value, has been constructed. S-RMWR is a technology demonstration reactor which also conducts material and fuel testing for commercial use of Reduced-Moderation Water Reactor (RMWR) in large-scale power plants. It has a very tight triangular fuel rod lattice and a high coolant void fraction. The RMWR core axially has two short and flat uranium plutonium mixed oxide (MOX) regions with an internal blanket region in between, in order to avoid a positive void reactivity coefficient. The MOX regions are sandwiched between upper and lower blanket regions, in order to increase a conversion ratio.

In this small reactor core, leakage of neutrons is expected to be larger than in a large core. Therefore, a core design concept different from that for a large core is necessary. Core burnup calculations and nuclear and thermal-hydraulic coupled calculations were performed in the present study with SRAC and MOSRA codes. MVP code was also used to obtain control rod worth. Because of its large neutron leakage, keeping the void reactivity coefficient negative is easier for S-RMWR than RMWR. Thus, the heights of MOX region can be taller and the plutonium enrichment can be lower than in RMWR. On the other hand, to achieve the conversion ratio of 1.0, radial blanket and stainless steel reflector assemblies are necessary, whereas they are not needed for RMWR.     

Experimental Study of Nuclear Characteristics of Large Axially Heterogeneous Core Using Fast Critical Assembly

The physics characteristics of large axially heterogeneous liquid-metal fast breeder reactors (LMFBRs), particularly the parameters for use in design and safety assessment, were examined using the JAERI fast critical assembly facility, arranged in Assembly XH-1, a partial mock-up of axially heterogeneous LMFBR. The properties measured were (1) criticality, (2) reaction rates and reaction rate ratios, (3) material sample worths, (4) sodium-void worths and (5) B4C control rod worths.

The results were compared with those of prior experiments with assemblies representing conventional homogeneous core. Confirmation was obtained of the typical nuclear characteristics attributed to axially heterogeneous LMFBRs, including flattening of the axial distribution of power and of the differential worth of control rod, as also lower sodium void worth.

Theoretical analyses paralleling the experiments, using JENDL-2 cross section library and JAERI standard calculation code system for fast reactor neutronics, resulted in some discrepancies, particularly for the internal blanket, in respect of plutonium sample worth, fission rate and fission rate ratio.     

Experimental investigation on bubble characteristics entrained by surface vortex

The cover gas entrainment at the free surface of sodium coolant becomes one of the significant issues according to the compact sizing of reactor vessel in the latest reactor design. In the present study, some basic experiments for the gas entrainment due to the surface vortex were performed in order to obtain the fundamental knowledge about the entrained bubble size. Distributions of entrained bubble diameters in several experimental conditions were obtained from bubble images using an image processing technique. Velocity fields around vortices and surface dimple shapes (gas cores) due to surface vortices were measured to grasp those influences on bubble shapes. The result showed that mean equivalent diameters of bubbles were varied from 1.3 to 2.1 mm in the range of present experimental conditions. The bubble sizes were influenced by the thickness of gas core.     

ARMA models in the analysis of temperature noise measurements in a Super-Phenix subassembly mock-up

The noise analysis methods play an important role in the early, reliable detection of local cooling disturbances in a fast reactor subassembly such as sodium boiling or blockage, which are considered among the initiating events of major disruptive core accidents. In this paper we apply the Box and Jenkins auto-regressive moving-average ARMA models to the analysis of several temperature time-series measured by the Commissariat à l'Energie Atomique in the course of the CFNa experiment carried out at the nuclear center of Grenoble. The source of data is a thermocouple placed at the outlet of a Super-Phenix subassembly mock-up. The analysis shows that a simple ARMA (3,2) model adequately accounts for the observed fluctuations. This model provides methods for a continuous, in situ estimate of the thermocouple time constant, for the identification of a suitable boiling and blockage indicator and for the detection in real time of suddenly occurring disturbances.     

The design and experimental validation of an emergency core cooling system for a pool type research reactor

This paper presents the design of the Emergency Core Cooling System (ECCS) for the IEA-R1m pool type research reactor. This system with passive features, uses sprays installed above the core. The experimental program performed to define system parameters and to demonstrate to the licensing authorities, that the fuel elements limiting temperature is not exceeded, is also presented. Flow distribution experiments using a core mock-up in full-scale were performed to define the spray header geometry and spray nozzles specifications as well as the system total flow rate. Another set of experiments using electrically heated plates simulating heat fluxes corresponding to the decay heat curve after full power operation at 5 MW was conducted to measure the temperature distribution at the most critical position. The observed water flow pattern through the plates has a very peculiar behavior resulting in a temperature distribution which was modelled by a 2D energy equation numerical solution. In all tested conditions the measured temperatures were shown to be below the limiting value.     

Feasibility study of large MOX fueled FBR core aimed at the Self-Consistent Nuclear Energy System

The potential of a MOX fueled fast breeder reactor (FBR) is evaluated with regard to its ability to transmute radioactive nuclides and its safety when incorporated in the so-called self-consistent nuclear energy system (SCNES). The FBR's annual production amounts of selected long-lived fission products (LLFPs), Se-79, Tc-99 Pd-107, I-129, Cs-135 and Sm-151, can be transmuted by using a radial blanket region and a part of a lower axial blanket region without any significant impact on its nuclear and safety characteristics. The other LLFPs are confined in the system. The hazard index level of the LLFPs per one ton of spent fuel from the system after 1000 years is as small as that of a typical uranium ore. To realize self-controllability (passive safety), the proposed FBR core concept employs gas expansion modules and sodium plenum above the core. To realize self-terminability, even if MOX fuel melting should cause a core compaction, recriticality of the core can be avoided by a fuel dilution and relocation module. The results show the MOX fueled FBR core has potential applicability to the SCNES. With the final goal of the ideal SCNES, fundamental applicability of various coolants and fuels is evaluated based on neutron balance. It is shown that the harder the core spectra is, the larger the potential for transmuting LLFPs would be.     

Design and analysis on tritium system of multi-functional experimental fusion–fission hybrid reactor (FDS-MFX)

As early application of fusion technology, the fusion–fission hybrid systems/reactors could be used to transmute long-lived radioactive waste and produce fissile nuclear fuel. A fusion–fission hybrid reactor named FDS-MFX was designated for checking and validating the DEMO reactor blanket relevant technologies. The reactor design is based on easy-achieved plasma parameters extrapolated from the successful operation of tokamaks and the subcritical blanket is designed based on the well-developed technologies of fission reactors. In this contribution, a concept of the tritium system was designed for the FDS-MFX: the tritium was extracted from LiPb by the helium purge gas which contains a small amount of hydrogen gas, then the impurity gas was removed by cold trap, finally tritium was separated from hydrogen isotope by the cryogenic distillation and supply to reactor core. On the basis of data obtained by present design and experimental research, the system parameters were presented and discussed in detail. The results preliminarily demonstrated the engineering feasibility of the design.