Analysis of beryllium poisoning effect on liquid metal reactor with U-Be alloy fuel

Xiao-Liang Zou·Yun-Qing Bai·Ming-Huang Wang·Bing Hong

Abstract A liquid metal reactor(LMR)loaded with a fuel compound of uranium and beryllium(U-Be alloy fuel),which was cooled by a lead-bismuth eutectic alloy(PbBi),has been applied in Russian Alfa-class nuclear submarines.Because of the large amount of beryllium in the core,the reaction between the beryllium atoms and neutrons could result in the accumulation of3He and6Li,which are called the ‘‘poisoned elements''owing to their large thermal neutron capture cross section.The accumulation of neutron absorber can affect the performance of a reactor.In this study,the Super Multi-functional Calculation Program(SuperMC)code,which was developed by Institute of Nuclear Energy Safety Technology of the Chinese Academy of Sciences(INEST,CAS),was adopted to illustrate the in fluence of beryllium on an LMR.

Keywords LMR ·U-Be alloy fuel·Beryllium poisoning·Super Multi-functional Calculation Program

1 Introduction

Lead-based reactors have attracted much attention because they are intrinsically safe and have high generating ef ficiency,mainly because of the inherent properties of lead coolant[1-5].The first application of lead-based reactors in the world was Russian Alfa-class nuclear submarines in the 1960s,and the reactor,cooled by PbBi,was loaded with U-Be alloy fuel,in which the content of beryllium was as much as 90%[6].With large amount of beryllium in the reactor,it was easier to achieve a reactor core with a smaller size and lower235U loading[7,8].In addition,beryllium had another type of reaction with gamma rays as a(γ,n)reaction to produce photoneutrons with the advantage of overcoming the instrumentation dif ficulty of the ‘‘blind zone''during reactor startup after any shutdown[9].

However,because of the large amount of beryllium in liquid metal reactor(LMR)core,except for the(n,n′)reaction of beryllium atoms with neutrons as re flector and moderator,the reaction(n,α)was also signi ficant with a cross section threshold energy only above 0.74 MeV.The nuclides3He and6Li,which were called ‘‘poisoned elements''because their large thermal neutrons capture a cross section of 5327×10-24cm2and 940×10-24cm2,were generated by the interaction between the beryllium and neutron.The reactivity of the reactor continues to decrease with the increased neutron poisoning of poisoned element(3He and6Li)accumulation.Many studies have been carried out on the beryllium poisoning effect.Omar et al.studied the effect of beryllium re flector poisoning on the Syrian MNSR[10].Kalcheva et al.[11]illustrated the impact of the poisoning of the beryllium re flector on reactivity variationsofthe Belgian MTR BR2 in SCK·CEN.Andrzejewski et al.[12]discussed beryllium poisoning in the MARIA reactor.However,all the research focused on the beryllium re flector,and no research was carried out on U-Be alloy fuel.

In this study,the analysis of the beryllium poisoning effect on an LMR with U-Be alloy fuel was investigated using the Super Multi-functional Calculation Program(SuperMC)code and MATLAB program.The LMR is introduced in Sect.2.The computational tools and calculations are presented in Sect.3.The results and discussion are given in Sect.4.The conclusions are given in Sect.5.

2 LMR description

The LMR was a prototype for the Alfa-class submarine,which had held the world record of nuclear submarine speed of 42 knots.The operation of the LMR started in September 1967 at 10%full power and lasted about 250 days[13].The overall structural diagram of the LMR is shown in Fig.1.The height and diameter of the LMR core were approximately 800 and 780 mm,respectively.A radial re flector was used,and it was composed of two layers of SS and a layer of BeO.Radial thermal shields of nine concentric cylinders were positioned around each reactor core between the outer surface of the re flector region and the inner surface of the reactor pressure vessel(RPV)wall.The nine concentric cylinders and the inner surface of the RPV wall were separated by,from the re flector region surface outward,nine annular PbBi coolant channels of 3-mm thickness each.Above the core was a special shield plug through which the 3 emergency protection rods,10 control or compensation rods,and 24 emergency cooling tubes passed.There was approximately 300 mm of PbBi coolant between the top of the reactorcore and the reactor shield plug.Below the reactor were two cylindrical plates of the bottom thermal shield made of stainless steel(SS),from the lower core surface downward,with a 200 and 100 mm thickness,respectively.The lower cylindrical plate was adjacent to the inner surface of the RPV bottom and was separated from the upper cylindrical plate by a 100-mm-thick zone of the PbBi coolant.

Fig.1(Color online)LMR isometric view

