diff --git a/ProcessLib/HT/Tests.cmake b/ProcessLib/HT/Tests.cmake
index f25e2368e2e96123b2ca8297885c66ac232fa251..cc2b9c39e4339fa6ba16da62f37def1675be5c76 100644
--- a/ProcessLib/HT/Tests.cmake
+++ b/ProcessLib/HT/Tests.cmake
@@ -25,6 +25,20 @@ AddTest(
VIS ConstViscosityThermalConvection_pcs_1_ts_149_t_50000000000.000000.vtu
)
+AddTest(
+ NAME LARGE_2D_Adaptive_dt_ThermalConvection_constviscosityStaggeredScheme
+ PATH Parabolic/HT/StaggeredCoupling/ConstViscosity
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS square_5500x5500_staggered_scheme_adaptive_dt.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT (OGS_USE_LIS OR OGS_USE_MPI)
+ DIFF_DATA
+ square_5500x5500.vtu ConstViscosityThermalConvectionStaggeredAdaptive_dt_pcs_1_ts_135_t_50000000000.000000.vtu T_ref T 1e-3 1.e-3
+ square_5500x5500.vtu ConstViscosityThermalConvectionStaggeredAdaptive_dt_pcs_1_ts_135_t_50000000000.000000.vtu p_ref p 1e-3 1.e-3
+ VIS ConstViscosityThermalConvectionStaggeredAdaptive_dt_pcs_1_ts_135_t_50000000000.000000.vtu
+)
+
# MPI/PETSc tests
AddTest(
NAME Parallel_LARGE_2D_ThermalConvection_constviscosity
@@ -58,3 +72,17 @@ AddTest(
square_5500x5500.vtu ConstViscosityThermalConvection_pcs_1_ts_149_t_50000000000_000000_0.vtu p_ref p 1e-14 1.e-14
VIS ConstViscosityThermalConvection_pcs_1_ts_149_t_50000000000.000000_0.vtu
)
+
+AddTest(
+ NAME LARGE_2D_Adaptive_dt_ThermalConvection_constviscosityStaggeredScheme
+ PATH Parabolic/HT/StaggeredCoupling/ConstViscosity
+ EXECUTABLE_ARGS square_5500x5500_staggered_scheme_adaptive_dt.prj
+ WRAPPER mpirun
+ WRAPPER_ARGS -np 1
+ TESTER vtkdiff
+ REQUIREMENTS OGS_USE_MPI
+ DIFF_DATA
+ square_5500x5500.vtu ConstViscosityThermalConvectionStaggeredAdaptive_dt_pcs_1_ts_135_t_50000000000_000000_0.vtu T_ref T 1e-3 1.e-3
+ square_5500x5500.vtu ConstViscosityThermalConvectionStaggeredAdaptive_dt_pcs_1_ts_135_t_50000000000_000000_0.vtu p_ref p 1e-3 1.e-3
+ VIS ConstViscosityThermalConvectionStaggeredAdaptive_dt_pcs_1_ts_135_t_50000000000_000000_0.vtu
+)
diff --git a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/quad_5500x5500_adaptive_dt_constant_viscosity.prj b/Tests/Data/Parabolic/HT/StaggeredCoupling/ConstViscosity/square_5500x5500_staggered_scheme_adaptive_dt.prj
similarity index 58%
rename from Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/quad_5500x5500_adaptive_dt_constant_viscosity.prj
rename to Tests/Data/Parabolic/HT/StaggeredCoupling/ConstViscosity/square_5500x5500_staggered_scheme_adaptive_dt.prj
index bfbeaad34c2a3a769c835edcd56b9801ee06f281..0eb5b012447bf38ea90743e637299c5f7b915637 100644
--- a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/quad_5500x5500_adaptive_dt_constant_viscosity.prj
+++ b/Tests/Data/Parabolic/HT/StaggeredCoupling/ConstViscosity/square_5500x5500_staggered_scheme_adaptive_dt.prj
@@ -1,103 +1,95 @@
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<mesh>square_5500x5500.vtu</mesh>
- <geometry>quad_5500x5500.gml</geometry>
+ <geometry>square_5500x5500.gml</geometry>
<processes>
<process>
- <name>HeatConduction</name>
- <type>HEAT_CONDUCTION</type>
+ <name>ConstViscosityThermalConvection</name>
+ <type>HT</type>
+ <coupling_scheme>staggered</coupling_scheme>
<integration_order>2</integration_order>
- <thermal_conductivity>lambda_solid</thermal_conductivity>
- <heat_capacity>c_solid</heat_capacity>
- <density>rho_solid</density>
<process_variables>
- <process_variable>temperature</process_variable>
+ <temperature>T</temperature>
+ <pressure>p</pressure>
</process_variables>
- <secondary_variables>
- <secondary_variable type="static" internal_name="heat_flux_x" output_name="heat_flux_x"/>
- <secondary_variable type="static" internal_name="heat_flux_y" output_name="heat_flux_y"/>
- </secondary_variables>
- </process>
- <process>
- <name>LiquidFlow</name>
- <type>LIQUID_FLOW</type>
- <integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>9.81</g>
- </darcy_gravity>
- <process_variables>
- <process_variable>pressure</process_variable>
- </process_variables>
- <secondary_variables>
- <secondary_variable type="static" internal_name="darcy_velocity" output_name="v"/>
- </secondary_variables>
- <material_property>
- <fluid>
- <density>
- <type>TemperatureDependent</type>
- <rho0>1000</rho0>
- <temperature0>20</temperature0>
- <beta>4.3e-4</beta>
- </density>
- <viscosity>
+ <fluid>
+ <density>
+ <type>TemperatureDependent</type>
+ <rho0>1000</rho0>
+ <temperature0>20</temperature0>
+ <beta>4.3e-4</beta>
+ </density>
+ <viscosity>
+ <type>Constant</type>
+ <value>1.0e-3</value>
+ </viscosity>
+ <specific_heat_capacity>
+ <type>Constant</type>
+ <value>4200</value>
+ </specific_heat_capacity>
+ </fluid>
+ <porous_medium>
+ <porous_medium id="0">
+ <permeability>
+ <permeability_tensor_entries>kappa1</permeability_tensor_entries>
<type>Constant</type>
- <value>1.0e-3</value>
- </viscosity>
- <specific_heat_capacity>
+ </permeability>
+ <porosity>
<type>Constant</type>
- <value> 4200. </value>
- </specific_heat_capacity>
- <thermal_conductivity>
+ <porosity_parameter>constant_porosity_parameter</porosity_parameter>
+ </porosity>
+ <storage>
<type>Constant</type>
- <value> .65 </value>
- </thermal_conductivity>
- </fluid>
- <porous_medium>
- <porous_medium id="0">
- <permeability>
- <permeability_tensor_entries>kappa1</permeability_tensor_entries>
- <type>Constant</type>
- </permeability>
- <porosity>
- <type>Constant</type>
- <porosity_parameter>constant_porosity_parameter</porosity_parameter>
- </porosity>
- <storage>
- <type>Constant</type>
- <value>0.0</value>
- </storage>
- </porous_medium>
+ <value>0.0</value>
+ </storage>
</porous_medium>
- <solid>
- <thermal_expansion>solid_thermal_expansion</thermal_expansion>
- <biot_constant>biot_const</biot_constant>
- </solid>
- </material_property>
+ </porous_medium>
+ <density_solid>rho_solid</density_solid>
+ <specific_heat_capacity_solid>c_solid</specific_heat_capacity_solid>
+ <thermal_conductivity_solid>lambda_solid</thermal_conductivity_solid>
+ <thermal_conductivity_fluid>lambda_fluid</thermal_conductivity_fluid>
+ <specific_body_force>0 -9.81</specific_body_force>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="darcy_velocity" output_name="darcy_velocity"/>
+ </secondary_variables>
</process>
</processes>
<time_loop>
<global_process_coupling>
- <max_iter> 5 </max_iter>
- <convergence_criterion>
- <type>Residual</type>
- <norm_type>NORM2</norm_type>
- <reltol>1.e-8</reltol>
- </convergence_criterion>
+ <max_iter> 6 </max_iter>
+ <convergence_criteria>
+ <!-- convergence criterion for the first process -->
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1.e-14</reltol>
+ </convergence_criterion>
+ <!-- convergence criterion for the second process -->
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1.e-14</reltol>
+ </convergence_criterion>
+ </convergence_criteria>
</global_process_coupling>
+
<processes>
- <process ref="HeatConduction">
- <nonlinear_solver>basic_picard_T</nonlinear_solver>
+ <process ref="ConstViscosityThermalConvection">
+ <nonlinear_solver>picard_T</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
- <abstol>1.e-6</abstol>
+ <abstol>1.e-3</abstol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
+ <output>
+ <variables>
+ <variable>T</variable>
+ <variable>darcy_velocity</variable>
+ </variables>
+ </output>
<time_stepping>
<type>EvolutionaryPIDcontroller</type>
<t_initial> 0.0 </t_initial>
@@ -107,30 +99,25 @@
<dt_max> 4e8 </dt_max>
<rel_dt_min> 0.1 </rel_dt_min>
<rel_dt_max> 10 </rel_dt_max>
+ <tol> 10.0 </tol>
<specific_times>5.0e9 2.e10</specific_times>
- <tol> 1.e-2 </tol>
</time_stepping>
- <coupled_processes>
- <coupled_process> LiquidFlow </coupled_process>
- </coupled_processes>
- <output>
- <variables>
- <variable> temperature </variable>
- <variable> heat_flux_x </variable>
- <variable> heat_flux_y </variable>
- </variables>
- </output>
</process>
- <process ref="LiquidFlow">
- <nonlinear_solver>basic_picard_H</nonlinear_solver>
+ <process ref="ConstViscosityThermalConvection">
+ <nonlinear_solver>picard_H</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
- <abstol>1.e-6</abstol>
+ <abstol>1.e-3</abstol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
+ <output>
+ <variables>
+ <variable>p</variable>
+ </variables>
+ </output>
<time_stepping>
<type>EvolutionaryPIDcontroller</type>
<t_initial> 0.0 </t_initial>
@@ -143,39 +130,20 @@
<tol> 10.0 </tol>
<specific_times>5.0e9 2.e10</specific_times>
</time_stepping>
- <coupled_processes>
- <coupled_process> HeatConduction </coupled_process>
- </coupled_processes>
- <output>
- <variables>
- <variable> pressure </variable>
- <variable> v </variable>
- </variables>
- </output>
</process>
</processes>
<output>
<type>VTK</type>
- <prefix>ConstViscosityThermalConvection</prefix>
+ <prefix>ConstViscosityThermalConvectionStaggeredAdaptive_dt</prefix>
<timesteps>
<pair>
<repeat> 1 </repeat>
- <each_steps> 137 </each_steps>
+ <each_steps> 149 </each_steps>
</pair>
</timesteps>
</output>
</time_loop>
<parameters>
- <parameter>
- <name>poro</name>
- <type>Constant</type>
- <value>0.001</value>
- </parameter>
- <parameter>
- <name>kappa</name>
- <type>Constant</type>
- <value>1e-14</value>
- </parameter>
<parameter>
<name>rho_solid</name>
<type>Constant</type>
@@ -197,24 +165,19 @@
<value>0.65</value>
</parameter>
<parameter>
- <name>c_solid</name>
- <type>Constant</type>
- <value>0</value>
- </parameter>
- <parameter>
- <name>c_fluid</name>
+ <name>alpha_l</name>
<type>Constant</type>
- <value>4200</value>
+ <value>0.0</value>
</parameter>
<parameter>
- <name>solid_thermal_expansion</name>
+ <name>alpha_t</name>
<type>Constant</type>
- <value>0e-5</value>
+ <value>0.0</value>
</parameter>
<parameter>
- <name>biot_const</name>
+ <name>c_solid</name>
<type>Constant</type>
- <value>1.0</value>
+ <value>0</value>
</parameter>
<parameter>
<name>T0</name>
@@ -242,19 +205,19 @@
<value>20</value>
</parameter>
<parameter>
- <name>kappa1</name>
+ <name>constant_porosity_parameter</name>
<type>Constant</type>
- <values>1.0e-14</values>
+ <value>0.001</value>
</parameter>
<parameter>
- <name>constant_porosity_parameter</name>
+ <name>kappa1</name>
<type>Constant</type>
- <value>0.001</value>
+ <values>1.e-14 0 0 1.e-14</values>
</parameter>
</parameters>
<process_variables>
<process_variable>
- <name>temperature</name>
+ <name>T</name>
<components>1</components>
<order>1</order>
<initial_condition>T0</initial_condition>
@@ -274,7 +237,7 @@
</boundary_conditions>
</process_variable>
<process_variable>
- <name>pressure</name>
+ <name>p</name>
<components>1</components>
<order>1</order>
<initial_condition>P0</initial_condition>
@@ -290,27 +253,28 @@
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
- <name>basic_picard_T</name>
+ <name>picard_T</name>
<type>Picard</type>
- <max_iter>100</max_iter>
+ <max_iter>20</max_iter>
<linear_solver>linear_solver_T</linear_solver>
</nonlinear_solver>
<nonlinear_solver>
- <name>basic_picard_H</name>
+ <name>picard_H</name>
<type>Picard</type>
- <max_iter>100</max_iter>
+ <max_iter>20</max_iter>
<linear_solver>linear_solver_H</linear_solver>
</nonlinear_solver>
</nonlinear_solvers>
+
<linear_solvers>
<linear_solver>
<name>linear_solver_T</name>
- <lis>-i bicgstab -p jacobi -tol 1e-14 -maxiter 3000</lis>
+ <lis>-i bicgstab -p ilu -tol 1e-20 -maxiter 10000</lis>
<eigen>
<solver_type>BiCGSTAB</solver_type>
- <precon_type>DIAGONAL</precon_type>
+ <precon_type>ILUT</precon_type>
<max_iteration_step>10000</max_iteration_step>
- <error_tolerance>1e-14</error_tolerance>
+ <error_tolerance>1e-20</error_tolerance>
</eigen>
<petsc>
<prefix>T</prefix>
@@ -319,16 +283,16 @@
</linear_solver>
<linear_solver>
<name>linear_solver_H</name>
- <lis>-i bicgstab -p jacobi -tol 1e-14 -maxiter 3000</lis>
+ <lis>-i bicgstab -p jacobi -tol 1e-14 -maxiter 10000</lis>
<eigen>
<solver_type>BiCGSTAB</solver_type>
<precon_type>DIAGONAL</precon_type>
<max_iteration_step>10000</max_iteration_step>
- <error_tolerance>1e-16</error_tolerance>
+ <error_tolerance>1e-14</error_tolerance>
</eigen>
<petsc>
- <prefix>gw</prefix>
- <parameters>-gw_ksp_type bcgs -gw_pc_type hypre -gw_pc_hypre_type boomeramg -gw_ksp_rtol 1e-16 -gw_ksp_max_it 3000</parameters>
+ <prefix>H</prefix>
+ <parameters>-H_ksp_type bcgs -H_pc_type hypre -H_pc_hypre_type boomeramg -H_ksp_rtol 1e-16 -H_ksp_max_it 3000</parameters>
</petsc>
</linear_solver>
</linear_solvers>
diff --git a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/ConstViscosityThermalConvection_pcs_1_ts_137_t_50000000000.000000.vtu b/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/ConstViscosityThermalConvection_pcs_1_ts_137_t_50000000000.000000.vtu
deleted file mode 100644
index 4d91bcf99c23d2dd4b56147cad24951419c065b1..0000000000000000000000000000000000000000
--- a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/ConstViscosityThermalConvection_pcs_1_ts_137_t_50000000000.000000.vtu
+++ /dev/null
@@ -1,3 +0,0 @@
-version https://git-lfs.github.com/spec/v1
-oid sha256:257dbd5da28f0613e39dc0d1431ad0e3def245cf4f0bb744e083c3895a9555be
-size 69833
diff --git a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/ThermalConvection_pcs_1_ts_232_t_50000000000.000000_non_const_mu.vtu b/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/ThermalConvection_pcs_1_ts_232_t_50000000000.000000_non_const_mu.vtu
deleted file mode 100644
index c524f6e18adeabfe0af11073db819d8abaab5d59..0000000000000000000000000000000000000000
--- a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/ThermalConvection_pcs_1_ts_232_t_50000000000.000000_non_const_mu.vtu
+++ /dev/null
@@ -1,3 +0,0 @@
-version https://git-lfs.github.com/spec/v1
-oid sha256:203c376f2e5dafb6252abde3148d7ff981cb135dfa4dd8fb57bd109c1b8c92b5
-size 69893
diff --git a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/quad_5500x5500.gml b/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/quad_5500x5500.gml
deleted file mode 100644
index 032286fc8f97784b6f71509a15840125436e14fe..0000000000000000000000000000000000000000
--- a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/quad_5500x5500.gml
+++ /dev/null
@@ -1,3 +0,0 @@
-version https://git-lfs.github.com/spec/v1
-oid sha256:8b35708cfa327cd080308516153e5d97cfca995b75378c47462f9b64a63b21d5
-size 1044
diff --git a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/quad_5500x5500_adaptive_dt.prj b/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/quad_5500x5500_adaptive_dt.prj
deleted file mode 100644
index 2fdcff9483f881962f5e00fa62c72cf32bcfd3b2..0000000000000000000000000000000000000000
--- a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/quad_5500x5500_adaptive_dt.prj
+++ /dev/null
@@ -1,337 +0,0 @@
-<?xml version="1.0" encoding="ISO-8859-1"?>
-<OpenGeoSysProject>
- <mesh>square_5500x5500.vtu</mesh>
- <geometry>quad_5500x5500.gml</geometry>
- <processes>
- <process>
- <name>HeatConduction</name>
- <type>HEAT_CONDUCTION</type>
- <integration_order>2</integration_order>
- <thermal_conductivity>lambda_solid</thermal_conductivity>
- <heat_capacity>c_solid</heat_capacity>
- <density>rho_solid</density>
- <process_variables>
- <process_variable>temperature</process_variable>
- </process_variables>
- <secondary_variables>
- <secondary_variable type="static" internal_name="heat_flux_x" output_name="heat_flux_x"/>
- <secondary_variable type="static" internal_name="heat_flux_y" output_name="heat_flux_y"/>
- </secondary_variables>
- </process>
- <process>
- <name>LiquidFlow</name>
- <type>LIQUID_FLOW</type>
- <integration_order>2</integration_order>
- <darcy_gravity>
- <!-- axis_id: 0, 1, or the dimension of space minus one -->
- <axis_id>1</axis_id>
- <!-- g>=0. g=0: non gravity term -->
- <g>9.81</g>
- </darcy_gravity>
- <process_variables>
- <process_variable>pressure</process_variable>
- </process_variables>
- <secondary_variables>
- <secondary_variable type="static" internal_name="darcy_velocity" output_name="v"/>
- </secondary_variables>
- <material_property>
- <fluid>
- <density>
- <type>TemperatureDependent</type>
- <rho0>1000</rho0>
- <temperature0>20</temperature0>
- <beta>4.3e-4</beta>
- </density>
- <viscosity>
- <type>TemperatureDependent</type>
- <mu0>1.0e-3</mu0>
- <tc>20</tc>
- <tv>75</tv>
- </viscosity>
- <specific_heat_capacity>
- <type>Constant</type>
- <value> 4200. </value>
- </specific_heat_capacity>
- <thermal_conductivity>
- <type>Constant</type>
- <value> .65 </value>
- </thermal_conductivity>
- </fluid>
- <porous_medium>
- <porous_medium id="0">
- <permeability>
- <permeability_tensor_entries>kappa1</permeability_tensor_entries>
- <type>Constant</type>
- </permeability>
- <porosity>
- <type>Constant</type>
- <porosity_parameter>constant_porosity_parameter</porosity_parameter>
- </porosity>
- <storage>
- <type>Constant</type>
- <value>0.0</value>
- </storage>
- </porous_medium>
- </porous_medium>
- <solid>
- <thermal_expansion>solid_thermal_expansion</thermal_expansion>
- <biot_constant>biot_const</biot_constant>
- </solid>
- </material_property>
- </process>
- </processes>
- <time_loop>
- <global_process_coupling>
- <max_iter> 6 </max_iter>
- <convergence_criterion>
- <type>Residual</type>
- <norm_type>NORM2</norm_type>
- <reltol>1.e-6</reltol>
- </convergence_criterion>
- </global_process_coupling>
- <processes>
- <process ref="HeatConduction">
- <nonlinear_solver>basic_picard_T</nonlinear_solver>
- <convergence_criterion>
- <type>DeltaX</type>
- <norm_type>NORM2</norm_type>
- <abstol>1.e-6</abstol>
- </convergence_criterion>
- <time_discretization>
- <type>BackwardEuler</type>
- </time_discretization>
- <time_stepping>
- <type>EvolutionaryPIDcontroller</type>
- <t_initial> 0.0 </t_initial>
- <t_end> 5e10 </t_end>
- <dt_guess> 1 </dt_guess>
- <dt_min> 1 </dt_min>
- <dt_max> 4e8 </dt_max>
- <rel_dt_min> 0.1 </rel_dt_min>
- <rel_dt_max> 2 </rel_dt_max>
- <specific_times>5.0e9 2.e10</specific_times>
- <tol> 1.e-2 </tol>
- </time_stepping>
- <coupled_processes>
- <coupled_process> LiquidFlow </coupled_process>
- </coupled_processes>
- <output>
- <variables>
- <variable> temperature </variable>
- <variable> heat_flux_x </variable>
- <variable> heat_flux_y </variable>
- </variables>
- </output>
- </process>
- <process ref="LiquidFlow">
- <nonlinear_solver>basic_picard_H</nonlinear_solver>
- <convergence_criterion>
- <type>DeltaX</type>
- <norm_type>NORM2</norm_type>
- <abstol>1.e-6</abstol>
- </convergence_criterion>
- <time_discretization>
- <type>BackwardEuler</type>
- </time_discretization>
- <time_stepping>
- <type>EvolutionaryPIDcontroller</type>
- <t_initial> 0.0 </t_initial>
- <t_end> 5e10 </t_end>
- <dt_guess> 1 </dt_guess>
- <dt_min> 1 </dt_min>
- <dt_max> 4e8 </dt_max>
- <rel_dt_min> 0.1 </rel_dt_min>
- <rel_dt_max> 10 </rel_dt_max>
- <tol> 10.0 </tol>
- <specific_times>5.0e9 2.e10</specific_times>
- </time_stepping>
- <coupled_processes>
- <coupled_process> HeatConduction </coupled_process>
- </coupled_processes>
- <output>
- <variables>
- <variable> pressure </variable>
- <variable> v </variable>
- </variables>
- </output>
- </process>
- </processes>
- <output>
- <type>VTK</type>
- <prefix>ThermalConvection</prefix>
- <timesteps>
- <pair>
- <repeat> 1 </repeat>
- <each_steps> 231 </each_steps>
- </pair>
- </timesteps>
- </output>
- </time_loop>
- <parameters>
- <parameter>
- <name>poro</name>
- <type>Constant</type>
- <value>0.001</value>
- </parameter>
- <parameter>
- <name>kappa</name>
- <type>Constant</type>
- <value>1e-14</value>
- </parameter>
- <parameter>
- <name>rho_solid</name>
- <type>Constant</type>
- <value>0.0</value>
- </parameter>
- <parameter>
- <name>rho_fluid</name>
- <type>Constant</type>
- <value>1000</value>
- </parameter>
- <parameter>
- <name>lambda_solid</name>
- <type>Constant</type>
- <value>3.0</value>
- </parameter>
- <parameter>
- <name>lambda_fluid</name>
- <type>Constant</type>
- <value>0.65</value>
- </parameter>
- <parameter>
- <name>c_solid</name>
- <type>Constant</type>
- <value>0</value>
- </parameter>
- <parameter>
- <name>c_fluid</name>
- <type>Constant</type>
- <value>4200</value>
- </parameter>
- <parameter>
- <name>solid_thermal_expansion</name>
- <type>Constant</type>
- <value>0e-5</value>
- </parameter>
- <parameter>
- <name>biot_const</name>
- <type>Constant</type>
- <value>1.0</value>
- </parameter>
- <parameter>
- <name>T0</name>
- <type>MeshNode</type>
- <field_name>initial_temperature</field_name>
- </parameter>
- <parameter>
- <name>P0</name>
- <type>MeshNode</type>
- <field_name>initial_pressure</field_name>
- </parameter>
- <parameter>
- <name>p_Dirichlet_topleft</name>
- <type>Constant</type>
- <value>0</value>
- </parameter>
- <parameter>
- <name>t_Dirichlet_bottom</name>
- <type>Constant</type>
- <value>170</value>
- </parameter>
- <parameter>
- <name>t_Dirichlet_top</name>
- <type>Constant</type>
- <value>20</value>
- </parameter>
- <parameter>
- <name>kappa1</name>
- <type>Constant</type>
- <values>1.0e-14</values>
- </parameter>
- <parameter>
- <name>constant_porosity_parameter</name>
- <type>Constant</type>
- <value>0.001</value>
- </parameter>
- </parameters>
- <process_variables>
- <process_variable>
- <name>temperature</name>
- <components>1</components>
- <order>1</order>
- <initial_condition>T0</initial_condition>
- <boundary_conditions>
- <boundary_condition>
- <geometrical_set>square_1x1_geometry</geometrical_set>
- <geometry>bottom</geometry>
- <type>Dirichlet</type>
- <parameter>t_Dirichlet_bottom</parameter>
- </boundary_condition>
- <boundary_condition>
- <geometrical_set>square_1x1_geometry</geometrical_set>
- <geometry>top</geometry>
- <type>Dirichlet</type>
- <parameter>t_Dirichlet_top</parameter>
- </boundary_condition>
- </boundary_conditions>
- </process_variable>
- <process_variable>
- <name>pressure</name>
- <components>1</components>
- <order>1</order>
- <initial_condition>P0</initial_condition>
- <boundary_conditions>
- <boundary_condition>
- <geometrical_set>square_1x1_geometry</geometrical_set>
- <geometry>topleft</geometry>
- <type>Dirichlet</type>
- <parameter>p_Dirichlet_topleft</parameter>
- </boundary_condition>
- </boundary_conditions>
- </process_variable>
- </process_variables>
- <nonlinear_solvers>
- <nonlinear_solver>
- <name>basic_picard_T</name>
- <type>Picard</type>
- <max_iter>100</max_iter>
- <linear_solver>linear_solver_T</linear_solver>
- </nonlinear_solver>
- <nonlinear_solver>
- <name>basic_picard_H</name>
- <type>Picard</type>
- <max_iter>100</max_iter>
- <linear_solver>linear_solver_H</linear_solver>
- </nonlinear_solver>
- </nonlinear_solvers>
- <linear_solvers>
- <linear_solver>
- <name>linear_solver_T</name>
- <lis>-i bicgstab -p jacobi -tol 1e-10 -maxiter 3000</lis>
- <eigen>
- <solver_type>BiCGSTAB</solver_type>
- <precon_type>DIAGONAL</precon_type>
- <max_iteration_step>10000</max_iteration_step>
- <error_tolerance>1e-14</error_tolerance>
- </eigen>
- <petsc>
- <prefix>T</prefix>
- <parameters>-T_ksp_type bcgs -T_pc_type bjacobi -T_ksp_rtol 1e-14 -T_ksp_max_it 3000</parameters>
- </petsc>
- </linear_solver>
- <linear_solver>
- <name>linear_solver_H</name>
- <lis>-i bicgstab -p jacobi -tol 1e-10 -maxiter 3000</lis>
- <eigen>
- <solver_type>BiCGSTAB</solver_type>
- <precon_type>DIAGONAL</precon_type>
- <max_iteration_step>10000</max_iteration_step>
- <error_tolerance>1e-16</error_tolerance>
- </eigen>
- <petsc>
- <prefix>gw</prefix>
- <parameters>-gw_ksp_type bcgs -gw_pc_type bjacobi -gw_pc_hypre_type boomeramg -gw_ksp_rtol 1e-14 -gw_ksp_max_it 3000</parameters>
- </petsc>
- </linear_solver>
- </linear_solvers>
-</OpenGeoSysProject>
diff --git a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/square_5500x5500.vtu b/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/square_5500x5500.vtu
deleted file mode 100644
index 5fbe991650a85b238bc70e653e61bf2765b2c8a1..0000000000000000000000000000000000000000
--- a/Tests/Data/StaggeredCoupledProcesses/TH/ThermalConvection2D/square_5500x5500.vtu
+++ /dev/null
@@ -1,3 +0,0 @@
-version https://git-lfs.github.com/spec/v1
-oid sha256:d1cf26cf2835f5509e88b6e6156589305bb5bd8038970c92c554f1fa09585f13
-size 115302