A numerical model was developed to simulate BHE induced soil freezing.
Latent heat from freezing slows down the BHE outlet temperature drop.
COP corrected boundary condition produces more realistic estimation.
Longer BHE do not always lead to a better financial performance over 30 years.
Total cost of GSHPS depends heavily on electricity price and interest rate.