Pith. sign in

REVIEW 3 major objections 6 minor 37 references

Novel discretization method to calculate g-functions of vertical geothermal boreholes with improved accuracy and efficiency

T0 review · 3 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Geothermal borehole g-functions can be computed 20–200× faster by solving the underlying Fredholm integral equation with Gauss-Legendre quadrature.

desk verdict Genuinely faster g-function computation via Gauss-Legendre quadrature, but the accuracy improvement is measured against the method's own regularized high-order solution and needs independent validation. read the letter →

arxiv 2508.11154 v1 pith:KDA33RM5 submitted 2025-08-15 physics.comp-ph

Yue Yang , Xiaodong Yang , Chenhui Lin , Luo Xu , Qi Wang , Shuwei Xu , Wenchuan Wu This is my paper · ORCID
classification physics.comp-ph MSC 65R2045B0580A19
keywords geothermalenergyboreholeheatexchangerg-functiondiscretizationmethodGauss-LegendrequadratureFredholmintegralequationregularizationnumericalsimulation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that the dimensionless thermal response function of a geothermal borehole field, the g-function used to design ground-source heat pumps and seasonal thermal storage, should be computed by solving the spatio-temporal integral equations that enforce a uniform borehole wall temperature rather than by the usual stacked finite line source (SFLS) approximation. The SFLS model is shown to be one special discretization of those equations, in which the heat extraction rate is taken constant within each depth segment. The paper replaces that stepwise approximation with Gauss-Legendre quadrature, so depth integrals become weighted sums at prescribed points and the costly segment-to-segment response integrals become analytic point-to-point factors; it also shows the discretized system is a Fredholm equation of the first kind that can become ill-conditioned enough at high order for the g-function to diverge, and that a small diagonal regularization restores convergence. If the numerical tests are representative, the method computes g-functions with comparable or better accuracy at 20 to 200 times the speed of optimized SFLS, reducing the tested 3×3 bore-field case from about 31 minutes to about 14 seconds.

What carries the argument

The load-bearing object is the Nyström discretization (a quadrature-based way of solving integral equations) applied to the governing equations: integrals over depth are replaced by quadrature sums, so the heat-extraction profile enters only through values at Gauss-Legendre nodes $z^{\rm heat}_{ij}$ with weights $w^{\rm heat}_{ij}$, and the wall-temperature constraint is evaluated at temperature quadrature points $z^{\rm temp}_{mnl}$ with weights $w^{\rm temp}_{mnl}$. The coefficient matrix is built from analytic point-to-point response factors $J(d,t_1,t_2)=\frac{1}{4\pi k d}[\operatorname{erfc}(d/2\sqrt{\alpha t_2})-\operatorname{erfc}(d/2\sqrt{\alpha t_1})]$, which replace the double-integral segment-to-segment factors of SFLS. The second load-bearing piece is the diagnosis that the governing equations form a Fredholm integral equation of the first kind, meaning the unknown heat extraction rate appears only inside the integral and the inverse problem is ill-posed; the paper's remedy is to add $\lambda I$ to the coefficient block (Eq. 41), turning the unstable first-kind system into a stable regularized one.

What would settle it

Take the two-borehole test case at $N_{\rm heat}^{\rm quad}=512$ and compute the regularized g-function for $\lambda$ ranging from $10^{-12}$ to $10^{-4}$; if the value moves by more than the quoted relative error of roughly $10^{-5}$ as $\lambda$ varies, the accuracy claim depends on an unexamined parameter choice. Cross-check the regularized value against an independent, finely meshed finite-element simulation with the same uniform borehole wall temperature boundary condition.

Watch

Extended reading notes

Core claim

