Geothermal energy can be promising to help the energy transition. Warm water is pumped up from a geothermal well and after it is cooled down, the cold water is reinjected into the hot aquifer underground. The temperature change could cause chemical balances to be disturbed and result in precipitation on the porous medium of the aquifer. This process of clogging of an aquifer is not desirable since it re- duces the efficiency of the production of warm water. This research is a stepping stone for modelling these reactions and clogging of the aquifer. The reactions rates depend on the temperature and the concentration of the reactants in the aquifer. For this reason the temperature distribution of the aquifer was numerically derived with finite difference methods. The velocity field is important for modelling both the temperature and the concentrations. Therefore the pressure and velocity field were modelled with use of mass conservation and the Darcy model. The models were solved numerically with finite difference methods and some results were compared to the analytical solution of the model. It was concluded that the numerical model for the velocity field is accurate when compared to analytical solu- tions of the model equations. In the model used, an increase in porosity resulted in a delay of cooling the aquifer. This is because there was no relation between porosity and permeability implemented in the model. A lower flow rate in the aquifer delays the cooling of the aquifer as well. For further research it is recommended to implement the chemical reactions into the model and use this model to look into how to prevent the clogging in the aquifer.