Institute of Metals Division - Mathematics of the Thermal Diffusion of Hydrogen in Zircaloy-2

By means of mathematical solutions to the appropriate diffusion equations, we describe the kinetics of the thermal diffusion of hydrogen in Zircaloy-2 for the various temperatures and concentrations encountered in a heavy water moderated reactor. When the hydrogen concentration is below terminal solid solubility only the a Phase is present. Redistributions are then described in terms of the characteristic functions of the difhsion equation. For higher concentrations both the a and 6 phases are present. We assume the two phases to be always in equilibrium. For moderately small hydrogen concentrations exact solutions of the two-phase equation approach the the approximate solutions derived by Sawatzky and for all concentrations the exact solutions exhibit the qualitative features of his result: the two-phase concentration increases with time, everywhere; in the absence of a hydrogen current at the hot end of the sample an a-phase region always exists there; the interface of the a + 6 , a-phase boundary moves toward the cold md of the sample and the hydrogen concentration is discontinuous at the interface. Simultaneous solulions of the a and a + 6 hydrogen distributims and of the concomitant interface motion are obtained and compared to the observations of Sawatzky and Markowitz. The kinetics of the hydrogen diffusion process are shown to lead to m apparent heat of transport of the a phase which is lower than the actual value (even for samples with long anneals) thus resolving at least partially the disparity between experimental measurements of this quantity. A number of recent papers1"4 have reported measurements on the diffusion and redistribution of hydrogen in Zircaloy-2 under temperature and concentration gradients. These studies were instigated by problems arising from increasing use of Zircaloy-2 as a fuel element cladding material in pressurized-water power-reactors. The Zircaloy-2 picks up hydrogen during the operation of the reactor: the consequent precipitation of zirconium hydride in the Zircaloy-2 has pronounced effects on its mechanical properties. The Purpose of the present Paper is to describe the kinetics of the thermal diffusion of hydrogen in Zircaloy-2 for the various hydrogen concentrations and temperatures likely to be encountered in reactors. When the hydrogen in the Zircaloy-2 is entirely in the solid solution phase (a phase), the differential equation for thermal diffusion is well known and the redistribution can be described by standard mathematical methods some of which are given in Section 11 below. The treatment of the a + 6 (hydride) region differs from the earlier treatment of sawatzky3 by taking into account several modifications of the two-phase diffusion equations as suggested to us by Markowitz4 and Kidson.5 Even with the modifications the results obtained in our more complete treatment are substantially the same as those found by Sawatzky. As we show, the difference between the results of sawatzky3 and kIarkowitz4 is largely due to the difference in the geometry of their experiments. In the next section we derive the modified two-phase equation for arbitrary geometry. The exact solution of this equation is given in Section III for linear and cylindrical geometry. It is shown that Sawatzky's approximate solution is quite accurate for almost all the temperatures and hydrogen concentrations which are actually encountered. In Sawatzky's approximate theory the hydrogen concentration everywhere in the two-phase region increases continuously. As his paper pointed out, this result, together with hydrogen conservation, implies that in a sample with no hydrogen current flowing into the hot end a two-phase region is always accompanied by a single-phase region at the hot end of the sample. The net hydrogen gain in the two-phase region is supplied by the decrease of hydrogen in the single-phase region and by the movement of the (a, a + 6) boundary toward the cold end of the sample. Hydrogen conservation at the (a, a + 6) boundary leads to a discontinuity in the concentration and its derivative there. In Section IV it is shown how these qualitative features of the thermal diffusion kinetics arise from the simultaneous solution of the a-phase differential equation and a + 6 phase equation. Methods derived by us previously for the solution of the redistribution in both phases and the accompanying motion of the boundary are used to obtain approximate solutions valid for large times. The solutions account well for the observed redistributions of Sawatzky and Markowitz. It is shown how the continuing motion of the (a, a + 6) boundary leads to measured values of the heat of transport lower than the actual value, even for specimens with relatively long anneals.