Abstract Based on the one-dimensional diffusion theory, the diffusion parameters were obtained from numerical fitting by the Cyclic Initial Value (CIV) program written by MATLAB. Taking the through-diffusion experimental of cesium (stable isotope 133Cs) as an example, on the premise of ensuring accuracy, fitting calculation was used to obtain the diffusion equilibrium time of Cs in different lengths bentonite column. The fitting results of diffusion test for tritium water (HTO) and Cs are both very well. The calculation results of the equilibrium time for Cs diffusion show that the equilibrium time obtained by CIV is less than the experimental period in both groundwater (GW) and seawater (SW). In GW environment, when the sensitivity is at the maximum setting value, the diffusion coefficient of Cs in 1.5 cm bentonite column could be calculated in a shorter period of time. Compared with the experimental period, the time was shortened by 110 days. The main purpose is to verify the feasibility of CIV through the experimental data of Cs in different column lengths. The CIV program can also be used to fit and calculate the experimental data of other ongoing diffusion experiments and get the equilibrium time required for diffusion experiments. It shows that the proposed models offer the advantages of saving experimental time and reducing experimental waste.
Read full abstract