Commit 1a064ecb authored by Dmitry Yu. Naumov's avatar Dmitry Yu. Naumov
Browse files

[T/RM] Restart for MFront internal variables.

parent 6a2fbcfc
......@@ -23,6 +23,11 @@ if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE RichardsMechanics/LiakopoulosHM/liakopoulos_QN.prj RUNTIME 50)
endif()
if (NOT OGS_USE_MPI AND OGS_USE_MFRONT)
OgsTest(PROJECTFILE RichardsMechanics/mfront_restart_part1.prj RUNTIME 1)
OgsTest(PROJECTFILE RichardsMechanics/mfront_restart_part2.prj RUNTIME 1)
endif()
AddTest(
NAME RichardsMechanics_square_1e2_confined_compression_restart
PATH RichardsMechanics
......
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<meshes>
<mesh axially_symmetric="true">square_1x1_quad_1e0.vtu</mesh>
<mesh axially_symmetric="true">square_1x1_quad_1e0_left.vtu</mesh>
<mesh axially_symmetric="true">square_1x1_quad_1e0_right.vtu</mesh>
<mesh axially_symmetric="true">square_1x1_quad_1e0_top.vtu</mesh>
<mesh axially_symmetric="true">square_1x1_quad_1e0_bottom.vtu</mesh>
</meshes>
<processes>
<process>
<name>RM</name>
<type>RICHARDS_MECHANICS</type>
<integration_order>2</integration_order>
<dimension>2</dimension>
<jacobian_assembler>
<type>Analytical</type>
</jacobian_assembler>
<constitutive_relation>
<type>MFront</type>
<behaviour>ModCamClay_semiExplParaInit</behaviour>
<material_properties>
<material_property name="YoungModulus" parameter="YoungModulus"/>
<material_property name="PoissonRatio" parameter="PoissonRatio"/>
<material_property name="CriticalStateLineSlope" parameter="CriticalStateLineSlope"/>
<material_property name="SwellingLineSlope" parameter="SwellingLineSlope"/>
<material_property name="VirginConsolidationLineSlope" parameter="VirginConsolidationLineSlope"/>
<material_property name="InitialPreConsolidationPressure" parameter="InitialPreConsolidationPressure"/>
<material_property name="InitialPorosity" parameter="InitialPorosity"/>
</material_properties>
</constitutive_relation>
<process_variables>
<pressure>pressure</pressure>
<displacement>displacement</displacement>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="ElasticStrain" output_name="ElasticStrain"/>
<secondary_variable internal_name="EquivalentPlasticStrain" output_name="EquivalentPlasticStrain"/>
<secondary_variable internal_name="sigma" output_name="sigma"/>
<secondary_variable internal_name="swelling_stress" output_name="swelling_stress"/>
<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_variable internal_name="porosity" output_name="porosity"/>
<secondary_variable internal_name="transport_porosity" output_name="transport_porosity"/>
<secondary_variable internal_name="dry_density_solid" output_name="dry_density_solid"/>
</secondary_variables>
<specific_body_force>0 0</specific_body_force>
<initial_stress>sigma0</initial_stress>
<explicit_hm_coupling_in_unsaturated_zone>false</explicit_hm_coupling_in_unsaturated_zone>
<mass_lumping>true</mass_lumping>
</process>
</processes>
<media>
<medium>
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1e-3</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>1e3</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>density</name>
<type>Constant</type>
<value>2780.0</value>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>permeability</name>
<type>Constant</type>
<value>1e-12</value>
</property>
<property>
<name>porosity</name>
<type>PorosityFromMassBalance</type>
<initial_porosity>phi0</initial_porosity>
<minimal_porosity>0</minimal_porosity>
<maximal_porosity>1</maximal_porosity>
</property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
<value>293.15</value>
</property>
<property>
<name>relative_permeability</name>
<type>RelativePermeabilityVanGenuchten</type>
<residual_liquid_saturation>0.0</residual_liquid_saturation>
<residual_gas_saturation>0</residual_gas_saturation>
<exponent>0.01</exponent>
<minimum_relative_permeability_liquid>1e-2</minimum_relative_permeability_liquid>
</property>
<property>
<name>saturation</name>
<type>SaturationVanGenuchten</type>
<residual_liquid_saturation>0.0</residual_liquid_saturation>
<residual_gas_saturation>0</residual_gas_saturation>
<exponent>0.4</exponent>
<p_b>1e4</p_b>
</property>
<property>
<name>bishops_effective_stress</name>
<type>BishopsPowerLaw</type>
<exponent>1</exponent>
</property>
<property>
<name>biot_coefficient</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</medium>
</media>
<time_loop>
<processes>
<process ref="RM">
<nonlinear_solver>basic_newton</nonlinear_solver>
<convergence_criterion>
<type>PerComponentDeltaX</type>
<norm_type>NORM2</norm_type>
<reltols>1e-10 1e0 1e-9</reltols>
</convergence_criterion>
<compensate_non_equilibrium_initial_residuum>false</compensate_non_equilibrium_initial_residuum>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>IterationNumberBasedTimeStepping</type>
<t_initial>0</t_initial>
<t_end>1000</t_end>
<initial_dt>1000</initial_dt>
<minimum_dt>100</minimum_dt>
<maximum_dt>1e5</maximum_dt>
<number_iterations>1 5 10 15</number_iterations>
<multiplier>1.2 1.0 0.8 0.6</multiplier>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>mfront_restart_part1</prefix>
<suffix>_ts_{:timestep}_t_{:time}</suffix>
<timesteps>
<pair>
<repeat>1</repeat>
<each_steps>1</each_steps>
</pair>
</timesteps>
<variables>
<variable>displacement</variable>
<variable>pressure</variable>
<variable>sigma</variable>
<variable>epsilon</variable>
<variable>velocity</variable>
<variable>saturation</variable>
<variable>porosity</variable>
</variables>
</output>
</time_loop>
<parameters>
<!--Modified Cam clay parameters-->
<parameter>
<name>sigma0</name>
<type>Function</type>
<expression>-5e3</expression>
<expression>-5e3</expression>
<expression>-5e3</expression>
<expression>0</expression>
</parameter>
<parameter>
<name>YoungModulus</name>
<type>Constant</type>
<value>52e6</value><!--Pa-->
</parameter>
<parameter>
<name>PoissonRatio</name>
<type>Constant</type>
<value>0.3</value>
</parameter>
<parameter>
<name>CriticalStateLineSlope</name>
<type>Constant</type>
<value>1.2</value>
</parameter>
<parameter>
<name>SwellingLineSlope</name>
<type>Constant</type>
<value>6.6e-3</value>
</parameter>
<parameter>
<name>VirginConsolidationLineSlope</name>
<type>Constant</type>
<value>7.7e-2</value>
</parameter>
<parameter>
<name>InitialPreConsolidationPressure</name>
<type>Constant</type>
<value>2e5</value>
</parameter>
<parameter>
<name>InitialPorosity</name>
<type>Constant</type>
<value>0.44</value>
</parameter>
<!-- Initial and boundary values -->
<parameter>
<name>rho_sr</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>displacement0</name>
<type>Constant</type>
<values>0 0</values>
</parameter>
<parameter>
<name>zero</name>
<type>Constant</type>
<value>0.0</value>
</parameter>
<parameter>
<name>phi0</name>
<type>Constant</type>
<value>0.432</value>
</parameter>
<parameter>
<name>top_pressure</name>
<type>Constant</type>
<value>-1e5</value>
</parameter>
<parameter>
<name>pressure_ic</name>
<type>Constant</type>
<value>-5e3</value>
</parameter>
</parameters>
<process_variables>
<process_variable>
<name>displacement</name>
<components>2</components>
<order>1</order>
<initial_condition>displacement0</initial_condition>
<boundary_conditions>
<!--fix left in radial direction-->
<boundary_condition>
<mesh>square_1x1_quad_1e0_left</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>zero</parameter>
</boundary_condition>
<!--fix bottom in axial direction-->
<boundary_condition>
<mesh>square_1x1_quad_1e0_bottom</mesh>
<type>Dirichlet</type>
<component>1</component>
<parameter>zero</parameter>
</boundary_condition>
<!--compression in axial direction -->
<boundary_condition>
<mesh>square_1x1_quad_1e0_top</mesh>
<type>Neumann</type>
<component>1</component>
<parameter>top_pressure</parameter>
</boundary_condition>
<!--compression in -radial direction-->
<boundary_condition>
<mesh>square_1x1_quad_1e0_right</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>zero</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>square_1x1_quad_1e0_top</mesh>
<type>Dirichlet</type>
<parameter>pressure_ic</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_newton</name>
<type>Newton</type>
<max_iter>60</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>
<test_definition>
<vtkdiff>
<regex>mfront_restart_part1_ts_.*_t_.*.vtu</regex>
<field>displacement</field>
<absolute_tolerance>2e-14</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
<vtkdiff>
<regex>mfront_restart_part1_ts_.*_t_.*.vtu</regex>
<field>sigma</field>
<absolute_tolerance>1e-9</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
<vtkdiff>
<regex>mfront_restart_part1_ts_.*_t_.*.vtu</regex>
<field>epsilon</field>
<absolute_tolerance>1e-15</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
<vtkdiff>
<regex>mfront_restart_part1_ts_.*_t_.*.vtu</regex>
<field>pressure</field>
<absolute_tolerance>1e-15</absolute_tolerance>
<relative_tolerance>1e-15</relative_tolerance>
</vtkdiff>
<vtkdiff>
<regex>mfront_restart_part1_ts_.*_t_.*.vtu</regex>
<field>velocity</field>
<absolute_tolerance>1e-15</absolute_tolerance>
<relative_tolerance>1e-15</relative_tolerance>
</vtkdiff>
<vtkdiff>
<regex>mfront_restart_part1_ts_.*_t_.*.vtu</regex>
<field>HydraulicFlow</field>
<absolute_tolerance>1e-15</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
<vtkdiff>
<regex>mfront_restart_part1_ts_.*_t_.*.vtu</regex>
<field>NodalForces</field>
<absolute_tolerance>1e-10</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
</test_definition>
</OpenGeoSysProject>
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
<UnstructuredGrid>
<FieldData>
<DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="1173" format="appended" RangeMin="34" RangeMax="125" offset="0" />
<DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45" RangeMax="103" offset="352" />
<DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.0014206039571" RangeMax="0.0014206039571" offset="436" />
<DataArray type="Float64" Name="material_state_variable_ElasticStrain_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.0014600790048" RangeMax="0.0014600790048" offset="552" />
<DataArray type="Float64" Name="material_state_variable_EquivalentPlasticStrain_ip" NumberOfTuples="4" format="appended" RangeMin="0" RangeMax="0" offset="660" />
<DataArray type="Float64" Name="material_state_variable_PlasticVolumetricStrain_ip" NumberOfTuples="4" format="appended" RangeMin="0" RangeMax="0" offset="720" />
<DataArray type="Float64" Name="material_state_variable_Porosity_ip" NumberOfTuples="4" format="appended" RangeMin="0.43920389644" RangeMax="0.43920389644" offset="780" />
<DataArray type="Float64" Name="material_state_variable_PreConsolidationPressure_ip" NumberOfTuples="4" format="appended" RangeMin="200000" RangeMax="200000" offset="852" />
<DataArray type="Float64" Name="material_state_variable_Time_ip" NumberOfTuples="4" format="appended" RangeMin="1000" RangeMax="1000" offset="920" />
<DataArray type="Float64" Name="material_state_variable_VolumeRatio_ip" NumberOfTuples="4" format="appended" RangeMin="1.783179294" RangeMax="1.783179294" offset="988" />
<DataArray type="Float64" Name="porosity_ip" NumberOfTuples="4" format="appended" RangeMin="0.43119194903" RangeMax="0.43119194903" offset="1060" />
<DataArray type="Float64" Name="saturation_ip" NumberOfTuples="4" format="appended" RangeMin="0.89688116795" RangeMax="0.89858628419" offset="1132" />
<DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="124270.51012" RangeMax="124270.51012" offset="1216" />
<DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0" RangeMax="0" offset="1344" />
<DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="4" format="appended" RangeMin="0.432" RangeMax="0.432" offset="1404" />
</FieldData>
<Piece NumberOfPoints="4" NumberOfCells="1" >
<PointData>
<DataArray type="Float64" Name="ElasticStrain" NumberOfComponents="4" format="appended" RangeMin="0.0014600790048" RangeMax="0.0014600790048" offset="1476" />
<DataArray type="Float64" Name="EquivalentPlasticStrain" format="appended" RangeMin="0" RangeMax="0" offset="1620" />
<DataArray type="Float64" Name="HydraulicFlow" format="appended" RangeMin="-0.001336947427" RangeMax="8.3700407716e-17" offset="1680" />
<DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="104719.75512" RangeMax="249548.51621" offset="1768" />
<DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0.0014206039571" offset="1888" />
<DataArray type="Float64" Name="dry_density_solid" format="appended" RangeMin="1581.2863817" RangeMax="1581.2863817" offset="1968" />
<DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.0014206039571" RangeMax="0.0014206039571" offset="2040" />
<DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.43119194903" RangeMax="0.43119194903" offset="2164" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="-5000" RangeMax="-4896.7038722" offset="2240" />
<DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="-5000" RangeMax="-4896.7038722" offset="2316" />
<DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.89625705209" RangeMax="0.89921040004" offset="2392" />
<DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="124270.51012" RangeMax="124270.51012" offset="2476" />
<DataArray type="Float64" Name="swelling_stress" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="2620" />
<DataArray type="Float64" Name="transport_porosity" format="appended" RangeMin="0.432" RangeMax="0.432" offset="2680" />
<DataArray type="Float64" Name="velocity" NumberOfComponents="2" format="appended" RangeMin="1.0329612775e-09" RangeMax="1.0329612775e-09" offset="2760" />
</PointData>
<CellData>
<DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="2840" />
<DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0.43119194903" RangeMax="0.43119194903" offset="2900" />
<DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.89773372607" RangeMax="0.89773372607" offset="2968" />
<DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="124270.51012" RangeMax="124270.51012" offset="3036" />
</CellData>
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="3128" />
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="3200" />
<DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="3272" />
<DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="3332" />
</Cells>
</Piece>
</UnstructuredGrid>
<AppendedData encoding="base64">
_AQAAAAAAAAAAgAAAAAAAAJUEAAAAAAAA5gAAAAAAAAA=eF61kjFPAzEMhf9L5i4gpltR9xMgFoQstzUnS0kcbKeoqvrfSSkLQxuk9rY4lr/n95J94Ow0KTpLhiKtAlTFnYXhbf+nKbohDcP9ImRMFIZgPCUELqHd1LQiBfmAtaQimbI3wMNh0UOg19/eBc5dl1NExdh311FcMVtDOdyGZ18UI+cJzJXMrouKinG8nFMfktBJGWPbqJ1gi61YRYJlRHNeP7cEeC6Jz8pbjA0z/lesn/A5sVHpUbJJ5M3P5HjMvyrNpHYy9CqxJnKd2dpNvuY5+snD03FqHoEXTp1XeD98AzuXtOw=AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPRM9C1MNVNT0sxMje1sDSwMAUAL1kEmw==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAANAAAAAAAAAA=eF5jYICAlSyhjwrdwvdDuQz7ueQZjnn82Q3jo8ufsa5Pm7dAbA+Mv5qAfnR5mH4ABuUeGw==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAALgAAAAAAAAA=eF5r3SuaZqHFsn/qptalh5+H72+B8hnAgGVPK155Bbj8dAL6scsr7AEAXCAvDQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAACwAAAAAAAAA=eF5jYMAPAAAgAAE=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAACwAAAAAAAAA=eF5jYMAPAAAgAAE=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF7L1Hqw4pX0HftMHDQAFtcQ7Q==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5jYACCLA5HBhw0ACfsAs0=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5jYAACh34HBhw0ADtoBD0=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF57JPCC6Xn7H/tHOGgAKqcSFQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF4rY7wpu2zGbfsyHDQA89wPFQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAHAAAAAAAAAA=eF5bcvYZq8XBN/Zf9k1Xctj8xn4JGh8AMA8SEw==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPgAAAAAAAAA=eF5bwNd4xtnh+QH+idPyltT/PLAAyl/YUBlx7faEvQtwyM+Z9fVz2M9VexdB+UJQ+UVo+nHJw/QDAOGFSL8=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF67ee578OOlt+1v4qABeMsWTQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAASwAAAAAAAAA=eF5r2iuaZqHFsn/iptalh5+H72+E8hdsMWHZoXF/dwOUPwkqXw/lF0Dl66D8Pqg8jB8Rdc9+F5vKnlYofwZUvgXKj4LKAwDZXjplAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAACwAAAAAAAAA=eF5jYMAPAAAgAAE=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAIAAAAAAAAAA=eF5jYAAChV4bBjA9w6bhQU7chOeu+1vBdOh+AHeiCz0=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAOAAAAAAAAAA=eF5jYACBhbYmUT9E/kz46SBhGV3FN4HhoCmYz+kIlmYosDWGyB/4yenH2dDLcNAMIn8QAK1AF9g=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGQAAAAAAAAA=eF5jYCAOrGAJfVToFr4fxl8D5QMAWu4HZQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF6z/W7rqLp1hoMVGu0EpQHgjA2QAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOgAAAAAAAAA=eF5jYICApSyhjwrdwvdDuQwP2i6Hp3Bc2w3jL0OTP4AmvwBNnvsY50LPIJk9MP4aNHkeqDwAXIcdjw==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF4rYLwpu2zGbfs8KJ0LpUuhNADxxA79AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF7z3PXDYIvC5gOeUJoBBDoQNADwLAzbAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF7z3PXDYIvC5gOeUJoBBDoQNADwLAzbAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAHgAAAAAAAAA=eF5b4dX7OOTYG/slUPqR40cX5XVv7F9DaQAvQxH8AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAASQAAAAAAAAA=eF6bzdd4xtnh+QH2idPyltT/PDAbyr9U7MZ23zlr71wonwcqD+NvgMrPgPLZoPIw/tWfz8SbHm7YuwjKF4LKw/jXoPIA7UZFxQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGQAAAAAAAAA=eF67eu578OOlt+0vQekLUPoWlAYAdxsWOg==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGgAAAAAAAAA=eF5jYICAO4WTGJn3C9oxEMl/COUDAFsbC8o=AQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF4rY7wpu2zGbXsAD7sDxg==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF7zcbMXsdn1xh4AC2UDBw==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF5bwNd4xtnh+QGBidPyltT/PLAAyi/TlRRYuHnhXgALqBBKAQAAAAAAAAAAgAAAAAAAAGAAAAAAAAAAFQAAAAAAAAA=eF5jYMAHPtjjlcaQh/ER4gCW5AS9AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYIAARijNDKWZoDQAAHgABw==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYYAAAAAoAAU=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jBAAACgAK
</AppendedData>
</VTKFile>
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<meshes>
<mesh axially_symmetric="true">mfront_restart_part1_ts_1_t_1000.000000.vtu</mesh>
<mesh axially_symmetric="true">square_1x1_quad_1e0_left.vtu</mesh>
<mesh axially_symmetric="true">square_1x1_quad_1e0_right.vtu</mesh>
<mesh axially_symmetric="true">square_1x1_quad_1e0_top.vtu</mesh>
<mesh axially_symmetric="true">square_1x1_quad_1e0_bottom.vtu</mesh>
</meshes>
<processes>
<process>
<name>RM</name>
<type>RICHARDS_MECHANICS</type>
<integration_order>2</integration_order>
<dimension>2</dimension>
<jacobian_assembler>
<type>Analytical</type>
</jacobian_assembler>
<constitutive_relation>
<type>MFront</type>
<behaviour>ModCamClay_semiExplParaInit</behaviour>
<material_properties>
<material_property name="YoungModulus" parameter="YoungModulus"/>
<material_property name="PoissonRatio" parameter="PoissonRatio"/>
<material_property name="CriticalStateLineSlope" parameter="CriticalStateLineSlope"/>
<material_property name="SwellingLineSlope" parameter="SwellingLineSlope"/>
<material_property name="VirginConsolidationLineSlope" parameter="VirginConsolidationLineSlope"/>
<material_property name="InitialPreConsolidationPressure" parameter="InitialPreConsolidationPressure"/>
<material_property name="InitialPorosity" parameter="InitialPorosity"/>
</material_properties>
</constitutive_relation>
<process_variables>
<pressure>pressure</pressure>
<displacement>displacement</displacement>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="ElasticStrain" output_name="ElasticStrain"/>
<secondary_variable internal_name="EquivalentPlasticStrain" output_name="EquivalentPlasticStrain"/>
<secondary_variable internal_name="sigma" output_name="sigma"/>
<secondary_variable internal_name="swelling_stress" output_name="swelling_stress"/>
<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_variable internal_name="porosity" output_name="porosity"/>
<secondary_variable internal_name="transport_porosity" output_name="transport_porosity"/>
<secondary_variable internal_name="dry_density_solid" output_name="dry_density_solid"/>
</secondary_variables>
<specific_body_force>0 0</specific_body_force>
<explicit_hm_coupling_in_unsaturated_zone>false</explicit_hm_coupling_in_unsaturated_zone>
<mass_lumping>true</mass_lumping>
</process>
</processes>
<media>
<medium>
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1e-3</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>1e3</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>density</name>
<type>Constant</type>
<value>2780.0</value>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>permeability</name>
<type>Constant</type>
<value>1e-12</value>
</property>
<property>
<name>porosity</name>
<type>PorosityFromMassBalance</type>
<initial_porosity>phi0</initial_porosity>
<minimal_porosity>0</minimal_porosity>
<maximal_porosity>1</maximal_porosity>
</property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
<value>293.15</value>
</property>
<property>
<name>relative_permeability</name>
<type>RelativePermeabilityVanGenuchten</type>
<residual_liquid_saturation>0.0</residual_liquid_saturation>
<residual_gas_saturation>0</residual_gas_saturation>
<exponent>0.01</exponent>
<minimum_relative_permeability_liquid>1e-2</minimum_relative_permeability_liquid>
</property>
<property>
<name>saturation</name>
<type>SaturationVanGenuchten</type>
<residual_liquid_saturation>0.0</residual_liquid_saturation>
<residual_gas_saturation>0</residual_gas_saturation>
<exponent>0.4</exponent>
<p_b>1e4</p_b>
</property>
<property>
<name>bishops_effective_stress</name>
<type>BishopsPowerLaw</type>
<exponent>1</exponent>
</property>
<property>
<name>biot_coefficient</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</medium>
</media>
<time_loop>
<processes>
<process ref="RM">
<nonlinear_solver>basic_newton</nonlinear_solver>
<convergence_criterion>
<type>PerComponentDeltaX</type>
<norm_type>NORM2</norm_type>
<reltols>1e-10 1e0 1e-9</reltols>
</convergence_criterion>
<compensate_non_equilibrium_initial_residuum>false</compensate_non_equilibrium_initial_residuum>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>IterationNumberBasedTimeStepping</type>
<t_initial>1000</t_initial>
<t_end>2000</t_end>
<initial_dt>1000</initial_dt>
<minimum_dt>100</minimum_dt>
<maximum_dt>1e5</maximum_dt>
<number_iterations>1 5 10 15</number_iterations>
<multiplier>1.2 1.0 0.8 0.6</multiplier>
</time_stepping>
</process>