From 860e11d4720fee6ff57b71029751d95d000ccc80 Mon Sep 17 00:00:00 2001
From: Simon Scholz <simon.scholz@iws.uni-stuttgart.de>
Date: Mon, 20 Nov 2017 08:27:01 +0100
Subject: [PATCH] [2pncmin][next]add vtkoutputfields.hh do first changes,
unclear about permeability
---
.../2pncmin/implicit/vtkoutputfields.hh | 11 +++++++----
1 file changed, 7 insertions(+), 4 deletions(-)
diff --git a/dumux/porousmediumflow/2pncmin/implicit/vtkoutputfields.hh b/dumux/porousmediumflow/2pncmin/implicit/vtkoutputfields.hh
index b4b61ba576..6ea93e7f3a 100644
--- a/dumux/porousmediumflow/2pncmin/implicit/vtkoutputfields.hh
+++ b/dumux/porousmediumflow/2pncmin/implicit/vtkoutputfields.hh
@@ -41,7 +41,7 @@ class TwoPNCMinVtkOutputFields
using GridView = typename GET_PROP_TYPE(TypeTag, GridView);
static constexpr int numPhases = GET_PROP_VALUE(TypeTag, NumPhases);
- static constexpr int numSPhases = GET_PROP_VALUE(TypeTag, NumSPhases);
+ static constexpr int numSPhases = GET_PROP_VALE(TypteTag, NumSPhases);
static constexpr int numComponents = GET_PROP_VALUE(TypeTag, NumComponents);
static constexpr int dim = GridView::dimension;
@@ -52,11 +52,14 @@ public:
// use default fields from the 2pnc model
TwoPNCVtkOutputFields<TypeTag>::init(vtk);
- //output additional to TwoPNCMin output:
+ //output additional to TwoPNC output:
for (int i = 0; i < numSPhases; ++i)
- {
vtk.addVolumeVariable([i](const VolumeVariables& v){ return v.precipitateVolumeFraction(numPhases + i); },"precipVolFrac_"+ FluidSystem::phaseName(numPhases + i));
- }
+ vtk.addVolumeVariable([](const VolumeVariables& v){ return this->perm_(v.permeability())[0][0]; }, "Kxx"); //TODO: get correct permeability from where? add perm_ function in private?
+ if (dim >= 2)
+ vtk.addVolumeVariable([](const VolumeVariables& v){ return this->perm_(v.permeability())[1][1]; }, "Kyy"); //TODO: get correct permeability from where? add perm_ function in private?
+ if (dim >= 3)
+ vtk.addVolumeVariable([](const VolumeVariables& v){ return this->perm_(v.permeability())[2][2]; }, "Kzz"); //TODO: get correct permeability from where? add perm_ function in private?
}
};
--
GitLab