The central claim is that the g-function under the uniform borehole wall temperature boundary condition is governed by a closed system of integral equations (Eqs. 11–14) coupling the depth- and time-dependent heat extraction rate $Q_i(z,t_k)$ to the uniform wall temperature $T(t_k)$, with the total extraction rate normalized to 1. In this formulation the SFLS model is exactly the discretization that approximates $Q_i$ by a step function constant on each segment, and its segment-to-segment response factors are double integrals of the kernel $J(d,t_1,t_2)$. The proposed method instead applies a Nyström discretization: the depth integrals of heat extraction are replaced by weighted sums at Gauss-Legendre points, the average wall temperature is evaluated with a second Gauss-Legendre rule over temperature points, and the resulting linear system couples only point-to-point response factors $J(r(\cdot),t_k-t_p,t_k-t_{p-1})$, which are analytic. Because the infinite-dimensional problem is a Fredholm integral equation of the first kind, the discretized matrix is ill-conditioned and its condition number grows rapidly with the number of unknowns, reaching about $10^{18}$ at 1024 points in the tested case; adding $\lambda I$ with $\lambda=10^{-8}$ to the coefficient matrix stabilizes the system. With regularization, the method converges at $N_{\rm heat}^{\rm quad}=512$ to a relative error of $2.8\times 10^{-5}$, about an order of magnitude below the best nonuniform SFLS scheme tested, while taking 20–200 times less computation.

Load-bearing premise

All numerical results rest on a single fixed regularization parameter $\lambda=10^{-8}$ being simultaneously small enough to leave the physical g-function unchanged and large enough to suppress the ill-conditioning, but the paper gives no sensitivity study or error bound connecting the regularized solution to the true solution of the Fredholm equation.

Editorial extensions

If this is right

  • If the claim is right, fine-grained g-functions for large bore fields can be produced in seconds instead of minutes, making simulation-based design optimization and uncertainty quantification practical.
  • The divergence at high discretization order is explained by matrix conditioning rather than by a failure of the physical model, so users of fine SFLS discretizations now know when to mistrust the output and what to fix.
  • Gauss-Legendre point placement removes the need to hand-tune or empirically optimize segment lengths, so the accuracy claims transfer across borehole configurations without re-tuning.
  • At the highest tested resolution the quadrature method's relative error is about $2.8\times 10^{-5}$, roughly ten times smaller than the optimized nonuniform SFLS value of $2.6\times 10^{-4}$, so the accuracy advantage coexists with the speed advantage.
  • The linear system is built by pairing $N_b N_{\rm heat}^{\rm quad}$ heat points with $N_b N_{\rm heat}^{\rm quad} N_{\rm temp}^{\rm quad}$ temperature points, a structure whose cost grows quadratically with the number of unknowns but is still cheaper than SFLS because each factor is analytic.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: a sweep of $\lambda$ from $10^{-12}$ to $10^{-4}$ at $N_{\rm heat}^{\rm quad}=512$ would reveal whether the fixed $\lambda=10^{-8}$ is load-bearing; if the solution is insensitive across that range, the accuracy claim is robust, and if not, the durable contribution is the ill-conditioning diagnosis rather than the specific speed numbers.
  • Editorial inference: the Fredholm-first-kind formulation invites applying standard inverse-problem parameter selection rules, such as the discrepancy principle or L-curve, to choose $\lambda$ per time step and per bore field, which could make the regularization adaptive rather than fixed.
  • Editorial inference: the same Nyström-plus-regularization view may extend to other borehole boundary conditions, inclined boreholes, or boreholes with internal thermal resistance, and whether the first-kind ill-posedness survives those generalizations is a testable question.
  • Editorial inference: if the reported speed and accuracy hold in practice, on-the-fly g-function generation in optimization loops becomes feasible, replacing precomputed libraries as the standard way to size ground-source heat pump fields.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper reformulates the computation of g-functions for vertical geothermal boreholes under the uniform borehole wall temperature (UBWT) boundary condition as a system of spatio-temporal integral equations. It shows that the stacked finite line source (SFLS) model is a stepwise-constant discretization of these integral equations, and proposes an alternative Nyström-type discretization using Gauss-Legendre quadrature in space and point-to-point analytical response factors. The authors further identify the governing equations as Fredholm integral equations of the first kind, explain the observed divergence at high discretization orders as ill-conditioning, and add a diagonal regularization term lambda*I to stabilize the linear systems. Numerical experiments on 2, 2x2, and 3x3 borehole fields are used to claim that the proposed method achieves comparable or better accuracy than optimized SFLS discretizations while being 20-200 times faster.