The LMR with a thermal power of 70 MW was loaded with 90 kg of235U.The fuel rod pellets,which were constructed of U-Be alloy in a BeO ceramic matrix,were approximately 10 mm in diameter.Therefore,this results in as much as 90%beryllium in the fuel.Lead-bismuth eutectic alloy,in which the content of lead was 44.5 wt%and of bismuth was 55.5 wt%,was used as the coolant in the LMR.The LMR had several advantages with a leadbismuth eutectic alloy with a low melting point(approximately 125°C)and high boiling point(approximately 1670°C).It was safer and more compact than pressurized water reactors,because it operated at a low pressure with a less heavy pressure vessel.In addition,it could achieve a higher thermoelectric ef ficiency at a higher operational temperature.Typically,the distribution of LMR core materials was approximately 54%fuel,36%PbBi,and 10%SS.

3 Computational tools and calculations

3.1 Basic equations

For the operation of the LMR,the neutron absorbers6Li and3He were accumulated in the reactor and had negative effects on reactivity[14].The(n,α)reaction is described by the following equations:

The(n,α)reaction in Eq.(1)was a fast neutron reaction,and its threshold energy was approximately 0.74 MeV.The(n,α)reaction in Eq.(3)and the(n,p)reaction in Eq.(5)were thermal neutron reactions with cross sections of 945 barn and 5400 barn,respectively.In Eq.(2),the half time of6He was so short that6Li could be considered produced immediately in Eq.(1).Therefore,the number densities for each element could be obtained by solving the following equations:

The notation N in above equations was used to simplify to represent the number densities.The subscripts Be,Li,T,and He stand for,respectively,9Be,6Li,3H,and3He.The constant λTwas the tritium decay constant equal to 1.78×10-09s-1.R stood for different reaction rates of the isotopes:(n,T)for6Li,(n,p)for3He,and(n,α)reaction for beryllium:

The energy of the neutrons in the reactor was in the range of 0-20 MeV,and σi(E)was the relevant neutroninduced cross section for isotope i.Moreover,a MATLAB program[15]was used to solve the system of Eqs.(6)-(9)to get the number densities for lithium,tritium,and helium.

3.2 SuperMC calculation

In this research,the core analysis was performed by SuperMC with the Hybrid Evaluated Nuclear Data Library(HENDL).SuperMC,a general,intelligent,accurate,and precise simulation tool for the nuclear design and safety evaluation of nuclear systems,was designed by the FDS team in China[16,17]to perform comprehensive neutronics calculation,taking the radiation transport as the core and including the depletion,radiation source term/dose/biohazard,material activation and transmutation,etc.,as shown in Fig.2.SuperMC has been veri fied and validated by more than 2000 benchmark models and experiments,including ICSBEP and SINBAD[18,19].HENDL,a series of working nuclear data libraries developed by the FDS Team in China,has been applied in various applications pertaining to analysis of nuclear reactors.The accuracy of the library has been assessed and validated against various criticality safety and shielding benchmark experimental data[20].

The LMR core was modeled in three dimensions by the SuperMC code,and two different cross sections of the LMR are shown in Fig.3.The SuperMC model was built according to the real dimensions of the LMR mentioned in the International Atomic Energy Agency(IAEA)report[18],and some simpli fications were made.The SuperMC model consisted of four main parts:reactor core,radial re flector region,bottom thermal shield,and upper re flector region.The reactor core with an 800 mm height and 780 mm diameter was made of 54%U-Be,36%PbBi,and 10%SS.The parameters used in the SuperMC simulation are listed in Table 1.The input file for SuperMC included 500 cycles made of 50 inactive and 450 active cycles with 20,000 histories per cycle.In this work,the flux,the reaction rate,and keffwere calculated by the SuperMC calculation code.To normalize a criticality calculation by the steady-state power level of the reactor,the following conversion was used:P(70×106)×3.467×1010×ν(2.443)=5.9 × 1018.

4 Results and discussion

4.1 Time effects

Two calculation models were built for 250 days of fullpower operation:(1)the one-step calculation model for one step of 250 days and(2)the five-step calculation model for five steps of 250 days and each step of 50 days.For the one-step calculation model,the initial reaction rate was obtained from the SuperMC calculation.Then,the reaction rates were assumed to be constant in the operation.For the five-step calculation model, first,the initial reaction rate was the same as for the one-step model.Then,the number density was calculated for 50 days(a step)by the MATLAB program.The number density in the SuperMC input file was updated in the next-step calculation.

The results of the calculations are presented in Fig.4.For the one-step calculation model,constant reaction rates were used in the MATLAB program to make a relationship between the number density and the operation time,while reaction rates were changed for the next-step calculation in the five-step calculation model.As the figures show,forboth the one-step calculation model and the five-step calculation model,the number densities of the poisoned elements,including6Li,3He,and3H,continued increasing to the end of the period.For the one-step calculation model,the maximum values were 9.100×1018,1.575×1016,and 1.913×1018atoms/cm3,and they were 9.122×1018,1.564×1016and 1.869×1018atoms/cm3for the fivestep calculation model.Moreover,in the figures,the results of the one-step calculation model show good agreement with the five-step calculation model.This indicates that the poisoned elements accumulated in the LMR for 250 days of full-power operation showed almost no in fluence on the reaction rates.Figure 5 represents the fluxes for each step of the five-step calculation model,and they were almost the same for the poisoned element accumulation.The reaction rates were dependent on the flux in the reactor,which could explain why there was almost no difference among the results of the two calculation models.As a consequence,constant reaction rates were used for further study.

Fig.2(Color online)Functional architecture of SuperMC

Fig.3(Color online)Vertical and horizontal cross sections of LMR using the SuperMC code

Table 1 Parameters of LMR SuperMC model

Figure 6 shows the poisoned element accumulation effect on kefffor 250 days'operation of the LMR,in which clean for operation without the effect of the poisoned element accumulation,10%power for actual operation,and full power for maximum-power operation.With the accumulation of the poisoned elements,keffcontinued to decrease for both 10%power and full-power operation except clean.It can be derived that the reactivity loss caused by the accumulation of the poisoned elements6Li and3He was approximately 98 pcm,which represented approximately 1.4% ofthe excessreactivity,while 1697 pcm represented 23.0%for full-power operation.

Fig.4 Number density versus time:6Li(a),3He(b),and3H(c)

Fig.5 Each step flux for the five steps calculation model

4.2 Spatial effects

The magnitude and spectrum of the neutron flux have a direct in fluence on the poisoned element accumulation by the reaction rate.Therefore,the reaction rates depended on the position in the reactor.A detailed model was made to show the poisoned element accumulation for the different positions of the LMR core by dividing the LMR core from the radial direction into eight sections as follows:

Section 1:radius from 0 to 5 cm

Section 2:radius from 5 to 10 cm

Section 3:radius from 10 to 15 cm

Section 4:radius from 15 to 20 cm

Section 5:radius from 20 to 25 cm

Section 6:radius from 25 to 30 cm

Section 7:radius from 30 to 35 cm

Section 8:radius from 35 to 39 cm

The results of the calculations are presented in Fig.7.The poisoned elements,including6Li,3He,and3H,kept increasing in all sections,and the inside sections increased faster than the outside sections.Equation(10)indicates that the reaction rates that lead to the poisoned element accumulation depended on the flux distribution in the reactor.Figure 8 shows the flux distribution in the radial direction.As can be seen,the flux of the inside sections was bigger than that of the outside sections,and this explains the result above.

Fig.6 Poisoned element accumulation effect on keff

In Fig.7,one can also see that the difference between the inside sections was smaller than that of the outside sections.This was caused by the trend of flux variation with the radial direction,which decreased slowly first and then decreased sharply.

5 Conclusion

In this work,we focused on analysis of the beryllium poisoning effect on an LMR with U-Be alloy fuel.The beryllium poisoning effect was calculated by the SuperMC code and MATLAB program.The conclusions drawn from this study are as follows.

1. For actual operation for 250 days and 10%power,the reactivity loss caused by the accumulation of the poisoned elements6Li and3He was only 98 pcm,which represented approximately 1.4%of the excess reactivity,while 1697 pcm represented 23.0%for fullpower operation.

Fig.7(Color online)Number density versus time in different sections:6Li(a),3He(b),and3H(c)

Fig.8 Core radial distribution of flux

2. In particular,the one-step calculation model and the five-step calculation model were used to perform the poisoned element accumulation effect on the generation of the poisoned elements later.The results showed that the effect of the accumulated poison elements could be neglected,because the flux and the power of the reactor were unchanging during 250 days of fullpower operation.

3. The spatial effects were researched by simulating with models dividing the reactor core into eight sections.The results showed that the poisoned elements accumulated in inside sections increased faster than those of outside sections,and the differences of inside sections were much smaller than those of outside sections because of the radial flux distribution.

The results show that the beryllium poisoning effect on an LMR with U-Be alloy fuel cannot be neglected.In the future,the effect should be taken into consideration when designing lead-based reactors with beryllium as a moderator or re flector.

AcknowledgementsThe authors would like to show their great appreciation to other members of the FDS Team in this research.