Commit 1be7d8fa authored by Boyan Meng's avatar Boyan Meng
Browse files

Add test files of BHE with groundwater advection feature

replace orginal benchmark files with new ones using smaller mesh and shorter simulation time (+2 squashed commit)
use four spaces instead of tab in all prj files
reformat prj files using xmlstarlet
add ctest

Squashed commit:

[adb3f6bfa] add ctest

[f711c1f8d] reformat prj files using xmlstarlet (+1 squashed commits)

Squashed commits:

[a89b72dc1] use four spaces instead of tab in all prj files (+1 squashed commits)

Squashed commits:

[352748bb4] add the benchmark test (+1 squashed commits)

Squashed commits:

[a22a66353] add the benchmark test
parent 71543c2e
......@@ -64,3 +64,17 @@ AddTest(
3D_2U_BHE_pcs_0_ts_10_t_600.000000.vtu 3D_2U_BHE_pcs_0_ts_10_t_600.000000.vtu temperature_BHE1 temperature_BHE1 1e-12 1e-14
3D_2U_BHE_pcs_0_ts_10_t_600.000000.vtu 3D_2U_BHE_pcs_0_ts_10_t_600.000000.vtu temperature_soil temperature_soil 1e-12 1e-13
)
AddTest(
NAME HeatTransportBHE_3D_BHE_groundwater_advection
PATH Parabolic/T/3D_BHE_GW_advection
RUNTIME 30
EXECUTABLE ogs
EXECUTABLE_ARGS BHE_GW_advection.prj
WRAPPER time
TESTER vtkdiff
REQUIREMENTS NOT OGS_USE_MPI
DIFF_DATA
BHE_GW_advection_pcs_0_ts_10_t_500.000000.vtu BHE_GW_advection_pcs_0_ts_10_t_500.000000.vtu temperature_BHE1 temperature_BHE1 1e-12 1e-14
BHE_GW_advection_pcs_0_ts_10_t_500.000000.vtu BHE_GW_advection_pcs_0_ts_10_t_500.000000.vtu temperature_soil temperature_soil 1e-12 1e-13
)
......@@ -60,7 +60,7 @@
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<property>
<name>phase_velocity</name>
<type>Constant</type>
<value>0 0 0</value>
......@@ -109,15 +109,15 @@
</phase>
</phases>
<properties>
<property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>2.78018</value>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>2.78018</value>
</property>
<property>
<name>thermal_longitudinal_dispersivity</name>
......
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<meshes>
<mesh>BHE_GW_advection.vtu</mesh>
<mesh>front.vtu</mesh>
<mesh>back.vtu</mesh>
<mesh>left.vtu</mesh>
<mesh>right.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_variables>
<borehole_heat_exchangers>
<borehole_heat_exchanger>
<type>1U</type>
<flow_and_temperature_control>
<type>FixedPowerConstantFlow</type>
<flow_rate>2.0e-4</flow_rate>
<power>300</power>
</flow_and_temperature_control>
<borehole>
<length>15.0</length>
<diameter>0.13</diameter>
</borehole>
<grout>
<density>2190.0</density>
<porosity>0.0</porosity>
<heat_capacity>1735.160</heat_capacity>
<thermal_conductivity>0.73</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>4068</specific_heat_capacity>
<thermal_conductivity>0.62863</thermal_conductivity>
<reference_temperature>25</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 1e-7 0</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>4200</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>1000</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>880</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>2650</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.26</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>2.5</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>1e-6</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial> 0.0 </t_initial>
<t_end> 500 </t_end>
<timesteps>
<pair>
<repeat>10</repeat>
<delta_t>50</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>BHE_GW_advection</prefix>
<timesteps>
<pair>
<repeat> 1</repeat>
<each_steps> 10 </each_steps>
</pair>
</timesteps>
<variables>
<variable>temperature_soil</variable>
<variable>temperature_BHE1</variable>
</variables>
</output>
</time_loop>
<parameters>
<parameter>
<name>T0</name>
<type>Constant</type>
<value>25</value>
</parameter>
<parameter>
<name>T_DBC</name>
<type>Constant</type>
<value>25</value>
</parameter>
<parameter>
<name>T0_BHE1</name>
<type>Constant</type>
<values>25.36 25.13 25.23 25.115</values>
</parameter>
</parameters>
<process_variables>
<process_variable>
<name>temperature_soil</name>
<components>1</components>
<order>1</order>
<initial_condition>T0</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>front</mesh>
<type>Dirichlet</type>
<parameter>T_DBC</parameter>
</boundary_condition>
<boundary_condition>
<mesh>back</mesh>
<type>Dirichlet</type>
<parameter>T_DBC</parameter>
</boundary_condition>
<boundary_condition>
<mesh>left</mesh>
<type>Dirichlet</type>
<parameter>T_DBC</parameter>
</boundary_condition>
<boundary_condition>
<mesh>right</mesh>
<type>Dirichlet</type>
<parameter>T_DBC</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>temperature_BHE1</name>
<components>4</components>
<order>1</order>
<initial_condition>T0_BHE1</initial_condition>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_picard</name>
<type>Picard</type>
<max_iter>50</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>
</OpenGeoSysProject>
......@@ -60,7 +60,7 @@
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<property>
<name>phase_velocity</name>
<type>Constant</type>
<value>0 0 0</value>
......@@ -109,12 +109,12 @@
</phase>
</phases>
<properties>
<property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>2.78018</value>
......
......@@ -54,13 +54,13 @@
</borehole_heat_exchangers>
</process>
</processes>
<media>
<media>
<medium id="0">
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<property>
<name>phase_velocity</name>
<type>Constant</type>
<value>0 0 0</value>
......@@ -92,7 +92,7 @@
</property>
</properties>
</phase>
<phase>
<phase>
<type>Gas</type>
<properties>
<property>
......@@ -109,12 +109,12 @@
</phase>
</phases>
<properties>
<property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>2.78018</value>
......
......@@ -59,7 +59,7 @@
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<property>
<name>phase_velocity</name>
<type>Constant</type>
<value>0 0 0</value>
......@@ -108,12 +108,12 @@
</phase>
</phases>
<properties>
<property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>2.5</value>
......
......@@ -59,7 +59,7 @@
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<property>
<name>phase_velocity</name>
<type>Constant</type>
<value>0 0 0</value>
......@@ -108,12 +108,12 @@
</phase>
</phases>
<properties>
<property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>2.5</value>
......
Markdown is supported
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