135
represented the outer wall of the lysimeter and equations (75) and (80)
were written as
k+1/2 k+1/2
- AR h + BR h = DR (95)
i,m-l i,m i,m
where
and
A V h
BR = 1 +
AR
(96)
k+1
k+1
k+1
+
BV h
- CV h = DV
(97)
i-1 ,m
i,m
i+l,m i,m
where
k+1/2 k+1/2
1 At K K 1 At
k i-l/2,m i+l/2,m k
DV = h + ( ) S
i,m i,m
4 A z C 2 C
i,m i,m
k+1/2
1 At r K
i.m-1/2 i,m-l/2 k+1/2 k+1/2
+ ( ) ( h h )
2 i,m-l i,m
2 a r r
i ,m
i ,m
(98)
The boundary condition for the inner radius was also represented as
a no flux boundary. For the inner radius r(i,j) = 0, therefore its
inclusion into equation (8) would result in a division by zero.
Therefore, the second term in equation (8) was rewritten as described in
equations (50) and (51) such that the finite difference equations for
the ADI solution of equation (51) were expressed as