Philipp

Die Infiltrationsmenge nach einer Zeit wird nach Philipp berechnet mit:

from pylab import * def PhilippInt(t): S = 85.0 K = 5.0 Ip = S*t**(0.5)+K*t return Ip t = arange(0.0, 60.0, 0.5) plot(t, PhilippInt(t)) ytext = ylabel('Infiltrationsmenge (mm)') xtext = xlabel('Zeit (Min.)') show()

mit der Sorptivität in mm mm/und der hydraulischen Leitfähigkeit in mm/h.

Aus dieser Gleichung kann die Gleichung für die Infiltrationsrate durch Ableiten erstellt werden:

from pylab import * def Philipp(t): S = 85.0 K = 5.0 ip = 0.5*S*t**(-0.5)+K return ip t = arange(0.0, 60.0, 0.5) plot(t, Philipp(t)) ytext = ylabel('Infiltrationsmenge (mm)') xtext = xlabel('Zeit (Min.)') show()