Merge pull request #2221 from HBShaoUFZ/bhe_impl

Implementation of HeatTransportBHE process

Applications/ApplicationsLib/ProjectData.cpp

@@ -53,6 +53,9 @@
 #ifdef OGS_BUILD_PROCESS_HEATCONDUCTION
 #include "ProcessLib/HeatConduction/CreateHeatConductionProcess.h"
 #endif
+#ifdef OGS_BUILD_PROCESS_HEATTRANSPORTBHE
+#include "ProcessLib/HeatTransportBHE/CreateHeatTransportBHEProcess.h"
+#endif
 #ifdef OGS_BUILD_PROCESS_HYDROMECHANICS
 #include "ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.h"
 #endif
@@ -371,6 +374,24 @@ void ProjectData::parseProcesses(BaseLib::ConfigTree const& processes_config,
         }
         else
 #endif
+#ifdef OGS_BUILD_PROCESS_HEATTRANSPORTBHE
+            if (type == "HEAT_TRANSPORT_BHE")
+        {
+            if (_mesh_vec[0]->getDimension() != 3)
+            {
+                OGS_FATAL(
+                    "HEAT_TRANSPORT_BHE can only work with a 3-dimentional "
+                    "mesh! ");
+            }
+
+            process =
+                ProcessLib::HeatTransportBHE::createHeatTransportBHEProcess(
+                    *_mesh_vec[0], std::move(jacobian_assembler),
+                    _process_variables, _parameters, integration_order,
+                    process_config, _curves);
+        }
+        else
+#endif
 #ifdef OGS_BUILD_PROCESS_HYDROMECHANICS
             if (type == "HYDRO_MECHANICS")
         {
@@ -455,17 +476,17 @@ void ProjectData::parseProcesses(BaseLib::ConfigTree const& processes_config,
             {
                 case 2:
                     process =
-                        ProcessLib::PhaseField::createPhaseFieldProcess<
-                            2>(*_mesh_vec[0], std::move(jacobian_assembler),
-                               _process_variables, _parameters,
-                               integration_order, process_config);
+                        ProcessLib::PhaseField::createPhaseFieldProcess<2>(
+                            *_mesh_vec[0], std::move(jacobian_assembler),
+                            _process_variables, _parameters, integration_order,
+                            process_config);
                     break;
                 case 3:
                     process =
-                        ProcessLib::PhaseField::createPhaseFieldProcess<
-                            3>(*_mesh_vec[0], std::move(jacobian_assembler),
-                               _process_variables, _parameters,
-                               integration_order, process_config);
+                        ProcessLib::PhaseField::createPhaseFieldProcess<3>(
+                            *_mesh_vec[0], std::move(jacobian_assembler),
+                            _process_variables, _parameters, integration_order,
+                            process_config);
                     break;
             }
         }

CMakeLists.txt

@@ -90,6 +90,7 @@ set(ProcessesList
     GroundwaterFlow
     HT
     HeatConduction
+    HeatTransportBHE
     HydroMechanics
     LiquidFlow
     LIE

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/borehole/i_borehole.md

+It represents the well used to install pipes for heat transport.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/borehole/t_diameter.md

+It is the borehole diameter. One of the most important geometric parameters of borehole and needs to be specified here.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/borehole/t_length.md

+It is the borehole length. It is another one of the most important geometric parameters of borehole and should be provided.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/flow_and_temperature_control/FixedPowerConstantFlow/c_FixedPowerConstantFlow.md

+It means the borehole heat exchanger has a fixed heat load and the flow rate in the pipes is constant.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/flow_and_temperature_control/FixedPowerConstantFlow/t_flow_rate.md

+It is the flow rate for the FixedPowerConstantFlow condition. Here it needs to be specified while using FixedPowerConstantFlow as BHE boundary condition with the unit of m3/s.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/flow_and_temperature_control/FixedPowerConstantFlow/t_power.md

+It is the fixed power for the FixedPowerConstantFlow condition. It should be provided and its unit is W. Positve value means extracting energy from the soil.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/flow_and_temperature_control/FixedPowerFlowCurve/c_FixedPowerFlowCurve.md

+It represents that the BHE has the fixed power and the flow rate in pipes is various with time given as a curve. It is kind of another BHE boundary condition to calculate the temperature difference between inlet and outlet fluid. With this BHE boundary condition, the flow rate of circulating fluid is various and inflow temperature is calculated by the heat load, flow rate at particular time and the outflow temperature got from the previous time step.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/flow_and_temperature_control/FixedPowerFlowCurve/t_flow_rate_curve.md

+It is the flow rate curve. Here the flow rate in the pipes is not a constant, the flow rate curve should be presented and set as input data.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/flow_and_temperature_control/FixedPowerFlowCurve/t_power.md

+It is the fixed power on the BHE in the FixedPowerFlowCurve condition. It needs to be provided and the unit is W.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/flow_and_temperature_control/TemperatureCurveConstantFlow/c_TemperatureCurveConstantFlow.md

+This is another type of controlling temperature difference with specific inflow temperature curve and constant flow rate.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/flow_and_temperature_control/TemperatureCurveConstantFlow/t_flow_rate.md

+It is the flow rate in the BHE pipes. Flow rate needs to be offered by using m3/s as its unit.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/flow_and_temperature_control/TemperatureCurveConstantFlow/t_temperature_curve.md

+It means the inflow temperature curve varying with time.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/flow_and_temperature_control/i_flow_and_temperature_control.md

+It means the way of controlling flow and temperature in the BHE. In the BHE, the inflow temperature at this time step will be calculated by the heat load on the BHE, flow rate of circulating fluid and previous outflow temperature. Thus, There is a need to specify the way of getting the heat load and flow rate values.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/flow_and_temperature_control/t_type.md

+The type of controlling flow rate and heat load on the BHE.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/grout/i_grout.md

+It is the grout in the borehole. Grout is usually used to keep the pipes in the place. It is like a bridge of exchanging heat energy between soil and pipe.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/grout/t_density.md

+It is the density of the grout.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/grout/t_heat_capacity.md

+Here means specific heat capacity of the grout.

Documentation/ProjectFile/prj/processes/process/HEAT_TRANSPORT_BHE/borehole_heat_exchangers/borehole_heat_exchanger/grout/t_porosity.md

+Porosity of the grout. We always would like to have a kind of grout with low porosity because the efficiency of the heat conduction will be higher. In grout, there is an assume that no heat advection exists.