From fc46149a0eb9ad4f0357f06dd8159124cee8c43f Mon Sep 17 00:00:00 2001
From: Lars von Wolff <lars.von-wolff@ians.uni-stuttgart.de>
Date: Tue, 26 Oct 2021 14:54:14 +0200
Subject: [PATCH] Cleaned up for shared use
---
dune/phasefield/porenetwork/pn_2pflow.hh | 16 ++++++++--------
1 file changed, 8 insertions(+), 8 deletions(-)
diff --git a/dune/phasefield/porenetwork/pn_2pflow.hh b/dune/phasefield/porenetwork/pn_2pflow.hh
index 9af50e0..c1dfca8 100644
--- a/dune/phasefield/porenetwork/pn_2pflow.hh
+++ b/dune/phasefield/porenetwork/pn_2pflow.hh
@@ -349,7 +349,7 @@ private:
void findStationaryState()
{
- writeVTK(0);
+ //writeVTK(0);
Dune::FieldVector<double,1> integral;
param.reaction.reactionRate = 0;
int iter = 0;
@@ -358,6 +358,7 @@ private:
do
{
+ std::cout << "Find Stationary State, Iteration " << iter << " of (maximum) " << param.PCS_Search_MaxCHNSIter << std::endl;
iter++;
applyChPressureSaturation(1);
applyNsPressureSaturation(1);
@@ -440,14 +441,12 @@ public:
initializeGridOperators();
-/* applyNs();
- updatePhiSolid();
+ findStationaryState();
+ calculateKf();
+ updateSaturation();
initialPhiSolid = phiSolid;
- updateKf();
-*/
- //findStationaryState();
- updateKf();
- std::cout << kf_11 << ", "
+
+ std::cout << "Pn_2PFlow Initialized, initial Kf (viscosity ratio possibly not right)" << kf_11 << ", "
<< kf_12 << ", "
<< kf_21 << ", "
<< kf_22 << std::endl;
@@ -640,6 +639,7 @@ void calculateKf(RF ViscosityRatio = 0)
void findSaturation(RF Curvature)
{
+ std::cout << "Pn_2PFlow, function find Saturation not yet implemented" << std::endl;
//TODO!!!!
}
--
GitLab