From b65a63987de62413406655479111e4fddd730e75 Mon Sep 17 00:00:00 2001
From: Lars von Wolff <lars.von-wolff@ians.uni-stuttgart.de>
Date: Fri, 11 Jun 2021 10:50:06 +0200
Subject: [PATCH] Increased VelDissipationRate to model a No-Slip boundary
condition
---
src/pntest/pn1p.ini | 2 +-
src/pntest/pntest.cc | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/pntest/pn1p.ini b/src/pntest/pn1p.ini
index e60b946..449a68a 100644
--- a/src/pntest/pn1p.ini
+++ b/src/pntest/pn1p.ini
@@ -7,7 +7,7 @@ rho2 = 1.0
ReactionRate=1#1.5
#ReactionLimiter = 1000
SurfaceTensionScaling=0.1
-VelDissipationRate = 100
+VelDissipationRate = 4000
CurvatureAlpha=0#0.001
[domain]
diff --git a/src/pntest/pntest.cc b/src/pntest/pntest.cc
index 30751c7..42fb795 100644
--- a/src/pntest/pntest.cc
+++ b/src/pntest/pntest.cc
@@ -64,7 +64,7 @@ int main(int argc, char** argv)
//Example of usecase for Pn_1PFlow
std::cout << "Kf = " << pn_1pFlow.getKf() << " phiSolid = " << pn_1pFlow.getPhiSolid() << std::endl;
pn_1pFlow.writeVTK(0.0);
- pn_1pFlow.grow(0.2e-6); //Growth is measured in m
+ pn_1pFlow.grow(2e-6); //Growth is measured in m
pn_1pFlow.writeVTK(1.0);
std::cout << "Kf = " << pn_1pFlow.getKf() << " phiSolid = " << pn_1pFlow.getPhiSolid() << std::endl;
std::cout << "Total Flux = - Kf/viscosity * gradient p" << std::endl;
--
GitLab