Skip to content
Snippets Groups Projects
Commit 31fe8476 authored by Gabi Seitz's avatar Gabi Seitz
Browse files

[1pncmin] add test and reference

parent e2de77eb
No related branches found
No related tags found
2 merge requests!617[WIP] Next,!565Transfer 1pncmin to next
Expand all Hide whitespace changes
Inline Side-by-side
Showing
with 573 additions and 130 deletions
test/porousmediumflow/1pnc/implicit/test_box1pnc.input
+ 1
− 1
View file @ 31fe8476
...@@ -8,7 +8,7 @@ UpperRight = 1 1 ...@@ -8,7 +8,7 @@ UpperRight = 1 1
Cells = 20 2 Cells = 20 2
[Problem] [Problem]
Name = outflowbox # name passed to the output routines Name = 1pnctestbox # name passed to the output routines
EnableGravity = 0 # disable gravity EnableGravity = 0 # disable gravity
[Vtk] [Vtk]
......
test/porousmediumflow/1pncmin/implicit/CMakeLists.txt
+ 24
− 45
View file @ 31fe8476
...@@ -2,50 +2,29 @@ ...@@ -2,50 +2,29 @@
add_input_file_links() add_input_file_links()
add_executable(test_box1pncmin test_box1pncmin.cc) # add_executable(test_box1pncmin test_box1pncmin.cc)
add_executable(test_cc1pncmin test_cc1pncmin.cc) # add_executable(test_cc1pncmin test_cc1pncmin.cc)
# add_executable_all(test_box1p2c test_box1p2c.cc) add_dumux_test(test_box1pncmin test_box1pncmin test_box1pncmin.cc
# python ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
# add_executable_all(test_cc1p2c test_cc1p2c.cc) --script fuzzy
--files ${CMAKE_SOURCE_DIR}/test/references/box1pncmin-00042.vtu
# dune_symlink_to_source_files(FILES test_2pnc.input) ${CMAKE_CURRENT_BINARY_DIR}/box1pncmin-00042.vtu
# # dune_symlink_to_source_files(FILES test_box2pnc.input) --command "${CMAKE_CURRENT_BINARY_DIR}/test_box1pncmin"
# dune_symlink_to_source_files(FILES test_box1p2c.input) --zeroThreshold {"temperature":5e-6})
# dune_symlink_to_source_files(FILES test_cc1p2c.input)
add_dumux_test(test_cc1pncmin test_cc1pncmin test_cc1pncmin.cc
# set(BUILD_TYPE Debug) python ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
--script fuzzy
#install sources --files ${CMAKE_SOURCE_DIR}/test/references/cc1pncmin-00042.vtu
# install(FILES ${CMAKE_CURRENT_BINARY_DIR}/cc1pncmin-00042.vtu
# thermochemproblem.hh --command "${CMAKE_CURRENT_BINARY_DIR}/test_cc1pncmin"
# thermochemspatialparams.hh --zeroThreshold {"temperature":5e-6})
# 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)
# #
# #install sources # #install sources
# install(FILES install(FILES
# thermochemproblem.hh thermochempatialparams.hh
# thermochemspatialparams.hh thermochemproblem.hh
# test_box2pnc.cc test_box1pncmin.cc
# test_cc2pnc.cc test_cc1pncmin.cc
# DESTINATION ${CMAKE_INSTALL_INCLUDEDIR}/dumux/test/implicit/2pnc) DESTINATION ${CMAKE_INSTALL_INCLUDEDIR}/dumux/test/implicit/1pncmin)
test/porousmediumflow/1pncmin/implicit/thermochemproblem.hh
+ 53
− 54
View file @ 31fe8476
...@@ -71,8 +71,7 @@ SET_PROP(ThermoChemProblem, FluidSystem) ...@@ -71,8 +71,7 @@ SET_PROP(ThermoChemProblem, FluidSystem)
// Enable velocity output // Enable velocity output
SET_BOOL_PROP(ThermoChemProblem, VtkAddVelocity, false); 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 the spatial parameters
SET_TYPE_PROP(ThermoChemProblem, SpatialParams, ThermoChemSpatialParams<TypeTag>); SET_TYPE_PROP(ThermoChemProblem, SpatialParams, ThermoChemSpatialParams<TypeTag>);
...@@ -211,7 +210,7 @@ public: ...@@ -211,7 +210,7 @@ public:
if(globalPos[0] < eps_ ) if(globalPos[0] < eps_ )
{ {
// values.setDirichlet(pressureIdx); values.setDirichlet(pressureIdx);
values.setDirichlet(firstMoleFracIdx); values.setDirichlet(firstMoleFracIdx);
values.setDirichlet(temperatureIdx); values.setDirichlet(temperatureIdx);
} }
...@@ -234,43 +233,44 @@ public: ...@@ -234,43 +233,44 @@ public:
*/ */
PrimaryVariables dirichletAtPos(const GlobalPosition &globalPos) const PrimaryVariables dirichletAtPos(const GlobalPosition &globalPos) const
{ {
PrimaryVariables priVars(0.0); PrimaryVariables priVars(0.0);
//input parameters //input parameters
// Scalar pIn; Scalar pIn;
Scalar pOut; Scalar pOut;
Scalar tIn; Scalar tIn;
Scalar tOut; Scalar tOut;
Scalar vaporIn; Scalar vaporIn;
// read input parameters // read input parameters
if (isCharge_ == true){ if (isCharge_ == true){
// pIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureIn); pIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureIn);
pOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureOut); pOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureOut);
tIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureIn); tIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureIn);
tOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureOut); tOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureOut);
vaporIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, VaporIn); vaporIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, VaporIn);
}
}
else{ else{
// pIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureIn); pIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureIn);
pOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureOut); pOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureOut);
tIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureIn); tIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureIn);
tOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureOut); tOut = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureOut);
vaporIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, VaporIn); vaporIn = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, VaporIn);
} }
if(globalPos[0] < eps_) if(globalPos[0] < eps_)
{ {
priVars[pressureIdx] = pIn;
priVars[firstMoleFracIdx] = vaporIn; // Saturation outer boundary priVars[firstMoleFracIdx] = vaporIn; // Saturation outer boundary
priVars[temperatureIdx] = tIn; priVars[temperatureIdx] = tIn;
} }
if(globalPos[0] > this->bBoxMax()[0] - eps_) if(globalPos[0] > this->bBoxMax()[0] - eps_)
{ {
priVars[pressureIdx] = pOut; priVars[pressureIdx] = pOut;
priVars[firstMoleFracIdx] = 0.01; // Saturation inner boundary priVars[firstMoleFracIdx] = 0.01; // Saturation inner boundary
priVars[temperatureIdx] = tOut; priVars[temperatureIdx] = tOut;
} }
return priVars; return priVars;
...@@ -299,16 +299,16 @@ public: ...@@ -299,16 +299,16 @@ public:
{ {
PrimaryVariables priVars(0.0); 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; return priVars;
} }
...@@ -334,22 +334,21 @@ public: ...@@ -334,22 +334,21 @@ public:
Scalar CaO2H2Init; Scalar CaO2H2Init;
if (isCharge_ == true){ if (isCharge_ == true){
pInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureInitial); pInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, PressureInitial);
tInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureInitial); tInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, TemperatureInitial);
h2oInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, VaporInitial); h2oInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, VaporInitial);
CaOInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, CaOInitial); CaOInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, CaOInitial);
CaO2H2Init = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, CaO2H2Initial); CaO2H2Init = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Charge, CaO2H2Initial);
} }
else { else {
std::cout << "false " << "\n"; pInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureInitial);
pInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, PressureInitial); tInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureInitial);
tInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, TemperatureInitial); h2oInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, VaporInitial);
h2oInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, VaporInitial); CaOInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, CaOInitial);
CaOInit = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, CaOInitial); CaO2H2Init = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, CaO2H2Initial);
CaO2H2Init = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, Scalar, Discharge, CaO2H2Initial);
} }
std::cout << "CaO2H2Init = " << CaO2H2Init << "\n";
priVars[pressureIdx] = pInit; priVars[pressureIdx] = pInit;
priVars[firstMoleFracIdx] = h2oInit; priVars[firstMoleFracIdx] = h2oInit;
#if NONISOTHERMAL #if NONISOTHERMAL
...@@ -385,10 +384,10 @@ public: ...@@ -385,10 +384,10 @@ public:
if(volVars.moleFraction(phaseIdx, firstMoleFracIdx) > 1e-3) if(volVars.moleFraction(phaseIdx, firstMoleFracIdx) > 1e-3)
moleFractionVapor = volVars.moleFraction(phaseIdx, firstMoleFracIdx) ; moleFractionVapor = volVars.moleFraction(phaseIdx, firstMoleFracIdx) ;
if(volVars.moleFraction(phaseIdx, firstMoleFracIdx) >= 1.0){
moleFractionVapor = 1; if(volVars.moleFraction(phaseIdx, firstMoleFracIdx) >= 1.0) moleFractionVapor = 1;
// std::cout << " test vapor = " << "\n";
}
Scalar vaporPressure = volVars.pressure(phaseIdx) *moleFractionVapor ; Scalar vaporPressure = volVars.pressure(phaseIdx) *moleFractionVapor ;
vaporPressure *= 1.0e-5; vaporPressure *= 1.0e-5;
Scalar pFactor = log(vaporPressure); Scalar pFactor = log(vaporPressure);
...@@ -423,9 +422,9 @@ public: ...@@ -423,9 +422,9 @@ public:
// make sure not more CaO reacts than present // 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(); q = moleFracCaO_sPhase/this->timeManager().timeStepSize();
// std::cout << "q_discharge = " << q << "\n";
} }
source[conti0EqIdx+CaO2H2Idx] = q; source[conti0EqIdx+CaO2H2Idx] = q;
...@@ -449,9 +448,9 @@ public: ...@@ -449,9 +448,9 @@ public:
Scalar q = -rDehydration; 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(); q = moleFracCaO2H2_sPhase/this->timeManager().timeStepSize();
} }
source[conti0EqIdx+CaO2H2Idx] = -q; source[conti0EqIdx+CaO2H2Idx] = -q;
......
test/porousmediumflow/1pncmin/implicit/thermochemspatialparams.hh
+ 3
− 30
View file @ 31fe8476
...@@ -70,13 +70,10 @@ class ThermoChemSpatialParams : public ImplicitSpatialParamsOneP<TypeTag> ...@@ -70,13 +70,10 @@ class ThermoChemSpatialParams : public ImplicitSpatialParamsOneP<TypeTag>
enum { enum {
dim=GridView::dimension, dim=GridView::dimension,
dimWorld=GridView::dimensionworld, dimWorld=GridView::dimensionworld,
// wPhaseIdx = FluidSystem::wPhaseIdx
}; };
using GlobalPosition = Dune::FieldVector<CoordScalar, dimWorld>; using GlobalPosition = Dune::FieldVector<CoordScalar, dimWorld>;
using Tensor = Dune::FieldMatrix<CoordScalar, dimWorld, 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 SubControlVolume = typename GET_PROP_TYPE(TypeTag, SubControlVolume);
using Element = typename GridView::template Codim<0>::Entity; using Element = typename GridView::template Codim<0>::Entity;
...@@ -94,12 +91,6 @@ public: ...@@ -94,12 +91,6 @@ public:
ThermoChemSpatialParams(const Problem& problem, const GridView &gridView) ThermoChemSpatialParams(const Problem& problem, const GridView &gridView)
: ParentType(problem, 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 //thermal conductivity of CaO
lambdaSolid_ = 0.4; //[W/(m*K)] Nagel et al [2013b] lambdaSolid_ = 0.4; //[W/(m*K)] Nagel et al [2013b]
...@@ -119,19 +110,6 @@ public: ...@@ -119,19 +110,6 @@ public:
permLaw_.init(*this); 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 * \brief Define the initial permeability \f$[m^2]\f$ distribution
* *
...@@ -237,9 +215,6 @@ public: ...@@ -237,9 +215,6 @@ public:
const ElementSolutionVector& elemSol) const const ElementSolutionVector& elemSol) const
{ {
return 790; 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: ...@@ -272,11 +247,9 @@ public:
{ return lambdaSolid_; } { return lambdaSolid_; }
private: private:
// DimMatrix K_;
// Scalar porosity_; Scalar eps_;
Scalar eps_; Scalar lambdaSolid_;
// MaterialLawParams materialParams_;
Scalar lambdaSolid_;
bool isCharge_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, bool, Problem, IsCharge); bool isCharge_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, bool, Problem, IsCharge);
PorosityLaw poroLaw_; PorosityLaw poroLaw_;
PermeabilityLaw permLaw_; PermeabilityLaw permLaw_;
......
test/references/box1pncmin-00042.vtu 0 → 100644
+ 273
− 0
View file @ 31fe8476
This diff is collapsed.
test/references/cc1pncmin-00042.vtu 0 → 100644
+ 219
− 0
View file @ 31fe8476
<?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>
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment