Die Infiltration kann mit dem Modell nach Horton berechnet werden mit:
in mm/h mit: Infiltrationsrate,
konstante Rate nach der Anfangszeit,
der anfänglichen Infiltrationsrate und
einem Abklingfaktor.
from pylab import * def Horton(t): fu = 3.0 f0 = 12.0 k = 0.5 H = fu+(f0-fu)*exp(-k*t) return H t = arange(0.0, 60.0, 0.5) plot(t, Horton(t)) ytext = ylabel('Infiltrationsrate (mm/min)') xtext = xlabel('Zeit (Min.)') show()
Die kumulative Infiltration beträgt:
in mm.
from pylab import * def HortonInt(t): fu = 3.0 f0 = 12.0 k = 0.05 Hi = fu*t+(f0-fu)/k*(1-exp(-k*t)) return Hi t = arange(0.0, 60.0, 0.5) plot(t, HortonInt(t)) ytext = ylabel('Infiltration (mm)') xtext = xlabel('Zeit (Min.)') show()