A mathematical model for mass and heat flow and a computer program have been developed to demonstrate the effect of heat released from a hypothetical radioactive waste repository on the groundwater flow regime. The model, based on the continuum approach, conceptualizes the fracture pattern and the solid blocks as two overlapping continua and consists of a set of coupled nonlinear partial differential equations. The general form of the model is three-dimensional and can treat the fluid and rock either as two separate media with a quasi-steady exchange of heat between them or as a single equivalent medium with instantaneous thermal equilibrium. Numerical solutions have been obtained by the Galerkin finite element method. Examples have been presented for topographically different locations of the repository: below a horizontal ground surface, below a hill crest, below a hillside, and close to major fractures. The effects of constant permeability and porosity or downward decreasing with depth as well as the effect of anisotropic permeability have been investigated. Solutions include the velocity field, path lines, and traveling times of water particles passing the repository and the temperature distribution. The examples have been worked out for a two-dimensional flow domain, assuming that instantaneous thermal equilibrium takes place. This assumption was found to be justified by the relatively low flow velocities that occurred in the examples. Except for the location close to a major draining fracture, heat released from the radioactive waste repository may have a significant influence on the flow regime around the repository.