Commit 64199c96 authored by wenqing's avatar wenqing
Browse files

[TRM/Test] Added Liakopoulos HM test

parent a2603167
......@@ -4,6 +4,27 @@ if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE ThermoRichardsMechanics/RichardsFlow2D/RichardsFlow_2d_small.prj)
endif()
AddTest(
NAME ThermoRichardsMechanics_liakopoulosHM
PATH ThermoRichardsMechanics/LiakopoulosHM
EXECUTABLE ogs
EXECUTABLE_ARGS liakopoulos.prj
WRAPPER time
TESTER vtkdiff
REQUIREMENTS NOT (OGS_USE_LIS OR OGS_USE_MPI)
RUNTIME 17
DIFF_DATA
liakopoulos_t_300.vtu liakopoulos_t_300.vtu sigma sigma 1e-9 1e-12
liakopoulos_t_300.vtu liakopoulos_t_300.vtu displacement displacement 1e-10 1e-12
liakopoulos_t_300.vtu liakopoulos_t_300.vtu saturation saturation 1e-10 1e-12
liakopoulos_t_600.vtu liakopoulos_t_600.vtu sigma sigma 1e-9 1e-12
liakopoulos_t_600.vtu liakopoulos_t_600.vtu displacement displacement 1e-10 1e-12
liakopoulos_t_600.vtu liakopoulos_t_600.vtu saturation saturation 1e-10 1e-12
liakopoulos_t_7200.vtu liakopoulos_t_7200.vtu sigma sigma 1e-9 1e-12
liakopoulos_t_7200.vtu liakopoulos_t_7200.vtu displacement displacement 1e-10 1e-12
liakopoulos_t_7200.vtu liakopoulos_t_7200.vtu saturation saturation 1e-10 1e-12
)
AddTest(
NAME ThermoRichardsMechanics_3D_ThermoElastic_Stress_Analysis
PATH ThermoRichardsMechanics/Simple3DThermoMechanicsFromTM
......
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<meshes>
<mesh axially_symmetric="true">liakopoulos.vtu</mesh>
<mesh axially_symmetric="true">liakopoulos_left.vtu</mesh>
<mesh axially_symmetric="true">liakopoulos_right.vtu</mesh>
<mesh axially_symmetric="true">liakopoulos_top.vtu</mesh>
<mesh axially_symmetric="true">liakopoulos_bottom.vtu</mesh>
</meshes>
<processes>
<process>
<name>Liakopoulos</name>
<type>THERMO_RICHARDS_MECHANICS</type>
<integration_order>3</integration_order>
<mass_lumping>true</mass_lumping>
<constitutive_relation>
<type>LinearElasticIsotropic</type>
<youngs_modulus>E</youngs_modulus>
<poissons_ratio>nu</poissons_ratio>
</constitutive_relation>
<process_variables>
<temperature>temperature</temperature>
<displacement>displacement</displacement>
<pressure>pressure</pressure>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="sigma" output_name="sigma"/>
<secondary_variable internal_name="epsilon" output_name="epsilon"/>
<secondary_variable internal_name="velocity" output_name="velocity"/>
<secondary_variable internal_name="saturation" output_name="saturation"/>
</secondary_variables>
<specific_body_force>0 -9.81</specific_body_force>
<initial_stress>Initial_stress</initial_stress>
</process>
</processes>
<media>
<medium>
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<name>bulk_modulus</name>
<type>Constant</type>
<value>2e9</value>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1e-3</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>1e3</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>density</name>
<type>Constant</type>
<value>2e3</value>
</property>
<property>
<name>bulk_modulus</name>
<type>Constant</type>
<value>1e12</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>1.0</value>
</property>
<property>
<name>thermal_expansivity</name>
<type>Constant</type>
<value>0.0</value>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0.2975</value>
</property>
<property>
<name>biot_coefficient</name>
<type>Constant</type>
<value>1.0</value>
</property>
<property>
<name>saturation</name>
<type>SaturationLiakopoulos</type>
</property>
<property>
<name>bishops_effective_stress</name>
<type>BishopsPowerLaw</type>
<exponent>1</exponent>
</property>
<property>
<name>permeability</name>
<type>Constant</type>
<value>4.5e-13</value>
</property>
<property>
<name>relative_permeability</name>
<type>Curve</type>
<independent_variable>liquid_saturation</independent_variable>
<curve>k_rel</curve>
</property>
</properties>
</medium>
</media>
<time_loop>
<processes>
<process ref="Liakopoulos">
<nonlinear_solver>basic_newton</nonlinear_solver>
<compensate_non_equilibrium_initial_residuum>false
</compensate_non_equilibrium_initial_residuum>
<convergence_criterion>
<type>PerComponentDeltaX</type>
<norm_type>NORM2</norm_type>
<abstols>1e-5 1e-6 1e-10 1e-10</abstols>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial> 0.0 </t_initial>
<t_end>7200</t_end>
<timesteps>
<pair><repeat>10</repeat>
<delta_t>1</delta_t></pair>
<pair><repeat>9</repeat>
<delta_t>10</delta_t></pair>
<pair><repeat>11</repeat>
<delta_t>100</delta_t></pair>
<pair><repeat>3</repeat>
<delta_t>200</delta_t></pair>
<pair><repeat>3</repeat>
<delta_t>400</delta_t></pair>
<pair><repeat>7</repeat>
<delta_t>600</delta_t></pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>{:meshname}</prefix>
<suffix>_t_{:gtime}</suffix>
<timesteps>
<pair>
<repeat>500</repeat>
<each_steps>10</each_steps>
</pair>
</timesteps>
<variables>
<variable>displacement</variable>
<variable>pressure</variable>
<variable>sigma</variable>
<variable>epsilon</variable>
<variable>velocity</variable>
<variable>saturation</variable>
</variables>
<fixed_output_times>
0.06
60.
120.
300.0
600.0
1200.0
2400.0
4800.0
6000.0
7200.0
</fixed_output_times>
</output>
</time_loop>
<parameters>
<parameter>
<mesh>liakopoulos</mesh>
<name>Initial_stress</name>
<type>Function</type>
<expression>0</expression>
<expression>-((1-0.2975)*2e3+0.2975*1e3)*9.81*(1-y)
</expression>
<expression>0</expression>
<expression>0</expression>
</parameter>
<!-- Mechanics -->
<parameter>
<name>E</name>
<type>Constant</type>
<value>1.3e6</value>
</parameter>
<parameter>
<name>nu</name>
<type>Constant</type>
<value>0.4</value>
</parameter>
<!-- Model parameters -->
<parameter>
<name>displacement_ic</name>
<type>Constant</type>
<values>0 0</values>
</parameter>
<parameter>
<name>pressure_ic</name>
<type>Constant</type>
<value>0.0</value>
</parameter>
<parameter>
<name>pressure_bc</name>
<type>Constant</type>
<value>0.0</value>
</parameter>
<parameter>
<name>dirichlet</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>T_ref</name>
<type>Constant</type>
<value>293.15</value>
</parameter>
</parameters>
<curves>
<curve>
<name>k_rel</name>
<coords>
0.2
0.528
0.536
0.544
0.552
0.56
0.568
0.576
0.584
0.592
0.6
0.608
0.616
0.624
0.632
0.64
0.648
0.656
0.664
0.672
0.68
0.688
0.696
0.704
0.712
0.72
0.728
0.736
0.744
0.752
0.76
0.768
0.776
0.784
0.792
0.8
0.808
0.816
0.824
0.832
0.84
0.848
0.856
0.864
0.872
0.88
0.888
0.896
0.904
0.912
0.92
0.928000000000001
0.936000000000001
0.944000000000001
0.952000000000001
0.960000000000001
0.968000000000001
0.976000000000001
0.984000000000001
0.992000000000001
1
</coords>
<values>
0.0
0.0
0.0
0.003125124673949
0.020823885169721
0.038518821680155
0.056209865519575
0.073896945473639
0.091579987656767
0.109258915358633
0.126933648878653
0.144604105347267
0.162270198532675
0.179931838631514
0.197588932041764
0.215241381115967
0.232889083892571
0.250531933802915
0.268169819351031
0.285802623763035
0.303430224602408
0.3210524933469
0.338669294922165
0.356280487186422
0.373885920359564
0.391485436388984
0.409078868243118
0.426666039122081
0.444246761572863
0.461820836494194
0.479388052013308
0.496948182213303
0.514500985685365
0.532046203874627
0.549583559181513
0.56711275277159
0.584633462035706
0.602145337627643
0.619647999987513
0.637141035234136
0.654623990276238
0.672096366947284
0.689557614907064
0.707007122967474
0.724444208378889
0.741868103439419
0.759277938533726
0.776672720324285
0.79405130322659
0.811412351361665
0.828754286640385
0.846075216009299
0.863372826201009
0.880644225491632
0.897885694137663
0.915092266056643
0.932256968686883
0.949369276790233
0.966411379494378
0.983345955951454
1
</values>
</curve>
</curves>
<process_variables>
<process_variable>
<name>temperature</name>
<components>1</components>
<order>1</order>
<initial_condition>T_ref</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>liakopoulos_bottom</mesh>
<type>Dirichlet</type>
<parameter>T_ref</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>
<mesh>liakopoulos_bottom</mesh>
<type>Dirichlet</type>
<parameter>pressure_bc</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>displacement</name>
<components>2</components>
<order>2</order>
<initial_condition>displacement_ic</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>liakopoulos_left</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>dirichlet</parameter>
</boundary_condition>
<boundary_condition>
<mesh>liakopoulos_right</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>dirichlet</parameter>
</boundary_condition>
<boundary_condition>
<mesh>liakopoulos_bottom</mesh>
<type>Dirichlet</type>
<component>1</component>
<parameter>dirichlet</parameter>
</boundary_condition>
</boundary_conditions>
</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>true</scaling>
</eigen>
</linear_solver>
</linear_solvers>
</OpenGeoSysProject>
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