Proper generalized decomposition for geothermal applications