- Add test in optimisation methods to use fixed Rb when constant fluid and flow rate data

Author: wouterpeere

GHEtool/Methods/optimise_load_profile.py

@@ -2,7 +2,8 @@
 import numpy as np
 
 from typing import Union
-from GHEtool.VariableClasses import HourlyBuildingLoad, MonthlyBuildingLoadMultiYear, HourlyBuildingLoadMultiYear
+from GHEtool.VariableClasses import HourlyBuildingLoad, MonthlyBuildingLoadMultiYear, HourlyBuildingLoadMultiYear, \
+    ConstantFluidData, ConstantFlowRate
 
 
 def optimise_load_profile_power(
@@ -59,6 +60,12 @@ def optimise_load_profile_power(
     if temperature_threshold < 0:
         raise ValueError(f"The temperature threshold is {temperature_threshold}, but it cannot be below 0!")
 
+    # since the depth does not change, the Rb* value is constant, if there is no temperature dependent fluid data
+    if isinstance(borefield.borehole.fluid_data, ConstantFluidData) \
+            and isinstance(borefield.borehole.flow_data, ConstantFlowRate):
+        borefield.Rb = borefield.borehole.get_Rb(borefield.H, borefield.D, borefield.r_b,
+                                                 borefield.ground_data.k_s(borefield.depth, borefield.D))
+
     # set load
     borefield.load = copy.deepcopy(building_load)
 
@@ -194,6 +201,12 @@ def optimise_load_profile_energy(
     if temperature_threshold < 0:
         raise ValueError(f"The temperature threshold is {temperature_threshold}, but it cannot be below 0!")
 
+    # since the depth does not change, the Rb* value is constant, if there is no temperature dependent fluid data
+    if isinstance(borefield.borehole.fluid_data, ConstantFluidData) \
+            and isinstance(borefield.borehole.flow_data, ConstantFlowRate):
+        borefield.Rb = borefield.borehole.get_Rb(borefield.H, borefield.D, borefield.r_b,
+                                                 borefield.ground_data.k_s(borefield.depth, borefield.D))
+
     building_load_copy = copy.deepcopy(building_load)
 
     # if building load is not a multi-year load, convert to multiyear
@@ -398,6 +411,12 @@ def optimise_load_profile_balance(
     if imbalance_factor > 1 or imbalance_factor < 0:
         raise ValueError(f"The imbalance factor is {imbalance_factor}, but it should be between 0-1!")
 
+    # since the depth does not change, the Rb* value is constant, if there is no temperature dependent fluid data
+    if isinstance(borefield.borehole.fluid_data, ConstantFluidData) \
+            and isinstance(borefield.borehole.flow_data, ConstantFlowRate):
+        borefield.Rb = borefield.borehole.get_Rb(borefield.H, borefield.D, borefield.r_b,
+                                                 borefield.ground_data.k_s(borefield.depth, borefield.D))
+
     # set load
     borefield.load = copy.deepcopy(building_load)