Skip to content
Snippets Groups Projects
Commit 282f1d21 authored by wenqing's avatar wenqing Committed by Lars Bilke
Browse files

[Test/TRM/PETSc] Run the point heat source benchmark with parallel computing

parent 8e016ce8
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
Showing
with 695 additions and 0 deletions
ProcessLib/ThermoRichardsMechanics/Tests.cmake
+ 60
− 0
View file @ 282f1d21
......@@ -117,3 +117,63 @@ AddTest(
flow_fully_saturated_ts_2_t_2_000000_0.vtu flow_fully_saturated_ts_2_t_2_000000_0.vtu epsilon epsilon 1e-10 1e-9
flow_fully_saturated_ts_2_t_2_000000_1.vtu flow_fully_saturated_ts_2_t_2_000000_1.vtu epsilon epsilon 1e-10 1e-9
)
AddTest(
NAME ParallelFEM_ThermoRichardsMechanics_point_heat_injection
PATH ThermoRichardsMechanics/PointHeatSource
RUNTIME 45
EXECUTABLE ogs
EXECUTABLE_ARGS point_heat_source_2D.prj
WRAPPER mpirun
WRAPPER_ARGS -np 3
TESTER vtkdiff
REQUIREMENTS OGS_USE_MPI
DIFF_DATA
PointHeatSource_ts_10_t_50000_000000_0.vtu PointHeatSource_ts_10_t_50000_000000_0.vtu displacement displacement 1e-10 1.0e-9
PointHeatSource_ts_10_t_50000_000000_0.vtu PointHeatSource_ts_10_t_50000_000000_0.vtu pressure pressure 1e-10 1.0e-6
PointHeatSource_ts_10_t_50000_000000_0.vtu PointHeatSource_ts_10_t_50000_000000_0.vtu temperature temperature 1e-10 1.0e-9
PointHeatSource_ts_10_t_50000_000000_0.vtu PointHeatSource_ts_10_t_50000_000000_0.vtu epsilon epsilon 1e-10 1.0e-9
PointHeatSource_ts_10_t_50000_000000_0.vtu PointHeatSource_ts_10_t_50000_000000_0.vtu sigma sigma 1e-10 1.0e-6
#
PointHeatSource_ts_10_t_50000_000000_1.vtu PointHeatSource_ts_10_t_50000_000000_1.vtu displacement displacement 1e-10 1.0e-9
PointHeatSource_ts_10_t_50000_000000_1.vtu PointHeatSource_ts_10_t_50000_000000_1.vtu pressure pressure 1e-10 1.0e-6
PointHeatSource_ts_10_t_50000_000000_1.vtu PointHeatSource_ts_10_t_50000_000000_1.vtu temperature temperature 1e-10 1.0e-9
PointHeatSource_ts_10_t_50000_000000_1.vtu PointHeatSource_ts_10_t_50000_000000_1.vtu epsilon epsilon 1e-10 1.0e-9
PointHeatSource_ts_10_t_50000_000000_1.vtu PointHeatSource_ts_10_t_50000_000000_1.vtu sigma sigma 1e-10 1.0e-6
#
PointHeatSource_ts_10_t_50000_000000_2.vtu PointHeatSource_ts_10_t_50000_000000_2.vtu displacement displacement 1e-10 1.0e-9
PointHeatSource_ts_10_t_50000_000000_2.vtu PointHeatSource_ts_10_t_50000_000000_2.vtu pressure pressure 1e-10 1.0e-6
PointHeatSource_ts_10_t_50000_000000_2.vtu PointHeatSource_ts_10_t_50000_000000_2.vtu temperature temperature 1e-10 1.0e-9
PointHeatSource_ts_10_t_50000_000000_2.vtu PointHeatSource_ts_10_t_50000_000000_2.vtu epsilon epsilon 1e-10 1.0e-9
PointHeatSource_ts_10_t_50000_000000_2.vtu PointHeatSource_ts_10_t_50000_000000_2.vtu sigma sigma 1e-10 1.0e-6
)
AddTest(
NAME ParallelFEM_ThermoRichardsMechanics_point_heat_injection_gml
PATH ThermoRichardsMechanics/PointHeatSource
RUNTIME 45
EXECUTABLE ogs
EXECUTABLE_ARGS point_heat_source_2D_gml.prj
WRAPPER mpirun
WRAPPER_ARGS -np 3
TESTER vtkdiff
REQUIREMENTS OGS_USE_MPI
DIFF_DATA
PointHeatSource_gml_ts_10_t_50000_000000_0.vtu PointHeatSource_gml_ts_10_t_50000_000000_0.vtu displacement displacement 1e-10 1.0e-9
PointHeatSource_gml_ts_10_t_50000_000000_0.vtu PointHeatSource_gml_ts_10_t_50000_000000_0.vtu pressure pressure 1e-10 1.0e-6
PointHeatSource_gml_ts_10_t_50000_000000_0.vtu PointHeatSource_gml_ts_10_t_50000_000000_0.vtu temperature temperature 1e-10 1.0e-9
PointHeatSource_gml_ts_10_t_50000_000000_0.vtu PointHeatSource_gml_ts_10_t_50000_000000_0.vtu epsilon epsilon 1e-10 1.0e-9
PointHeatSource_gml_ts_10_t_50000_000000_0.vtu PointHeatSource_gml_ts_10_t_50000_000000_0.vtu sigma sigma 1e-10 1.0e-6
#
PointHeatSource_gml_ts_10_t_50000_000000_1.vtu PointHeatSource_gml_ts_10_t_50000_000000_1.vtu displacement displacement 1e-10 1.0e-9
PointHeatSource_gml_ts_10_t_50000_000000_1.vtu PointHeatSource_gml_ts_10_t_50000_000000_1.vtu pressure pressure 1e-10 1.0e-6
PointHeatSource_gml_ts_10_t_50000_000000_1.vtu PointHeatSource_gml_ts_10_t_50000_000000_1.vtu temperature temperature 1e-10 1.0e-9
PointHeatSource_gml_ts_10_t_50000_000000_1.vtu PointHeatSource_gml_ts_10_t_50000_000000_1.vtu epsilon epsilon 1e-10 1.0e-9
PointHeatSource_gml_ts_10_t_50000_000000_1.vtu PointHeatSource_gml_ts_10_t_50000_000000_1.vtu sigma sigma 1e-10 1.0e-6
#
PointHeatSource_gml_ts_10_t_50000_000000_2.vtu PointHeatSource_gml_ts_10_t_50000_000000_2.vtu displacement displacement 1e-10 1.0e-9
PointHeatSource_gml_ts_10_t_50000_000000_2.vtu PointHeatSource_gml_ts_10_t_50000_000000_2.vtu pressure pressure 1e-10 1.0e-6
PointHeatSource_gml_ts_10_t_50000_000000_2.vtu PointHeatSource_gml_ts_10_t_50000_000000_2.vtu temperature temperature 1e-10 1.0e-9
PointHeatSource_gml_ts_10_t_50000_000000_2.vtu PointHeatSource_gml_ts_10_t_50000_000000_2.vtu epsilon epsilon 1e-10 1.0e-9
PointHeatSource_gml_ts_10_t_50000_000000_2.vtu PointHeatSource_gml_ts_10_t_50000_000000_2.vtu sigma sigma 1e-10 1.0e-6
)
Tests/Data/ThermoRichardsMechanics/PointHeatSource/point_heat_source_2D.prj
+ 7
− 0
View file @ 282f1d21
......@@ -306,6 +306,13 @@
<solver_type>SparseLU</solver_type>
<scaling>1</scaling>
</eigen>
<petsc>
<parameters>-ksp_type tfqmr
-pc_type jacobi
-ksp_rtol 1.e-10 -ksp_atol 1.e-7
-ksp_max_it 4000
</parameters>
</petsc>
</linear_solver>
</linear_solvers>
</OpenGeoSysProject>
Tests/Data/ThermoRichardsMechanics/PointHeatSource/point_heat_source_2D_gml.prj 0 → 100644
+ 323
− 0
View file @ 282f1d21
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<mesh axially_symmetric="true">quarter_002_2nd.vtu</mesh>
<geometry>quarter_circle_geometry.gml</geometry>
<search_length_algorithm>
<type>fixed</type>
<value>0.001</value>
</search_length_algorithm>
<processes>
<process>
<name>PointHeatSource</name>
<type>THERMO_RICHARDS_MECHANICS</type>
<integration_order>4</integration_order>
<constitutive_relation>
<type>LinearElasticIsotropic</type>
<youngs_modulus>E</youngs_modulus>
<poissons_ratio>nu</poissons_ratio>
</constitutive_relation>
<process_variables>
<temperature>temperature</temperature>
<pressure>pressure</pressure>
<displacement>displacement</displacement>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="sigma" output_name="sigma"/>
<secondary_variable internal_name="epsilon" output_name="epsilon"/>
</secondary_variables>
<specific_body_force>0 0</specific_body_force>
</process>
</processes>
<media>
<medium id="0">
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>4280.0</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>0.6</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>999.1</value>
<!--type>Linear</type>
<reference_value>999.1</reference_value>
<independent_variable>
<variable_name>temperature</variable_name>
<reference_condition>273.15</reference_condition>
<slope>-4e-4</slope>
</independent_variable-->
</property>
<property>
<name>thermal_expansivity</name>
<type>Constant</type>
<value>4.0e-4</value>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1.e-3</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>density</name>
<type>Constant</type>
<value>2290</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>1.838</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>917.654</value>
</property>
<property>
<name>thermal_expansivity</name>
<type>Constant</type>
<value>1.5e-5</value>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>saturation</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>relative_permeability</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>permeability</name>
<type>Constant</type>
<value>2e-20</value>
</property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0.16</value>
</property>
<property>
<name>bishops_effective_stress</name>
<type>BishopsSaturationCutoff</type>
<cutoff_value>1.0</cutoff_value>
</property>
<property>
<name>biot_coefficient</name>
<type>Constant</type>
<value>1.0</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>EffectiveThermalConductivityPorosityMixing</type>
</property>
</properties>
</medium>
</media>
<time_loop>
<processes>
<process ref="PointHeatSource">
<nonlinear_solver>basic_newton</nonlinear_solver>
<convergence_criterion>
<type>PerComponentDeltaX</type>
<norm_type>NORM2</norm_type>
<reltols>1e-14 1e-12 1e-13 1e-13</reltols>
<abstols>1e-6 1e-5 1e-13 1e-13</abstols>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0</t_initial>
<t_end>50000</t_end>
<timesteps>
<pair>
<repeat>10</repeat>
<delta_t>5000</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>PointHeatSource_gml</prefix>
<timesteps>
<pair>
<repeat>1</repeat>
<each_steps>10</each_steps>
</pair>
</timesteps>
<variables>
<variable>displacement</variable>
<variable>pressure</variable>
<variable>temperature</variable>
<variable>sigma</variable>
<variable>epsilon</variable>
</variables>
<suffix>_ts_{:timestep}_t_{:time}</suffix>
</output>
</time_loop>
<parameters>
<parameter>
<name>E</name>
<type>Constant</type>
<value>5000000000</value>
</parameter>
<parameter>
<name>nu</name>
<type>Constant</type>
<value>0.3</value>
</parameter>
<parameter>
<name>T0</name>
<type>Constant</type>
<value>273.15</value>
</parameter>
<parameter>
<name>displacement0</name>
<type>Constant</type>
<values>0 0</values>
</parameter>
<parameter>
<name>pressure_ic</name>
<type>Constant</type>
<values>0</values>
</parameter>
<parameter>
<name>dirichlet0</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>Neumann0</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>temperature_ic</name>
<type>Constant</type>
<value>273.15</value>
</parameter>
<parameter>
<name>pressure_bc_left</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>temperature_bc_left</name>
<type>Constant</type>
<value>273.15</value>
</parameter>
<parameter>
<name>temperature_source_term</name>
<type>Constant</type>
<value>150</value>
</parameter>
</parameters>
<process_variables>
<process_variable>
<name>displacement</name>
<components>2</components>
<order>2</order>
<initial_condition>displacement0</initial_condition>
<boundary_conditions>
<boundary_condition>
<geometrical_set>quarter_circle_geometry</geometrical_set>
<geometry>bottom</geometry>
<type>Dirichlet</type>
<component>1</component>
<parameter>dirichlet0</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>quarter_circle_geometry</geometrical_set>
<geometry>left</geometry>
<type>Dirichlet</type>
<component>0</component>
<parameter>dirichlet0</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>pressure</name>
<components>1</components>
<order>1</order>
<initial_condition>pressure_ic</initial_condition>
<boundary_conditions>
<boundary_condition>
<geometrical_set>quarter_circle_geometry</geometrical_set>
<geometry>out</geometry>
<type>Dirichlet</type>
<component>0</component>
<parameter>pressure_bc_left</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>temperature</name>
<components>1</components>
<order>1</order>
<initial_condition>temperature_ic</initial_condition>
<boundary_conditions>
<boundary_condition>
<geometrical_set>quarter_circle_geometry</geometrical_set>
<geometry>out</geometry>
<type>Dirichlet</type>
<component>0</component>
<parameter>temperature_bc_left</parameter>
</boundary_condition>
</boundary_conditions>
<source_terms>
<source_term>
<geometrical_set>quarter_circle_geometry</geometrical_set>
<geometry>center</geometry>
<type>Nodal</type>
<parameter>temperature_source_term</parameter>
</source_term>
</source_terms>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_newton</name>
<type>Newton</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>
<eigen>
<solver_type>SparseLU</solver_type>
<scaling>1</scaling>
</eigen>
<petsc>
<parameters>-ksp_type tfqmr
-pc_type jacobi
-ksp_rtol 1.e-10 -ksp_atol 1.e-7
-ksp_max_it 4000
</parameters>
</petsc>
</linear_solver>
</linear_solvers>
</OpenGeoSysProject>
Tests/Data/ThermoRichardsMechanics/PointHeatSource/quarter_002.vtu 0 → 100644
+ 23
− 0
View file @ 282f1d21
This diff is collapsed.
Tests/Data/ThermoRichardsMechanics/PointHeatSource/quarter_002_2nd_partitioned_msh_cfg3.bin 0 → 100644
+ 0
− 0
View file @ 282f1d21
File added
Tests/Data/ThermoRichardsMechanics/PointHeatSource/quarter_002_2nd_partitioned_msh_ele3.bin 0 → 100644
+ 0
− 0
View file @ 282f1d21
File added
Tests/Data/ThermoRichardsMechanics/PointHeatSource/quarter_002_2nd_partitioned_msh_nod3.bin 0 → 100644
+ 0
− 0
View file @ 282f1d21
File added
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