Globally there is an increasing need to reduce the greenhouse gas emissions and increase the use of renewable sources of energy. The storage of solar thermal energy is a crucial aspect for implementing the solar energy for space heating in high latitudes, where solar insolation is high in summer and almost negligible in winter when the domestic heating demand is high. To use the solar heating during winter thermal energy storage is required. In this paper, equations representing the single U-tube heat exchanger are implemented in weak form edge elements in COMSOL Multiphysics (R) to speed up the calculation process for modelling of a borehole storage layout. Multiple borehole seasonal solar thermal energy storage scenarios are successfully simulated. After 5 years of operation, the most efficient simulated borehole pattern containing 168 borehole heat exchangers recovers 69% of the stored seasonal thermal energy and provides 971 MWh of thermal energy for heating in winter.