Merge branch 'GS_model' into 'master'

Guenther Salzer model with damage See merge request ogs/ogs!4516
ufz · Jul 20, 2023 · 6fea0a8 · 6fea0a8
2 parents bca07ca + 3018ea1
commit 6fea0a8
Show file tree

Hide file tree

Showing 14 changed files with 7,670 additions and 132 deletions.
diff --git a/MaterialLib/SolidModels/MFront/CMakeLists.txt b/MaterialLib/SolidModels/MFront/CMakeLists.txt
@@ -29,6 +29,7 @@ set(_mfront_behaviours
  BDT
  DruckerPrager
  GuentherSalzer
+ GuentherSalzerDilatancy_semi_expl
  Lubby2
  Lubby2mod
  ModCamClay_semiExpl

diff --git a/MaterialLib/SolidModels/MFront/GuentherSalzerDilatancy_semi_expl.mfront b/MaterialLib/SolidModels/MFront/GuentherSalzerDilatancy_semi_expl.mfront
@@ -0,0 +1,193 @@
+
+@DSL Implicit;
+@Behaviour GuentherSalzerDilatancy_semi_expl;
+@Author Thomas Nagel;
+@Description
+{
+Günther Salzer model for rock salt.
+Semi-explicit treatment of damage/dilatancy. Requires very small time
+steps, otherwise neither robust nor accurate.
+Fully implicit treatment WIP.
+
+Günther, R., & Salzer, K. (2012). Advanced strain-hardening approach.
+In Mechanical Behaviour of Salt VII (Issue December 2016).
+CRC Press. https://doi.org/10.1201/b12041-4
+
+Günther, R.-M., Salzer, K., Popp, T., & Lüdeling, C. (2015).
+Steady-State Creep of Rock Salt: Improved Approaches for Lab Determination
+and Modelling. Rock Mechanics and Rock Engineering,
+48(6), 2603–2613. https://doi.org/10.1007/s00603-015-0839-2
+
+Fixed parameter set: Asse.
+Units: MPa d
+};
+
+@Algorithm NewtonRaphson;
+@MaximumNumberOfIterations 100;
+//@CompareToNumericalJacobian true;
+//@PerturbationValueForNumericalJacobianComputation 1.e-8;
+//@JacobianComparisonCriterion 1e-7;
+
+@Epsilon 1.e-14;
+@Theta 1.0;
+
+@ModellingHypotheses{".+"};
+@Brick StandardElasticity;
+
+// Activation energies in input files in J/mol
+@Parameter real Ru = 8.314472; // J/(Kmol)
+Ru.setEntryName("UniversalGasConstant");
+
+// Units: MPa, d
+
+// Reference stress
+@Parameter real sig0 = 1.;
+
+// initial hardenging -- can be critical for initial convergence!
+@Parameter real epsV0 = 1.e-8;
+
+// secondary power law
+@Parameter real As1 = 0.001;
+@Parameter real ns1 = 4.0;
+@Parameter real Q1 = 40.e3;
+@Parameter real As2 = 0.001;
+@Parameter real ns2 = 1.0;
+@Parameter real Q2 = 65.e3;
+
+// primary power law and stress/dilatancy dependence
+@Parameter real Ap = 3.98e-29;
+@Parameter real mup = 3.0;
+@Parameter real np0 = 15.18;
+@Parameter real n1 = 0.8;
+@Parameter real n2 = -0.33;
+@Parameter real n3 = 2000.;
+@Parameter real n4 = 0.5;
+
+// Dilatancy boundary
+@Parameter real D1 = 8.;
+@Parameter real D2 = 0.25;
+@Parameter real D3 = 4.8;
+
+// dilatancy strain evolution
+@Parameter real a1 = 5.e-3;
+@Parameter real a2 = 0.0135;
+@Parameter real a3 = -0.4;
+
+@Parameter real a4 = 1.e-6;
+@Parameter real a5 = 2.5e-4;
+@Parameter real a6 = -6.25;
+
+@Parameter real a7 = 0.;
+@Parameter real a8 = 5.13;
+@Parameter real a9 = -0.3;
+
+// Hardening strain
+@StateVariable real epsCrV;
+epsCrV.setEntryName("HardeningStrain");
+
+// Dilatancy strain
+@StateVariable real eps_vol;
+eps_vol.setEntryName("DilatancyStrain");
+
+@AuxiliaryStateVariable stress Udil;
+Udil.setEntryName("DilatancyWork");
+
+// Recovery strain increment
+@LocalVariable real deps_cr_E;
+// Creep strain increment
+@LocalVariable Stensor deps_cr;
+
+@LocalVariable real np; // corrected stress exponents
+
+//! Second Lamé coefficient
+@LocalVariable stress mu;
+//! First Lamé coefficient
+@LocalVariable stress lambda;
+
+@LocalVariable real As1T, As2T, s3_prev, seq_prev_Dil;
+
+@InitLocalVariables
+{
+mu = computeMu(young, nu);
+lambda = computeLambda(young, nu);
+// Compute initial elastic strain
+eel = 1. / (2. * mu) * sig - nu / young * trace(sig) * Stensor::Id();
+As1T = As1 * std::exp(-Q1 / (Ru * (T + dT)));
+As2T = As2 * std::exp(-Q2 / (Ru * (T + dT)));
+
+const Stensor eps_prev = eel - deel;
+const Stensor sig_prev =
+ lambda * trace(eps_prev) * Stensor::Id() + 2. * mu * eps_prev;
+const auto sig_p = sig_prev.computeEigenValues(Stensor::DESCENDING);
+s3_prev = -sig_p(0); // sign change for compression convention
+seq_prev_Dil = s3_prev * (D1 / (D2 + s3_prev) + D3);
+
+np = (np0 + n1 * std::exp(n2 * s3_prev) * (1. - std::exp(-n3 * eps_vol))) /
+ std::pow(1. - eps_vol, n4);
+}
+
+@Integrator
+{
+const auto eeps = real(1e-16);
+const auto seps = eeps * young;
+const auto s = deviator(sig);
+const auto seq = std::max(sigmaeq(s), seps);
+const auto norm_s = seq / std::sqrt(3. / 2.);
+
+constexpr auto Pdev = Stensor4::K();
+const auto A1sig = a1 + a2 * std::exp(a3 * s3_prev);
+const auto A2sig = a4 + a5 * std::exp(a6 * s3_prev);
+const auto A3sig = a7 + a8 * std::exp(a9 * s3_prev);
+
+// secondary (recovery)
+const auto b =
+ As1T * std::pow(seq / sig0, ns1) + As2T * std::pow(seq / sig0, ns2);
+const auto db_dseq = ns1 / sig0 * As1T * std::pow(seq / sig0, ns1 - 1) +
+ ns2 / sig0 * As2T * std::pow(seq / sig0, ns2 - 1);
+
+deps_cr_E = b * dt;
+const auto epsCrV_new = epsCrV + depsCrV;
+const auto seq_to_np = std::pow(seq / sig0, np);
+const auto dseq_ddeel = 3. * s / (2. * seq) * 2. * mu;
+
+const auto deps_cr_tot = depsCrV + deps_cr_E + deps_vol;
+
+fepsCrV =
+ deps_cr_tot - dt * (Ap * seq_to_np / std::pow(epsV0 + epsCrV_new, mup));
+dfepsCrV_ddepsCrV =
+ 1. + dt * mup * Ap * seq_to_np / std::pow(epsV0 + epsCrV_new, mup + 1);
+const auto dfepsCrV_dseq = -dt * (Ap * np * std::pow(seq / sig0, np - 1) /
+ std::pow(epsV0 + epsCrV_new, mup)) +
+ db_dseq * dt;
+dfepsCrV_ddeel = dfepsCrV_dseq * dseq_ddeel;
+dfepsCrV_ddeps_vol = 1.;
+deps_cr = eval(std::sqrt(3. / 2.) * deps_cr_tot * s / norm_s);
+feel += deps_cr;
+dfeel_ddepsCrV = std::sqrt(3. / 2.) / norm_s * s;
+dfeel_ddeps_vol = std::sqrt(3. / 2.) / norm_s * s;
+dfeel_ddeel += 2. * mu * std::sqrt(3. / 2.) / norm_s *
+ (deps_cr_tot * (Pdev - (s ^ s) / norm_s / norm_s) +
+ db_dseq * dt * 3. / (2. * seq) * (s ^ s));
+
+const auto s_m_sdil = seq - seq_prev_Dil;
+if (s_m_sdil > 0.)
+{
+ const auto Ai_Udil = A1sig + A2sig * std::exp(A3sig * Udil);
+ const auto s_m_sdil_Ai_Udil = Ai_Udil * s_m_sdil;
+ feps_vol = deps_vol - s_m_sdil_Ai_Udil * deps_cr_tot;
+ dfeps_vol_ddeps_vol = 1. - s_m_sdil_Ai_Udil;
+ dfeps_vol_ddepsCrV = -s_m_sdil_Ai_Udil;
+ dfeps_vol_ddeel =
+ -(s_m_sdil_Ai_Udil * dt * db_dseq + Ai_Udil * deps_cr_tot) * dseq_ddeel;
+}
+}
+
+@UpdateAuxiliaryStateVariables
+{
+const auto seq = sigmaeq(sig);
+if (seq - seq_prev_Dil > 0.)
+{
+ const auto dUdil = (seq - seq_prev_Dil) * (depsCrV + deps_cr_E + deps_vol);
+ Udil += dUdil;
+}
+}
diff --git a/ProcessLib/SmallDeformation/Tests.cmake b/ProcessLib/SmallDeformation/Tests.cmake
@@ -100,7 +100,8 @@ if (OGS_USE_MFRONT)
  OgsTest(PROJECTFILE Mechanics/Ehlers/MFront/square_1e1_2_matIDs_restart.prj RUNTIME 4)
  OgsTest(PROJECTFILE Mechanics/Ehlers/MFront/two_material_ids_single_solid.prj RUNTIME 1)
  if(NOT OGS_USE_MPI)
- OgsTest(PROJECTFILE Mechanics/GuentherSalzer/model_triaxtest.prj)
+ OgsTest(PROJECTFILE Mechanics/GuentherSalzer/model_uniaxial_test.prj)
+ OgsTest(PROJECTFILE Mechanics/GuentherSalzer/model_triaxial_test.prj RUNTIME 100)
  endif()
  OgsTest(PROJECTFILE Mechanics/PLLC/uniax_compression.prj)
 

diff --git a/Tests/Data/Mechanics/GuentherSalzer/geometry_bottom.vtu b/Tests/Data/Mechanics/GuentherSalzer/geometry_bottom.vtu
@@ -1,34 +1,25 @@
 <?xml version="1.0"?>
-<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">
-<!--This file was created by meshio v4.3.5-->
-<UnstructuredGrid>
-<Piece NumberOfPoints="5" NumberOfCells="4">
-<Points>
-<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary">
-AQAAAACAAAB4AAAAJgAAAA==eJxjYMAHLB3QRf61/AcCCQxxnvkgcQ0M8YxbQOH9RhjiABHwEM8=
-</DataArray>
-</Points>
-<Cells>
-<DataArray type="Int64" Name="connectivity" format="binary">
-AQAAAACAAABAAAAAFgAAAA==eJxjYIAAJjSaGY1mQaMZoTQAApAAFA==
-</DataArray>
-<DataArray type="Int64" Name="offsets" format="binary">
-AQAAAACAAAAgAAAAEwAAAA==eJxjYoAAFijNBqU5oDQAAWAAFQ==
-</DataArray>
-<DataArray type="Int64" Name="types" format="binary">
-AQAAAACAAAAgAAAADgAAAA==eJxjZoAAZhw0AAEQAA0=
-</DataArray>
-</Cells>
-<PointData>
-<DataArray type="UInt64" Name="bulk_node_ids" format="binary">
-AQAAAACAAAAoAAAAFgAAAA==eJxjYIAARijNAqVZoTQblAYAASgAEQ==
-</DataArray>
-</PointData>
-<CellData>
-<DataArray type="UInt64" Name="bulk_elem_ids" format="binary">
-AQAAAACAAAAgAAAAEwAAAA==eJxjYIAAYSitBqUtoTQABhAAcw==
-</DataArray>
-</CellData>
-</Piece>
-</UnstructuredGrid>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="5" NumberOfCells="4" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="6" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_elem_ids" format="appended" RangeMin="0" RangeMax="57" offset="64" />
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="57" offset="120" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="25" offset="176" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="348" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="444" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="500" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _KAAAAAAAAAAAAAAAAAAAAAEAAAAAAAAABAAAAAAAAAAFAAAAAAAAAAYAAAAAAAAAIAAAAAAAAAAAAAAAAAAAABMAAAAAAAAAJgAAAAAAAAA5AAAAAAAAAA==IAAAAAAAAAAAAAAAAAAAABMAAAAAAAAAJgAAAAAAAAA5AAAAAAAAAA==eAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA5QAAAAAAAAAAAAAAAAAAAAAD+hP////8YQAAAAAAAAAAAAAAAAAAAAAAMn/////8oQAAAAAAAAAAAAAAAAAAAAABo2v///78yQAAAAAAAAAAAAAAAAAAAAAA=QAAAAAAAAAAAAAAAAAAAAAIAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAEAAAAAAAAAIAAAAAAAAAACAAAAAAAAAAQAAAAAAAAABgAAAAAAAAAIAAAAAAAAAA==BAAAAAAAAAADAwMD
+ </AppendedData>
 </VTKFile>
diff --git a/Tests/Data/Mechanics/GuentherSalzer/geometry_left.vtu b/Tests/Data/Mechanics/GuentherSalzer/geometry_left.vtu
@@ -1,34 +1,25 @@
 <?xml version="1.0"?>
-<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">
-<!--This file was created by meshio v4.3.5-->
-<UnstructuredGrid>
-<Piece NumberOfPoints="20" NumberOfCells="19">
-<Points>
-<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary">
-AQAAAACAAADgAQAAoQAAAA==eJxjYCAGRDqgi1ipnsrWWB+OIV7iNeV6QFwYhvi6/DjnSt5QDPEXUzTWL9wTjCGutOuT1KnsIAzxmPt72j5JBWKIa7PzPTWd4o8hLmi7CKjBB0P8S5EZ0KFeGOK3V4D85Y4hvv8eyP2uGOJLRT4D3emEIf47H+QeewzxJZtA9lphiPt/BZlviiH+AWyOPoa4O1i9Koa4A1hcFC4OACQZPiE=
-</DataArray>
-</Points>
-<Cells>
-<DataArray type="Int64" Name="connectivity" format="binary">
-AQAAAACAAAAwAQAANQAAAA==eJxdxTcSgCAAADDEgiBY8P9/ZfBckiVT+ESeeeGVN068c+bCB1dufPLFNz/c+eXfAKwoAXw=
-</DataArray>
-<DataArray type="Int64" Name="offsets" format="binary">
-AQAAAACAAACYAAAANQAAAA==eJwtxdEGgDAAAMDJTDLJZCbJTPr/T+xhdy+3hCk6efXm7N2Hi09XN1++/bh7+PXnH1O4AX0=
-</DataArray>
-<DataArray type="Int64" Name="types" format="binary">
-AQAAAACAAACYAAAADwAAAA==eJxjZoAA5kFGAwASaAA6
-</DataArray>
-</Cells>
-<PointData>
-<DataArray type="UInt64" Name="bulk_node_ids" format="binary">
-AQAAAACAAACgAAAAMwAAAA==eJwtxckNgCAAADAS10AFQcFr//V40H4awrQ4evXm3cnZh4urT19u7r79+PXn3wO2UAKV
-</DataArray>
-</PointData>
-<CellData>
-<DataArray type="UInt64" Name="bulk_elem_ids" format="binary">
-AQAAAACAAACYAAAAMQAAAA==eJxVxUkCgBAAAECpiKLl/4/tYC7mMl8YXj++3d18+XR18eHs5N2bV0cvYfYDR9gArA==
-</DataArray>
-</CellData>
-</Piece>
-</UnstructuredGrid>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="20" NumberOfCells="19" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="45" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_elem_ids" format="appended" RangeMin="0" RangeMax="18" offset="224" />
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="18" offset="440" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="100" offset="656" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1308" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="1724" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1940" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _oAAAAAAAAAAAAAAAAAAAAAMAAAAAAAAAHAAAAAAAAAAdAAAAAAAAAB4AAAAAAAAAHwAAAAAAAAAgAAAAAAAAACEAAAAAAAAAIgAAAAAAAAAjAAAAAAAAACQAAAAAAAAAJQAAAAAAAAAmAAAAAAAAACcAAAAAAAAAKAAAAAAAAAApAAAAAAAAACoAAAAAAAAAKwAAAAAAAAAsAAAAAAAAAC0AAAAAAAAAmAAAAAAAAAASAAAAAAAAABEAAAAAAAAAEAAAAAAAAAAPAAAAAAAAAA4AAAAAAAAADQAAAAAAAAAMAAAAAAAAAAsAAAAAAAAACgAAAAAAAAAJAAAAAAAAAAgAAAAAAAAABwAAAAAAAAAGAAAAAAAAAAUAAAAAAAAABAAAAAAAAAADAAAAAAAAAAIAAAAAAAAAAQAAAAAAAAAAAAAAAAAAAA==mAAAAAAAAAASAAAAAAAAABEAAAAAAAAAEAAAAAAAAAAPAAAAAAAAAA4AAAAAAAAADQAAAAAAAAAMAAAAAAAAAAsAAAAAAAAACgAAAAAAAAAJAAAAAAAAAAgAAAAAAAAABwAAAAAAAAAGAAAAAAAAAAUAAAAAAAAABAAAAAAAAAADAAAAAAAAAAIAAAAAAAAAAQAAAAAAAAAAAAAAAAAAAA==4AEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAFlAAAAAAAAAAAAAAAAAAAAAADolymsor1dAAAAAAAAAAAAAAAAAAAAAAHRKlNdQXlZAAAAAAAAAAAAAAAAAAAAAAK5vXkN5DVVAAAAAAAAAAAAAAAAAAAAAAOiUKK+hvFNAAAAAAAAAAAAAAAAAAAAAACK68hrKa1JAAAAAAAAAAAAAAAAAAAAAAFzfvIbyGlFAAAAAAAAAAAAAAAAAAAAAACsHDuU1lE9AAAAAAAAAAAAAAAAAAAAAABE9oryG8kxAAAAAAAAAAAAAAAAAAAAAAPRyNpTXUEpAAAAAAAAAAAAAAAAAAAAAANuoymsor0dAAAAAAAAAAAAAAAAAAAAAAL/eXkN5DUVAAAAAAAAAAAAAAAAAAAAAAKUU8xrKa0JAAAAAAAAAAAAAAAAAAAAAAPtvDuU1lD9AAAAAAAAAAAAAAAAAAAAAAKSyNpTXUDpAAAAAAAAAAAAAAAAAAAAAAE/1XkN5DTVAAAAAAAAAAAAAAAAAAAAAAPBvDuU1lC9AAAAAAAAAAAAAAAAAAAAAAEf1XkN5DSVAAAAAAAAAAAAAAAAAAAAAAED1XkN5DRVAAAAAAAAAAAA=MAEAAAAAAAABAAAAAAAAAAIAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAHAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAkAAAAAAAAACgAAAAAAAAAKAAAAAAAAAAsAAAAAAAAACwAAAAAAAAAMAAAAAAAAAAwAAAAAAAAADQAAAAAAAAANAAAAAAAAAA4AAAAAAAAADgAAAAAAAAAPAAAAAAAAAA8AAAAAAAAAEAAAAAAAAAAQAAAAAAAAABEAAAAAAAAAEQAAAAAAAAASAAAAAAAAABIAAAAAAAAAEwAAAAAAAAATAAAAAAAAAAAAAAAAAAAAmAAAAAAAAAACAAAAAAAAAAQAAAAAAAAABgAAAAAAAAAIAAAAAAAAAAoAAAAAAAAADAAAAAAAAAAOAAAAAAAAABAAAAAAAAAAEgAAAAAAAAAUAAAAAAAAABYAAAAAAAAAGAAAAAAAAAAaAAAAAAAAABwAAAAAAAAAHgAAAAAAAAAgAAAAAAAAACIAAAAAAAAAJAAAAAAAAAAmAAAAAAAAAA==EwAAAAAAAAADAwMDAwMDAwMDAwMDAwMDAwMD
+ </AppendedData>
 </VTKFile>
diff --git a/Tests/Data/Mechanics/GuentherSalzer/geometry_right.vtu b/Tests/Data/Mechanics/GuentherSalzer/geometry_right.vtu
@@ -1,34 +1,25 @@
 <?xml version="1.0"?>
-<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">
-<!--This file was created by meshio v4.3.5-->
-<UnstructuredGrid>
-<Piece NumberOfPoints="20" NumberOfCells="19">
-<Points>
-<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary">
-AQAAAACAAADgAQAAqwAAAA==eJxjYAABSwcGrAAmHumALn7iU6xzJa8ohnjcD5C4KoZ4+huep6ZT9DHE49+A1JtiiHclmU65HmCFqf4myBx7DPF5hp+kTmU7YYj/+QYy3xVDPHbPyWyN9e4Y4nsaQfZ6YYhLuy/c0/bJB0O8kocX6B5/DHGTIqByqUAMce14kDuDMMSVfDTWL9wTjCEuZRkHdH8ohriQGtCZcWEY4hzCp4D+CoeLAwCUPkQ7
-</DataArray>
-</Points>
-<Cells>
-<DataArray type="Int64" Name="connectivity" format="binary">
-AQAAAACAAAAwAQAANAAAAA==eJxdxTECQCAAAMAQRSL5/18Nme6WC2GYeeHIK2+cOPPOBxc+ufLFNzd+uPPL0/8HqwABfA==
-</DataArray>
-<DataArray type="Int64" Name="offsets" format="binary">
-AQAAAACAAACYAAAANQAAAA==eJwtxdEGgDAAAMDJTDLJZCbJTPr/T+xhdy+3hCk6efXm7N2Hi09XN1++/bh7+PXnH1O4AX0=
-</DataArray>
-<DataArray type="Int64" Name="types" format="binary">
-AQAAAACAAACYAAAADwAAAA==eJxjZoAA5kFGAwASaAA6
-</DataArray>
-</Cells>
-<PointData>
-<DataArray type="UInt64" Name="bulk_node_ids" format="binary">
-AQAAAACAAACgAAAAMgAAAA==eJwtxbcBgCAAADBksSAIFuT/Sx1MlkzhFz178erNybuziw9XN5++fPtx9+vhD0YgARs=
-</DataArray>
-</PointData>
-<CellData>
-<DataArray type="UInt64" Name="bulk_elem_ids" format="binary">
-AQAAAACAAACYAAAAMAAAAA==eJwtxUECQCAAALC+pYgUiv//x6HtsiVM0cmrN2fvPlx8urr58u3H3cOvP/92twTn
-</DataArray>
-</CellData>
-</Piece>
-</UnstructuredGrid>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="20" NumberOfCells="19" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="1" RangeMax="24" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_elem_ids" format="appended" RangeMin="57" RangeMax="75" offset="224" />
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="57" RangeMax="75" offset="440" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="25" RangeMax="103.07764064" offset="656" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1308" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="1724" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1940" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _oAAAAAAAAAABAAAAAAAAAAIAAAAAAAAABwAAAAAAAAAIAAAAAAAAAAkAAAAAAAAACgAAAAAAAAALAAAAAAAAAAwAAAAAAAAADQAAAAAAAAAOAAAAAAAAAA8AAAAAAAAAEAAAAAAAAAARAAAAAAAAABIAAAAAAAAAEwAAAAAAAAAUAAAAAAAAABUAAAAAAAAAFgAAAAAAAAAXAAAAAAAAABgAAAAAAAAAmAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAEAAAAAAAAAAQQAAAAAAAABCAAAAAAAAAEMAAAAAAAAARAAAAAAAAABFAAAAAAAAAEYAAAAAAAAARwAAAAAAAABIAAAAAAAAAEkAAAAAAAAASgAAAAAAAABLAAAAAAAAAA==mAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAEAAAAAAAAAAQQAAAAAAAABCAAAAAAAAAEMAAAAAAAAARAAAAAAAAABFAAAAAAAAAEYAAAAAAAAARwAAAAAAAABIAAAAAAAAAEkAAAAAAAAASgAAAAAAAABLAAAAAAAAAA==4AEAAAAAAAAAAAAAAAA5QAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA5QAAAAAAAAFlAAAAAAAAAAAAAAAAAAAA5QMjyXUN5DRVAAAAAAAAAAAAAAAAAAAA5QF74XUN5DSVAAAAAAAAAAAAAAAAAAAA5QGfsDOU1lC9AAAAAAAAAAAAAAAAAAAA5QF/sXUN5DTVAAAAAAAAAAAAAAAAAAAA5QIpiNZTXUDpAAAAAAAAAAAAAAAAAAAA5QF/ZDOU1lD9AAAAAAAAAAAAAAAAAAAA5QJ4x8hrKa0JAAAAAAAAAAAAAAAAAAAA5QPz2XUN5DUVAAAAAAAAAAAAAAAAAAAA5QF28yWsor0dAAAAAAAAAAAAAAAAAAAA5QLyBNZTXUEpAAAAAAAAAAAAAAAAAAAA5QBtHobyG8kxAAAAAAAAAAAAAAAAAAAA5QHkMDeU1lE9AAAAAAAAAAAAAAAAAAAA5QDRyvIbyGlFAAAAAAAAAAAAAAAAAAAA5QCtf8hrKa1JAAAAAAAAAAAAAAAAAAAA5QCJMKK+hvFNAAAAAAAAAAAAAAAAAAAA5QBo5XkN5DVVAAAAAAAAAAAAAAAAAAAA5QBImlNdQXlZAAAAAAAAAAAAAAAAAAAA5QAgTymsor1dAAAAAAAAAAAA=MAEAAAAAAAAAAAAAAAAAAAIAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAHAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAkAAAAAAAAACgAAAAAAAAAKAAAAAAAAAAsAAAAAAAAACwAAAAAAAAAMAAAAAAAAAAwAAAAAAAAADQAAAAAAAAANAAAAAAAAAA4AAAAAAAAADgAAAAAAAAAPAAAAAAAAAA8AAAAAAAAAEAAAAAAAAAAQAAAAAAAAABEAAAAAAAAAEQAAAAAAAAASAAAAAAAAABIAAAAAAAAAEwAAAAAAAAATAAAAAAAAAAEAAAAAAAAAmAAAAAAAAAACAAAAAAAAAAQAAAAAAAAABgAAAAAAAAAIAAAAAAAAAAoAAAAAAAAADAAAAAAAAAAOAAAAAAAAABAAAAAAAAAAEgAAAAAAAAAUAAAAAAAAABYAAAAAAAAAGAAAAAAAAAAaAAAAAAAAABwAAAAAAAAAHgAAAAAAAAAgAAAAAAAAACIAAAAAAAAAJAAAAAAAAAAmAAAAAAAAAA==EwAAAAAAAAADAwMDAwMDAwMDAwMDAwMDAwMD
+ </AppendedData>
 </VTKFile>