Stochastic finite element modelling of water flow in variably saturated heterogeneous soils |
| |
| Authors: | M Mousavi Nezhad A A Javadi F Abbasi |
| |
| Institution: | 1. Computational Geomechanics Group, College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter EX4 4QF, U.K.;2. Agricultural Engineering Research Institute, AERI, P.O. Box 31585‐845, Karaj, Iran |
| |
| Abstract: | Water flow is greatly influenced by the characteristics of the domain through which the process occurs. It is generally accepted that earth materials have extreme variations from point to point in space. Consequently, this heterogeneity results in high variation in hydraulic properties of soil. In order to develop an accurate predictive model for transport processes in soil, the effects of this variability should be considered. In this study a two‐dimensional stochastic finite element flow model was developed for simulation of water flow through unsaturated soils. In this model, the stochastic partial differential governing equation of water flow, obtained from implementation of the perturbation‐spectral stochastic method on classical Richard's equation, was solved using a finite element method in the space domain and a finite difference scheme in the time domain. The effective hydrological parameters embedded in the mathematical model depend on time derivatives of capillary tension head; this makes possible to consider the hysteresis due to large‐scale variability of soil hydrological properties. The model is also capable of simulating infiltration and evaporation events and rapid change in the land surface boundary condition from one type event to another, based on a scheme used in the model for implementation of land surface boundary condition. The model was validated with the data obtained from a layered lysimeter test. The model was also used to simulate water flow under a long irrigation furrow. The results obtained with this model show better agreement with experimental measurements in comparison with a deterministic model. The possible reason for this agreement is that in the developed model, the influence of the variability of the properties of soil and effects of parameter hysteresis on water flow and water content redistribution are considered. Copyright © 2010 John Wiley & Sons, Ltd. |
| |
| Keywords: | stochastic finite element method heterogeneous soil water flow unsaturated soil |
|
|
