From 30a01f21bbc4838b25ca2c1125ba1fbe9102e317 Mon Sep 17 00:00:00 2001
From: FZill <florian.zill@ufz.de>
Date: Wed, 3 Mar 2021 21:51:15 +0100
Subject: [PATCH] [MPL] EmbeddedFracturePermeability fix for 2D
---
.../EmbeddedFracturePermeability.cpp | 5 +-
ProcessLib/HydroMechanics/Tests.cmake | 16 +-
.../EmbeddedFracturePermeability/square.prj | 242 ++++++++++++++++++
.../EmbeddedFracturePermeability/square.vtu | 42 +++
.../square_1x1.gml | 31 +++
.../square_ts_1_t_1.000000.vtu | 41 +++
6 files changed, 375 insertions(+), 2 deletions(-)
create mode 100644 Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square.prj
create mode 100644 Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square.vtu
create mode 100644 Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square_1x1.gml
create mode 100644 Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square_ts_1_t_1.000000.vtu
diff --git a/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.cpp b/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.cpp
index d4edd046a14..8c01cb4ec6c 100644
--- a/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.cpp
+++ b/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.cpp
@@ -72,7 +72,9 @@ PropertyDataType EmbeddedFracturePermeability<DisplacementDim>::value(
double const coeff = H_de * (b_f / _a) * ((b_f * b_f / 12.0) - _k);
Eigen::Matrix3d I = Eigen::Matrix3d::Identity();
- return (coeff * (I - n * n.transpose()) + _k * I).eval();
+ return (coeff * (I - n * n.transpose()) + _k * I)
+ .template topLeftCorner<DisplacementDim, DisplacementDim>()
+ .eval();
}
template <int DisplacementDim>
@@ -106,6 +108,7 @@ PropertyDataType EmbeddedFracturePermeability<DisplacementDim>::dValue(
Eigen::Matrix3d const M = n * n.transpose();
return (H_de * (b_f * b_f / 4 - _k) * (Eigen::Matrix3d::Identity() - M) * M)
+ .template topLeftCorner<DisplacementDim, DisplacementDim>()
.eval();
}
diff --git a/ProcessLib/HydroMechanics/Tests.cmake b/ProcessLib/HydroMechanics/Tests.cmake
index a06bf438605..a6a663bf812 100644
--- a/ProcessLib/HydroMechanics/Tests.cmake
+++ b/ProcessLib/HydroMechanics/Tests.cmake
@@ -283,7 +283,21 @@ AddTest(
# Permeability models
AddTest(
- NAME HydroMechanics_Permeability_EmbeddedFracture
+ NAME HydroMechanics_Permeability_EmbeddedFracture_square
+ PATH HydroMechanics/EmbeddedFracturePermeability
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS square.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT OGS_USE_MPI
+ DIFF_DATA
+ square_ts_1_t_1.000000.vtu square_ts_1_t_1.000000.vtu pressure pressure 0 1e-14
+ square_ts_1_t_1.000000.vtu square_ts_1_t_1.000000.vtu velocity velocity 1e-15 0
+ square_ts_1_t_1.000000.vtu square_ts_1_t_1.000000.vtu displacement displacement 1e-15 0
+)
+
+AddTest(
+ NAME HydroMechanics_Permeability_EmbeddedFracture_cube
PATH HydroMechanics/EmbeddedFracturePermeability
EXECUTABLE ogs
EXECUTABLE_ARGS cube.prj
diff --git a/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square.prj b/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square.prj
new file mode 100644
index 00000000000..dd48cabe3fe
--- /dev/null
+++ b/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square.prj
@@ -0,0 +1,242 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>square.vtu</mesh>
+ <geometry>square_1x1.gml</geometry>
+ <processes>
+ <process>
+ <name>HM</name>
+ <type>HYDRO_MECHANICS</type>
+ <integration_order>3</integration_order>
+ <dimension>2</dimension>
+ <constitutive_relation>
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <process_variables>
+ <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_variables>
+ <specific_body_force>0 0</specific_body_force>
+ <mass_lumping>false</mass_lumping>
+ </process>
+ </processes>
+ <media>
+ <medium>
+ <phases>
+ <phase>
+ <type>Gas</type>
+ <properties>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>2.15e3</value>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.0 </value>
+ </property>
+ <property>
+ <name>reference_temperature</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>EmbeddedFracturePermeability</type>
+ <intrinsic_permeability>1e-20</intrinsic_permeability>
+ <initial_aperture>0</initial_aperture>
+ <mean_frac_distance>0.01</mean_frac_distance>
+ <threshold_strain>1e-5</threshold_strain>
+ <fracture_normal>1 0 0</fracture_normal>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <processes>
+ <process ref="HM">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltols>1e-15 1e1 1e-15</reltols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>1</t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>square</prefix>
+ <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>
+ </variables>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ </output>
+ </time_loop>
+ <parameters>
+ <!-- Mechanics -->
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1e9</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.3</value>
+ </parameter>
+ <!-- Model parameters -->
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0</values>
+ </parameter>
+ <parameter>
+ <name>pressure_ic</name>
+ <type>Constant</type>
+ <values>1e5</values>
+ </parameter>
+ <parameter>
+ <name>neumann</name>
+ <type>Constant</type>
+ <values>1e-10</values>
+ </parameter>
+ <parameter>
+ <name>dirichlet0</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>displacementTop</name>
+ <type>Constant</type>
+ <value>1e-4</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>square_1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>displacementTop</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>square_1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Neumann</type>
+ <component>0</component>
+ <parameter>neumann</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <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>50</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <lis>-i bicgstab -p ilu -tol 1e-10 -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>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square.vtu b/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square.vtu
new file mode 100644
index 00000000000..9cb905c45b3
--- /dev/null
+++ b/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square.vtu
@@ -0,0 +1,42 @@
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="9" NumberOfCells="1">
+ <PointData>
+ </PointData>
+ <CellData>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="1.4142135623730951">
+ 2AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADwPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADwPwAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAD+///////fPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADwP/7//////98/AAAAAAAAAAAAAAAAAADgPwAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOA/AAAAAAAAAAAAAAAAAADgPwAAAAAAAOA/AAAAAAAAAAA=
+ <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+ <Value index="0">
+ 0
+ </Value>
+ <Value index="1">
+ 1.4142135624
+ </Value>
+ </InformationKey>
+ <InformationKey name="L2_NORM_FINITE_RANGE" location="vtkDataArray" length="2">
+ <Value index="0">
+ 0
+ </Value>
+ <Value index="1">
+ 1.4142135624
+ </Value>
+ </InformationKey>
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="8">
+ SAAAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAA=
+ </DataArray>
+ <DataArray type="Int64" Name="offsets" format="binary" RangeMin="9" RangeMax="9">
+ CAAAAAAAAAAJAAAAAAAAAA==
+ </DataArray>
+ <DataArray type="UInt8" Name="types" format="binary" RangeMin="28" RangeMax="28">
+ AQAAAAAAAAAc
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square_1x1.gml b/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square_1x1.gml
new file mode 100644
index 00000000000..2f637aa9f95
--- /dev/null
+++ b/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square_1x1.gml
@@ -0,0 +1,31 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml-stylesheet type="text/xsl" href="OpenGeoSysGLI.xsl"?>
+
+<OpenGeoSysGLI xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ogs="http://www.opengeosys.org">
+ <name>square_1x1_geometry</name>
+ <points>
+ <point id="0" x="0" y="0" z="0"/>
+ <point id="1" x="0" y="1" z="0"/>
+ <point id="2" x="1" y="0" z="0"/>
+ <point id="3" x="1" y="1" z="0"/>
+ </points>
+
+ <polylines>
+ <polyline id="0" name="left">
+ <pnt>0</pnt>
+ <pnt>1</pnt>
+ </polyline>
+ <polyline id="1" name="right">
+ <pnt>2</pnt>
+ <pnt>3</pnt>
+ </polyline>
+ <polyline id="2" name="bottom">
+ <pnt>0</pnt>
+ <pnt>2</pnt>
+ </polyline>
+ <polyline id="3" name="top">
+ <pnt>1</pnt>
+ <pnt>3</pnt>
+ </polyline>
+ </polylines>
+</OpenGeoSysGLI>
diff --git a/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square_ts_1_t_1.000000.vtu b/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square_ts_1_t_1.000000.vtu
new file mode 100644
index 00000000000..32e2e24b273
--- /dev/null
+++ b/Tests/Data/HydroMechanics/EmbeddedFracturePermeability/square_ts_1_t_1.000000.vtu
@@ -0,0 +1,41 @@
+<?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="166" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45" RangeMax="103" offset="176" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.00010879675866" RangeMax="0.00010879675866" offset="260" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="114728.64296" RangeMax="114728.64296" offset="456" />
+ </FieldData>
+ <Piece NumberOfPoints="9" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="HydraulicFlow" format="appended" RangeMin="-5e-11" RangeMax="5e-11" offset="664" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="2.9515836866e-12" RangeMax="73260.07326" offset="752" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0.00010879675866" offset="960" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.00010879675866" RangeMax="0.00010879675866" offset="1088" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="0" RangeMax="10100000" offset="1328" />
+ <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="100000" RangeMax="10100000" offset="1416" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="114728.64296" RangeMax="114728.64296" offset="1516" />
+ <DataArray type="Float64" Name="velocity" NumberOfComponents="2" format="appended" RangeMin="1e-10" RangeMax="1e-10" offset="1848" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="permeability" NumberOfComponents="4" format="appended" RangeMin="1.4142135624e-20" RangeMax="1.4142135624e-20" offset="2000" />
+ <DataArray type="Float64" Name="principal_stress_values" NumberOfComponents="3" format="appended" RangeMin="114728.64296" RangeMax="114728.64296" offset="2072" />
+ <DataArray type="Float64" Name="principal_stress_vector_1" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="2156" />
+ <DataArray type="Float64" Name="principal_stress_vector_2" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="2224" />
+ <DataArray type="Float64" Name="principal_stress_vector_3" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1" offset="2288" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="2352" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="2448" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2536" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="2596" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAKYAAAAAAAAAYgAAAAAAAAA=eF6FzDEKgDAMheG7ZO6mU68iEqJGCdikpHUQ8e52ddHxvR++C0Qrb05VTDFbW0judBaIw/WK5gs7xC6AUmKIUGRLhJKhPUea2NFWnC1lU9bagP4OPwTnInsLn8h4PyXvPSI=AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPRM9A1NDTUTTcxMzI2T7UwtwAAK4QEXQ==AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAcAAAAAAAAAA=eF67ddrG3KCCbb+Os8zrR2ZS9gxo4A5UXg9N/ox1fdq8Bces70HltbHKb7O+j8P8ssmyql3z9lvfhcrr4rD/Lg77TTfG3Y9YfGj3Lai8Fpr8zZz7Z8p+PIC7D10/DNzB4b7gkPUzu848tgYAFA9UEg==AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAegAAAAAAAAA=eF5jYACBHbaHFGUOKV755cD3UIDvocQDBwYUsMP2MFReACp/yGzOp+umC2yh8nv3Q+V5UOQnwOUPoMkbbqjc8XHZNFuYDQfx2r8Abj8/VP7RXoWt6Vtm7YXJo5uvHPTq9eGAHXDzDxMwH93+K/ETX1jG7bEFAHipYec=AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAIQAAAAAAAAA=eF7bs/f6zfs1p213Qui9u6H0Hqg4Axj8qCdEAwCkCh6AAQAAAAAAAAAAgAAAAAAAAJAAAAAAAAAAegAAAAAAAAA=eF4rWnBfeefdnr2qB4V0Dj66eMD9Wvm9Fbdm22pD+fGz/1jzPCjZqwLhOxQzVan1s823VYPy1/IX7uH/OnGvPpj/8YB1uo/xS+sZtgxgsMHWQdDP/1Rgma0eRN7hyFE11Uus7VD5ir1mExYdcPNZvxfCP2QLANHmPkQ=AQAAAAAAAAAAgAAAAAAAAJAAAAAAAAAAPgAAAAAAAAA=eF5jYEAF907bmBtUsO1H5+s6y7x+ZCZlDxOH8W+D5b/tg4nfharXAstz2d+DyqPr14bK34HK60D5AMbUKTg=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAkwAAAAAAAAA=eF67cNrG3KCCbb+ys8zrR2ZS9gxQUKp7p/lfx9bd56Hy8mjyvTe/crRFfLa+DZXXRJMPfxOhYSvxy/oBVF4dTX6L/dQ7c+ZctX4BlTdDk//+68z2fw8OW9+FymujybftyK5yYFgIN98ITV7jUk3Sqekndj+FylugyVc/LZyttvuC9R0c5vNZvT2zMHTHDgCUb2JIAQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAHwAAAAAAAAA=eF5jYGBgWOKc7MgAAlk/HGD0fyBYDBMnAgAAE2QKVg==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAKQAAAAAAAAA=eF5jYGBgWOKc7MgAAlk/HGD0fyBYDBWvKAl2RBYvh/KXIMkDAPy1Ff4=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAA1wAAAAAAAAA=eF5bmS9l7bXyuO0+RZlDild+OfA8FOB7KPHA4SnXZAuu2317nUvXmXbPOGa7AyrPDpVP91e89G/7adugG401zNwzbPdA5fmg8r4l8s17Nl6y9Xja6Prh28W9u6DyrFD5zwsDndS+r7P9cSR/wvx/B/aegMoLQuUDi4RsO6rn2IYkJc514phkexAqzwuVZyyaV/EmusL2Q83lG4Fec23R9X+W3C8dOG/h3gOXg9d833fQ9gJUXhwqv998dtSrH8ttGUAgYcqu/Wjmh5tv2/I6b+IuAFF6h8Q=AQAAAAAAAAAAgAAAAAAAAJAAAAAAAAAAUAAAAAAAAAA=eF7bsvf6zfs1t225NFcy+opMttoC5d/i4fF4uynHag+Uz8aLyheCqj8E5TMxgECD1TYo/yNU/14onwEqfxTKF4Dq3wHl//0PAvVWANAxSoo=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF5TcprEs3DOcWslKM2ABgCblgaDAQAAAAAAAAAAgAAAAAAAABgAAAAAAAAAHgAAAAAAAAA=eF4LDQWBVlu+hwJ8DyUeOBxUlDmkeOWXAwBvIgnWAQAAAAAAAAAAgAAAAAAAABgAAAAAAAAAEAAAAAAAAAA=eF5jYACBD/YMaAAAFScBMA==AQAAAAAAAAAAgAAAAAAAABgAAAAAAAAADgAAAAAAAAA=eF5jYMAGPtgDAAI3ATA=AQAAAAAAAAAAgAAAAAAAABgAAAAAAAAADwAAAAAAAAA=eF5jYEAGH+xhLAALrwEwAQAAAAAAAAAAgAAAAAAAANgAAAAAAAAAJgAAAAAAAAA=eF5jYMAHPtjjFyckz8Dw7z8I3MepDrv8AwLmP8ChHiEOAMCFGcU=AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAIAAAAAAAAAA=eF5jYIAARijNBKWZoTQLlGaF0mxQmh1Kc0BpAAQIACU=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jZIAAAABQAAo=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF6TAQAAHQAd
+ </AppendedData>
+</VTKFile>
--
GitLab