From 75b20ee85894c6cdd3cb06d6413ec7ded632a095 Mon Sep 17 00:00:00 2001 From: Elisa Heim Date: Thu, 13 Nov 2025 15:43:03 +0100 Subject: [PATCH] Modify cylindrical heat source distance fixes #392 --- GHEtool/VariableClasses/Cylindrical_correction.py | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/GHEtool/VariableClasses/Cylindrical_correction.py b/GHEtool/VariableClasses/Cylindrical_correction.py index fe80447d..e516e6a6 100644 --- a/GHEtool/VariableClasses/Cylindrical_correction.py +++ b/GHEtool/VariableClasses/Cylindrical_correction.py @@ -221,7 +221,8 @@ def thermal_response_factors(self, time, alpha, kind='linear'): r_b = self.boreholes[i].r_b ii_segment = j_segment[j_segment == i_segment] h_ils = infinite_line_source(time, alpha, dis) - h_chs = cylindrical_heat_source(time, alpha, r_b, r_b) + # evaluate CHS 1 cm away from borehole wall to ensure numeric stability + h_chs = cylindrical_heat_source(time, alpha, r_b+0.01, r_b) h_ij[ii_segment, ii_segment, 1:] = ( h_ij[ii_segment, ii_segment, 1:] + 2 * np.pi * h_chs - 0.5 * h_ils)