Transient thermal response in nuclear waste repositories

In this paper, it is shown that a previously reported non-linear, one-dimensional, theoretical approximation simplifies - from a computational point of view - the calculation of the time-decay temperature field in nuclear waste repositories (NWR). This conclusion has been reached after solving, by using the control volume numerical method, the full three dimensional, transient, non-linear heat diffusion equation. The transient thermal field in a rock salt repository, is analytically solved and numerically predicted, along 100 years, after the disposal of a high-level waste (HLW). The nuclear waste, with a half-life of 32.9 years, releases an exponentially time dependent heat flux with 12 Wm (-2) as the initial thermal load. Two cases are studied, in the first one it is assumed that the conductivity (k) and the volumetric heat capacity (rho)Cp of the host rock (diffusion domain) remain constant (linear case), whereas in the second one, a more realistic situation is analysed. In this last case, the conductivity of the rock salt varies as a function of the temperature field and the product (rho) x Cp remains constant (non-linear case). In order to observe the effect of the salt conductivity (constant or variable) on the repository temperature distribution, a comparison of both cases is performed. It is concluded, that the theoretical model, which provides an analytical solution of the thermal fields may be a powerful low cost method for design purposes.