Laserfiche WebLink
<br />where Lp is the L value at t; ty' Because L in Eq, 47 <br />is not expressed as an explicit function of t, the value <br />of L for given t must be determined numerically, for <br />instance, by the Newton-Raphson method, After L is <br />computed, the value of L will be substituted either <br />into Eq, 38 or 44 for the computation of the f value, <br />Note that use of Eq. 44 has a slight advantage over <br />that of Eq, 38 because any enor involved in the <br />evaluation of dL/dt in Eq, 38 can be avoided, <br />However, use of either Eq, 38 or 44 requires the prior <br />computation of the values of tp and Lp' <br /> <br />Determination of the tp and Lp values, Before <br />ponding (t<O;tp), the cumulative inmtration, F(t), is rt <br />for constant rainfall intensity, r. Thus, from Eq, 38, <br />after integration and rearrangement, one obtains at t <br />;tp <br /> <br />(r - K ) <br />L 0 t <br />P - as - 00 p <br /> <br />, , , , ' , , .(48) <br /> <br />and from Eq. 44, at t; tp <br /> <br />~KS(h - ~o) <br />L - <br />P (r - Kg) <br /> <br />, , ' , , . .(49) <br /> <br />Equating Eqs, 48 and 49 yields <br /> <br />t = <br />p <br /> <br />(8 - 8 )(h - ~ ) <br />K s 0 0 <br />~ s (r - K )(r - K ) <br />o s <br /> <br />, , , , ,(50) <br /> <br />The values of Lp and tp can be determined from Eqs, <br />49 and 50 respectively, Like Eq, 47, the fl value must <br />be evaluated or given first in order for Eqs, 49 and 50 <br />to be useful, <br /> <br />EvalWltion of the fl value, As defined in Eq, 42, <br />the fl value depends on the K(z} and ",(z} <br />distributions, which are, of course, unknown. How- <br />ever, if one knows the K-"'relationship for the soil, <br />the fl value can be readily evaluated from Eq, 42 and <br />the K. '" relation, It is understood that the K.w <br />relation is hysteretic, but will be assumed unique <br />herein as long as the wetting process is only con. <br />sidered, Several forms of the K.",relation have been <br />assumed by different investigators in an attempt to <br />best fit the field or experimental data, Because the <br />soil after ponding (for h>O) has both the saturated <br />and unsaturated zones, all the forms of the K.", <br />relation must be modified to include K; K, for ",;;.0. <br /> <br />The simplest of the available K-'" relations is a <br />linear relationship between K and "', proposed by <br />Richards (1931), <br /> <br />K .. a;' + Ks <br /> <br />for ;. :OS 0 <br /> <br />(5la) <br /> <br />K = K <br />s <br /> <br />for ;. ~ 0 <br /> <br />(Sib) <br /> <br />where "an is simply a constant. <br /> <br />Raats (1971) and Braester (1973) among others <br />used an exponential form in K and "'. <br /> <br />K" K ealjJ <br />s <br /> <br />for 1/1 :S 0 <br /> <br />(52a) <br />(52b) <br /> <br />K ,., K <br />s <br /> <br />for I/J 2:: 0 <br /> <br />It should be noted that Eqs, 51 and 52 have only one <br />parameter, a, needed to be determined from the K-", <br />data, <br /> <br />Several researchers proposed multi-parameter <br />models which should fit the K-'" data better than the <br />one-parameter model. For example, Gardner's (1958) <br />original three.parameter model for "'<0 may be <br />extended to include "'; 0 so that <br /> <br />K - <br /> <br />a <br /> <br />(53a) <br /> <br />for I/J :S 0 <br /> <br />(_~)n + b <br /> <br />K - K <br />s <br /> <br />for tjI ~ 0 <br /> <br />(53b) <br /> <br />where a, b, and n are all constants. However, the <br />necessity of having K ; K, at "'; 0 results in the loss <br />of freedom for Eq, 53a to choose the a value other <br />than bI(,. Conversely, if K ; K, at '" ; 0 is merely <br />considered as a data point in the K-", relationship, the <br />values of a, b, and n in Eq, 53a can be determined by <br />using a least squares optimization procedure. Mter <br />that, one should take the value of K, equal to a/b, In <br />view of the difficulty in integrating Eq, 53a with <br />respe<:t to '" upon substitution of Eq, 53a into Eq, 42 <br />for the evaluation of the fl value, the following <br />compatible form similar to Eq, 53a is proposed <br />herein: <br /> <br />K- <br /> <br />a <br />(_". + b)n <br /> <br />(54a) <br /> <br />for 1/J:5 0 . . . . . <br /> <br />K - K <br />s <br /> <br />for ~~ 0 . , , , , (54b) <br /> <br />By the same token, the values of a, b, and n in Eq, <br />54a can be estimated by using a least squares <br />optimization technique, and hence the value of K, <br />may be taken to be equal to a/b". <br /> <br />The values of a, b, and n determined from Eqs. <br />53 and 54 are not all dimensionless, In order to make <br />these dimensionally different parameters consistent <br />(i.e" dimensionless), some investigators proposed <br />dimensionless forms of Eq, 53 with the same number <br />of parameters used in Eq, 53, Of them, Wei (1971) <br />has used the most general one that is nevertheless not <br />adopted in this study due to the same reason as given <br />for Eq, 53, <br /> <br />Substituting Eqs, 51, 52, and 54 into Eq, 42 <br />yields <br /> <br />34 <br />