Skip to content
Snippets Groups Projects
Commit b65a6398 authored by Lars von Wolff's avatar Lars von Wolff
Browse files

Increased VelDissipationRate to model a No-Slip boundary condition

parent 6f495aef
No related branches found
No related tags found
No related merge requests found
......@@ -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]
......
......@@ -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;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment