Nuclear fuel response to reactor accidents
From Wikipedia, the free encyclopedia
This page is devoted to a discussion of how uranium dioxide nuclear fuel behaves during both normal nuclear reactor operation and under reactor accident conditions such as overheating. Work in this area is often very expensive to conduct, and so has often been performed on a collaborative basis between groups of countries, usually under the aegis of the CSNI.
[edit] Swelling
[edit] Cladding
It is important to note that both the fuel can swell and the cladding which covers the fuel to form a fuel pin can be deformed. It is normal to fill the gap between the fuel and the cladding with helium gas to permit better thermal contact between the fuel and the cladding. During use the amount of gas inside the fuel pin can increase because of the formation of noble gases (krypton and xenon) by the fission process. If a Loss of coolant accident (LOCA) (e.g. Three Mile Island) or a Reactivity Initiated Accident (RIA) (e.g. Chernobyl or SL-1) occurs then the temperature of this gas can increase. As the fuel pin is sealed the pressure of the gas will increase (PV = nRT) and it is possible to deform and burst the cladding. It has been noticed that both corrosion and irradiation can alter the properties of the zirconium alloy commonly used as cladding, making it brittle. As a result the experiments using unirradated zirconium alloy tubes can be misleading.
According to one paper[1][5] the following difference between the cladding failure mode of unused and used fuel was seen.
Unirradated fuel rods were pressurized before being placed in a special reactor at the Japanese Nuclear Safety Reasearch Reactor (NSRR) where they were subjected to a simulated RIA transient. These rods failed after ballooning late in the transient when the cladding temperature was high. The failure of the cladding in these tests was ductile, and it was a burst opening.
The used fuel (61 GW days / ton of Uranium) failed early in the transient with a brittle fracture which was a longitudinal crack.
If has been found that hydrided zirconium tube is weaker and the bursting pressure is lower.[2]
[edit] Fuel
The fuel can swell during use, this is because of effects such as bubble formation in the fuel and the damage which occurs to the lattice of the solid. The swelling of the fuel can impose mechanical stresses upon the cladding which covers the fuel. A document on the subject of the swelling of the fuel can be downloaded from the NASA web site.[6].
[edit] Fission gas release
As the fuel is degraded or heated the more volatile fission products which are trapped within the uranium dioxide may become free. For example see [3]
A report on the release of 85Kr, 106Ru and 137Cs from uranium when air is present has been written. It was found that uranium dioxide was converted to U3O8 between about 300 and 500 °C in air. They report that this process requires some time to start, after the induction time the sample gains mass. The authors report that a layer of U3O7 was present on the uranium dioxide surface during this induction time. They report that 3 to 8% of the krypton-85 was released, and that much less of the ruthenium (0.5%) and caesium (2.6 x 10-3%) occurred during the oxidation of the uranium dioxide.[4]
[edit] Heat transfer between the cladding and the water
In a water cooled power reactor (or in a water filled used fuel store cooling pond), if a power surge occurs as a result of a reactivity initiated accident, an understanding of the transfer of heat from the surface of the cladding to the water is very useful. In a French study, metal pipe immersed in water (both under typical PWR and pond conditions), was electrically heated to simulate the generation of heat within a fuel pin by nuclear processes. The temperature of the pipe was monitored by thermocouples and for the tests conducted under PWR conditions the water entering the larger pipe (14.2 mm diameter) holding the test metal pipe (9.5 mm outside diameter and 600 mm long) was at 280 oC and 15 MPa. The water was flowing past the inner pipe at circa 4 ms-1 and the cladding was subjected to heating at 2200 to 4900 oC s-1 to simulate an RIA. It was found that as the temperature of the cladding increased the rate of heat transfer from the surface of the cladding increased at first as the water boiled at nucleation sites. When the heat flux is greater than the critical heat flux a boiling crisis occurs. This occurs as the temperature of the fuel cladding surface increases so that the surface of the metal was too hot (surface dries out) for nucleation boiling. When the surface dries out the rate of heat transfer decreases, after a further increase in the temperature of the metal surface the boiling resumes but it is now film boiling.[5]
[edit] Corrosion and other changes to materials in the reactor
[edit] Corrosion on the inside of the cladding
Zirconium alloys can undergo stress corrosion cracking when exposed to iodine,[6] the iodine is formed as a fission product which depending on the nature of the fuel can escape from the pellet.[7] It has been shown that iodine causes the rate of cracking in preasurised zircaloy-4 tubing to increase.[7]
[edit] Graphite moderated reactors
In the cases of carbon dioxide cooled graphite moderated reactors such as magnox and AGR power reactors an important corrosion reaction is the reaction of a molecule of carbon dioxide with graphite (carbon) to form two molecules of carbon monoxide. This is one of the processes which limits the working life of this type of reactor.
[edit] Water cooled reactors
[edit] Corrosion
In a water cooled reactor the action of radiation on the water forms hydrogen peroxide and oxygen. These can cause stress corrosion cracking of metal parts which include fuel cladding and other pipework. To mitigate this hydrazine and hydrogen are injected into a BWR or PWR primary cooling circuit to adjust the redox properties of the system. A review of recent developments on this topic has been published[8]
[edit] Thermal stresses upon quenching
In a loss of coolant accident (LOCA) it is thought that the surface of the cladding could reach a temperature between 800 and 1400 K, and the cladding will be exposed to steam for some time before water is reintroduced into the reactor to cool the fuel. During this time when the hot cladding is exposed to steam some oxidation of the zirconium will occur to form a zirconium oxide which is more zirconium rich than zirconia. This Zr(O) phase is the α-phase, further oxidation forms zirconia. The longer the cladding is exposed to steam the less ductile it will be. One measure of the ductility is to compress a ring along a diameter (at a constant rate of displacement, in this case 2 mm min-1) until the first crack occurs, then the ring will start to fail. The elongation which occurs between when the maximium force is applied and when the mechanical load is declined to 80% of the load required to induce the first crack is the L0.8 value in mm. The more ductile a sample is the greater this L0.8 value will be.
In one experiment the zirconium is heated in steam to 1473 K, the sample is slowly cooled in steam to 1173 K before being quenched in water. As the heating time at 1473 K is increased the zirconium becomes more brittle and the L0.8 value declines.[9]
[edit] Aging of steels
Irradiation causes the properties of steels to become poorer, for instance SS316 becomes less ductile and less tough. Also creep and stress corrosion cracking become worse. Papers on this effect continue to be published[10]
[edit] Cracking and overheating of the fuel
This is due to the fact that as the fuel expands on heating, the core of the pellet expands more than the rim. Because of the thermal stress thus formed the fuel cracks, the cracks tend to go from the centre to the edge in a star shaped pattern. A PhD thesis on the subject has been published [8] by a student at the Royal Institute of Technology in Stockholm (Sweden).
The cracking of the fuel has an effect on the release of radioactivity from fuel both under accident conditions and also when the spent fuel is used as the final disposal form. The cracking increases the surface area of the fuel which increases the rate at which fission products can leave the fuel.
The temperature of the fuel varies as a function of the distance from the centre to the rim. At distance x from the centre the temperature (Tx) is described by the equation where ρ is the power desnity (W m-3) and Kf is the thermal conductivity.
Tx = TRim + ρ (rpellet² - x²) (4 Kf)-1
To explain this a for a series of fuel pellets being used with a rim temperature of 200 °C (typical for a BWR) with different diameters and power densities of 250 Wm-3 have been modeled using the above equation. Note that these fuel pellets are rather large; it is normal to use oxide pellets which are about 10 mm in diameter.
To show the effects of different power densitys on the centreline temperatures two graphs for 20 mm pellets at different power levels are shown below. It is clear that for all pellets (and most true of uranium dioxide) that for a given sized pellet that a limit must be set on the power density. It is likely that the maths used for these calculations would be used to explain how electrical fuses function and also it could be used to predict the centreline temperature in any system where heat is released throughout a cylinder shaped object.[11]
[edit] Loss of volatile fission products from pellets
It is important to note that the heating of pellets can result in some of the fission products being lost from the core of the pellet. If the xenon can rapidly leave the pellet then the amount of 134Cs and 137Cs which is present in the gap between the cladding and the fuel will increase. As a result if the zirconium tubes holding the pellet are broken then a greater release of radioactive caesium from the fuel will occur. It is important to understand that the 134Cs and 137Cs are formed in different ways, and hence as a result the two caesium isotopes can be found at different pats of a fuel pin.
It is clear that the volatile iodine and xenon isotopes have minutes in which they can diffuse out of the pellet and into the gap between the fuel and the cladding. Here the xenon can decay to the long lived caesium isotope.
[edit] Genesis of 137Cs
Element | Isotope | decay mode | half life | direct fission yield |
---|---|---|---|---|
Sn | 137 | β | very short (<1 s) | 0.00% |
Sb | 137 | β | very short (<1 s) | 0.03% |
Te | 137 | β | 2.5 seconds | 0.19% |
I | 137 | β | 24.5 seconds | 1.40% |
Xe | 137 | β | 3.8 minutes | 1.44% |
Cs | 137 | β | 30 years | 0.08% |
Note that these fission yields were calculated for 235U assuming thermal neutrons (0.0253 eV) using data from the chart of the nuclides.[9]
[edit] Genesis of 134Cs
In the case of 134Cs the precursor to this isotope is stable 133Cs which is formed by the decay of much longer lived xenon and iodine isotopes. No 134Cs is formed without neutron activation as 134Xe is a stable isotope. As a result of this different mode of formation the physical location of 134Cs can differ from that of 137Cs.
Element | Isotope | decay mode | half life | direct fission yield |
---|---|---|---|---|
In | 133 | β | 0.18 seconds | 0.00% |
Sn | 133 | β | 1.45 seconds | 0.07% |
Sb | 133 | β | 2.5 minutes | 1.11% |
Te | 133m | β (82.5%) | 55.4 minutes | 0.49% |
Te | 133 | β | 12.5 minutes | 0.15% |
I | 133 | β | 20.8 hours | 1.22% |
Xe | 133 | β | 5.2 days | 0.00% |
Cs | 133 | - | stable (undergoes neutron activation in the core) | 0.00% |
Cs | 134 | β | 2.1 years | 6.4 x 10-6% |
Note that these fission yields were calculated for 235U assuming thermal neutrons (0.0253 eV) using data from the chart of the nuclides.[10]
[edit] An example of a recent PIE study
In a recent study used 20% enriched uranium dispersed in a range of different matrixes was examined to determine the physical locations of different isotopes and chemical elements.
- A solid solution of urania in yttrium stablised zirconia (YSZ) {Y:Zr atom ratio of 1:4}).
The fuels varied in their abiulity to retain the fission xenon, the first of the three fuels retained 97% of the 133Xe the second retained 94% while the last fuel only retained 76% of this xenon isotope. The 133Xe is a long lived radioactive isotope which can diffuse slowly out of the pellet before being neutron activated to form 134Cs. The more short lived 137Xe was less able to leach out of the pellets, 99%, 98% and 95% of the 137Xe was retained within the pellets. It was also found that the 137Cs concentration in the core of the pellet was much lower than the concentration in the rim of the pellet, while the less volatile 106Ru was spread more evenly through out the pellets.[12]
The following fuel is particles of solid solution of urania in yttrium stablised zirconia dispersed in alumina which had burnt up to 105 GW days per cubic meter.[13] The SEM is of the interface between the alumina and a fuel particle. It can be seen that the fission products are well confined to within the fuel, little of the fission products have entered the alumina matrix. The neodymium is spread throughout the fuel in a uniform manner, while the caesium is almost homogenously spread out throughout the fuel. The caesium concentration is slightly higher at two points where xenon bubbles are present. The much of the xenon is present in bubbles, while almost all of the ruthenium is present in the form of nanoparticles. It is important to note that the ruthenium nanoparticles are not always colocated with the xenon bubbles.
[edit] The release of fission products into coolant water in a Three Mile Island type accident
At Three Mile Island a recently SCRAMed core was staved of cooling water, as a result of the decay heat the core dried out and the fuel was damaged. Attempts were made to recool the core using water. According to the International Atomic Energy Agency for a 3000 MW(t) PWR the normal coolant radioactivity levels are shown below in the table, and the coolant activities for reactors which have been allowed to dry out (and over heat) before being recovered with water. In a gap release the activity in the fuel/cladding gap has been released while in the core melt release the core was melted before being recovered by water.[14]
Isotope | Normal | >20% Gap release | >10% Core melt |
---|---|---|---|
131I | 2 | 200000 | 700000 |
134Cs | 0.3 | 10000 | 60000 |
137Cs | 0.3 | 6000 | 30000 |
140Ba | 0.5 | - | 100000 |
[edit] The Chernobyl release
The release of radioactivity from the used fuel is greatly controlled by the volitility of the elements. At Chernobyl much of the xenon and iodine was released while much less of the zirconium was released. The fact that only the more volatile fission products are released with ease will greatly retard the release of radioactivity in the event of an accident which causes serious damage to the core. Using two sources of data it is possible to see that the elements which were in the form of gases, volatile compounds or semi-volatile compounds (such as CsI) were released at Chernobyl while the less volatile elements which form solid solutions with the fuel remained inside the reactor fuel.
According to the OECD NEA report on Chernobyl (ten years on)[11], the following proportions of the core inventory were released. The physical and chemical forms of the release included gases, aerosols and finely fragmented solid fuel. According to some research the ruthenium is very mobile when the nuclear fuel is heated with air.[12]
Some work has been done on TRISO fuel under similar conditions.[13]
[edit] Table of chemical data
Element | Gas | Metal | Oxide | Solid solution | Radioisotopes | Release at Chernobyl[15] | T required for 10% release from UO2 | T required for 10% release from U3O8 |
---|---|---|---|---|---|---|---|---|
Br | Yes | - | - | - | - | - | - | - |
Kr | Yes | - | - | - | 85Kr | 100% | - | - |
Rb | Yes | - | Yes | - | - | - | - | - |
Sr | - | - | Yes | Yes | 89Sr and 90Sr | 4-6% | 1950 K | - |
Y | - | - | - | Yes | - | 3.5% | - | - |
Zr | - | - | Yes | Yes | 95Zr | 3.5% | 2600 K | - |
Nb | - | - | Yes | - | - | - | - | - |
Mo | - | Yes | Yes | - | 99Mo | >3.5% | - | 1200 K |
Tc | - | Yes | - | - | - | - | - | 1300 K |
Ru | - | Yes | - | - | 103Ru and 106Ru | >3.5% | - | - |
Rh | - | Yes | - | - | - | - | - | - |
Pd | - | Yes | - | - | - | - | - | - |
Ag | - | Yes | - | - | - | - | - | - |
Cd | - | Yes | - | - | - | - | - | - |
In | - | Yes | - | - | - | - | - | - |
Sn | - | Yes | - | - | - | - | - | - |
Sb | - | Yes | - | - | - | - | - | - |
Te | Yes | Yes | Yes | Yes | 132Te | 25-60% | 1400 K | 1200 K |
I | Yes | - | - | - | 131I | 50-60% | 1300 K | 1100 K |
Xe | Yes | - | - | - | 133Xe | 100% | 1450 K | - |
Cs | Yes | - | Yes | - | 134Cs and 137Cs | 20-40% | 1300 K | 1200 to 1300 K |
Ba | - | - | Yes | Yes | 140Ba | 4-6% | 1850 K | 1300 K |
La | - | - | - | Yes | - | 3.5% | 2300 K | - |
Ce | - | - | - | Yes | 141Ce and 144Ce | 3.5% | 2300 K | - |
Pr | - | - | - | Yes | - | 3.5% | 2300 K | - |
Nd | - | - | - | Yes | - | 3.5% | 2300 K | - |
Pm | - | - | - | Yes | - | 3.5% | 2300 K | - |
Sm | - | - | - | Yes | - | 3.5% | 2300 K | - |
Eu | - | - | - | Yes | - | 3.5% | 2300 K | - |
The releases of fission products and uranium from uranium dioxide (from spent BWR fuel, burn-up was 65 GWd t-1) which was heated in a Knudsen cell has been repeated.[15] Fuel was heated in the Knudsen cell both with and without preoxidation in oxygen at c 650 K. It was found even for the noble gases that a high temperature was required to liberate them from the uranium oxide solid. For unoxidized fuel 2300 K was required to release 10% of the uranium while oxidized fuel only requires 1700 K to release 10% of the uranium.
According to the report on Chernobyl used in the above table 3.5% of the following isotopes in the core were released 239Np, 238Pu, 239Pu, 240Pu, 241Pu and 242Cm.
[edit] Degradation of the whole fuel element
It is important to note that water and zirconium can react violently at 1200 °C, at the same temperature the zirconium cladding can react with uranium dioxide to form zirconium oxide and a uranium/zirconium alloy melt.[16]
[edit] PHEBUS
In France a facility exists in which a fuel melting incident can be made to happen under strictly controlled conditions.[17][18] In the PHEBUS research program fuels have been allowed to heat up to temperatures in excess of the normal operating temperatures, the fuel in question is in a special channel which is in a toroidal nuclear reactor. The nuclear reactor is used as a driver core to irradate the test fuel. While the reactor is cooled as normal by its own cooling system the test fuel has its own cooling system, which is fitted with filters and equipment to study the release of radioactivity from the damaged fuel. Already the release of radioisotopes from fuel under different conditions has been studied. After the fuel has been used in the experiment it is subject to a detailed examination (PIE), In the 2004 annual report from the ITU some results of the PIE on PHEBUS (FPT2) fuel are reported in section 3.6.[19][20]
[edit] LOFT
The Loss of Fluid Tests (LOFT) were an early attempt to scope the response of real nuclear fuel to conditions under a Loss of Coolant Accident, funded by USNRC. The facility was built at Idaho National Laboratory, and was essentially a scale-model of a commercial PWR. ('Power/volume scaling' was used between the LOFT model, with a 50MWth core, and a commercial plant of 3000MWth).
The original intention (1963-1975) was to study only one or two major (large break) LOCA, since these had been the main concern of US 'rule-making' hearings in the late 1960s and early 1970s. These rules had focussed around a rather stylised large-break accident, and a set of criteria (eg for extent of fuel-clad oxidation) set out in 'Appendix K' of 10CFR50 (Code of Federal Regulations). However, following the accident at Three Mile Island, detailed modelling of much smaller LOCA became of equal concern.
38 LOFT tests were eventually performed and their scope was broadened to study a wide spectrum of breach sizes. These tests were used to help validate a series of computer codes (such as RELAP-4, RELAP-5 and TRAC) then being developed to calculate the thermal-hydraulics of LOCA.
- Some details of the tests can be read on-line.
[edit] See also
[edit] Contact of molten fuel with water and concrete
[edit] Water
Extensive work was done from 1970 to 1990 on the possibility of a steam explosion or FCI when molten 'corium' contacted water. Many experiments suggested quite low conversion of thermal to mechanical energy, whereas the theoretical models available appeared to suggest that much higher efficiencies were possible. A NEA/OECD report was written on the subject in 2000 which states that a steam explosion caused by contact of corium with water has four stages.[31]
- Premixing
As the jet of corium enters the water, it breaks up into droplets. During this stage the thermal contact between the corium and the water is not good because a vapour film surrounds the droplets of corium and this insulates the two from each other. It is possible for this meta-stable state to quench without an explosion or it can trigger in the next step
- Triggering
A externally or internally generated trigger (such as a pressure wave) causes a collapse of the vapour film between the corium and the water.
- Propagation
The local increase in pressure due to the increased heating of the water can generate enhanced heat transfer (usually due to rapid fragmentation of the hot fluid within the colder more volatile one) and a greater pressure wave, this process can be self-sustained. (The mechanics of this stage would then be similar to those in a classical ZND detonation wave).
- Expansion
This process leads to the whole of the water being suddenly heated to boiling. This causes an increase in pressure which can result in damage to the plant.
[edit] Recent work
Some work has been done in Japan where uranium dioxide and zirconium dioxide was melted in a crucible before being added to water. The fragmentation of the fuel which results is reported in the paper [32] which is in Journal of Nuclear Science and Technology[33]
[edit] Concrete
A review of the subject can be read at [34] and work on the subject continues to this day; in Germany at the FZK some work has been done on the effect of thermite on concrete, this is a simulation of the effect of the moltern core of a reactor breaking through the bottom of the pressure vessel into the containment.[35][36][37]
[edit] Lava flows from corium
It is possible to see in the photo shown below that the corium (molten core) will cool and change to a solid with time. It is thought that the solid is weathering with time. The solid can be described as Fuel Containing Mass, it is a mixture of sand, zirconium and uranium dixoide which had been heated at a very high temperature[38] until it has melted. The chemical nature of this FCM has been the subject of some research.[39] The amount of fuel left in this form within the plant has been considered[40]. A silicone polymer has been used to fix the contamination.
The Chernobyl melt was a silicate melt which did contain inclusions of Zr/U phases, molten steel and high uranium zircon. The lava flow consists of more than one type of material -- a brown lava and a porous ceramic material have been found. The uranium to zirconium for different parts of the solid differs a lot, in the brown lava a uranium rich phase with a U:Zr ratio of 19:3 to about 38:10 is found. The uranium poor phase in the brown lava has a U:Zr ratio of about 1:10. [16] It is possible from the examination of the Zr/U phases to know the thermal history of the mixture, it can be shown that before the explosion that in part of the core the temperature was higher than 2000 °C. While in some areas the temperature was over 2400-2600 °C.
[edit] Spent fuel corrosion
[edit] Uranium dioxide films
Uranium dioxide films can be deposited by reactive sputtering using an argon and oxygen mixture at a low pressure. This has been used to make a layer of the uranium oxide on a gold surface which was then studied with AC impedance spectrscopy.[17]
[edit] Noble metal nanoparticles and hydrogen
According to the work of the corrosion electrochemist Shoesmith[41][42] the nanoparticles of Mo-Tc-Ru-Pd have a strong effect on the corrosion of uranium dioxide fuel. For instance his work suggests that when the hydrogen (H2) concentration is high (due to the anaerobic corrosion of the steel waste can) the oxidation of hydrogen at the nanoparticles will exert a protective effect on the uranium dioxide. This effect can be thought of as an example of protection by a sacrificial anode where instead of a metal anode reacting and dissolving it is the hydrogen gas which is consumed.
[edit] References
- ^ T. Nakamura, T. Fuketa, T. Sugiyama and H. Sasajima, Journal of Nuclear Science and Technology, 2004, 41, 37.
- ^ F. Nagase and T. Fuketa, Journal of Nuclear Science and Technology, 2005, 42, 58-65
- ^ J.Y. Colle, J.P. Hiernaut, D. Papaioannou, C. Ronchi, A. Sasahara, Journal of Nuclear Materials, 2006, 348, 229.
- ^ P. Wood and G.H. Bannister, CEGB report [1]
- ^ V. Bessiron, Journal of Nuclear Science and Technology, 2007, 44, 211-221.
- ^ Gladkov, V.P. ; Petrov, V.I. ; Svetlov, A.V. ; Smirnov, E.A. ; Tenishev, V.I. ; Bibilashvili, Yu. K. ; Novikov, V.V, Atomic Energy, 1994, 75, 97-103
- ^ S.Y. Park, J.H. Kim, M.H. Lee and Y.H. Jeong, Journal of Nuclear Materials, 2007, in press (doi:10.1016/j.jnucmat.2007.03.258)
- ^ K. Ishida, Y. Wada, M. Tachibana, M. Aizawa, M. Fuse and E. Kadoi, Journal of Nuclear Science and Technology, 2006, 43, 65-76.[2]
- ^ Y. Udagawa, F. Nagase and T. Fuketa, Journal of Nuclear Science and Technology, 2006, 43, 844
- ^ K. Fukuya, K. Fujii, H. Nishioka and Y. Kitsunai, Journal of Nuclear Science and Technology, 2006, 43, 159-173.[3]
- ^ Radiochemistry and Nuclear Chemistry, G. Choppin, J-O Liljenzin and J. Rydberg, 3rd Ed, 2002, Butterworth-Heinemann, ISBN 0-7506-7463-6
- ^ N. Nitani, K. Kuramoto, T. Yamashita, K. Ichise, K. Ono and Y. Nihei, Journal of Nuclear Materials, 2006, 352, 365-371
- ^ N. Nitani, K. Kuramoto, T. Yamashita, Y. Nihel and Y. Kimura, Journal of Nuclear Materials, 2003, 319, 102.
- ^ Page 60, Generic assessment procedures for determining protective actions during a reactor accident, International Atomic Energy Agency technical document 955, published in Austria August 1997 ISSN 1011-4289
- ^ J.Y. Colle, J.-P. Hiernaut, D. Papaioannou, C. Ronchi and A. Sasahara, Journal of Nuclear Materials, 2006, 348, 229-242
- ^ S.V. Ushakov, B.E. Burakov, S.I. Shabalev and E.B. Anderson, Mater. Res. Soc. Symp. Proc., 1997, 465, 1313-1318.[4]
- ^ F. Miserque, T. Gouder, D.H. Wegen and P.D.W. Bottomley, Journal of Nuclear Materials, 2001, 298, 280-290.