Time‐dependent shape functions for modeling highly transient geothermal systems