Commit 2398b52a authored by jouke's avatar jouke
Browse files

Correction of the Hagen Pouiseulle calculation of the viscosity 

coeficient.
parent 067a11c0

OpenSimRoot/src/modules/WaterModule/DoussanModel.cpp

+1 −1
Original line number Diff line number Diff line
@@ -144,7 +144,7 @@ void HydraulicConductance::calculate(const Time &t, double &c) { 
			double numberOfXylemVessels, diameterOfXylemVessels;

			pNumberOfXylemVessels->get(t - rootSegmentStartTime, numberOfXylemVessels);

			pDiameterOfXylemVessels->get(t - rootSegmentStartTime, diameterOfXylemVessels);

			c = (M_PI*60*60*24*numberOfXylemVessels*pow(diameterOfXylemVessels/2, 4))/(8*8.9*pow(10, -6));

			c = (M_PI*numberOfXylemVessels*pow(diameterOfXylemVessels/2, 4))/(8*0.769);

		}

		//by length, the longer the root the lower the lateral con

		if (s > 1e-9) {