Research Article  Open Access
The Similarity Hypothesis and New Analytical Support on the Estimation of Horizontal Infiltration into Sand
Abstract
A method based on a specific powerlaw relationship between the hydraulic head and the Boltzmann variable, presented using a similarity hypothesis, was recently generalized to a range of powers to satisfy the Bruce and Klute equation exactly. Here, considerations are presented on the proposed similarity assumption, and new analytical support is given to estimate the water density flux into and inside the soil, based on the concept of sorptivity and on BuckinghamDarcy's law. Results show that the new analytical solution satisfies both theories in the calculation of water density fluxes and is in agreement with experimental results of water infiltrating horizontally into sand. However, the utility of this analysis still needs to be verified for a variety of different textured soils having a diverse range of initial soil water contents.
1. Introduction
Based on physical laws of similarity applied to the rate of work required for water to wet and move through a soil, Prevedello et al. [1] presented an analytic solution of a Boltzmann transformed equation of continuity for horizontal infiltration which is derived without invoking the concept or use of the soil water diffusivity function. The derivation assumes that a similarity exists between the shapes of the soil water retention function and the Boltzmann transformation , and the solution successfully described soil water content profiles experimentally measured for different infiltration times into a homogeneous sand. More recently, an extension of this theory generalized the solution to a range of powers to include the saturated zone [2], to satisfy the Bruce and Klute equation exactly.
With a similar assumption, but not exactly as expressed by Prevedello et al. [1], Prevedello et al. [3] obtained a new analytic solution of the Richards equation for the infiltration into the same sand that holds for all infiltration times from zero to infinity, including vertical directions without making use of empirical constants. In this case, the derivation assumed that a similarity exists between the soil water retention function and the soil water content distribution within the soil profile during infiltration.
Although the similarity assumptions used by [1, 3] seem different, we now show that the mathematical development of both leads to the same equation for the description of the soil water content profiles for horizontal infiltration. We also present an analysis to obtain analytical equations to estimate the water density flow into and inside the soil, according to the sorptivity concept and the BuckinghamDarcy law. These analytical equations are compared with experimental results of horizontal infiltration into the same homogeneous marine sand.
2. Theoretical Considerations
2.1. Extended Similarity Hypothesis
Two identical solutions for the horizontal water infiltration into a soil were presented by [1, 3], coincidently (39) in both papers, but using two apparently different initial postulates. The explanation for that is given by a generalization of the similarity theory proposed in [1]. A generalization of the similarity between the Boltzmann transformation and the soil water retention function can be written as follows: in which is a positive integer, the soil water content (m^{3}m^{3}), and the matric water potential head (). With (1) being the fundamental hypothetical equation, it follows that or and introducing the horizontal position coordinate and the time (s), we have note that .
When , From (2) it follows that
Substituting (7) into the BuckinghamDarcy equation written in the direction in terms of the soil hydraulic conductivity as a function of , we have where the subscript refers to the model hypothesis of (4). Knowing (8), changing the variables so that the dependent variable becomes the distance , and remembering that because the partial derivative of is taken with respect to , we have which leads to For a fixed time, based on (7), and multiplying by , we have or Holding constant and integrating, the above equation becomes Evaluating at when , Recalling (3) in the form of and equating the two above equations that define , we have Returning to (4), we know that Substituting into the above equation, we obtain Recalling that , the solution is
For the cases in which the soil is initially very dry (i.e., ), the above equation reduces to which is (39) presented in [1, 3].
From (5) we have According to (21), the generalized solution of is But the fundamental hypothesis that gave rise to (22) is The derivative of (22) with respect to is Dividing (22) by ,
Equating (24) and (25), since both follow the postulate, the result is . This shows that any value of different from 2 cannot be accepted in this theory of similarity, as originally postulated and confirmed in [1], since it is with the square root of the time that the advancement of the wetting front remains linear in the process of horizontal infiltration.
2.2. Explanation Notes
Considering and since , [1, equation (19)], for the cases in which at and when the soil is initially very dry , can be written as Therefore, we know that and we can write (26) as follows: which shows that the gradient formulated as an initial postulate in [1] is simplified during the theoretical development to . This statement can also be verified starting with (11). In this equation, holding constant and making , we can write where we can see that the initially formulated gradient becomes . The integration of the equation above leads to (22) in [1, 3].
The solution obtained by [3] for horizontal infiltration case, although identical to that obtained by [1], presupposes a relation of the type with , that is, or or , and, therefore, it cannot be taken as formal as that presented in [1]. As explained above, this is due to the fact that independently of the chosen value of in the generalised similarity theory here introduced, the gradient of (7), as initially postulated in [1], is always simplified during the theoretical development from to , so that the BD equation is always converted to , and that, therefore, leads to a final solution that is identical to that obtained by [1].
In this way it is shown that the same soil water profile , expressed by (22) in both papers [1, 3], can be obtained from both similarity approaches, that is, or as images of the soil water retention curve . However, only as an image of translates formally what was postulated into a solution, such as shown in [1] or in the extended similarity hypothesis shown above, in which only carries out formally this postulate.
2.3. The Similarity Hypothesis and New Analytical Support on the Estimation of Horizontal Infiltration
Since Prevedello et al. [1] found a solution such that is proportional to , so that is proportional to , it follows that for , from the definition of , at both and are equal to zero. At a first glance, one could assume that to have a finite flux density at , should be infinite at . According to this interpretation, Barry et al. [2] introduced a nonzero airentry value (e.g., Haverkamp et al. [4]), such that for , the soil is unsaturated and saturated for , and generalized the Prevedello et al. [1] solution to a range of powers. We here use the original interpretation of Prevedello et al. [1] to calculate the flux density at and at , without running into the shortcomings mentioned above. From Phillip’s analysis, we know that for the case of horizontal infiltration,
If in (30) we substitute by the Bruce and Klute equation or (1) of Barry et al. [2] are obtained. There are many analytical forms to prove the equality in (30). Here, our option was to obtain an analytical expression for both the members and to compare the results obtained for each one. According to the similarity assumption proposed by [1] we know that From van Genuchten we know that so that (31) becomes so that
The water flux density can also be evaluated through the BuckinghamDarcy law directly, using the Prevedello et al. [1, 3] similarity assumption: But, according to the similarity assumption, each value of of the wetting profile is associated to its position value . Therefore, does not represent directly the total or representative matric gradient at the wetting front. In this way, (35) should be interpreted as the water density flow that enters at in response to a “representative” gradient that occurs at the wetting front. Thus, if is taken as the average matric potential of all the active potentials, then we can find analytically this representative matric potential, as well as the position where it acts at the wetting front, applying the similarity assumption. In this way, the previous equation can be written as where is the representative potential active at the wetting front and the position where this potential acts, defined by or
Notice that this interpretation provides a constant value of independent of time, and it follows that
From the soil surface, where the soil is saturated (), the wetting front advances as a response to an average matric potential , and with a force (average matric gradient), which should yield the same results as (34).
The analytical equation to estimate water density flow inside the soil can be determined simply by substituting for or into (40), because the average gradient remains fixed, but inside the soil we have the soil water content or matric potential for each position :
3. Material and Methods
To validate experimentally the equations proposed above to quantify the water density into and inside the soil during the horizontal infiltration process for homogeneous sand, we used the same data presented in [1, 3].
The van Genuchten [6] model was used to describe the soil water retention curve where , , , (residual soil water content), and are the classical parameters of the soil water retention curve. Using a nonlinear regression program as the one suggested by Boratto [5], (42) was fitted to experimental values to obtain the independent parameters , , , , and .
According to van Genuchten [6], introducing (42) into the Mualem [7] theory, the function is obtained as
Data for the saturated hydraulic conductivity and for the parameters of the water retention curve for the marine sand are presented in Table 1, as previously shown in [1, 3]. With a coefficient of determination greater than 0.99, the van Genuchten [6] model provides an excellent representation of the observed values.

