update the benchmark input and output files

Tests/Data/Parabolic/T/3D_3BHEs_array/3bhes_1U.prj

@@ -145,14 +145,14 @@
                     <type>AqueousLiquid</type>
                     <properties>
                         <property>
-                            <name>specific_heat_capacity</name>
+                            <name>phase_velocity</name>
                             <type>Constant</type>
-                            <value>4068</value>
+                            <value>0 0 0</value>
                         </property>
                         <property>
-                            <name>thermal_conductivity</name>
+                            <name>specific_heat_capacity</name>
                             <type>Constant</type>
-                            <value>0.1</value>
+                            <value>4068</value>
                         </property>
                         <property>
                             <name>density</name>
@@ -169,11 +169,6 @@
                             <type>Constant</type>
                             <value>1780</value>
                         </property>
-                        <property>
-                            <name>thermal_conductivity</name>
-                            <type>Constant</type>
-                            <value>2.4</value>
-                        </property>
                         <property>
                             <name>density</name>
                             <type>Constant</type>
@@ -189,11 +184,6 @@
                             <type>Constant</type>
                             <value>1000</value>
                         </property>
-                        <property>
-                            <name>thermal_conductivity</name>
-                            <type>Constant</type>
-                            <value>3.2</value>
-                        </property>
                         <property>
                             <name>density</name>
                             <type>Constant</type>
@@ -203,6 +193,16 @@
                 </phase>
             </phases>
             <properties>
+                <property>
+                    <name>porosity</name>
+                    <type>Constant</type>
+                    <value>0</value>
+                </property>
+                <property>
+                    <name>thermal_conductivity</name>
+                    <type>Constant</type>
+                    <value>2.4</value>
+                </property>
                 <property>
                     <name>thermal_longitudinal_dispersivity</name>
                     <type>Constant</type>
@@ -268,7 +268,6 @@
                 </variables>
         </output>
     </time_loop>
-
     <parameters>
         <parameter>
             <name>T0</name>
@@ -281,7 +280,6 @@
             <values>305.36 305.13 305.23 305.115</values>
         </parameter>
     </parameters>
-
     <process_variables>
         <process_variable>
             <name>temperature_soil</name>

Tests/Data/Parabolic/T/3D_3BHEs_array/bcs_tespy.py

@@ -39,7 +39,7 @@ def timerange(t):
 #month demand
 def consumer_demand(t):#dynamic thermal demand from consumer
     #thermal demand in each month, assumed specific heat extraction rate * lenth of BHE * number of BHE
-    month_demand = [-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3,-10*50*3]
+    month_demand = [-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3,-25*50*3]
     return month_demand[t -1]
 
 # %% dynamic hydraulic flow rate