  This paper is devoted to numerical simulations of the short-term behavior of
the spatial temperature distribution in a geothermal energy storage. Such
simulations are needed for the optimal control and management of residential
heating systems equipped with an underground thermal storage. We apply
numerical methods derived in our companion paper [8]. There we studied the
governing initial boundary value problem for a linear heat equation with
convection. Further, we perform extensive numerical experiments in which we
investigate properties of the spatio-temporal temperature distribution and of
its aggregated characteristics.
