Merge pull request #2042 from ibpsa/issue2041_borefieldSteadyState

mwetter · web-flow · commit 18747cd6e584 · 2025-09-09T06:16:41.000-07:00
Guarded against division by zero in borefield data
diff --git a/IBPSA/Fluid/Geothermal/Borefields/Data/Filling/Template.mo b/IBPSA/Fluid/Geothermal/Borefields/Data/Filling/Template.mo
@@ -8,10 +8,11 @@ record Template
     "Specific heat capacity of the borehole filling material";
   parameter Modelica.Units.SI.Density dFil(displayUnit="kg/m3")
     "Density of the borehole filling material";
-  parameter Boolean steadyState = (cFil < Modelica.Constants.eps or dFil < Modelica.Constants.eps)
+  parameter Boolean steadyState = dFil*cFil < Modelica.Constants.eps
     "Flag, if true, then material is computed using steady-state heat conduction"
     annotation(Evaluate=true);
-  final parameter Modelica.Units.SI.ThermalDiffusivity aFil=kFil/(dFil*cFil)
+  final parameter Modelica.Units.SI.ThermalDiffusivity aFil=
+    if (dFil*cFil < Modelica.Constants.eps) then 1E20 else kFil/(dFil*cFil)
     "Heat diffusion coefficient of the borehole filling material";
   annotation (
   defaultComponentPrefixes="parameter",
@@ -25,6 +26,12 @@ IBPSA.Fluid.Geothermal.Borefields.Data.Filling</a>.</p>
 revisions="<html>
 <ul>
 <li>
+September 3, 2025, by Michael Wetter:<br/>
+Guarded against division by zero.<br/>
+This is for
+<a href=\"https://github.com/ibpsa/modelica-ibpsa/issues/2041\">IBPSA, issue 2041</a>.
+</li>
+<li>
 July 15, 2018, by Michael Wetter:<br/>
 Revised implementation, added <code>defaultComponentPrefixes</code> and
 <code>defaultComponentName</code>.