Significance. If the accuracy and efficiency claims are substantiated, the paper makes a useful contribution: it provides a principled quadrature-based discretization that avoids empirical segment selection, replaces segment-to-segment integrals with cheaper point-to-point response factors, and gives a plausible explanation for the previously observed divergence of g-functions at high discretization orders. The derivation from the Green's function to the discretized linear systems is coherent, and the complexity reduction is plausible from the stated O(N_b^2 N_heat_quad^2 N_temp_quad) structure. However, the central accuracy claims currently rest on a self-referential benchmark, and the regularization parameter is fixed without a sensitivity study. These issues must be addressed before the reported error reductions can be considered established.

major comments (3)
  1. [Section 3.2, Eq. (42), Tables 10 and 12-14] The 'exact' g-function used as the reference in the error metric is the proposed method's own regularized solution at Nheat_quad=1024 and Ntemp_quad=50. This reference is not independent of the method being evaluated, and it inherits any bias introduced by the fixed regularization parameter lambda=1e-8. The convergence tables do not resolve the ambiguity: at t=20a and N=1024, the proposed method gives g=7.2889 (Table 6) while the optimal nonuniform SFLS method gives g=7.2904 (Table 9), a relative difference of about 2e-4, which is the same order as several of the errors reported in Table 12. Without comparison to an external reference (for example, a high-resolution finite-element or super-position-borehole-model solution, or an independently converged unregularized solution), the claim of 'nearly 10 times improvement in accuracy' is not established.
  2. [Section 2.4, Eqs. (37) and (41)] The regularization parameter lambda is fixed at 1e-8 for every numerical result, with no sensitivity study and no error bound relating the solution of the regularized system to the solution of the original Fredholm first-kind equation. Since Eq. (37) modifies the physical boundary condition to T_reg = T - lambda*Q, and since the condition number at N=1024 remains as large as 1.48e10 (Table 6), one cannot exclude a regularization bias at the 1e-4 level. The paper should report g-function values for a range of lambda (for example, 1e-12 to 1e-4) at each discretization order and demonstrate that the chosen lambda leaves the results unchanged within the claimed tolerance.
  3. [Section 3.2, Eq. (42)] The scalar error defined in Eq. (42) measures the relative difference between sums of g-function values over all ten time points, rather than a pointwise or norm-based error. Errors with different signs at different time points can cancel in this aggregate, so Tables 10 and 12-14 may understate local inaccuracies. Reporting per-time-point relative errors or an L2/max norm over the time points would make the accuracy comparison more robust and would let the reader verify the claimed accuracy improvement at the scale of individual response times.
minor comments (6)
  1. [Eq. (5)] The first equality in Eq. (5) states that J(d,0,t2) equals the integral from t1 to t2, but the definition and the erfc expression correspond to an integral from 0 to t2; the integration limits should be corrected.
  2. [Eqs. (32)-(33)] The segment bounds use z_{(i-1)j}, which is not a well-defined quantity for the i-th borehole's quadrature points; this appears to be a typo for z_{i(j-1)} or a similar index, and should be clarified so the segmentation scheme is unambiguous.
  3. [Abstract and Section 3.2] The abstract and the introduction claim speedups of 20-200 times, but the reported numerical results show 20-160 times in Table 12 and 20-150 times in Table 14; the magnitudes should be reconciled or the claim should be stated as 'up to 20-160 times' based on the actual results.
  4. [Throughout] There are several typos, including 'thay' in Section 3.1.1, 'seperate' in Section 2.3, 'segementation' in Section 3, and 'discretimzation' in Section 3.2, which should be corrected in a revised version.
  5. [Table 7] The column header in Table 7 reads 'Nheat_quad' for the SFLS model, but the row variable is the number of segments Ns; using a consistent header would avoid confusion.
  6. [Reproducibility] The paper does not include a data or code availability statement; since the timing comparisons depend on implementation details, a link to the implementation or a detailed description of the numerical libraries and hardware would strengthen reproducibility.

