Numerical Evaluation of Fractional Vertical Soil Water Flow Equations

Significant deviations from standard Boltzmann scaling, which corresponds to normal or Fickian diffusion, have been observed in the literature for water movement in porous media. However, as demonstrated by various researchers, the widely used conventional Richards equation cannot mimic anomalous diffusion and ignores the features of natural soils which are heterogeneous. Within this framework, governing equations of transient water flow in porous media in fractional time and multi-dimensional fractional soil space in anisotropic media were recently introduced by the authors by coupling Brooks–Corey constitutive relationships with the fractional continuity and motion equations. In this study, instead of utilizing Brooks–Corey relationships, empirical expressions, obtained by least square fits through hydraulic measurements, were utilized to show the suitability of the proposed fractional approach with other constitutive hydraulic relations in the literature. Next, a finite difference numerical method was proposed to solve the fractional governing equations. The applicability of the proposed fractional governing equations was investigated numerically in comparison to their conventional counterparts. In practice, cumulative infiltration values are observed to deviate from conventional infiltration approximation, or the wetting front through time may not be consistent with the traditional estimates of Richards equation. In such cases, fractional governing equations may be a better alternative for mimicking the physical process as they can capture sub-, super-, and normal-diffusive soil water flow processes during infiltration.

[1]  D. Benson,et al.  Hydraulic conductivity, velocity, and the order of the fractional dispersion derivative in a highly heterogeneous system , 2002 .

[2]  B. Berkowitz,et al.  Measurement and analysis of non-Fickian dispersion in heterogeneous porous media. , 2003, Journal of contaminant hydrology.

[3]  Diego A. Murio,et al.  Implicit finite difference approximation for time fractional diffusion equations , 2008, Comput. Math. Appl..

[4]  R. Sidle,et al.  Spatially varying hydraulic and solute transport characteristics of a fractured till determined by field tracer tests, Funen, Denmark , 1998 .

[5]  Zaid M. Odibat,et al.  Computational algorithms for computing the fractional derivatives of functions , 2009, Math. Comput. Simul..

[6]  K. K. Watson.  An instantaneous profile method for determining the hydraulic conductivity of unsaturated porous materials , 1966 .

[7]  Luciano Lopez,et al.  The 1D Richards’ equation in two layered soils: a Filippov approach to treat discontinuities , 2017 .

[8]  Georges Vachaud,et al.  Hysteresis During Infiltration and Redistribution in a Soil Column at Different Initial Water Contents , 1971 .

[9]  W. Green.  Studies in soil physics : I. The flow of air and water through soils , 1911 .

[10]  Sean Andrew McKenna,et al.  On the late‐time behavior of tracer test breakthrough curves , 2000 .

[11]  John A. Cherry,et al.  Migration of contaminants in groundwater at a landfill: A case study: 4. A natural-gradient dispersion test , 1983 .

[12]  N. Su.  Mass-time and space-time fractional partial differential equations of water movement in soils: theoretical framework and application to infiltration , 2014 .

[13]  Keith Beven,et al.  Macropores and water flow in soils revisited , 2013 .

[14]  L. A. Richards.  Capillary conduction of liquids through porous mediums , 1931 .

[15]  Y. Pachepsky,et al.  Generalized Richards' equation to simulate water transport in unsaturated soils , 2003 .

[16]  Randel Haverkamp,et al.  A Comparison of Numerical Simulation Models For One-Dimensional Infiltration1 , 1977 .

[17]  J. Philip,et al.  THE THEORY OF INFILTRATION: 2. THE PROFILE OF INFINITY , 1957 .

[18]  W. Green,et al.  Studies on Soil Phyics. , 1911, The Journal of Agricultural Science.

[19]  S. L. Rawlins,et al.  A Test of the Validity of the Diffusion Equation for Unsaturated Flow of Soil Water 1 , 1963 .

[20]  D. Benson,et al.  Time and space nonlocalities underlying fractional-derivative models: Distinction and literature review of field applications , 2009 .

[21]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[22]  J. Philip,et al.  Theory of Infiltration , 1969 .

[23]  J. Parlange,et al.  Anomalous water absorption in porous materials , 2003 .

[24]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[25]  H. Gerke,et al.  Dual‐permeability model for flow in shrinking soil with dominant horizontal deformation , 2012 .

[26]  R. Kanivetsky,et al.  On the propagation of nonlinear transients of temperature and pore pressure in a thin porous boundary layer between two rocks , 2019, Journal of Hydrology.

[27]  M. Caputo.  Linear Models of Dissipation whose Q is almost Frequency Independent-II , 1967 .

[28]  N. Su.  Theory of infiltration: infiltration into swelling soils in a material coordinate. , 2010 .

[29]  Hongguang Sun,et al.  A fractal Richards' equation to capture the non-Boltzmann scaling of water transport in unsaturated media. , 2013, Advances in water resources.

[30]  Brian Berkowitz,et al.  Theory of anomalous chemical transport in random fracture networks , 1998 .

[31]  D. R. Nielsen,et al.  Experimental Consideration of Diffusion Analysis in Unsaturated Flow Problems 1 , 1962 .

[32]  M. Kavvas,et al.  Governing equations of transient soil water flow and soil water flux in multi-dimensional fractional anisotropic media and fractional time , 2016 .

[33]  David A. Benson,et al.  Space‐fractional advection‐dispersion equations with variable parameters: Diverse formulas, numerical solutions, and application to the Macrodispersion Experiment site data , 2007 .

[34]  D. R. Nielsen,et al.  Spatial variability of field-measured soil-water properties , 1973 .

[35]  J. Philip.  THE THEORY OF INFILTRATION: 1. THE INFILTRATION EQUATION AND ITS SOLUTION , 1957 .

[36]  E. Simpson,et al.  Laboratory evidence of the scale effect in dispersion of solutes in porous media , 1987 .

[37]  Vaughan R Voller,et al.  On a fractional derivative form of the Green–Ampt infiltration model , 2011 .

[38]  I. Vardoulakis,et al.  Modelling infiltration by means of a nonlinear fractional diffusion model , 2006 .

[39]  E. Park,et al.  Numerical solution of the Kirchhoff-transformed Richards equation for simulating variably saturated flow in heterogeneous layered porous media , 2019 .

[40]  J. Milczarek,et al.  Neutron radiography study of water absorption in porous building materials: anomalous diffusion analysis , 2004 .

[41]  T. Talsma,et al.  Infiltration and water movement in an in situ swelling soil during prolonged ponding , 1976 .

[42]  Ji-Huan He.  Approximate analytical solution for seepage flow with fractional derivatives in porous media , 1998 .

[43]  M. Küntz,et al.  Experimental evidence and theoretical analysis of anomalous diffusion during water infiltration in porous building materials , 2001 .

[44]  R. Horton.  An Approach Toward a Physical Interpretation of Infiltration-Capacity1 , 1941 .

[45]  W. R. Gardner.  Field measurement of soil water diffusivity. , 1970 .

[46]  G. Zaslavsky.  Chaos, fractional kinetics, and anomalous transport , 2002 .

[47]  M. V. Genuchten,et al.  A dual-porosity model for simulating the preferential movement of water and solutes in structured porous media , 1993 .

[48]  S. Logsdon.  TRANSIENT VARIATION IN THE INFILTRATION RATE DURING MEASUREMENT WITH TENSION INFILTROMETERS , 1997 .