Monte Carlo simulation of the radiolysis of water by fast neutrons at elevated temperatures up to 350°C

Abstract : Controlling the water chemistry in a water-cooled nuclear power reactor requires understanding and mitigating the effects of water radiolysis to limit the corrosion and degradation of materials by oxidizing radiolysis products. However, direct measurement of the chemistry in reactor cores is extremely difficult due to the extreme conditions of high temperature, pressure, and mixed neutron/γ-radiation fields, which are not compatible with normal chemical instrumentation. For these reasons, theoretical models and computer simulations are essential for predicting the detailed radiation chemistry of the cooling water in the core and the impact on materials. Monte Carlo simulations were used to calculate the yields for the primary species (e[superscript -][subscript aq], H[superscript •], H[subscript 2], [superscript •]OH, and H[subscript 2]O[subscript 2]) formed from the radiolysis of neutral liquid water by mono-energetic 2-MeV neutrons and the mechanisms involved at temperatures between 25 and 350 °C. In this work, we chose 2-MeV neutron as our energy of interest since it is known as representative of a fast neutron flux in a nuclear reactor. For light water, for that chosen energy, the most significant contribution to the radiolysis comes from the first four neutron collisions that generate mostly ejected protons with energies of ~1.264, 0.465, 0.171, and 0.063 MeV, which had, at 25 °C, mean linear energy transfers (LETs) of ~22, 43, 69, and 76 keV/[micro]m, respectively. In this work, we simply neglected the radiation effects due to oxygen ion recoils. The average final fast neutron yields could be estimated as the sum of the yields for these protons after allowance was made for the appropriate weightings (by using the Eq (2) in Chapter 4) according to their deposited energy. Yields were calculated at 10[superscript -7], 10[superscript -6] and 10[superscript -5] s. Our computed yield agreed reasonably well with the available experimental data. By comparing our results with data obtained for low-LET radiation ([superscript 60]Co γ-rays or fast electrons), our computed yields for fast neutron radiation showed essentially similar temperature dependences over the range of temperature studied, but with lower values for yields of free radicals and higher values for molecular yields. In this work, we also used our Monte Carlo simulation to investigate the existence of drop of hydrated electron (e[superscript -][subscript aq]) self-reaction rate constant at 150 °C. One of the main sources of H[subscript 2] formation is the self-reaction of hydrated electrons. The temperature dependence of the rate constant of this reaction (k[subscript 1]), measured under alkaline conditions, reveals that the rate constant drops abruptly above ~150 °C. However, when this drop in the e[superscript -][subscript aq] self-reaction rate constant is included in our code for low (isolated spurs) and high (cylindrical tracks) linear energy transfer (LET), g(H[subscript 2]) shows a marked downward discontinuity at ~150 °C which is not observed experimentally. The consequences of the presence of this discontinuity in H[subscript 2] yield for both low and high LET radiation are discussed. Another reaction that might explain the anomalous increasing of H[subscript 2] yield with temperature is the reaction of H[superscript •] atoms with water previously proposed by Swiatla-Wojcik and Buxton (2005) whose rate constant is still in controversial.

Name thesis:degree_name

M. Sc.

Level thesis:degree_level

Maîtrise

Discipline thesis:degree_discipline

Sciences des radiations et imagerie biomédicale

Grantor dc:publisher

Université de Sherbrooke

Year dc:date.issued

2014

Author dc:creator

• Butarbutar, Sofia Loren

Advisor dc:contributor.advisor

• Jay-Gerin, Jean-Paul

• Radiolyse de l'eau
• Neutrons rapides
• Protons de recul
• Transfert d'énergie linéaire (TEL)
• Haute température
• Rendements de radiolyse
• Auto-réaction de l'électron hydraté
• Constante de vitesse
• Rendement de H[indice inférieur 2]
• Simulation Monte Carlo
• Water radiolysis
• Fast neutrons
• Recoil protons
• Linear energy transfer (LET)
• High temperature
• Radiolytic yields
• Self-reaction of the hydrated electron
• Rate constant
• Yield of H[subscript 2]
• Monte Carlo

Language dc:language.iso

en

Handle dc:identifier.uri

http://hdl.handle.net/11143/5887

OAI identifier oai:identifier

oai:usherbrooke.scholaris.ca:11143/5887