From 0e0511bf212b13beb0bf5e234cf578b5de139104 Mon Sep 17 00:00:00 2001
From: tnagel <thomas.nagel@ufz.de>
Date: Mon, 2 Dec 2019 19:15:16 +0100
Subject: [PATCH] Added two test cohesive zone tests for new formulation.
---
ProcessLib/LIE/HydroMechanics/Tests.cmake | 46 +++
...nts_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu | 3 +
...w_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu | 3 +
...single_fracture_3compartments_flow_CHZ.prj | 347 ++++++++++++++++++
...fracture_3compartments_flow_CHZ_sigma0.prj | 347 ++++++++++++++++++
5 files changed, 746 insertions(+)
create mode 100644 Tests/Data/LIE/HydroMechanics/expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu
create mode 100644 Tests/Data/LIE/HydroMechanics/expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu
create mode 100644 Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ.prj
create mode 100644 Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ_sigma0.prj
diff --git a/ProcessLib/LIE/HydroMechanics/Tests.cmake b/ProcessLib/LIE/HydroMechanics/Tests.cmake
index be148f0e964..7f6a01a720d 100644
--- a/ProcessLib/LIE/HydroMechanics/Tests.cmake
+++ b/ProcessLib/LIE/HydroMechanics/Tests.cmake
@@ -174,3 +174,49 @@ AddTest(
expected_single_fracture_3compartments_flow_linear_aperture0_e_pcs_0_ts_1108_t_10000.000000.vtu single_fracture_3compartments_flow_linear_aperture0_e_pcs_0_ts_1108_t_10000.000000.vtu stress_yy stress_yy 1e-12 1e-15
expected_single_fracture_3compartments_flow_linear_aperture0_e_pcs_0_ts_1108_t_10000.000000.vtu single_fracture_3compartments_flow_linear_aperture0_e_pcs_0_ts_1108_t_10000.000000.vtu velocity velocity 1e-15 1e-15
)
+
+AddTest(
+ NAME LIE_HM_single_fracture_3compartments_flow_CHZ
+ PATH LIE/HydroMechanics
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS single_fracture_3compartments_flow_CHZ.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT (OGS_USE_LIS OR OGS_USE_MPI)
+ DIFF_DATA
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu pressure pressure 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu pressure_interpolated pressure_interpolated 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu displacement displacement 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu displacement_jump1 displacement_jump1 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu nodal_w nodal_w 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu nodal_aperture nodal_aperture 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu strain_xx strain_xx 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu strain_yy strain_yy 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu strain_xy strain_xy 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu stress_xx stress_xx 1e-10 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu stress_yy stress_yy 2e-10 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu velocity velocity 1e-15 1e-15
+)
+
+AddTest(
+ NAME LIE_HM_single_fracture_3compartments_flow_CHZ_sigma0
+ PATH LIE/HydroMechanics
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS single_fracture_3compartments_flow_CHZ_sigma0.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT (OGS_USE_LIS OR OGS_USE_MPI)
+ DIFF_DATA
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu pressure pressure 1e-10 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu pressure_interpolated pressure_interpolated 1e-10 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu displacement displacement 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu displacement_jump1 displacement_jump1 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu nodal_w nodal_w 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu nodal_aperture nodal_aperture 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu strain_xx strain_xx 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu strain_yy strain_yy 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu strain_xy strain_xy 1e-15 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu stress_xx stress_xx 3e-10 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu stress_yy stress_yy 1e-12 1e-15
+ expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu velocity velocity 1e-15 1e-15
+)
diff --git a/Tests/Data/LIE/HydroMechanics/expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu b/Tests/Data/LIE/HydroMechanics/expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu
new file mode 100644
index 00000000000..e505feb3200
--- /dev/null
+++ b/Tests/Data/LIE/HydroMechanics/expected_single_fracture_3compartments_flow_CHZ_pcs_0_ts_356_t_500.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:1a9fde7817d0014ccc63f17891d73de482267f1d40d65c0fd770635b194c5328
+size 19706
diff --git a/Tests/Data/LIE/HydroMechanics/expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu b/Tests/Data/LIE/HydroMechanics/expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu
new file mode 100644
index 00000000000..3e77d0cc3a5
--- /dev/null
+++ b/Tests/Data/LIE/HydroMechanics/expected_single_fracture_3compartments_flow_CHZ_sigma0_pcs_0_ts_356_t_500.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:88d67797d2e7619998c2b19fc11bf078326380671f5320005ebb2c62b2c10b7d
+size 19586
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ.prj
new file mode 100644
index 00000000000..563a05168ea
--- /dev/null
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ.prj
@@ -0,0 +1,347 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>single_fracture_2x2.vtu</mesh>
+ <geometry>single_fracture_2x2.gml</geometry>
+ <processes>
+ <process>
+ <name>HM</name>
+ <type>HYDRO_MECHANICS_WITH_LIE</type>
+ <integration_order>3</integration_order>
+ <dimension>2</dimension>
+ <process_variables>
+ <process_variable>pressure</process_variable>
+ <process_variable>displacement</process_variable>
+ <process_variable>displacement_jump1</process_variable>
+ </process_variables>
+ <constitutive_relation>
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <intrinsic_permeability>k</intrinsic_permeability>
+ <specific_storage>S</specific_storage>
+ <fluid_viscosity>mu</fluid_viscosity>
+ <biot_coefficient>biot_m</biot_coefficient>
+ <porosity>phi</porosity>
+ <solid_density>rho_sr</solid_density>
+ <fluid_density>rho_fr</fluid_density>
+ <specific_body_force>0 0</specific_body_force>
+ <initial_effective_stress>effective_stress0</initial_effective_stress>
+ <!--fracture_model>
+ <type>LinearElasticIsotropic</type>
+ <normal_stiffness>Kn</normal_stiffness>
+ <shear_stiffness>Ks</shear_stiffness>
+ <penalty_aperture_cutoff>0</penalty_aperture_cutoff>
+ <tension_cutoff>1</tension_cutoff>
+ </fracture_model-->
+ <fracture_model>
+ <type>CohesiveZoneModeI</type>
+ <normal_stiffness>Kn</normal_stiffness>
+ <shear_stiffness>Ks</shear_stiffness>
+ <fracture_toughness>Gc</fracture_toughness>
+ <peak_normal_traction>t_np</peak_normal_traction>
+ <penalty_aperture_cutoff>1e-8</penalty_aperture_cutoff>
+ <tension_cutoff>1</tension_cutoff>
+ </fracture_model>
+ <fracture_properties>
+ <material_id>0</material_id>
+ <initial_aperture>aperture0</initial_aperture>
+ <specific_storage>S_f</specific_storage>
+ <biot_coefficient>biot_f</biot_coefficient>
+ <permeability_model>
+ <type>CubicLaw</type>
+ </permeability_model>
+ </fracture_properties>
+ <initial_fracture_effective_stress>fracture_effective_stress0</initial_fracture_effective_stress>
+ <deactivate_matrix_in_flow>true</deactivate_matrix_in_flow>
+ <secondary_variables>
+ </secondary_variables>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="HM">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentResidual</type>
+ <norm_type>NORM2</norm_type>
+ <abstols>1e-10 1e-7 1e-5 1e-6 1e-6</abstols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>500</t_end>
+ <timesteps>
+ <pair>
+ <repeat>300</repeat>
+ <delta_t>0.1</delta_t>
+ </pair>
+ <pair>
+ <repeat>10</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>10</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>single_fracture_3compartments_flow_CHZ</prefix>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <each_steps>1000</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>pressure</variable>
+ <variable>displacement</variable>
+ <variable>displacement_jump1</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <!-- Fluid phase properties -->
+ <parameter>
+ <name>rho_fr</name>
+ <type>Constant</type>
+ <value>1e3</value>
+ </parameter>
+ <parameter>
+ <name>mu</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </parameter>
+ <!-- Solid phase properties -->
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <!-- Porous medium properties -->
+ <parameter>
+ <name>phi</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </parameter>
+ <parameter>
+ <name>k</name>
+ <type>Constant</type>
+ <value>1e-17</value>
+ </parameter>
+ <parameter>
+ <name>S</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>biot_m</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1e9</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.3</value>
+ </parameter>
+ <!-- Fracture properties -->
+ <parameter>
+ <name>aperture0</name>
+ <type>Constant</type>
+ <value>1e-5</value>
+ </parameter>
+ <parameter>
+ <name>S_f</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>biot_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <!--fracture model-->
+ <parameter>
+ <name>Kn</name>
+ <type>Constant</type>
+ <value>1e9</value>
+ </parameter>
+ <parameter>
+ <name>Ks</name>
+ <type>Constant</type>
+ <value>1e8</value>
+ </parameter>
+ <parameter>
+ <name>Gc</name>
+ <type>Constant</type>
+ <value>4</value>
+ </parameter>
+ <parameter>
+ <name>t_np</name>
+ <type>Constant</type>
+ <value>0.5e5</value>
+ </parameter>
+ <!-- Flow properties -->
+ <parameter>
+ <name>p0</name>
+ <type>Constant</type>
+ <values>0</values>
+ </parameter>
+ <parameter>
+ <name>p_in</name>
+ <type>Constant</type>
+ <values>1E5</values>
+ </parameter>
+ <parameter>
+ <name>p_out</name>
+ <type>Constant</type>
+ <values>0</values>
+ </parameter>
+ <!-- Mechanics properties -->
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0</values>
+ </parameter>
+ <parameter>
+ <name>effective_stress0</name>
+ <type>Constant</type>
+ <values>-0e5 -0e5 -0e5 0</values>
+ </parameter>
+ <parameter>
+ <name>fracture_effective_stress0</name>
+ <type>Constant</type>
+ <values>0 -0e5</values>
+ </parameter>
+ <parameter>
+ <name>zero_u</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>load_top</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <parameter>
+ <name>load_bottom</name>
+ <type>Constant</type>
+ <value>1e5</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>p0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>InjectionPoint</geometry>
+ <type>Dirichlet</type>
+ <parameter>p_in</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>OutletPoint</geometry>
+ <type>Dirichlet</type>
+ <parameter>p_out</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>displacement</name>
+ <components>2</components>
+ <order>2</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>PLY_SOUTH</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>PLY_WEST</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>PLY_EAST</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>PLY_NORTH</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>displacement_jump1</name>
+ <components>2</components>
+ <order>2</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>InjectionPoint</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>OutletPoint</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_newton</name>
+ <type>Newton</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>
+ <!--This is working but extremely slow
+ <lis>-i bicgstab -p ssor -ssor_w 1.5 -tol 1e-12 -maxiter 100000</lis>
+ -->
+ <eigen>
+ <!--
+ <solver_type>PardisoLU</solver_type>
+ -->
+ <solver_type>SparseLU</solver_type>
+ <scaling>1</scaling>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ_sigma0.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ_sigma0.prj
new file mode 100644
index 00000000000..66bd2d9f972
--- /dev/null
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ_sigma0.prj
@@ -0,0 +1,347 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>single_fracture_2x2.vtu</mesh>
+ <geometry>single_fracture_2x2.gml</geometry>
+ <processes>
+ <process>
+ <name>HM</name>
+ <type>HYDRO_MECHANICS_WITH_LIE</type>
+ <integration_order>3</integration_order>
+ <dimension>2</dimension>
+ <process_variables>
+ <process_variable>pressure</process_variable>
+ <process_variable>displacement</process_variable>
+ <process_variable>displacement_jump1</process_variable>
+ </process_variables>
+ <constitutive_relation>
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <intrinsic_permeability>k</intrinsic_permeability>
+ <specific_storage>S</specific_storage>
+ <fluid_viscosity>mu</fluid_viscosity>
+ <biot_coefficient>biot_m</biot_coefficient>
+ <porosity>phi</porosity>
+ <solid_density>rho_sr</solid_density>
+ <fluid_density>rho_fr</fluid_density>
+ <specific_body_force>0 0</specific_body_force>
+ <initial_effective_stress>effective_stress0</initial_effective_stress>
+ <!--fracture_model>
+ <type>LinearElasticIsotropic</type>
+ <normal_stiffness>Kn</normal_stiffness>
+ <shear_stiffness>Ks</shear_stiffness>
+ <penalty_aperture_cutoff>0</penalty_aperture_cutoff>
+ <tension_cutoff>1</tension_cutoff>
+ </fracture_model-->
+ <fracture_model>
+ <type>CohesiveZoneModeI</type>
+ <normal_stiffness>Kn</normal_stiffness>
+ <shear_stiffness>Ks</shear_stiffness>
+ <fracture_toughness>Gc</fracture_toughness>
+ <peak_normal_traction>t_np</peak_normal_traction>
+ <penalty_aperture_cutoff>1e-8</penalty_aperture_cutoff>
+ <tension_cutoff>1</tension_cutoff>
+ </fracture_model>
+ <fracture_properties>
+ <material_id>0</material_id>
+ <initial_aperture>aperture0</initial_aperture>
+ <specific_storage>S_f</specific_storage>
+ <biot_coefficient>biot_f</biot_coefficient>
+ <permeability_model>
+ <type>CubicLaw</type>
+ </permeability_model>
+ </fracture_properties>
+ <initial_fracture_effective_stress>fracture_effective_stress0</initial_fracture_effective_stress>
+ <deactivate_matrix_in_flow>true</deactivate_matrix_in_flow>
+ <secondary_variables>
+ </secondary_variables>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="HM">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentResidual</type>
+ <norm_type>NORM2</norm_type>
+ <abstols>1e-10 1e-7 1e-5 1e-6 1e-6</abstols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>500</t_end>
+ <timesteps>
+ <pair>
+ <repeat>300</repeat>
+ <delta_t>0.1</delta_t>
+ </pair>
+ <pair>
+ <repeat>10</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>10</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>single_fracture_3compartments_flow_CHZ_sigma0</prefix>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <each_steps>1000</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>pressure</variable>
+ <variable>displacement</variable>
+ <variable>displacement_jump1</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+ <!-- Fluid phase properties -->
+ <parameter>
+ <name>rho_fr</name>
+ <type>Constant</type>
+ <value>1e3</value>
+ </parameter>
+ <parameter>
+ <name>mu</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </parameter>
+ <!-- Solid phase properties -->
+ <parameter>
+ <name>rho_sr</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <!-- Porous medium properties -->
+ <parameter>
+ <name>phi</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </parameter>
+ <parameter>
+ <name>k</name>
+ <type>Constant</type>
+ <value>1e-17</value>
+ </parameter>
+ <parameter>
+ <name>S</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>biot_m</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1e9</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.3</value>
+ </parameter>
+ <!-- Fracture properties -->
+ <parameter>
+ <name>aperture0</name>
+ <type>Constant</type>
+ <value>1e-5</value>
+ </parameter>
+ <parameter>
+ <name>S_f</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>biot_f</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <!--fracture model-->
+ <parameter>
+ <name>Kn</name>
+ <type>Constant</type>
+ <value>1e9</value>
+ </parameter>
+ <parameter>
+ <name>Ks</name>
+ <type>Constant</type>
+ <value>1e8</value>
+ </parameter>
+ <parameter>
+ <name>Gc</name>
+ <type>Constant</type>
+ <value>4</value>
+ </parameter>
+ <parameter>
+ <name>t_np</name>
+ <type>Constant</type>
+ <value>0.5e5</value>
+ </parameter>
+ <!-- Flow properties -->
+ <parameter>
+ <name>p0</name>
+ <type>Constant</type>
+ <values>0</values>
+ </parameter>
+ <parameter>
+ <name>p_in</name>
+ <type>Constant</type>
+ <values>1E5</values>
+ </parameter>
+ <parameter>
+ <name>p_out</name>
+ <type>Constant</type>
+ <values>0</values>
+ </parameter>
+ <!-- Mechanics properties -->
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0</values>
+ </parameter>
+ <parameter>
+ <name>effective_stress0</name>
+ <type>Constant</type>
+ <values>-2e5 -2e5 -2e5 0</values>
+ </parameter>
+ <parameter>
+ <name>fracture_effective_stress0</name>
+ <type>Constant</type>
+ <values>0 -2e5</values>
+ </parameter>
+ <parameter>
+ <name>zero_u</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>load_top</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <parameter>
+ <name>load_bottom</name>
+ <type>Constant</type>
+ <value>1e5</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>p0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>InjectionPoint</geometry>
+ <type>Dirichlet</type>
+ <parameter>p_in</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>OutletPoint</geometry>
+ <type>Dirichlet</type>
+ <parameter>p_out</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>displacement</name>
+ <components>2</components>
+ <order>2</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>PLY_SOUTH</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>PLY_WEST</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>PLY_EAST</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>PLY_NORTH</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>displacement_jump1</name>
+ <components>2</components>
+ <order>2</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>InjectionPoint</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>single_fracture</geometrical_set>
+ <geometry>OutletPoint</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero_u</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_newton</name>
+ <type>Newton</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>
+ <!--This is working but extremely slow
+ <lis>-i bicgstab -p ssor -ssor_w 1.5 -tol 1e-12 -maxiter 100000</lis>
+ -->
+ <eigen>
+ <!--
+ <solver_type>PardisoLU</solver_type>
+ -->
+ <solver_type>SparseLU</solver_type>
+ <scaling>1</scaling>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
--
GitLab