Evaluation of Method of Line (MOL) for Solution of Soil Water Flow Equations and Comparison with Finite Element Method (FEM)
Main Article Content
Abstract
The objective of this research is to evaluate Method of Line (MOL) used for solution of water flow equations through porous media using MATLAB package functions for solution of ordinary differential equations ODE,s , instead of writing long programs codes. sink & source term to MOL model were included. Then Comparing MOL model with another model that uses finite element method in solving water flow equations(FEM) in one dimensional flow using computer program code in FORTRAN. Two cases were examined for evaluation and comparison of these two models. Firstly, infiltration phenomena using sandy soil was studied with the same parameter for both models. Results show that there is a divergence between the two models along time of 60 minutes of infiltration. Changes of moisture content with soil depth were sharp with FEM model. Second case, data of the volume of water content for wheat field where used taking irrigation and evaporation into account, along the growth period of wheat crop and different depths up to 100 cm. Results show that output of FEM model has high degree of agreement with the measured data for all depths and along all period of growth. Data given by MOL model were less in values than measured data for all depths and along all period of wheat growth time.
Metrics
Article Details

This work is licensed under a Creative Commons Attribution 4.0 International License.
THIS IS AN OPEN ACCESS ARTICLE UNDER THE CC BY LICENSE http://creativecommons.org/licenses/by/4.0/
Plaudit
References
Simunek , J., T. Vogel, and M. T. van Genuchten.1994, “The SWMS 2D Code for Simulating Water Flow and Solute Transport in Two-Dimensional Variably Saturated Media”, Version 1.21, 197 pp., U.S. Salinity Lab.,Agric. Res. Serv., U.S. Dep. of Agric., Riverside, Calif.
Schiesser, W.E., 1991, “The Numerical Method of Lines: Integration of Partial Differential Equations: ODEs, DAEs and PDEs”. Academic Press, San Diego.
Shampine, L.F., Reichelt, M.W., 1994. “The MATLAB ODE” Suite Report 94-6. Math. Dept SMU, Dallas.
Lee, H.S., Braddock, R.D., Sander, G.C., 1998, “Automating the Method of Lines for Modeling Moisture Flow in the Unsaturated Zone”, In: Noye, B.J., Teubner, M.D., Gill, A.W. (Eds.), Computational Techniques and Applications: CTAC-97. World Scientific Publishing Co, Singapore.
محمد, عصام محمود 2000 ,"نمذجة حركة
الماء في التربة تحت الظروف الحقلية باستخدام
طريقة العناصر المحدودة",رسالة دكتوراه, كلية
الزراعة والغابات, جامعة الموصل.
Richards, L. A. 1931, “Capillary Conduction of Liquids Through Porous Mediums,” Physics 1:318-333. DOI: https://doi.org/10.1063/1.1745010
Vogel, T., K. Huang, R. Zhang, and M. T. van Genuchten, 1996, “The HYDRUS Code for Simulating One-Dimensional Water Flow”, solute transport, and heat movement in variable-saturated media, version 5, U.S. Salinity Lab., Agric. Res. Serv., U.S. Dep. of Agric., Riverside, Calif.,.
Celia , M.A., and E.T.Boututas ,R.L.Zarba.Laao.1990, “A General Mass - Conservative Numerical Solution for the Unsaturated Flow Equation”, water Resour. Res. 26, 1483-1496. DOI: https://doi.org/10.1029/WR026i007p01483
Lee, H.S., C.J. Matthews a, R.D. Braddock a, G.C. Sander b, F. Gandola.2003, “A MATLAB Method of Lines Template for Transport Equations”, Environmental Modelling & Software 19 (2004) 603–614. DOI: https://doi.org/10.1016/j.envsoft.2003.08.017
Watson, K.K., Sardana, V.A., Sander, G.C., 1995, “Comparison of Analytical and Numerical Results for Constant Flux Infiltration”, Journal of Hydrology 165, 101–112. DOI: https://doi.org/10.1016/0022-1694(95)92768-9
Van Genuchten, M.Th., 1980, “A Closed - Form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils”, Soil Sei. Am., J., 44: 892-989. DOI: https://doi.org/10.2136/sssaj1980.03615995004400050002x