4. Results and Discussion
4.1. The Water Density Flow into the Soil (), according to Philip’s Concept of Sorptivity
From the boundary conditions and the soil parameters (Table 1) and assuming the initial soil water content equal to m^{3}·m^{3}, (34) reduces to which is shown in Figure 1 as a function of time, according to the analysis of Philip [8].
Measured soil water content profiles (Figure 5 of Prevedello et al. [1]) were integrated with respect to distance for times , , and equal to 15, 60, and 300 min, resulting in 165, 330, and 732 mm, respectively. We also know that the cumulative infiltration into the soil surface
is or where is the sorptivity.
In light of (46) and substituting the value mm·(min)^{1/2}, we have mm, mm, and mm, respectively, which are very close to those estimated by integration of experimental data and showed in Figure 2.
4.2. The Water Density Flow into the Soil (), According to BuckinghamDarcy Law
To verify the analytical validity of (40) to quantify the water density flow at , (38) and (39) were evaluated, resulting, respectively, in m and . These results show that the similarity theory proposed by [1], here represented by (36) or (40), does not imply in a water flux density zero at . Substituting these values and m/min into (40), we obtain which is exactly the same value obtained in the calculation through Philip’s concept for sorptivity (44) and what therefore also adjusts very well with the experimental results showed in the Figure 2.
Another apparently analogous way to obtain the results of (47) above is to consider the derivative from (20), that is, . But according to the similarity theory proposed in [1, 3], the value of is a constant for each position at which it is only dependent on the time of advancement, and the derivative of a constant is null. The same restriction is applicable to the Green and Ampt theory [9] when we want to obtain the matric gradient through derivative of (the pressure head at the wetting front) with respect to because is always a constant independent of . So, the derivative of with respect to cannot be used here, except if associated to the concept of the similarity, and as exposed previously in item b above, we know that the matric gradient is reduced to , and therefore we can also write (40) as
Comparing (30) with (48), we know that for , for , and for . It is now readily apparent from the above analysis that the distribution of was substituted by the constant in (48) to determine the water flow density at in response to the force that occurs at distance where . Inside of the soil the values of or are consistent with the values of found from the values of in (39) of Prevedello et al. [1, 3] because (39) provides the shape of the soil water content profile at the wetting profile while the average gradient / remains active on the time.
Equation (48) with values of from Table 1 and also yields values very close to the experimental results (Figure 2) and exactly the same results obtained in the calculation through Philip’s concept for sorptivity (44).
4.3. The Density of Flow inside of the Soil, for Any Position and Time
Considering the boundary conditions of soil water content varying between 0.0187 and 0.387 m^{3}·m^{3}, to use (41) we first estimate the corresponding matric potentials for any stipulated between 0.0187 and 0.387 m^{3}·m^{3}, to obtain their respective positions from (20) and their respective ) values from (43).
Through (41) we generate the water density flow for each position inside the soil, that is,
Analogously as before in (48), the water density flow inside the soil can also be determined through the following: in which , and the values of and are showed in the Table 1.
Both (49a) and (49b) yield exactly the same values, and Figure 3 shows the water density flow estimated from those equations for each reached at 15, 60 and 300 minutes of the horizontal infiltration into the homogeneous marine sand.
It can be noted that for the water density flow values coincide with those obtained from (44), (47), or (48), as it should be.
From the above, it follows that
or
indicating that the gradient decreases with the time, shown in Figure 4 for times of 15, 60, and 300 minutes.
To verify how accurate these estimated results are, the same boundary value problem was solved numerically from Richard’s equation without any similarity assumption, using the procedure of Philip [10] with appropriate initial and boundary conditions and using the measured soil water properties shown in Section 3 (calculations based on m^{3}·m^{3}). Results showed that both water density flows at and inside the soil were very similar in the transmission zone and wetting front for each of the three chosen times of horizontal infiltration. However, the utility of this analysis still needs to be verified for a variety of different textured soils having a diverse range of initial soil water contents.
5. Conclusion
The new analytical supports to estimate the water density into and inside of the soil during horizontal infiltration agree with the experimental results, with the concept of sorptivity proposed by Philip [8] as well as the BuckinghamDarcy law, when a “representative” matric gradient is considered acting at the wetting front.
Acknowledgment
The authors acknowledge Donald R. Nielsen and Klaus Reichardt for their support and guidance in preparing this paper.
References
 C. L. Prevedello, J. M. T. Loyola, K. Reichardt, and D. R. Nielsen, “New analytic solution of Boltzmann transform for horizontal water infiltration into sand,” Vadose Zone Journal, vol. 7, no. 4, pp. 1170–1177, 2008. View at: Publisher Site  Google Scholar
 D. A. Barry, J.Y. Parlange, C. L. Prevedello, J. M.T. Loyola, K. Reichardt, and D. R. Nielsen, “Extension of a recent method for obtaining exact solutions of the Bruce and Klute equation,” Vadose Zone Journal, vol. 9, no. 2, pp. 496–498, 2010. View at: Publisher Site  Google Scholar
 C. L. Prevedello, J. M. T. Loyola, K. Reichardt, and D. R. Nielsen, “New analytic solution related to the Richards, Philip, and GreenAmpt equations for infiltration,” Vadose Zone Journal, vol. 8, no. 1, pp. 127–135, 2009. View at: Publisher Site  Google Scholar
 R. Haverkamp, J. Y. Parlange, J. L. Starr, G. Schmitz, and C. Fuentes, “Infiltration under ponded conditions: 3. A predictive equation based on physical parameters,” Soil Science, vol. 149, no. 5, pp. 292–300, 1990. View at: Google Scholar
 F. Boratto, BASIC for Engineers and Scientists, Livros Técnicos e Científicos, Rio de Janeiro, Brazil, 2nd edition, 1984.
 M. TH. van Genuchten, “A closedform equation for predicting the hydraulic conductivity of unsaturated soils,” Soil Science Society of America Journal, vol. 44, no. 5, pp. 892–898, 1980. View at: Google Scholar
 Y. Mualem, “A new model for predicting the hydraulic conductivity of unsaturated porous media,” Water Resources Research, vol. 12, no. 3, pp. 513–522, 1976. View at: Google Scholar
 J. R. Philip, “The theory of infiltration: 1. The infiltration equation and its solution,” Soil Science, vol. 83, pp. 345–357, 1957. View at: Google Scholar
 W. H. Green and G. A. Ampt, “Studies in soil physics. Part I. The flow of air and water through soils,” The Journal of Agricultural Science, vol. 4, pp. 1–24, 1911. View at: Google Scholar
 J. R. Philip, “Numerical solution of equations of the diffusion type with diffusivity concentrationdependent,” Transactions of the Faraday Society, vol. 51, pp. 885–892, 1955. View at: Google Scholar
Copyright
Copyright © 2010 Celso L. Prevedello and Jocely M. T. Loyola. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.