SUMMARYThis paper presents a generalization of a computational method for the prediction of solute dispersion in fractured porous media. This method is specially useful for the prediction of subsurface flows in crystalline rocks: The model now includes a linear kinetics mechanism to represent the effects of sorption of radionuclides in the rock matrix. The method is improved in its accuracy and provides results useful for the assessment of radionuclide migration in high-level, radioactive waste repositories. Results including verification (analytical) and physical test simulations are given. These results provide a partial validation of the numerical model.