Circularity Check

1 steps flagged · score 4.0 of 10

Accuracy benchmark is anchored to the proposed method's own regularized high-order solution, so the claimed accuracy advantage over SFLS is partially self-referential; the core derivation itself is independent.

  1. self definitional [Section 3.2, Eq. (42) and the paragraph defining the error metric]
    "The accuracy is evaluated by the relative error between the approximate value and the assumed exact value of g-function when Nheat_quad = 1024 and Ntemp_quad = 50. The regularization parameter λ is set to 1×10−8. Error = (sum gapprox - sum gexact) / (sum gexact)"

    The reference 'exact' g-function is produced by the proposed method itself, at Nheat_quad=1024 and Ntemp_quad=50, using the same fixed regularization parameter λ=1e-8 that is used in all compared runs. Thus Eq. (42) measures the distance of each method to the proposed method's own regularized solution, not to an independent physical UBWT g-function. Any systematic bias introduced by the Gauss-Legendre discretization or by the λf term in Eq. (37) is shared by both the reference and the proposed candidates, so that bias cancels from the reported relative errors. The claimed 'nearly 10 times improvement in accuracy' over SFLS is therefore a statement of self-consistency with the chosen regularized reference, not an independently established accuracy comparison.

full rationale

The paper's central methodological contribution, reformulating g-function calculation as spatio-temporal integral equations and discretizing them with Gauss-Legendre quadrature to obtain analytical point-to-point response factors, is derived from the governing equations (Eqs. 11-14) and is not itself circular. The speed comparison is also independent: the measured wall-clock times and the complexity argument comparing segment-to-segment versus point-to-point factors do not rely on a fitted or self-defined target. No load-bearing self-citation chain was found; the paper cites prior work for standard results such as segment-to-segment response factors and regularization, but these are external. The identified circularity is confined to the validation benchmark: Eq. (42) defines g_exact as the proposed method's own high-order regularized solution, so the accuracy tables measure agreement with that self-chosen reference rather than with an external solution. Because the efficiency claim is unaffected and the mathematical derivation is self-contained, a moderate score of 4 is appropriate rather than a higher score.

Assumptions & free parameters 1 free parameters · 7 assumptions · 0 invented entities

The method introduces no new physical entities. Its dependence on prior literature is the line-source heat conduction model, which is standard in the field; the only hand-tuned numerical parameter is the regularization coefficient lambda.

free parameters (1)
  • regularization parameter lambda = 1e-8
    Chosen by hand to stabilize the ill-conditioned linear system; no sensitivity analysis or optimality criterion is provided. All accuracy results depend on this single value.
assumptions (7)
  • domain assumption Heat conduction in the ground is governed by the linear diffusion equation in an infinite isotropic medium with constant thermal properties, so point-source superposition (Green's function, Eq. 3) applies.
    Section 2.1, Eqs. (2)-(3). This is standard for FLS models but ignores groundwater, heterogeneity, and temperature-dependent properties.
  • domain assumption The mirror-borehole construction (opposite heat extraction above ground) exactly enforces a constant temperature ground surface.
    Introduced in Section 2.1 and used in Eq. (7). Standard method, but assumes semi-infinite ground and no surface temperature variation.
  • domain assumption The borehole can be represented as a line source of zero radius, with the borehole radius r_b used only to avoid the singularity in Eq. (8).
    Section 2.1, Eq. (8). This neglects borehole thermal resistance and internal pipe geometry.
  • domain assumption Under the uniform borehole wall temperature boundary condition, the temperature is uniform along each borehole wall (Eq. 14), and the average heat extraction rate per length is normalized to unity (Eq. 13).
    Section 2.1, Eqs. (13)-(14). This is the definition of the g-function problem being solved.
  • ad hoc to paper The g-function can be computed recursively in time by treating the heat extraction rate as piecewise constant on each time interval [t_{p-1}, t_p] (Eq. 9).
    Section 2.1, Eq. (9). Standard temporal discretization with stepwise Q(t); error from this approximation is not quantified.
  • standard math Gauss-Legendre quadrature converges for the spatial integrals and the temperature-average integrals, so replacing integrals by weighted point sums (Eqs. 29-31, 34) yields the true solution as quadrature order increases.
    Section 2.3. Valid for smooth integrands; the heat extraction rate near borehole ends may not be smooth, and convergence is demonstrated numerically only for one case.
  • ad hoc to paper The discretized Fredholm first-kind system can be regularized by adding lambda*I (Eq. 41), and the solution of the regularized system approaches the true solution as lambda goes to zero.
    Section 2.4, Eqs. (37)-(41). Standard theory, but the choice of lambda and the approximation error are not analyzed.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Novel discretization method to calculate g-functions of vertical geothermal boreholes with improved accuracy and efficiency." pith.science (2026). https://pith.science/paper/KDA33RM5

@misc{pith2026250811154,
  author       = {Pith},
  title        = {Pith review of: Novel discretization method to calculate g-functions of vertical geothermal boreholes with improved accuracy and efficiency},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KDA33RM5}},
  note         = {Machine review of arXiv:2508.11154}
}
read the original abstract

