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
und der hydraulischen Leitfähigkeit
in
.
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()