Shallow open-loop geothermal are being increasingly used to provide affordable and low-carbon heating and cooling of buildings. The design of these plants requires the transient simulation of flow and heat transport in the aquifer to assess the propagation of the so-called thermal plume, i.e. the aquifer volume where groundwater temperature is altered by the reinjection of colder/warmer water. The propagation of thermal plumes can also result in the thermal alteration of the abstracted water, a phenomenon that is known as thermal recycling. Thermal recycling has been addressed in the literature with analytical and numerical models. Recent developments allow to assess thermal recycling with constant flow rates: however, this is not a realistic assumption for a building heating/cooling system, where the flow rate follows the thermal needs of the building and hence it has a seasonal variation.