Efficient finite element formulation for geothermal heating systems. Part I: steady state | Synapse