Introduced distance between wells as parameter

ibpsa · Oct 5, 2023 · 1da8fe7 · 1da8fe7
1 parent 85ff288
commit 1da8fe7
Show file tree

Hide file tree

Showing 4 changed files with 12 additions and 9 deletions.
diff --git a/IBPSA/Fluid/Geothermal/Aquifer/Examples/CoolingOffice.mo b/IBPSA/Fluid/Geothermal/Aquifer/Examples/CoolingOffice.mo
@@ -23,6 +23,7 @@ model CoolingOffice
     redeclare package Medium = IBPSA.Media.Water,
     nVol=80,
     h=20,
+    d=200,
     length=40,
     rMax=500,
     TCol_start=285.15,

diff --git a/IBPSA/Fluid/Geothermal/Aquifer/MultiWell.mo b/IBPSA/Fluid/Geothermal/Aquifer/MultiWell.mo
@@ -5,12 +5,13 @@ model MultiWell "Model of a single well for aquifer thermal energy storage"
   parameter Integer nVol(min=1)=10 "Number of control volumes used in discretization" annotation (
       Dialog(group="Subsurface"));
   parameter Modelica.Units.SI.Height h "Aquifer thickness";
+  parameter Modelica.Units.SI.Height d "Distance between wells";
   parameter Modelica.Units.SI.Height length
     "Length of one well, used to compute pressure drop";
   parameter Real nPai=1 "Number of paired wells";
   parameter Modelica.Units.SI.Radius rWB=0.1 "Wellbore radius" annotation (
       Dialog(group="Subsurface"));
-  parameter Modelica.Units.SI.Radius rMax "Domain radius" annotation (
+  parameter Modelica.Units.SI.Radius rMax=d/2 "Domain radius" annotation (
       Dialog(group="Subsurface"));
   parameter Real griFac(min=1) = 1.15 "Grid factor for spacing" annotation (
       Dialog(group="Subsurface"));
@@ -31,7 +32,7 @@ model MultiWell "Model of a single well for aquifer thermal energy storage"
   parameter Modelica.Units.SI.MassFlowRate m_flow_nominal "Nominal mass flow rate" annotation (
       Dialog(group="Hydraulic circuit"));
   parameter Modelica.Units.SI.PressureDifference dpAquifer_nominal(displayUnit= "Pa")=
-    m_flow_nominal/nPai*Modelica.Constants.g_n/2/Modelica.Constants.pi/h/aquDat.K*log(rMax/rWB)
+    m_flow_nominal/nPai*Modelica.Constants.g_n/2/Modelica.Constants.pi/h/aquDat.K*log(d/rWB)
      "Pressure drop at nominal mass flow rate in the aquifer"  annotation (
       Dialog(group="Hydraulic circuit"));
   final parameter Modelica.Units.SI.PressureDifference dpWell_nominal(

diff --git a/IBPSA/Fluid/Geothermal/Aquifer/Validation/NumberWells.mo b/IBPSA/Fluid/Geothermal/Aquifer/Validation/NumberWells.mo
@@ -16,10 +16,11 @@ model NumberWells
     m_flow_nominal=0.1,
     dpExt_nominal=0) "Well model";
 
-  MyWell aquWel1
+  MyWell aquWel1(d=4800)
   "ATES with one pair of wells"
     annotation (Placement(transformation(extent={{-20,20},{0,40}})));
-  MyWell aquWel2(nPai=2, m_flow_nominal=0.2)
+  MyWell aquWel2(
+    d=4800,      nPai=2, m_flow_nominal=0.2)
             "ATES with two pairs of wells"
     annotation (Placement(transformation(extent={{-20,-60},{0,-40}})));
   Modelica.Blocks.Sources.Constant uPum(k=1) "Pump control signal"
@@ -36,10 +37,10 @@ model NumberWells
     "Assertion that checks for equality of results"
     annotation (Placement(transformation(extent={{60,60},{80,80}})));
 equation
-  connect(uPum.y, aquWel1.u) annotation (Line(points={{-59,0},{-40,0},{-40,35},{-22,35}},
-                     color={0,0,127}));
-  connect(aquWel2.u,uPum. y) annotation (Line(points={{-22,-45},{-40,-45},{-40,0},{-59,0}},
-                       color={0,0,127}));
+  connect(uPum.y, aquWel1.u) annotation (Line(points={{-59,0},{-40,0},{-40,30},
+          {-22,30}}, color={0,0,127}));
+  connect(aquWel2.u,uPum. y) annotation (Line(points={{-22,-50},{-40,-50},{-40,
+          0},{-59,0}}, color={0,0,127}));
   connect(aquWel2.port_Col, aquWel2.port_Hot) annotation (Line(points={{-16,-40},
           {-16,-20},{-4,-20},{-4,-40}}, color={0,127,255}));
   connect(aquWel1.port_Col, aquWel1.port_Hot) annotation (Line(points={{-16,40},

diff --git a/IBPSA/Fluid/Geothermal/Aquifer/Validation/SimulationTest.mo b/IBPSA/Fluid/Geothermal/Aquifer/Validation/SimulationTest.mo
@@ -6,8 +6,8 @@ model SimulationTest
     redeclare package Medium = IBPSA.Media.Water,
     nVol=232,
     h=200,
+    d=4800,
     length=40,
-    rMax=2400,
     griFac=1.1,
     TCol_start=307.15,
     THot_start=393.15,