From 31fe8476998fe3ae91d3c31cb5357c9a714111af Mon Sep 17 00:00:00 2001
From: Gabriele Seitz <gabriele.seitz@iws.uni-stuttgart.de>
Date: Wed, 18 Oct 2017 14:21:34 +0200
Subject: [PATCH] [1pncmin] add test and reference
---
.../1pnc/implicit/test_box1pnc.input | 2 +-
.../1pncmin/implicit/CMakeLists.txt | 69 ++---
.../1pncmin/implicit/thermochemproblem.hh | 107 ++++---
.../implicit/thermochemspatialparams.hh | 33 +--
test/references/box1pncmin-00042.vtu | 273 ++++++++++++++++++
test/references/cc1pncmin-00042.vtu | 219 ++++++++++++++
6 files changed, 573 insertions(+), 130 deletions(-)
create mode 100644 test/references/box1pncmin-00042.vtu
create mode 100644 test/references/cc1pncmin-00042.vtu
diff --git a/test/porousmediumflow/1pnc/implicit/test_box1pnc.input b/test/porousmediumflow/1pnc/implicit/test_box1pnc.input
index 85a6863ffc..56c07d2043 100644
--- a/test/porousmediumflow/1pnc/implicit/test_box1pnc.input
+++ b/test/porousmediumflow/1pnc/implicit/test_box1pnc.input
@@ -8,7 +8,7 @@ UpperRight = 1 1
Cells = 20 2
[Problem]
-Name = outflowbox # name passed to the output routines
+Name = 1pnctestbox # name passed to the output routines
EnableGravity = 0 # disable gravity
[Vtk]
diff --git a/test/porousmediumflow/1pncmin/implicit/CMakeLists.txt b/test/porousmediumflow/1pncmin/implicit/CMakeLists.txt
index 91d858c04b..813b2ec775 100644
--- a/test/porousmediumflow/1pncmin/implicit/CMakeLists.txt
+++ b/test/porousmediumflow/1pncmin/implicit/CMakeLists.txt
@@ -2,50 +2,29 @@
add_input_file_links()
-add_executable(test_box1pncmin test_box1pncmin.cc)
-add_executable(test_cc1pncmin test_cc1pncmin.cc)
-
-# add_executable_all(test_box1p2c test_box1p2c.cc)
-#
-# add_executable_all(test_cc1p2c test_cc1p2c.cc)
-
-# dune_symlink_to_source_files(FILES test_2pnc.input)
-# # dune_symlink_to_source_files(FILES test_box2pnc.input)
-# dune_symlink_to_source_files(FILES test_box1p2c.input)
-# dune_symlink_to_source_files(FILES test_cc1p2c.input)
-
-# set(BUILD_TYPE Debug)
-
-#install sources
-# install(FILES
-# thermochemproblem.hh
-# thermochemspatialparams.hh
-# test_box2pnc.cc
-# test_cc2pnc.cc
-# DESTINATION ${CMAKE_INSTALL_INCLUDEDIR}/dumux-devel/appl/thermochemistry)
-
-
-# # isothermal tests
-# add_dumux_test(test_box2pnc test_box2pnc test_box2pnc.cc
-# python ${dumux_INCLUDE_DIRS}/bin/testing/runtest.py
-# --script fuzzy
-# --files ${CMAKE_SOURCE_DIR}/test/references/fuelcell2pncbox-reference.vtu
-# ${CMAKE_CURRENT_BINARY_DIR}/fuelcell-box-00016.vtu
-# --command "${CMAKE_CURRENT_BINARY_DIR}/test_box2pnc -ParameterFile test_2pnc.input -Problem.Name fuelcell-box")
-#
-# add_dumux_test(test_cc2pnc test_cc2pnc test_cc2pnc.cc
-# python ${dumux_INCLUDE_DIRS}/bin/testing/runtest.py
-# --script fuzzy
-# --files ${CMAKE_SOURCE_DIR}/test/references/fuelcell2pnccc-reference.vtu
-# ${CMAKE_CURRENT_BINARY_DIR}/fuelcell-cc-00016.vtu
-# --command "${CMAKE_CURRENT_BINARY_DIR}/test_cc2pnc -ParameterFile test_2pnc.input -Problem.Name fuelcell-cc")
-#
-# dune_symlink_to_source_files(FILES test_2pnc.input)
+# add_executable(test_box1pncmin test_box1pncmin.cc)
+# add_executable(test_cc1pncmin test_cc1pncmin.cc)
+
+add_dumux_test(test_box1pncmin test_box1pncmin test_box1pncmin.cc
+ python ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+ --script fuzzy
+ --files ${CMAKE_SOURCE_DIR}/test/references/box1pncmin-00042.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/box1pncmin-00042.vtu
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_box1pncmin"
+ --zeroThreshold {"temperature":5e-6})
+
+add_dumux_test(test_cc1pncmin test_cc1pncmin test_cc1pncmin.cc
+ python ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+ --script fuzzy
+ --files ${CMAKE_SOURCE_DIR}/test/references/cc1pncmin-00042.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/cc1pncmin-00042.vtu
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_cc1pncmin"
+ --zeroThreshold {"temperature":5e-6})
#
# #install sources
-# install(FILES
-# thermochemproblem.hh
-# thermochemspatialparams.hh
-# test_box2pnc.cc
-# test_cc2pnc.cc
-# DESTINATION ${CMAKE_INSTALL_INCLUDEDIR}/dumux/test/implicit/2pnc)
+install(FILES
+thermochempatialparams.hh
+thermochemproblem.hh
+test_box1pncmin.cc
+test_cc1pncmin.cc
+DESTINATION ${CMAKE_INSTALL_INCLUDEDIR}/dumux/test/implicit/1pncmin)
diff --git a/test/porousmediumflow/1pncmin/implicit/thermochemproblem.hh b/test/porousmediumflow/1pncmin/implicit/thermochemproblem.hh
index 21847e375f..cdea15dd73 100644
--- a/test/porousmediumflow/1pncmin/implicit/thermochemproblem.hh
+++ b/test/porousmediumflow/1pncmin/implicit/thermochemproblem.hh
@@ -71,8 +71,7 @@ SET_PROP(ThermoChemProblem, FluidSystem)
// Enable velocity output
SET_BOOL_PROP(ThermoChemProblem, VtkAddVelocity, false);
-SET_TYPE_PROP(ThermoChemProblem, LinearSolver, UMFPackBackend<TypeTag>);
-
+// SET_TYPE_PROP(ThermoChemProblem, LinearSolver, UMFPackBackend<TypeTag>);
// Set the spatial parameters
SET_TYPE_PROP(ThermoChemProblem, SpatialParams, ThermoChemSpatialParams<TypeTag>);
@@ -211,7 +210,7 @@ public:
if(globalPos[0] < eps_ )
{
-// values.setDirichlet(pressureIdx);
+ values.setDirichlet(pressureIdx);
values.setDirichlet(firstMoleFracIdx);
values.setDirichlet(temperatureIdx);
}
@@ -234,43 +233,44 @@ public:
*/
PrimaryVariables dirichletAtPos(const GlobalPosition &globalPos) const
{
- PrimaryVariables priVars(0.0);
+ PrimaryVariables priVars(0.0);
- //input parameters
-// Scalar pIn;
- Scalar pOut;
- Scalar tIn;
- Scalar tOut;
- Scalar vaporIn;
+ //input parameters
+ Scalar pIn;
+ Scalar pOut;
+ Scalar tIn;
+ Scalar tOut;
+ Scalar vaporIn;
// read input parameters
- if (isCharge_ == true){
-// pIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureIn);
- pOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureOut);
- tIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureIn);
- tOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureOut);
- vaporIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, VaporIn);
-
- }
+ if (isCharge_ == true){
+ pIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureIn);
+ pOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureOut);
+ tIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureIn);
+ tOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureOut);
+ vaporIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, VaporIn);
+ }
+
else{
-// pIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureIn);
- pOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureOut);
- tIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureIn);
- tOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureOut);
- vaporIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, VaporIn);
- }
+ pIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureIn);
+ pOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureOut);
+ tIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureIn);
+ tOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureOut);
+ vaporIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, VaporIn);
+ }
if(globalPos[0] < eps_)
{
+ priVars[pressureIdx] = pIn;
priVars[firstMoleFracIdx] = vaporIn; // Saturation outer boundary
priVars[temperatureIdx] = tIn;
}
+
if(globalPos[0] > this->bBoxMax()[0] - eps_)
{
priVars[pressureIdx] = pOut;
priVars[firstMoleFracIdx] = 0.01; // Saturation inner boundary
priVars[temperatureIdx] = tOut;
-
}
return priVars;
@@ -299,16 +299,16 @@ public:
{
PrimaryVariables priVars(0.0);
- if(globalPos[0] < eps_)
- {
+// if(globalPos[0] < eps_)
+// {
+//
+// if (isCharge_ == true){
+// priVars[pressureIdx] = -GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, InFlow); //[mol/s] gas inflow; negative sign: inflow
+// }
+// else
+// priVars[pressureIdx] = -GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, InFlow); //[mol/s] gas inflow
+// }
- if (isCharge_ == true){
- priVars[pressureIdx] = -GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, InFlow); //[mol/s] gas inflow; negative sign: inflow
- }
- else
- priVars[pressureIdx] = -GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, InFlow); //[mol/s] gas inflow
- }
-// std::cout << " test neumann " << "\n";
return priVars;
}
@@ -334,22 +334,21 @@ public:
Scalar CaO2H2Init;
if (isCharge_ == true){
- pInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureInitial);
- tInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureInitial);
- h2oInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, VaporInitial);
- CaOInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, CaOInitial);
- CaO2H2Init = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, CaO2H2Initial);
+ pInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureInitial);
+ tInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureInitial);
+ h2oInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, VaporInitial);
+ CaOInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, CaOInitial);
+ CaO2H2Init = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, CaO2H2Initial);
}
else {
- std::cout << "false " << "\n";
- pInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureInitial);
- tInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureInitial);
- h2oInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, VaporInitial);
- CaOInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, CaOInitial);
- CaO2H2Init = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, CaO2H2Initial);
+ pInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureInitial);
+ tInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureInitial);
+ h2oInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, VaporInitial);
+ CaOInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, CaOInitial);
+ CaO2H2Init = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, CaO2H2Initial);
}
- std::cout << "CaO2H2Init = " << CaO2H2Init << "\n";
+
priVars[pressureIdx] = pInit;
priVars[firstMoleFracIdx] = h2oInit;
#if NONISOTHERMAL
@@ -385,10 +384,10 @@ public:
if(volVars.moleFraction(phaseIdx, firstMoleFracIdx) > 1e-3)
moleFractionVapor = volVars.moleFraction(phaseIdx, firstMoleFracIdx) ;
- if(volVars.moleFraction(phaseIdx, firstMoleFracIdx) >= 1.0){
- moleFractionVapor = 1;
-// std::cout << " test vapor = " << "\n";
- }
+
+ if(volVars.moleFraction(phaseIdx, firstMoleFracIdx) >= 1.0) moleFractionVapor = 1;
+
+
Scalar vaporPressure = volVars.pressure(phaseIdx) *moleFractionVapor ;
vaporPressure *= 1.0e-5;
Scalar pFactor = log(vaporPressure);
@@ -423,9 +422,9 @@ public:
// make sure not more CaO reacts than present
- if (- q*this->timeManager().timeStepSize() + moleFracCaO_sPhase* volVars.molarDensity(cPhaseIdx) < 0 + eps_){
+ if (- q*this->timeManager().timeStepSize() + moleFracCaO_sPhase* volVars.molarDensity(cPhaseIdx) < 0 + eps_)
+ {
q = moleFracCaO_sPhase/this->timeManager().timeStepSize();
-// std::cout << "q_discharge = " << q << "\n";
}
source[conti0EqIdx+CaO2H2Idx] = q;
@@ -449,9 +448,9 @@ public:
Scalar q = -rDehydration;
- if (- q*this->timeManager().timeStepSize() + moleFracCaO2H2_sPhase*volVars.molarDensity(hPhaseIdx) < 0){
+ if (- q*this->timeManager().timeStepSize() + moleFracCaO2H2_sPhase*volVars.molarDensity(hPhaseIdx) < 0)
+ {
q = moleFracCaO2H2_sPhase/this->timeManager().timeStepSize();
-
}
source[conti0EqIdx+CaO2H2Idx] = -q;
diff --git a/test/porousmediumflow/1pncmin/implicit/thermochemspatialparams.hh b/test/porousmediumflow/1pncmin/implicit/thermochemspatialparams.hh
index 36ee0d5715..80cc28c118 100644
--- a/test/porousmediumflow/1pncmin/implicit/thermochemspatialparams.hh
+++ b/test/porousmediumflow/1pncmin/implicit/thermochemspatialparams.hh
@@ -70,13 +70,10 @@ class ThermoChemSpatialParams : public ImplicitSpatialParamsOneP<TypeTag>
enum {
dim=GridView::dimension,
dimWorld=GridView::dimensionworld,
-
-// wPhaseIdx = FluidSystem::wPhaseIdx
};
using GlobalPosition = Dune::FieldVector<CoordScalar, dimWorld>;
using Tensor = Dune::FieldMatrix<CoordScalar, dimWorld, dimWorld>;
-// using FVElementGeometry = typename GET_PROP_TYPE(TypeTag, FVElementGeometry);
using SubControlVolume = typename GET_PROP_TYPE(TypeTag, SubControlVolume);
using Element = typename GridView::template Codim<0>::Entity;
@@ -94,12 +91,6 @@ public:
ThermoChemSpatialParams(const Problem& problem, const GridView &gridView)
: ParentType(problem, gridView)
{
- // intrinsic permeabilities
-// K_[0][0] = 5e-12;
-// K_[1][1] = 5e-12;
-// K_[0][0] = 2.23e-14;
-// K_[1][1] = 2.23e-14;
-
//thermal conductivity of CaO
lambdaSolid_ = 0.4; //[W/(m*K)] Nagel et al [2013b]
@@ -119,19 +110,6 @@ public:
permLaw_.init(*this);
}
- /*! Old
- * \brief Apply the intrinsic permeability tensor to a pressure
- * potential gradient.
- *
- * \param element The current finite element
- * \param fvGeometry The current finite volume geometry of the element
- * \param scvIdx The index of the sub-control volume
- */
-// const DimMatrix intrinsicPermeability(const Element &element,
-// const FVElementGeometry &fvGeometry,
-// const int scvIdx) const
-// { return K_; }
-//
/*!
* \brief Define the initial permeability \f$[m^2]\f$ distribution
*
@@ -237,9 +215,6 @@ public:
const ElementSolutionVector& elemSol) const
{
return 790;
-// 42 // specific heat capacity of CaO [J / (kg K)]
-// * 3370 // density of CaO [kg/m^3]
-// * (1 - porosity(element, fvGeometry, scvIdx)); // for CaO only!!
}
/*!
@@ -272,11 +247,9 @@ public:
{ return lambdaSolid_; }
private:
-// DimMatrix K_;
-// Scalar porosity_;
- Scalar eps_;
-// MaterialLawParams materialParams_;
- Scalar lambdaSolid_;
+
+ Scalar eps_;
+ Scalar lambdaSolid_;
bool isCharge_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, bool, Problem, IsCharge);
PorosityLaw poroLaw_;
PermeabilityLaw permLaw_;
diff --git a/test/references/box1pncmin-00042.vtu b/test/references/box1pncmin-00042.vtu
new file mode 100644
index 0000000000..972d3f0946
--- /dev/null
+++ b/test/references/box1pncmin-00042.vtu
@@ -0,0 +1,273 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="80" NumberOfPoints="123">
+ <PointData Scalars="temperature">
+ <DataArray type="Float32" Name="temperature" NumberOfComponents="1" format="ascii">
+ 773.15 748.138 773.15 748.138 732.976 732.976 731.372 731.372 731.432 731.432 731.682 731.682
+ 731.957 731.957 732.236 732.236 732.518 732.518 732.8 732.8 733.082 733.082 733.363 733.363
+ 733.639 733.639 733.908 733.908 734.167 734.167 734.411 734.411 734.632 734.632 734.823 734.823
+ 734.971 734.971 735.062 735.062 735.074 735.074 734.972 734.972 734.673 734.673 734.151 734.151
+ 733.369 733.369 732.295 732.295 730.895 730.895 729.151 729.151 727.06 727.06 724.624 724.624
+ 721.845 721.845 718.727 718.727 715.275 715.275 711.496 711.496 707.396 707.396 702.976 702.976
+ 698.225 698.225 693.11 693.11 687.567 687.567 681.467 681.467 673.15 673.15 773.15 748.138
+ 732.976 731.372 731.432 731.682 731.957 732.236 732.518 732.8 733.082 733.363 733.639 733.908
+ 734.167 734.411 734.632 734.823 734.971 735.062 735.074 734.972 734.673 734.151 733.369 732.295
+ 730.895 729.151 727.06 724.624 721.845 718.727 715.275 711.496 707.396 702.976 698.225 693.11
+ 687.567 681.467 673.15
+ </DataArray>
+ <DataArray type="Float32" Name="pressure" NumberOfComponents="1" format="ascii">
+ 102000 101999 102000 101999 101993 101993 101961 101961 101914 101914 101864 101864
+ 101813 101813 101762 101762 101710 101710 101658 101658 101606 101606 101553 101553
+ 101500 101500 101446 101446 101392 101392 101337 101337 101282 101282 101227 101227
+ 101171 101171 101115 101115 101059 101059 101002 101002 100945 100945 100889 100889
+ 100832 100832 100775 100775 100718 100718 100661 100661 100604 100604 100548 100548
+ 100491 100491 100436 100436 100380 100380 100325 100325 100270 100270 100217 100217
+ 100164 100164 100113 100113 100064 100064 100017 100017 100000 100000 102000 101999
+ 101993 101961 101914 101864 101813 101762 101710 101658 101606 101553 101500 101446
+ 101392 101337 101282 101227 101171 101115 101059 101002 100945 100889 100832 100775
+ 100718 100661 100604 100548 100491 100436 100380 100325 100270 100217 100164 100113
+ 100064 100017 100000
+ </DataArray>
+ <DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
+ 0.457777 0.441312 0.457777 0.441312 0.426818 0.426818 0.423712 0.423712 0.422608 0.422608 0.421757 0.421757
+ 0.420929 0.420929 0.420097 0.420097 0.419258 0.419258 0.418413 0.418413 0.417566 0.417566 0.416719 0.416719
+ 0.415876 0.415876 0.415044 0.415044 0.41423 0.41423 0.413446 0.413446 0.412705 0.412705 0.412028 0.412028
+ 0.411439 0.411439 0.410971 0.410971 0.410667 0.410667 0.41056 0.41056 0.41061 0.41061 0.410836 0.410836
+ 0.411252 0.411252 0.411857 0.411857 0.412646 0.412646 0.413615 0.413615 0.414779 0.414779 0.416167 0.416167
+ 0.417819 0.417819 0.419787 0.419787 0.422157 0.422157 0.425065 0.425065 0.428737 0.428737 0.433544 0.433544
+ 0.440106 0.440106 0.449467 0.449467 0.463401 0.463401 0.484953 0.484953 0.515472 0.515472 0.457777 0.441312
+ 0.426818 0.423712 0.422608 0.421757 0.420929 0.420097 0.419258 0.418413 0.417566 0.416719 0.415876 0.415044
+ 0.41423 0.413446 0.412705 0.412028 0.411439 0.410971 0.410667 0.41056 0.41061 0.410836 0.411252 0.411857
+ 0.412646 0.413615 0.414779 0.416167 0.417819 0.419787 0.422157 0.425065 0.428737 0.433544 0.440106 0.449467
+ 0.463401 0.484953 0.515472
+ </DataArray>
+ <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
+ 0.797746 0.760039 0.797746 0.760039 0.654209 0.654209 0.613832 0.613832 0.608052 0.608052 0.607362 0.607362
+ 0.607301 0.607301 0.607307 0.607307 0.607313 0.607313 0.607313 0.607313 0.607307 0.607307 0.607296 0.607296
+ 0.607281 0.607281 0.607261 0.607261 0.607237 0.607237 0.607208 0.607208 0.607174 0.607174 0.607132 0.607132
+ 0.607081 0.607081 0.607016 0.607016 0.606918 0.606918 0.606757 0.606757 0.606541 0.606541 0.606263 0.606263
+ 0.605907 0.605907 0.605473 0.605473 0.605005 0.605005 0.604598 0.604598 0.604336 0.604336 0.604205 0.604205
+ 0.604133 0.604133 0.604085 0.604085 0.604051 0.604051 0.60403 0.60403 0.604018 0.604018 0.604012 0.604012
+ 0.60401 0.60401 0.604014 0.604014 0.604039 0.604039 0.653058 0.653058 0.751006 0.751006 0.797746 0.760039
+ 0.654209 0.613832 0.608052 0.607362 0.607301 0.607307 0.607313 0.607313 0.607307 0.607296 0.607281 0.607261
+ 0.607237 0.607208 0.607174 0.607132 0.607081 0.607016 0.606918 0.606757 0.606541 0.606263 0.605907 0.605473
+ 0.605005 0.604598 0.604336 0.604205 0.604133 0.604085 0.604051 0.60403 0.604018 0.604012 0.60401 0.604014
+ 0.604039 0.653058 0.751006
+ </DataArray>
+ <DataArray type="Float32" Name="x_Air" NumberOfComponents="1" format="ascii">
+ 0.99 0.813249 0.99 0.813249 0.684603 0.684603 0.663337 0.663337 0.658583 0.658583 0.655867 0.655867
+ 0.653364 0.653364 0.650855 0.650855 0.648319 0.648319 0.645758 0.645758 0.643182 0.643182 0.640603 0.640603
+ 0.63804 0.63804 0.635516 0.635516 0.633065 0.633065 0.630732 0.630732 0.628576 0.628576 0.626676 0.626676
+ 0.625138 0.625138 0.624098 0.624098 0.623722 0.623722 0.624088 0.624088 0.624716 0.624716 0.625629 0.625629
+ 0.626787 0.626787 0.628078 0.628078 0.629361 0.629361 0.630547 0.630547 0.631684 0.631684 0.632929 0.632929
+ 0.634475 0.634475 0.63659 0.63659 0.639708 0.639708 0.644545 0.644545 0.652256 0.652256 0.664719 0.664719
+ 0.685029 0.685029 0.71833 0.71833 0.773238 0.773238 0.864316 0.864316 0.99 0.99 0.99 0.813249
+ 0.684603 0.663337 0.658583 0.655867 0.653364 0.650855 0.648319 0.645758 0.643182 0.640603 0.63804 0.635516
+ 0.633065 0.630732 0.628576 0.626676 0.625138 0.624098 0.623722 0.624088 0.624716 0.625629 0.626787 0.628078
+ 0.629361 0.630547 0.631684 0.632929 0.634475 0.63659 0.639708 0.644545 0.652256 0.664719 0.685029 0.71833
+ 0.773238 0.864316 0.99
+ </DataArray>
+ <DataArray type="Float32" Name="x_H2O" NumberOfComponents="1" format="ascii">
+ 0.01 0.186751 0.01 0.186751 0.315397 0.315397 0.336663 0.336663 0.341417 0.341417 0.344133 0.344133
+ 0.346636 0.346636 0.349145 0.349145 0.351681 0.351681 0.354242 0.354242 0.356818 0.356818 0.359397 0.359397
+ 0.36196 0.36196 0.364484 0.364484 0.366935 0.366935 0.369268 0.369268 0.371424 0.371424 0.373324 0.373324
+ 0.374862 0.374862 0.375902 0.375902 0.376278 0.376278 0.375912 0.375912 0.375284 0.375284 0.374371 0.374371
+ 0.373213 0.373213 0.371922 0.371922 0.370639 0.370639 0.369453 0.369453 0.368316 0.368316 0.367071 0.367071
+ 0.365525 0.365525 0.36341 0.36341 0.360292 0.360292 0.355455 0.355455 0.347744 0.347744 0.335281 0.335281
+ 0.314971 0.314971 0.28167 0.28167 0.226762 0.226762 0.135684 0.135684 0.01 0.01 0.01 0.186751
+ 0.315397 0.336663 0.341417 0.344133 0.346636 0.349145 0.351681 0.354242 0.356818 0.359397 0.36196 0.364484
+ 0.366935 0.369268 0.371424 0.373324 0.374862 0.375902 0.376278 0.375912 0.375284 0.374371 0.373213 0.371922
+ 0.370639 0.369453 0.368316 0.367071 0.365525 0.36341 0.360292 0.355455 0.347744 0.335281 0.314971 0.28167
+ 0.226762 0.135684 0.01
+ </DataArray>
+ <DataArray type="Float32" Name="X_Air" NumberOfComponents="1" format="ascii">
+ 0.993761 0.875098 0.993761 0.875098 0.777396 0.777396 0.760195 0.760195 0.756306 0.756306 0.754077 0.754077
+ 0.752018 0.752018 0.74995 0.74995 0.747854 0.747854 0.745734 0.745734 0.743597 0.743597 0.741452 0.741452
+ 0.739315 0.739315 0.737206 0.737206 0.735154 0.735154 0.733197 0.733197 0.731384 0.731384 0.729784 0.729784
+ 0.728487 0.728487 0.727608 0.727608 0.72729 0.72729 0.7276 0.7276 0.72813 0.72813 0.728901 0.728901
+ 0.729877 0.729877 0.730965 0.730965 0.732044 0.732044 0.733041 0.733041 0.733996 0.733996 0.73504 0.73504
+ 0.736335 0.736335 0.738104 0.738104 0.740706 0.740706 0.744728 0.744728 0.751105 0.751105 0.761321 0.761321
+ 0.777737 0.777737 0.80404 0.80404 0.845825 0.845825 0.911101 0.911101 0.993761 0.993761 0.993761 0.875098
+ 0.777396 0.760195 0.756306 0.754077 0.752018 0.74995 0.747854 0.745734 0.743597 0.741452 0.739315 0.737206
+ 0.735154 0.733197 0.731384 0.729784 0.728487 0.727608 0.72729 0.7276 0.72813 0.728901 0.729877 0.730965
+ 0.732044 0.733041 0.733996 0.73504 0.736335 0.738104 0.740706 0.744728 0.751105 0.761321 0.777737 0.80404
+ 0.845825 0.911101 0.993761
+ </DataArray>
+ <DataArray type="Float32" Name="X_H2O" NumberOfComponents="1" format="ascii">
+ 0.00623908 0.124902 0.00623908 0.124902 0.222604 0.222604 0.239805 0.239805 0.243694 0.243694 0.245923 0.245923
+ 0.247982 0.247982 0.25005 0.25005 0.252146 0.252146 0.254266 0.254266 0.256403 0.256403 0.258548 0.258548
+ 0.260685 0.260685 0.262794 0.262794 0.264846 0.264846 0.266803 0.266803 0.268616 0.268616 0.270216 0.270216
+ 0.271513 0.271513 0.272392 0.272392 0.27271 0.27271 0.2724 0.2724 0.27187 0.27187 0.271099 0.271099
+ 0.270123 0.270123 0.269035 0.269035 0.267956 0.267956 0.266959 0.266959 0.266004 0.266004 0.26496 0.26496
+ 0.263665 0.263665 0.261896 0.261896 0.259294 0.259294 0.255272 0.255272 0.248895 0.248895 0.238679 0.238679
+ 0.222263 0.222263 0.19596 0.19596 0.154175 0.154175 0.0888991 0.0888991 0.00623908 0.00623908 0.00623908 0.124902
+ 0.222604 0.239805 0.243694 0.245923 0.247982 0.25005 0.252146 0.254266 0.256403 0.258548 0.260685 0.262794
+ 0.264846 0.266803 0.268616 0.270216 0.271513 0.272392 0.27271 0.2724 0.27187 0.271099 0.270123 0.269035
+ 0.267956 0.266959 0.266004 0.26496 0.263665 0.261896 0.259294 0.255272 0.248895 0.238679 0.222263 0.19596
+ 0.154175 0.0888991 0.00623908
+ </DataArray>
+ <DataArray type="Float32" Name="Kxx" NumberOfComponents="1" format="ascii">
+ 4.848e-12 2.97842e-12 4.848e-12 2.97842e-12 9.14707e-13 9.14707e-13 6.05835e-13 6.05835e-13 5.71642e-13 5.71642e-13 5.677e-13 5.677e-13
+ 5.67351e-13 5.67351e-13 5.67386e-13 5.67386e-13 5.67422e-13 5.67422e-13 5.67421e-13 5.67421e-13 5.67387e-13 5.67387e-13 5.67325e-13 5.67325e-13
+ 5.67237e-13 5.67237e-13 5.67125e-13 5.67125e-13 5.66989e-13 5.66989e-13 5.66825e-13 5.66825e-13 5.66629e-13 5.66629e-13 5.66393e-13 5.66393e-13
+ 5.66103e-13 5.66103e-13 5.65731e-13 5.65731e-13 5.65176e-13 5.65176e-13 5.64265e-13 5.64265e-13 5.63041e-13 5.63041e-13 5.61476e-13 5.61476e-13
+ 5.59474e-13 5.59474e-13 5.57047e-13 5.57047e-13 5.54437e-13 5.54437e-13 5.5218e-13 5.5218e-13 5.50735e-13 5.50735e-13 5.50011e-13 5.50011e-13
+ 5.49616e-13 5.49616e-13 5.49349e-13 5.49349e-13 5.49164e-13 5.49164e-13 5.4905e-13 5.4905e-13 5.48984e-13 5.48984e-13 5.4895e-13 5.4895e-13
+ 5.48939e-13 5.48939e-13 5.48962e-13 5.48962e-13 5.49096e-13 5.49096e-13 9.03863e-13 9.03863e-13 2.66878e-12 2.66878e-12 4.848e-12 2.97842e-12
+ 9.14707e-13 6.05835e-13 5.71642e-13 5.677e-13 5.67351e-13 5.67386e-13 5.67422e-13 5.67421e-13 5.67387e-13 5.67325e-13 5.67237e-13 5.67125e-13
+ 5.66989e-13 5.66825e-13 5.66629e-13 5.66393e-13 5.66103e-13 5.65731e-13 5.65176e-13 5.64265e-13 5.63041e-13 5.61476e-13 5.59474e-13 5.57047e-13
+ 5.54437e-13 5.5218e-13 5.50735e-13 5.50011e-13 5.49616e-13 5.49349e-13 5.49164e-13 5.4905e-13 5.48984e-13 5.4895e-13 5.48939e-13 5.48962e-13
+ 5.49096e-13 9.03863e-13 2.66878e-12
+ </DataArray>
+ <DataArray type="Float32" Name="Kyy" NumberOfComponents="1" format="ascii">
+ 4.848e-12 2.97842e-12 4.848e-12 2.97842e-12 9.14707e-13 9.14707e-13 6.05835e-13 6.05835e-13 5.71642e-13 5.71642e-13 5.677e-13 5.677e-13
+ 5.67351e-13 5.67351e-13 5.67386e-13 5.67386e-13 5.67422e-13 5.67422e-13 5.67421e-13 5.67421e-13 5.67387e-13 5.67387e-13 5.67325e-13 5.67325e-13
+ 5.67237e-13 5.67237e-13 5.67125e-13 5.67125e-13 5.66989e-13 5.66989e-13 5.66825e-13 5.66825e-13 5.66629e-13 5.66629e-13 5.66393e-13 5.66393e-13
+ 5.66103e-13 5.66103e-13 5.65731e-13 5.65731e-13 5.65176e-13 5.65176e-13 5.64265e-13 5.64265e-13 5.63041e-13 5.63041e-13 5.61476e-13 5.61476e-13
+ 5.59474e-13 5.59474e-13 5.57047e-13 5.57047e-13 5.54437e-13 5.54437e-13 5.5218e-13 5.5218e-13 5.50735e-13 5.50735e-13 5.50011e-13 5.50011e-13
+ 5.49616e-13 5.49616e-13 5.49349e-13 5.49349e-13 5.49164e-13 5.49164e-13 5.4905e-13 5.4905e-13 5.48984e-13 5.48984e-13 5.4895e-13 5.4895e-13
+ 5.48939e-13 5.48939e-13 5.48962e-13 5.48962e-13 5.49096e-13 5.49096e-13 9.03863e-13 9.03863e-13 2.66878e-12 2.66878e-12 4.848e-12 2.97842e-12
+ 9.14707e-13 6.05835e-13 5.71642e-13 5.677e-13 5.67351e-13 5.67386e-13 5.67422e-13 5.67421e-13 5.67387e-13 5.67325e-13 5.67237e-13 5.67125e-13
+ 5.66989e-13 5.66825e-13 5.66629e-13 5.66393e-13 5.66103e-13 5.65731e-13 5.65176e-13 5.64265e-13 5.63041e-13 5.61476e-13 5.59474e-13 5.57047e-13
+ 5.54437e-13 5.5218e-13 5.50735e-13 5.50011e-13 5.49616e-13 5.49349e-13 5.49164e-13 5.4905e-13 5.48984e-13 5.4895e-13 5.48939e-13 5.48962e-13
+ 5.49096e-13 9.03863e-13 2.66878e-12
+ </DataArray>
+ <DataArray type="Float32" Name="permeabilityFactor" NumberOfComponents="1" format="ascii">
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0
+ </DataArray>
+ <DataArray type="Float32" Name="precipitateVolumeFraction_CaO" NumberOfComponents="1" format="ascii">
+ 0.200016 0.190883 0.200016 0.190883 0.0643481 0.0643481 0.0119243 0.0119243 0.0043627 0.0043627 0.00345262 0.00345262
+ 0.00336744 0.00336744 0.00337309 0.00337309 0.00338107 0.00338107 0.00338252 0.00338252 0.00337775 0.00337775 0.00336773 0.00336773
+ 0.0033531 0.0033531 0.00333415 0.00333415 0.00331085 0.00331085 0.00328283 0.00328283 0.00324932 0.00324932 0.00320904 0.00320904
+ 0.00316 0.00316 0.00309875 0.00309875 0.00301546 0.00301546 0.0028656 0.0028656 0.00265894 0.00265894 0.00239522 0.00239522
+ 0.00205291 0.00205291 0.00162532 0.00162532 0.00113847 0.00113847 0.000686056 0.000686056 0.000381022 0.000381022 0.000228894 0.000228894
+ 0.000149176 0.000149176 9.5797e-05 9.5797e-05 5.80078e-05 5.80078e-05 3.4229e-05 3.4229e-05 2.05082e-05 2.05082e-05 1.30132e-05 1.30132e-05
+ 9.96843e-06 9.96843e-06 1.1315e-05 1.1315e-05 2.37672e-05 2.37672e-05 8.59409e-05 8.59409e-05 0.200013 0.200013 0.200016 0.190883
+ 0.0643481 0.0119243 0.0043627 0.00345262 0.00336744 0.00337309 0.00338107 0.00338252 0.00337775 0.00336773 0.0033531 0.00333415
+ 0.00331085 0.00328283 0.00324932 0.00320904 0.00316 0.00309875 0.00301546 0.0028656 0.00265894 0.00239522 0.00205291 0.00162532
+ 0.00113847 0.000686056 0.000381022 0.000228894 0.000149176 9.5797e-05 5.80078e-05 3.4229e-05 2.05082e-05 1.30132e-05 9.96843e-06 1.1315e-05
+ 2.37672e-05 8.59409e-05 0.200013
+ </DataArray>
+ <DataArray type="Float32" Name="precipitateVolumeFraction_CaOH2" NumberOfComponents="1" format="ascii">
+ 4.35044e-11 0.0180816 4.35044e-11 0.0180816 0.268601 0.268601 0.372392 0.372392 0.387363 0.387363 0.389164 0.389164
+ 0.389333 0.389333 0.389322 0.389322 0.389306 0.389306 0.389303 0.389303 0.389313 0.389313 0.389332 0.389332
+ 0.389361 0.389361 0.389399 0.389399 0.389445 0.389445 0.389501 0.389501 0.389567 0.389567 0.389647 0.389647
+ 0.389744 0.389744 0.389865 0.389865 0.39003 0.39003 0.390327 0.390327 0.390736 0.390736 0.391258 0.391258
+ 0.391936 0.391936 0.392782 0.392782 0.393746 0.393746 0.394642 0.394642 0.395246 0.395246 0.395547 0.395547
+ 0.395705 0.395705 0.39581 0.39581 0.395885 0.395885 0.395932 0.395932 0.395959 0.395959 0.395974 0.395974
+ 0.39598 0.39598 0.395978 0.395978 0.395953 0.395953 0.39583 0.39583 7.18725e-06 7.18725e-06 4.35044e-11 0.0180816
+ 0.268601 0.372392 0.387363 0.389164 0.389333 0.389322 0.389306 0.389303 0.389313 0.389332 0.389361 0.389399
+ 0.389445 0.389501 0.389567 0.389647 0.389744 0.389865 0.39003 0.390327 0.390736 0.391258 0.391936 0.392782
+ 0.393746 0.394642 0.395246 0.395547 0.395705 0.39581 0.395885 0.395932 0.395959 0.395974 0.39598 0.395978
+ 0.395953 0.39583 7.18725e-06
+ </DataArray>
+ </PointData>
+ <CellData Scalars="process rank">
+ <DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+ 0 0 0 0.002 0 0 0 0.005 0 0.002 0.005 0
+ 0.004 0 0 0.004 0.005 0 0.006 0 0 0.006 0.005 0
+ 0.008 0 0 0.008 0.005 0 0.01 0 0 0.01 0.005 0
+ 0.012 0 0 0.012 0.005 0 0.014 0 0 0.014 0.005 0
+ 0.016 0 0 0.016 0.005 0 0.018 0 0 0.018 0.005 0
+ 0.02 0 0 0.02 0.005 0 0.022 0 0 0.022 0.005 0
+ 0.024 0 0 0.024 0.005 0 0.026 0 0 0.026 0.005 0
+ 0.028 0 0 0.028 0.005 0 0.03 0 0 0.03 0.005 0
+ 0.032 0 0 0.032 0.005 0 0.034 0 0 0.034 0.005 0
+ 0.036 0 0 0.036 0.005 0 0.038 0 0 0.038 0.005 0
+ 0.04 0 0 0.04 0.005 0 0.042 0 0 0.042 0.005 0
+ 0.044 0 0 0.044 0.005 0 0.046 0 0 0.046 0.005 0
+ 0.048 0 0 0.048 0.005 0 0.05 0 0 0.05 0.005 0
+ 0.052 0 0 0.052 0.005 0 0.054 0 0 0.054 0.005 0
+ 0.056 0 0 0.056 0.005 0 0.058 0 0 0.058 0.005 0
+ 0.06 0 0 0.06 0.005 0 0.062 0 0 0.062 0.005 0
+ 0.064 0 0 0.064 0.005 0 0.066 0 0 0.066 0.005 0
+ 0.068 0 0 0.068 0.005 0 0.07 0 0 0.07 0.005 0
+ 0.072 0 0 0.072 0.005 0 0.074 0 0 0.074 0.005 0
+ 0.076 0 0 0.076 0.005 0 0.078 0 0 0.078 0.005 0
+ 0.08 0 0 0.08 0.005 0 0 0.01 0 0.002 0.01 0
+ 0.004 0.01 0 0.006 0.01 0 0.008 0.01 0 0.01 0.01 0
+ 0.012 0.01 0 0.014 0.01 0 0.016 0.01 0 0.018 0.01 0
+ 0.02 0.01 0 0.022 0.01 0 0.024 0.01 0 0.026 0.01 0
+ 0.028 0.01 0 0.03 0.01 0 0.032 0.01 0 0.034 0.01 0
+ 0.036 0.01 0 0.038 0.01 0 0.04 0.01 0 0.042 0.01 0
+ 0.044 0.01 0 0.046 0.01 0 0.048 0.01 0 0.05 0.01 0
+ 0.052 0.01 0 0.054 0.01 0 0.056 0.01 0 0.058 0.01 0
+ 0.06 0.01 0 0.062 0.01 0 0.064 0.01 0 0.066 0.01 0
+ 0.068 0.01 0 0.07 0.01 0 0.072 0.01 0 0.074 0.01 0
+ 0.076 0.01 0 0.078 0.01 0 0.08 0.01 0
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 1 4 5 3 4 6 7 5
+ 6 8 9 7 8 10 11 9 10 12 13 11
+ 12 14 15 13 14 16 17 15 16 18 19 17
+ 18 20 21 19 20 22 23 21 22 24 25 23
+ 24 26 27 25 26 28 29 27 28 30 31 29
+ 30 32 33 31 32 34 35 33 34 36 37 35
+ 36 38 39 37 38 40 41 39 40 42 43 41
+ 42 44 45 43 44 46 47 45 46 48 49 47
+ 48 50 51 49 50 52 53 51 52 54 55 53
+ 54 56 57 55 56 58 59 57 58 60 61 59
+ 60 62 63 61 62 64 65 63 64 66 67 65
+ 66 68 69 67 68 70 71 69 70 72 73 71
+ 72 74 75 73 74 76 77 75 76 78 79 77
+ 78 80 81 79 2 3 83 82 3 5 84 83
+ 5 7 85 84 7 9 86 85 9 11 87 86
+ 11 13 88 87 13 15 89 88 15 17 90 89
+ 17 19 91 90 19 21 92 91 21 23 93 92
+ 23 25 94 93 25 27 95 94 27 29 96 95
+ 29 31 97 96 31 33 98 97 33 35 99 98
+ 35 37 100 99 37 39 101 100 39 41 102 101
+ 41 43 103 102 43 45 104 103 45 47 105 104
+ 47 49 106 105 49 51 107 106 51 53 108 107
+ 53 55 109 108 55 57 110 109 57 59 111 110
+ 59 61 112 111 61 63 113 112 63 65 114 113
+ 65 67 115 114 67 69 116 115 69 71 117 116
+ 71 73 118 117 73 75 119 118 75 77 120 119
+ 77 79 121 120 79 81 122 121
+ </DataArray>
+ <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+ 4 8 12 16 20 24 28 32 36 40 44 48
+ 52 56 60 64 68 72 76 80 84 88 92 96
+ 100 104 108 112 116 120 124 128 132 136 140 144
+ 148 152 156 160 164 168 172 176 180 184 188 192
+ 196 200 204 208 212 216 220 224 228 232 236 240
+ 244 248 252 256 260 264 268 272 276 280 284 288
+ 292 296 300 304 308 312 316 320
+ </DataArray>
+ <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
diff --git a/test/references/cc1pncmin-00042.vtu b/test/references/cc1pncmin-00042.vtu
new file mode 100644
index 0000000000..7dc3f2edbb
--- /dev/null
+++ b/test/references/cc1pncmin-00042.vtu
@@ -0,0 +1,219 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="80" NumberOfPoints="123">
+ <CellData Scalars="temperature">
+ <DataArray type="Float32" Name="temperature" NumberOfComponents="1" format="ascii">
+ 773.15 732.152 695.627 685.878 685.266 686.19 687.337 688.478 689.571 690.602 691.565 692.456
+ 693.268 693.998 694.642 695.195 695.653 696.01 696.262 696.403 696.425 696.321 696.082 695.697
+ 695.156 694.443 693.541 692.433 691.094 689.498 687.613 685.401 682.817 679.802 676.263 672.217
+ 667.688 662.84 660.326 673.15 773.15 732.152 695.627 685.878 685.266 686.19 687.337 688.478
+ 689.571 690.602 691.565 692.456 693.268 693.998 694.642 695.195 695.653 696.01 696.262 696.403
+ 696.425 696.321 696.082 695.697 695.156 694.443 693.541 692.433 691.094 689.498 687.613 685.401
+ 682.817 679.802 676.263 672.217 667.688 662.84 660.326 673.15
+ </DataArray>
+ <DataArray type="Float32" Name="pressure" NumberOfComponents="1" format="ascii">
+ 102000 102001 102002 101987 101950 101903 101854 101805 101755 101704 101652 101600
+ 101548 101494 101440 101385 101329 101273 101216 101158 101100 101041 100981 100921
+ 100860 100799 100738 100676 100614 100551 100489 100426 100364 100302 100240 100178
+ 100117 100061 100020 100000 102000 102001 102002 101987 101950 101903 101854 101805
+ 101755 101704 101652 101600 101548 101494 101440 101385 101329 101273 101216 101158
+ 101100 101041 100981 100921 100860 100799 100738 100676 100614 100551 100489 100426
+ 100364 100302 100240 100178 100117 100061 100020 100000
+ </DataArray>
+ <DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
+ 0.457777 0.477048 0.495121 0.498935 0.497984 0.49629 0.494508 0.492753 0.491054 0.489423 0.48787 0.486402
+ 0.48503 0.48376 0.4826 0.481559 0.480643 0.479863 0.479225 0.478741 0.478419 0.478271 0.478309 0.478546
+ 0.478996 0.479676 0.480604 0.4818 0.483286 0.485089 0.487236 0.489761 0.4927 0.496096 0.500008 0.504448
+ 0.509425 0.514897 0.519572 0.515472 0.457777 0.477048 0.495121 0.498935 0.497984 0.49629 0.494508 0.492753
+ 0.491054 0.489423 0.48787 0.486402 0.48503 0.48376 0.4826 0.481559 0.480643 0.479863 0.479225 0.478741
+ 0.478419 0.478271 0.478309 0.478546 0.478996 0.479676 0.480604 0.4818 0.483286 0.485089 0.487236 0.489761
+ 0.4927 0.496096 0.500008 0.504448 0.509425 0.514897 0.519572 0.515472
+ </DataArray>
+ <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
+ 0.799984 0.798476 0.734853 0.643529 0.617236 0.612379 0.611589 0.611481 0.611468 0.611459 0.611442 0.611421
+ 0.611398 0.611377 0.61136 0.611347 0.61134 0.611341 0.61135 0.611368 0.611396 0.611434 0.611482 0.611539
+ 0.611604 0.611676 0.611752 0.611826 0.61189 0.611932 0.611931 0.611849 0.611618 0.611067 0.610318 0.609744
+ 0.609625 0.622453 0.708308 0.799829 0.799984 0.798476 0.734853 0.643529 0.617236 0.612379 0.611589 0.611481
+ 0.611468 0.611459 0.611442 0.611421 0.611398 0.611377 0.61136 0.611347 0.61134 0.611341 0.61135 0.611368
+ 0.611396 0.611434 0.611482 0.611539 0.611604 0.611676 0.611752 0.611826 0.61189 0.611932 0.611931 0.611849
+ 0.611618 0.611067 0.610318 0.609744 0.609625 0.622453 0.708308 0.799829
+ </DataArray>
+ <DataArray type="Float32" Name="x_Air" NumberOfComponents="1" format="ascii">
+ 0.99 0.955329 0.919215 0.90314 0.896951 0.892888 0.889217 0.885655 0.882181 0.878817 0.875594 0.872542
+ 0.869692 0.867075 0.864719 0.862655 0.860913 0.859522 0.858512 0.857913 0.857755 0.858069 0.858884 0.860233
+ 0.862144 0.864646 0.867769 0.871538 0.875974 0.881094 0.886906 0.89341 0.900585 0.908385 0.916693 0.925546
+ 0.934981 0.945218 0.959856 0.99 0.99 0.955329 0.919215 0.90314 0.896951 0.892888 0.889217 0.885655
+ 0.882181 0.878817 0.875594 0.872542 0.869692 0.867075 0.864719 0.862655 0.860913 0.859522 0.858512 0.857913
+ 0.857755 0.858069 0.858884 0.860233 0.862144 0.864646 0.867769 0.871538 0.875974 0.881094 0.886906 0.89341
+ 0.900585 0.908385 0.916693 0.925546 0.934981 0.945218 0.959856 0.99
+ </DataArray>
+ <DataArray type="Float32" Name="x_H2O" NumberOfComponents="1" format="ascii">
+ 0.01 0.0446713 0.0807855 0.0968603 0.103049 0.107112 0.110783 0.114345 0.117819 0.121183 0.124406 0.127458
+ 0.130308 0.132925 0.135281 0.137345 0.139087 0.140478 0.141488 0.142087 0.142245 0.141931 0.141116 0.139767
+ 0.137856 0.135354 0.132231 0.128462 0.124026 0.118906 0.113094 0.10659 0.0994153 0.091615 0.0833068 0.0744539
+ 0.0650189 0.0547821 0.0401436 0.01 0.01 0.0446713 0.0807855 0.0968603 0.103049 0.107112 0.110783 0.114345
+ 0.117819 0.121183 0.124406 0.127458 0.130308 0.132925 0.135281 0.137345 0.139087 0.140478 0.141488 0.142087
+ 0.142245 0.141931 0.141116 0.139767 0.137856 0.135354 0.132231 0.128462 0.124026 0.118906 0.113094 0.10659
+ 0.0994153 0.091615 0.0833068 0.0744539 0.0650189 0.0547821 0.0401436 0.01
+ </DataArray>
+ <DataArray type="Float32" Name="X_Air" NumberOfComponents="1" format="ascii">
+ 0.993761 0.971757 0.948204 0.937506 0.933351 0.930612 0.92813 0.925714 0.923352 0.921059 0.918855 0.916764
+ 0.914806 0.913004 0.911379 0.909953 0.908747 0.907784 0.907083 0.906667 0.906558 0.906776 0.907342 0.908276
+ 0.909599 0.911329 0.913483 0.916074 0.919115 0.922612 0.926564 0.930964 0.935793 0.941012 0.946535 0.952382
+ 0.958568 0.96523 0.974664 0.993761 0.993761 0.971757 0.948204 0.937506 0.933351 0.930612 0.92813 0.925714
+ 0.923352 0.921059 0.918855 0.916764 0.914806 0.913004 0.911379 0.909953 0.908747 0.907784 0.907083 0.906667
+ 0.906558 0.906776 0.907342 0.908276 0.909599 0.911329 0.913483 0.916074 0.919115 0.922612 0.926564 0.930964
+ 0.935793 0.941012 0.946535 0.952382 0.958568 0.96523 0.974664 0.993761
+ </DataArray>
+ <DataArray type="Float32" Name="X_H2O" NumberOfComponents="1" format="ascii">
+ 0.00623908 0.0282428 0.0517955 0.0624941 0.0666492 0.0693879 0.0718702 0.0742855 0.0766478 0.0789412 0.0811445 0.0832359
+ 0.0851936 0.0869957 0.0886206 0.0900467 0.0912525 0.0922164 0.0929169 0.0933326 0.0934422 0.0932244 0.0926584 0.0917237
+ 0.0904007 0.0886707 0.0865172 0.0839256 0.0808848 0.0773884 0.0734361 0.0690358 0.064207 0.0589882 0.0534647 0.0476183
+ 0.0414318 0.0347705 0.025336 0.00623908 0.00623908 0.0282428 0.0517955 0.0624941 0.0666492 0.0693879 0.0718702 0.0742855
+ 0.0766478 0.0789412 0.0811445 0.0832359 0.0851936 0.0869957 0.0886206 0.0900467 0.0912525 0.0922164 0.0929169 0.0933326
+ 0.0934422 0.0932244 0.0926584 0.0917237 0.0904007 0.0886707 0.0865172 0.0839256 0.0808848 0.0773884 0.0734361 0.0690358
+ 0.064207 0.0589882 0.0534647 0.0476183 0.0414318 0.0347705 0.025336 0.00623908
+ </DataArray>
+ <DataArray type="Float32" Name="Kxx" NumberOfComponents="1" format="ascii">
+ 4.99888e-12 4.89654e-12 2.20489e-12 8.19248e-13 6.26977e-13 5.97043e-13 5.92321e-13 5.91677e-13 5.91602e-13 5.91544e-13 5.91445e-13 5.91318e-13
+ 5.91185e-13 5.9106e-13 5.90957e-13 5.90881e-13 5.90839e-13 5.90842e-13 5.90897e-13 5.91006e-13 5.91172e-13 5.91398e-13 5.91681e-13 5.92021e-13
+ 5.92411e-13 5.9284e-13 5.93291e-13 5.93732e-13 5.94117e-13 5.94368e-13 5.94357e-13 5.93869e-13 5.92492e-13 5.89217e-13 5.84799e-13 5.81439e-13
+ 5.80744e-13 6.60902e-13 1.63146e-12 4.98826e-12 4.99888e-12 4.89654e-12 2.20489e-12 8.19248e-13 6.26977e-13 5.97043e-13 5.92321e-13 5.91677e-13
+ 5.91602e-13 5.91544e-13 5.91445e-13 5.91318e-13 5.91185e-13 5.9106e-13 5.90957e-13 5.90881e-13 5.90839e-13 5.90842e-13 5.90897e-13 5.91006e-13
+ 5.91172e-13 5.91398e-13 5.91681e-13 5.92021e-13 5.92411e-13 5.9284e-13 5.93291e-13 5.93732e-13 5.94117e-13 5.94368e-13 5.94357e-13 5.93869e-13
+ 5.92492e-13 5.89217e-13 5.84799e-13 5.81439e-13 5.80744e-13 6.60902e-13 1.63146e-12 4.98826e-12
+ </DataArray>
+ <DataArray type="Float32" Name="Kyy" NumberOfComponents="1" format="ascii">
+ 4.99888e-12 4.89654e-12 2.20489e-12 8.19248e-13 6.26977e-13 5.97043e-13 5.92321e-13 5.91677e-13 5.91602e-13 5.91544e-13 5.91445e-13 5.91318e-13
+ 5.91185e-13 5.9106e-13 5.90957e-13 5.90881e-13 5.90839e-13 5.90842e-13 5.90897e-13 5.91006e-13 5.91172e-13 5.91398e-13 5.91681e-13 5.92021e-13
+ 5.92411e-13 5.9284e-13 5.93291e-13 5.93732e-13 5.94117e-13 5.94368e-13 5.94357e-13 5.93869e-13 5.92492e-13 5.89217e-13 5.84799e-13 5.81439e-13
+ 5.80744e-13 6.60902e-13 1.63146e-12 4.98826e-12 4.99888e-12 4.89654e-12 2.20489e-12 8.19248e-13 6.26977e-13 5.97043e-13 5.92321e-13 5.91677e-13
+ 5.91602e-13 5.91544e-13 5.91445e-13 5.91318e-13 5.91185e-13 5.9106e-13 5.90957e-13 5.90881e-13 5.90839e-13 5.90842e-13 5.90897e-13 5.91006e-13
+ 5.91172e-13 5.91398e-13 5.91681e-13 5.92021e-13 5.92411e-13 5.9284e-13 5.93291e-13 5.93732e-13 5.94117e-13 5.94368e-13 5.94357e-13 5.93869e-13
+ 5.92492e-13 5.89217e-13 5.84799e-13 5.81439e-13 5.80744e-13 6.60902e-13 1.63146e-12 4.98826e-12
+ </DataArray>
+ <DataArray type="Float32" Name="permeabilityFactor" NumberOfComponents="1" format="ascii">
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0
+ </DataArray>
+ <DataArray type="Float32" Name="precipitateVolumeFraction_CaO" NumberOfComponents="1" format="ascii">
+ 0.200016 0.198477 0.133545 0.0403422 0.0135083 0.0085513 0.0077455 0.00763502 0.00762206 0.00761215 0.0075952 0.00757343
+ 0.00755052 0.00752905 0.00751135 0.0074983 0.00749108 0.00749165 0.00750101 0.00751977 0.00754837 0.00758704 0.00763572 0.00769398
+ 0.00776079 0.00783433 0.00791144 0.00798699 0.00805276 0.00809552 0.0080937 0.00801041 0.00777469 0.00721207 0.00644775 0.00586238
+ 0.00574085 0.0188325 0.106454 0.199858 0.200016 0.198477 0.133545 0.0403422 0.0135083 0.0085513 0.0077455 0.00763502
+ 0.00762206 0.00761215 0.0075952 0.00757343 0.00755052 0.00752905 0.00751135 0.0074983 0.00749108 0.00749165 0.00750101 0.00751977
+ 0.00754837 0.00758704 0.00763572 0.00769398 0.00776079 0.00783433 0.00791144 0.00798699 0.00805276 0.00809552 0.0080937 0.00801041
+ 0.00777469 0.00721207 0.00644775 0.00586238 0.00574085 0.0188325 0.106454 0.199858
+ </DataArray>
+ <DataArray type="Float32" Name="precipitateVolumeFraction_CaOH2" NumberOfComponents="1" format="ascii">
+ 8.84454e-08 0.00304719 0.131602 0.316129 0.369256 0.37907 0.380665 0.380884 0.38091 0.380929 0.380963 0.381006
+ 0.381051 0.381094 0.381129 0.381155 0.381169 0.381168 0.381149 0.381112 0.381055 0.380979 0.380883 0.380767
+ 0.380635 0.380489 0.380337 0.380187 0.380057 0.379972 0.379976 0.380141 0.380607 0.381721 0.383235 0.384393
+ 0.384634 0.358715 0.185239 0.000312714 8.84454e-08 0.00304719 0.131602 0.316129 0.369256 0.37907 0.380665 0.380884
+ 0.38091 0.380929 0.380963 0.381006 0.381051 0.381094 0.381129 0.381155 0.381169 0.381168 0.381149 0.381112
+ 0.381055 0.380979 0.380883 0.380767 0.380635 0.380489 0.380337 0.380187 0.380057 0.379972 0.379976 0.380141
+ 0.380607 0.381721 0.383235 0.384393 0.384634 0.358715 0.185239 0.000312714
+ </DataArray>
+ <DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+ 0 0 0 0.002 0 0 0 0.005 0 0.002 0.005 0
+ 0.004 0 0 0.004 0.005 0 0.006 0 0 0.006 0.005 0
+ 0.008 0 0 0.008 0.005 0 0.01 0 0 0.01 0.005 0
+ 0.012 0 0 0.012 0.005 0 0.014 0 0 0.014 0.005 0
+ 0.016 0 0 0.016 0.005 0 0.018 0 0 0.018 0.005 0
+ 0.02 0 0 0.02 0.005 0 0.022 0 0 0.022 0.005 0
+ 0.024 0 0 0.024 0.005 0 0.026 0 0 0.026 0.005 0
+ 0.028 0 0 0.028 0.005 0 0.03 0 0 0.03 0.005 0
+ 0.032 0 0 0.032 0.005 0 0.034 0 0 0.034 0.005 0
+ 0.036 0 0 0.036 0.005 0 0.038 0 0 0.038 0.005 0
+ 0.04 0 0 0.04 0.005 0 0.042 0 0 0.042 0.005 0
+ 0.044 0 0 0.044 0.005 0 0.046 0 0 0.046 0.005 0
+ 0.048 0 0 0.048 0.005 0 0.05 0 0 0.05 0.005 0
+ 0.052 0 0 0.052 0.005 0 0.054 0 0 0.054 0.005 0
+ 0.056 0 0 0.056 0.005 0 0.058 0 0 0.058 0.005 0
+ 0.06 0 0 0.06 0.005 0 0.062 0 0 0.062 0.005 0
+ 0.064 0 0 0.064 0.005 0 0.066 0 0 0.066 0.005 0
+ 0.068 0 0 0.068 0.005 0 0.07 0 0 0.07 0.005 0
+ 0.072 0 0 0.072 0.005 0 0.074 0 0 0.074 0.005 0
+ 0.076 0 0 0.076 0.005 0 0.078 0 0 0.078 0.005 0
+ 0.08 0 0 0.08 0.005 0 0 0.01 0 0.002 0.01 0
+ 0.004 0.01 0 0.006 0.01 0 0.008 0.01 0 0.01 0.01 0
+ 0.012 0.01 0 0.014 0.01 0 0.016 0.01 0 0.018 0.01 0
+ 0.02 0.01 0 0.022 0.01 0 0.024 0.01 0 0.026 0.01 0
+ 0.028 0.01 0 0.03 0.01 0 0.032 0.01 0 0.034 0.01 0
+ 0.036 0.01 0 0.038 0.01 0 0.04 0.01 0 0.042 0.01 0
+ 0.044 0.01 0 0.046 0.01 0 0.048 0.01 0 0.05 0.01 0
+ 0.052 0.01 0 0.054 0.01 0 0.056 0.01 0 0.058 0.01 0
+ 0.06 0.01 0 0.062 0.01 0 0.064 0.01 0 0.066 0.01 0
+ 0.068 0.01 0 0.07 0.01 0 0.072 0.01 0 0.074 0.01 0
+ 0.076 0.01 0 0.078 0.01 0 0.08 0.01 0
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 1 4 5 3 4 6 7 5
+ 6 8 9 7 8 10 11 9 10 12 13 11
+ 12 14 15 13 14 16 17 15 16 18 19 17
+ 18 20 21 19 20 22 23 21 22 24 25 23
+ 24 26 27 25 26 28 29 27 28 30 31 29
+ 30 32 33 31 32 34 35 33 34 36 37 35
+ 36 38 39 37 38 40 41 39 40 42 43 41
+ 42 44 45 43 44 46 47 45 46 48 49 47
+ 48 50 51 49 50 52 53 51 52 54 55 53
+ 54 56 57 55 56 58 59 57 58 60 61 59
+ 60 62 63 61 62 64 65 63 64 66 67 65
+ 66 68 69 67 68 70 71 69 70 72 73 71
+ 72 74 75 73 74 76 77 75 76 78 79 77
+ 78 80 81 79 2 3 83 82 3 5 84 83
+ 5 7 85 84 7 9 86 85 9 11 87 86
+ 11 13 88 87 13 15 89 88 15 17 90 89
+ 17 19 91 90 19 21 92 91 21 23 93 92
+ 23 25 94 93 25 27 95 94 27 29 96 95
+ 29 31 97 96 31 33 98 97 33 35 99 98
+ 35 37 100 99 37 39 101 100 39 41 102 101
+ 41 43 103 102 43 45 104 103 45 47 105 104
+ 47 49 106 105 49 51 107 106 51 53 108 107
+ 53 55 109 108 55 57 110 109 57 59 111 110
+ 59 61 112 111 61 63 113 112 63 65 114 113
+ 65 67 115 114 67 69 116 115 69 71 117 116
+ 71 73 118 117 73 75 119 118 75 77 120 119
+ 77 79 121 120 79 81 122 121
+ </DataArray>
+ <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+ 4 8 12 16 20 24 28 32 36 40 44 48
+ 52 56 60 64 68 72 76 80 84 88 92 96
+ 100 104 108 112 116 120 124 128 132 136 140 144
+ 148 152 156 160 164 168 172 176 180 184 188 192
+ 196 200 204 208 212 216 220 224 228 232 236 240
+ 244 248 252 256 260 264 268 272 276 280 284 288
+ 292 296 300 304 308 312 316 320
+ </DataArray>
+ <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
--
GitLab