Computational Modeling of Shallow Geothermal Systems