Skip to content
Snippets Groups Projects
Commit 993024ab authored by Max Jaeschke's avatar Max Jaeschke Committed by Dmitri Naumov
Browse files

Add tests for BHE id feature

parent d19243b0
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 639 additions and 0 deletions
ProcessLib/HeatTransportBHE/Tests.cmake
+ 54
− 0
View file @ 993024ab
...@@ -212,6 +212,60 @@ AddTest( ...@@ -212,6 +212,60 @@ AddTest(
3bhes_1_ts_10_t_600.000000.vtu 3bhes_1U_BHE_1_mesh_ts_10_t_600.000000.vtu temperature_BHE1 temperature_BHE1 2e-13 1e-12 3bhes_1_ts_10_t_600.000000.vtu 3bhes_1U_BHE_1_mesh_ts_10_t_600.000000.vtu temperature_BHE1 temperature_BHE1 2e-13 1e-12
) )
AddTest(
NAME HeatTransportBHE_3D_3BHEs_BHE_id
PATH Parabolic/T/3D_3BHEs
RUNTIME 14
EXECUTABLE ogs
EXECUTABLE_ARGS 3bhes_id_1U_2U_1U.prj
WRAPPER time
TESTER vtkdiff
REQUIREMENTS NOT OGS_USE_MPI
DIFF_DATA
3bhes_ts_10_t_600.000000.vtu 3bhes_id_ts_10_t_600.000000.vtu temperature_BHE1 temperature_BHE1 0 1e-12
3bhes_ts_10_t_600.000000.vtu 3bhes_id_ts_10_t_600.000000.vtu temperature_BHE2 temperature_BHE2 0 1e-12
3bhes_ts_10_t_600.000000.vtu 3bhes_id_ts_10_t_600.000000.vtu temperature_BHE3 temperature_BHE3 0 1e-12
3bhes_ts_10_t_600.000000.vtu 3bhes_id_ts_10_t_600.000000.vtu temperature_soil temperature_soil 0 1e-12
)
AddTest(
NAME HeatTransportBHE_3D_3BHEs_BHE_id_failcase_id_out_of_range
PATH Parabolic/T/3D_3BHEs
RUNTIME 1
EXECUTABLE ogs
EXECUTABLE_ARGS id_out_of_range.xml
REQUIREMENTS NOT OGS_USE_MPI
PROPERTIES
PASS_REGULAR_EXPRESSION "The maximum given BHE id '100' did not match the number of given BHE definitions '3'. The BHE ids needs to be defined starting from 0, so the maximum BHE id needs to be number of BHE definitions minus 1. After all definitions there are no gaps allowed between the given ids."
)
AddTest(
NAME HeatTransportBHE_3D_3BHEs_BHE_id_failcase_duplicate_id
PATH Parabolic/T/3D_3BHEs
RUNTIME 1
EXECUTABLE ogs
EXECUTABLE_ARGS duplicate_id.xml
REQUIREMENTS NOT OGS_USE_MPI
PROPERTIES
PASS_REGULAR_EXPRESSION "BHE with id '1' is already present in the list! Check for duplicate definitions of BHE ids."
)
AddTest(
NAME HeatTransportBHE_3D_3BHEs_BHE_id_1U
PATH Parabolic/T/3D_3BHEs
RUNTIME 14
EXECUTABLE ogs
EXECUTABLE_ARGS 3bhes_id_1U.prj
WRAPPER time
TESTER vtkdiff
REQUIREMENTS NOT OGS_USE_MPI
DIFF_DATA
3bhes_1U_ts_10_t_600.000000.vtu 3bhes_1U_id_ts_10_t_600.000000.vtu temperature_BHE1 temperature_BHE1 0 1e-12
3bhes_1U_ts_10_t_600.000000.vtu 3bhes_1U_id_ts_10_t_600.000000.vtu temperature_BHE2 temperature_BHE2 0 1e-12
3bhes_1U_ts_10_t_600.000000.vtu 3bhes_1U_id_ts_10_t_600.000000.vtu temperature_BHE3 temperature_BHE3 0 1e-12
3bhes_1U_ts_10_t_600.000000.vtu 3bhes_1U_id_ts_10_t_600.000000.vtu temperature_soil temperature_soil 0 1e-12
)
AddTest( AddTest(
NAME HeatTransportBHE_3D_BHE_groundwater_advection NAME HeatTransportBHE_3D_BHE_groundwater_advection
PATH Parabolic/T/3D_BHE_GW_advection PATH Parabolic/T/3D_BHE_GW_advection
......
Tests/Data/Parabolic/T/3D_3BHEs/3bhes_1U_ts_10_t_600.000000.vtu 0 → 100644
+ 29
− 0
View file @ 993024ab
Source diff could not be displayed: it is too large. Options to address this: view the blob.
Tests/Data/Parabolic/T/3D_3BHEs/3bhes_id_1U.prj 0 → 100644
+ 252
− 0
View file @ 993024ab
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<meshes>
<mesh>3bhes.vtu</mesh>
</meshes>
<processes>
<process>
<name>HeatTransportBHE</name>
<type>HEAT_TRANSPORT_BHE</type>
<integration_order>2</integration_order>
<process_variables>
<process_variable>temperature_soil</process_variable>
<process_variable>temperature_BHE1</process_variable>
<process_variable>temperature_BHE2</process_variable>
<process_variable>temperature_BHE3</process_variable>
</process_variables>
<borehole_heat_exchangers>
<borehole_heat_exchanger id="*">
<type>1U</type>
<flow_and_temperature_control>
<type>TemperatureCurveConstantFlow</type>
<flow_rate>2.0e-4</flow_rate>
<temperature_curve>inflow_temperature</temperature_curve>
</flow_and_temperature_control>
<borehole>
<length>50.0</length>
<diameter>0.126</diameter>
</borehole>
<grout>
<density>2190.0</density>
<porosity>0.0</porosity>
<specific_heat_capacity>1735.160</specific_heat_capacity>
<thermal_conductivity>0.806</thermal_conductivity>
</grout>
<pipes>
<inlet>
<diameter> 0.013665</diameter>
<wall_thickness>0.003035</wall_thickness>
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
</inlet>
<outlet>
<diameter>0.013665</diameter>
<wall_thickness>0.003035</wall_thickness>
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
</outlet>
<distance_between_pipes>0.053</distance_between_pipes>
<longitudinal_dispersion_length>0.001</longitudinal_dispersion_length>
</pipes>
<refrigerant>
<density>992.92</density>
<viscosity>0.00067418</viscosity>
<specific_heat_capacity>4198</specific_heat_capacity>
<thermal_conductivity>0.62863</thermal_conductivity>
<reference_temperature>298.15</reference_temperature>
</refrigerant>
</borehole_heat_exchanger>
</borehole_heat_exchangers>
</process>
</processes>
<media>
<medium id="0">
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<name>phase_velocity</name>
<type>Constant</type>
<value>0 0 0</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>4068</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>992.92</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>1780</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>1120</value>
</property>
</properties>
</phase>
<phase>
<type>Gas</type>
<properties>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>1000</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>2500</value>
</property>
</properties>
</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>
<value>0</value>
</property>
<property>
<name>thermal_transversal_dispersivity</name>
<type>Constant</type>
<value>0</value>
</property>
</properties>
</medium>
</media>
<time_loop>
<processes>
<process ref="HeatTransportBHE">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1.e-10</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial> 0.0 </t_initial>
<t_end> 600 </t_end>
<timesteps>
<pair>
<repeat>10</repeat>
<delta_t>60</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>3bhes_1U_id</prefix>
<timesteps>
<pair>
<repeat> 10</repeat>
<each_steps> 1 </each_steps>
</pair>
</timesteps>
<variables>
<variable>temperature_soil</variable>
<variable>temperature_BHE1</variable>
<variable>temperature_BHE2</variable>
<variable>temperature_BHE3</variable>
</variables>
<suffix>_ts_{:timestep}_t_{:time}</suffix>
</output>
</time_loop>
<parameters>
<parameter>
<name>T0</name>
<type>Constant</type>
<value>305</value>
</parameter>
<parameter>
<name>T0_BHE_1</name>
<type>Constant</type>
<values>295 295 295 295</values>
</parameter>
</parameters>
<process_variables>
<process_variable>
<name>temperature_soil</name>
<components>1</components>
<order>1</order>
<initial_condition>T0</initial_condition>
</process_variable>
<process_variable>
<name>temperature_BHE1</name>
<components>4</components>
<order>1</order>
<initial_condition>T0_BHE_1</initial_condition>
</process_variable>
<process_variable>
<name>temperature_BHE2</name>
<components>4</components>
<order>1</order>
<initial_condition>T0_BHE_1</initial_condition>
</process_variable>
<process_variable>
<name>temperature_BHE3</name>
<components>4</components>
<order>1</order>
<initial_condition>T0_BHE_1</initial_condition>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_picard</name>
<type>Picard</type>
<max_iter>100</max_iter>
<linear_solver>general_linear_solver</linear_solver>
</nonlinear_solver>
</nonlinear_solvers>
<linear_solvers>
<linear_solver>
<name>general_linear_solver</name>
<lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
<eigen>
<solver_type>BiCGSTAB</solver_type>
<precon_type>ILUT</precon_type>
<max_iteration_step>1000</max_iteration_step>
<error_tolerance>1e-16</error_tolerance>
</eigen>
<petsc>
<prefix>gw</prefix>
<parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
</petsc>
</linear_solver>
</linear_solvers>
<curves>
<curve>
<name>inflow_temperature</name>
<coords>0 600
</coords>
<values>295 295
</values>
</curve>
</curves>
</OpenGeoSysProject>
Tests/Data/Parabolic/T/3D_3BHEs/3bhes_id_1U_2U_1U.prj 0 → 100644
+ 296
− 0
View file @ 993024ab
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<meshes>
<mesh>3bhes.vtu</mesh>
</meshes>
<processes>
<process>
<name>HeatTransportBHE</name>
<type>HEAT_TRANSPORT_BHE</type>
<integration_order>2</integration_order>
<process_variables>
<process_variable>temperature_soil</process_variable>
<process_variable>temperature_BHE1</process_variable>
<process_variable>temperature_BHE2</process_variable>
<process_variable>temperature_BHE3</process_variable>
</process_variables>
<borehole_heat_exchangers>
<borehole_heat_exchanger id="0,2">
<type>1U</type>
<flow_and_temperature_control>
<type>TemperatureCurveConstantFlow</type>
<flow_rate>2.0e-4</flow_rate>
<temperature_curve>inflow_temperature</temperature_curve>
</flow_and_temperature_control>
<borehole>
<length>50.0</length>
<diameter>0.126</diameter>
</borehole>
<grout>
<density>2190.0</density>
<porosity>0.0</porosity>
<specific_heat_capacity>1735.160</specific_heat_capacity>
<thermal_conductivity>0.806</thermal_conductivity>
</grout>
<pipes>
<inlet>
<diameter> 0.013665</diameter>
<wall_thickness>0.003035</wall_thickness>
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
</inlet>
<outlet>
<diameter>0.013665</diameter>
<wall_thickness>0.003035</wall_thickness>
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
</outlet>
<distance_between_pipes>0.053</distance_between_pipes>
<longitudinal_dispersion_length>0.001</longitudinal_dispersion_length>
</pipes>
<refrigerant>
<density>992.92</density>
<viscosity>0.00067418</viscosity>
<specific_heat_capacity>4198</specific_heat_capacity>
<thermal_conductivity>0.62863</thermal_conductivity>
<reference_temperature>298.15</reference_temperature>
</refrigerant>
</borehole_heat_exchanger>
<borehole_heat_exchanger id="1">
<type>2U</type>
<flow_and_temperature_control>
<type>TemperatureCurveConstantFlow</type>
<flow_rate>2.0e-4</flow_rate>
<temperature_curve>inflow_temperature</temperature_curve>
</flow_and_temperature_control>
<borehole>
<length>50.0</length>
<diameter>0.126</diameter>
</borehole>
<grout>
<density>2190.0</density>
<porosity>0.0</porosity>
<specific_heat_capacity>1735.160</specific_heat_capacity>
<thermal_conductivity>0.806</thermal_conductivity>
</grout>
<pipes>
<inlet>
<diameter> 0.013665</diameter>
<wall_thickness>0.003035</wall_thickness>
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
</inlet>
<outlet>
<diameter>0.013665</diameter>
<wall_thickness>0.003035</wall_thickness>
<wall_thermal_conductivity>0.39</wall_thermal_conductivity>
</outlet>
<distance_between_pipes>0.053</distance_between_pipes>
<longitudinal_dispersion_length>0.001</longitudinal_dispersion_length>
</pipes>
<refrigerant>
<density>992.92</density>
<viscosity>0.00067418</viscosity>
<specific_heat_capacity>4198</specific_heat_capacity>
<thermal_conductivity>0.62863</thermal_conductivity>
<reference_temperature>298.15</reference_temperature>
</refrigerant>
</borehole_heat_exchanger>
</borehole_heat_exchangers>
</process>
</processes>
<media>
<medium id="0">
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<name>phase_velocity</name>
<type>Constant</type>
<value>0 0 0</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>4068</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>992.92</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>1780</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>1120</value>
</property>
</properties>
</phase>
<phase>
<type>Gas</type>
<properties>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>1000</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>2500</value>
</property>
</properties>
</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>
<value>0</value>
</property>
<property>
<name>thermal_transversal_dispersivity</name>
<type>Constant</type>
<value>0</value>
</property>
</properties>
</medium>
</media>
<time_loop>
<processes>
<process ref="HeatTransportBHE">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1.e-10</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial> 0.0 </t_initial>
<t_end> 600 </t_end>
<timesteps>
<pair>
<repeat>10</repeat>
<delta_t>60</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>3bhes_id</prefix>
<timesteps>
<pair>
<repeat> 10</repeat>
<each_steps> 1 </each_steps>
</pair>
</timesteps>
<variables>
<variable>temperature_soil</variable>
<variable>temperature_BHE1</variable>
<variable>temperature_BHE2</variable>
<variable>temperature_BHE3</variable>
</variables>
<suffix>_ts_{:timestep}_t_{:time}</suffix>
</output>
</time_loop>
<parameters>
<parameter>
<name>T0</name>
<type>Constant</type>
<value>305</value>
</parameter>
<parameter>
<name>T0_BHE_1</name>
<type>Constant</type>
<values>295 295 295 295</values>
</parameter>
<parameter>
<name>T0_BHE_2</name>
<type>Constant</type>
<values>295 295 295 295 295 295 295 295</values>
</parameter>
</parameters>
<process_variables>
<process_variable>
<name>temperature_soil</name>
<components>1</components>
<order>1</order>
<initial_condition>T0</initial_condition>
</process_variable>
<process_variable>
<name>temperature_BHE1</name>
<components>4</components>
<order>1</order>
<initial_condition>T0_BHE_1</initial_condition>
</process_variable>
<process_variable>
<name>temperature_BHE2</name>
<components>8</components>
<order>1</order>
<initial_condition>T0_BHE_2</initial_condition>
</process_variable>
<process_variable>
<name>temperature_BHE3</name>
<components>4</components>
<order>1</order>
<initial_condition>T0_BHE_1</initial_condition>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_picard</name>
<type>Picard</type>
<max_iter>100</max_iter>
<linear_solver>general_linear_solver</linear_solver>
</nonlinear_solver>
</nonlinear_solvers>
<linear_solvers>
<linear_solver>
<name>general_linear_solver</name>
<lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
<eigen>
<solver_type>BiCGSTAB</solver_type>
<precon_type>ILUT</precon_type>
<max_iteration_step>1000</max_iteration_step>
<error_tolerance>1e-16</error_tolerance>
</eigen>
<petsc>
<prefix>gw</prefix>
<parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
</petsc>
</linear_solver>
</linear_solvers>
<curves>
<curve>
<name>inflow_temperature</name>
<coords>0 600
</coords>
<values>295 295
</values>
</curve>
</curves>
</OpenGeoSysProject>
Tests/Data/Parabolic/T/3D_3BHEs/duplicate_id.xml 0 → 100644
+ 4
− 0
View file @ 993024ab
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProjectDiff base_file="3bhes_id_1U_2U_1U.prj">
<replace sel="/*/processes/process/borehole_heat_exchangers/borehole_heat_exchanger[1]/@id">0,1</replace>
</OpenGeoSysProjectDiff>
Tests/Data/Parabolic/T/3D_3BHEs/id_out_of_range.xml 0 → 100644
+ 4
− 0
View file @ 993024ab
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProjectDiff base_file="3bhes_id_1U_2U_1U.prj">
<replace sel="/*/processes/process/borehole_heat_exchangers/borehole_heat_exchanger[1]/@id">0,100</replace>
</OpenGeoSysProjectDiff>
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment