The thermal response of a borehole field is often described by non?dimensional response factors called gfunctions.The g?function was firstly generated as a numerical solution based on SBM (Superposition BoreholeModel). An analytical approach, the FLS (Finite Line Source), is also accepted for generating the g?function. In thiswork the potential to numerically produce g?functions is studied for circular borehole fields using the commercialsoftware COMSOL.