Environmental Engineering Reference
In-Depth Information
% -----------------------------------------------------------------
function
[c,f,s] = unsatpde(x,t,u,DuDx,s1,s2,qr,f,a,n,ks)
[q,k,c] = sedprop(u,qr,f,a,n,ks);
f = k.*DuDx-k;
s = 0;
% -----------------------------------------------------------------
function
u0 = unsatic(x,s1,s2,qr,f,a,n,ks)
u0 = -200+x;
if x < 10 u0 = -0.5; end
% -----------------------------------------------------------------
function
[pl,ql,pr,qr] = unsatbc(xl,ul,xr,ur,t,s1,s2,qr,f,a,n,ks)
pl = s1;
ql = 1;
pr = ur(1)-s2;
qr = 0;
%------------------- soil hydraulic properties --------------------
function
[q,k,c] = sedprop(u,qr,f,a,n,ks)
m = 1-1/n;
if
u >= 0
c=1e-20;
k=ks;
q=f;
else
q=qr+(f-qr)*(1+(-a*u)^n)^-m;
c=((f-qr)*n*m*a*(-a*u)^(n-1))/((1+(-a*u)^n)^(m+1))+1.e-20;
k=ks*((q-qr)/(f-qr))^0.5*(1-(1-((q-qr)/(f-qr))^(1/m))^m)^2;
end
52
0
0
Pressure
Head [L]
Hydraulic
Head [L]
Water
Content [%]
-20
-20
50
-40
-40
48
-60
-60
46
-80
-80
44
-100
-100
42
-120
-120
40
-140
-140
38
-160
-160
36
-180
-180
- 200
- 200
34
0
100
200
0
100
200
0
100
200
Depth [L]
Depth [L]
Depth [L]
Fig. 11.9 Solution of Richards equation for infiltration within the soil compartment