The calculation of g-functions is essential for the design and simulation of geothermal boreholes. However, existing methods, such as the stacked finite line source (SFLS) model, face challenges regarding computational efficiency and accuracy, particularly with fine-grained discretization. This paper introduces a novel discretization method to address these limitations. We reformulate the g-function calculation under the uniform borehole wall temperature boundary condition as the solution to spatio-temporal integral equations. The SFLS model is identified as a special case using stepwise approximation of the heat extraction rate. Our proposed method employs the Gauss-Legendre quadrature to approximate the spatial integrals with a weighted sum of function values at strategically chosen points. This transforms the time-consuming segment-to-segment integral calculations in SFLS model into simpler and analytical point-to-point response factors. Furthermore, we identify that the governing integral equations are of the Fredholm first kind, leading to ill-conditioned linear systems that can cause g-function to diverge at high discretization orders. To address this, a regularization technique is implemented to ensure stable and convergent solutions. Numerical tests demonstrate that the proposed method is significantly more efficient, achieving comparable or improved accuracy at speeds 20 to 200 times faster than the SFLS model with optimized nonuniform discretization schemes.

Figures

Figures reproduced from arXiv: 2508.11154 by the authors.

Figure 1
Figure 1. Spatial superposition of boreholes and their mirrors [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Stacked finite line source model Stacked finite line source (SFLS) model divides each borehole into segments and approximates Qi(z, tk) as a stepwise function that remains constant in each segment. If each borehole is divided into Ns segments, the integral of heat extraction can be decomposed into summation of Ns discrete terms: T(X, Y, Z, tk) = T0(X, Y, Z, tk) − X Nb i=1 X Ns j=1 Qi,j (tk) Z UBi,j LBi,j J(r(X − xi … view at source ↗
Figure 3
Figure 3. Quadrature points of Gauss-Legendre rule with different order in an interval [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4 [PITH_FULL_IMAGE:figures/full_fig_p017_4.png]

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    Sarbu, C

    I. Sarbu, C. Sebarchievici, General review of ground-source heat pump systems for heating and cooling of buildings, Energy and buildings 70 (2014) 441–454

    work page 2014

  2. [2]

    J. Luo, Q. Zhang, C. Liang, H. Wang, X. Ma, An overview of the recent development of the ground source heat pump (gshp) system in china, Renewable Energy 210 (2023) 269–279

    work page 2023

  3. [3]

    X. Wang, H. Zhang, L. Cui, J. Wang, C. Lee, X. Zhu, Y. Dong, Borehole thermal energy storage for building heating application: A review, Renewable and Sustainable Energy Reviews 203 (2024) 114772

    work page 2024

  4. [4]

    Eskilson, Thermal analysis of heat extraction boreholes (1987)

    P. Eskilson, Thermal analysis of heat extraction boreholes (1987). 20

    work page 1987

  5. [5]

    Rees, Advances in ground-source heat pump systems, Woodhead Publishing, 2016

    S. Rees, Advances in ground-source heat pump systems, Woodhead Publishing, 2016

    work page 2016

  6. [6]

    T. N. West, J. D. Spitler, Ground heat exchanger design tool with rowwise placement of boreholes, Science and Technology for the Built Environment 30 (9) (2024) 1148–1166

    work page 2024

  7. [7]

    X. Liu, J. D. Spitler, M. Qu, L. Shi, Recent developments in the design of vertical borehole ground heat exchangers for cost reduction and thermal energy storage, Journal of Energy Resources Technology 143 (10) (2021) 100803

    work page 2021

  8. [8]

    Marcotte, P

    D. Marcotte, P. Pasquier, Fast fluid and ground temperature computation for geothermal ground-loop heat exchanger systems, Geothermics 37 (6) (2008) 651–665

    work page 2008

Show all 37 references
  1. [9]

    Y. Chen, B. Pan, X. Zhang, C. Du, Thermal response factors for fast parameterized design and long-term performance simulation of vertical gchp systems, Renewable energy 136 (2019) 793–804

    2019

  2. [10]

    M. A. Bernier, P. Pinel, R. Labib, R. Paillot, A multiple load aggregation algorithm for annual hourly simulations of gchp systems, Hvac&R Research 10 (4) (2004) 471–487

    2004

  3. [11]

    Claesson, S

    J. Claesson, S. Javed, A load-aggregation method to calculate extraction temperatures of borehole heat exchangers, ASHRAE transactions 118 (1) (2012) 530–540

    2012

  4. [12]

    Pahud, G

    D. Pahud, G. Hellstr¨ om, L. Mazzarella, Duct ground heat storage model for trnsys (trnvdst), Laboratory of Energy Systems, Lausanne (1997)

    1997

  5. [13]

    Y. L. E. Law, S. B. Dworkin, Characterization of the effects of borehole configuration and interference with long term ground temperature modelling of ground source heat pumps, Applied Energy 179 (2016) 1032–1047

    2016

  6. [14]

    X. Yang, W. Cai, Y. Li, M. Wang, Y. Kong, F. Wang, C. Chen, Numerical investigation on the influence of groundwater flow on long-term heat extraction performance of deep borehole heat exchanger array, Geothermal Energy 12 (1) (2024) 45

    2024

  7. [15]

    H. Y. Zeng, N. R. Diao, Z. H. Fang, A finite line-source model for boreholes in geothermal heat exchangers, Heat Transfer Asian Research 31 (7) (2002) 558–567

    2002

  8. [16]

    Lamarche, B

    L. Lamarche, B. Beauchamp, A new contribution to the finite line-source model for geothermal boreholes, Energy and buildings 39 (2) (2007) 188–198

    2007

  9. [17]

    P. Cui, H. Yang, Z. Fang, Heat transfer analysis of ground heat exchangers with inclined boreholes, Applied thermal engineering 26 (11-12) (2006) 1169–1175

    2006

  10. [18]

    Marcotte, P

    D. Marcotte, P. Pasquier, The effect of borehole inclination on fluid and ground temperature for glhe systems, Geothermics 38 (4) (2009) 392–398

    2009

  11. [19]

    Lamarche, A fast algorithm for the hourly simulations of ground-source heat pumps using arbitrary response factors, Renewable Energy 34 (10) (2009) 2252–2258

    L. Lamarche, A fast algorithm for the hourly simulations of ground-source heat pumps using arbitrary response factors, Renewable Energy 34 (10) (2009) 2252–2258

    2009

  12. [20]

    Fossa, The temperature penalty approach to the design of borehole heat exchangers for heat pump applications, Energy and Buildings 43 (6) (2011) 1473–1479

    M. Fossa, The temperature penalty approach to the design of borehole heat exchangers for heat pump applications, Energy and Buildings 43 (6) (2011) 1473–1479

    2011

  13. [21]

    Cimmino, M

    M. Cimmino, M. Bernier, A semi-analytical method to generate g-functions for geothermal bore fields, International Journal of Heat and Mass Transfer 70 (2014) 641–650

    2014

  14. [22]

    Monzo, P

    P. Monzo, P. Mogensen, J. Acu˜ na, F. Ruiz-Calvo, C. Montagud, A novel numerical approach for imposing a temperature boundary condition at the borehole wall in borehole fields, Geothermics 56 (2015) 35–44

    2015

  15. [23]

    Lazzarotto, A methodology for the calculation of response functions for geothermal fields with arbitrarily oriented boreholes–part 1, Renewable Energy 86 (2016) 1380–1393

    A. Lazzarotto, A methodology for the calculation of response functions for geothermal fields with arbitrarily oriented boreholes–part 1, Renewable Energy 86 (2016) 1380–1393

    2016

  16. [24]

    Lazzarotto, F

    A. Lazzarotto, F. Bj¨ ork, A methodology for the calculation of response functions for geothermal fields with arbitrarily oriented boreholes–part 2, Renewable energy 86 (2016) 1353–1361

    2016

  17. [25]

    M. Cimmino, The effects of borehole thermal resistances and fluid flow rate on the g-functions of geothermal bore fields, International Journal of Heat and Mass Transfer 91 (2015) 1119–1127

    2015

  18. [26]

    M. Cimmino, Semi-analytical method for g-function calculation of bore fields with series-and parallel-connected boreholes, Science and Technology for the Built Environment 25 (8) (2019) 1007–1022

    2019

  19. [27]

    Cimmino, g-functions for fields of series-and parallel-connected boreholes with variable fluid mass flow rate and re- versible flow direction, Renewable Energy 228 (2024) 120661

    M. Cimmino, g-functions for fields of series-and parallel-connected boreholes with variable fluid mass flow rate and re- versible flow direction, Renewable Energy 228 (2024) 120661

    2024

  20. [28]

    Cimmino, D

    M. Cimmino, D. Thorne, C. Langevin, M. Sukop, Pygfunction: an open-source toolbox for the evaluation of thermal response factors for geothermal borehole fields, Proceedings of eSim 2018 (2018)

    2018

  21. [29]

    Cimmino, J

    M. Cimmino, J. C. Cook, pygfunction 2.2: New features and improvements in accuracy and computational efficiency (2022)

    2022

  22. [30]

    J. C. Cook, Development of computer programs for fast computation of g-functions and automated ground heat exchanger design, Master’s thesis, Oklahoma State University (2021)

    2021

  23. [31]

    Lamarche, G-function generation using a piecewise-linear profile applied to ground heat exchangers, International Journal of Heat and Mass Transfer 115 (2017) 354–360

    L. Lamarche, G-function generation using a piecewise-linear profile applied to ground heat exchangers, International Journal of Heat and Mass Transfer 115 (2017) 354–360

    2017

  24. [32]

    Cimmino, J

    M. Cimmino, J. C. Cook, J. A. Isiordia Farrera, Optimal discretization of geothermal boreholes for the calculation of g-functions, Science and Technology for the Built Environment 30 (3) (2024) 234–249

    2024

  25. [33]

    H. S. Carslaw, J. C. Jaeger, Conduction of Heat in Solids, 2nd Edition, Oxford University Press, Oxford, 1986

    1986

  26. [34]

    Claesson, S

    J. Claesson, S. Javed, An analytical method to calculate borehole fluid temperatures for time-scales from minutes to decades., ASHRAE Transactions 117 (2) (2011)

    2011

  27. [35]

    Hackbusch, Integral equations: theory and numerical treatment, Vol

    W. Hackbusch, Integral equations: theory and numerical treatment, Vol. 120, Birkh¨ auser, 2012

    2012

  28. [36]

    Wazwaz, Linear and nonlinear integral equations, Vol

    A.-M. Wazwaz, Linear and nonlinear integral equations, Vol. 639, Springer, 2011

    2011

  29. [37]

    Wazwaz, The regularization method for fredholm integral equations of the first kind, Computers & Mathematics with Applications 61 (10) (2011) 2981–2986

    A.-M. Wazwaz, The regularization method for fredholm integral equations of the first kind, Computers & Mathematics with Applications 61 (10) (2011) 2981–2986. 21

    2011

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.