A new methodology for numerical simulation of geothermal down-hole heat exchangers