Abstract
In this paper, the numerical solutions for groundwater flow in unsaturated layered soil using the Richards equation are presented. A linearisation process for the nonlinear Richards equation to deal with groundwater flow in unsaturated layered soil is derived. To solve one-dimensional flow in the unsaturated zone of layered soil profiles, flux conservation and the continuity of pressure potential at the interface between two consecutive layers are considered in the numerical model. In addition, a novel method, named the dynamical Jacobian-inverse free method, incorporated with a two-side equilibration algorithm for solving ill-conditioned systems with extreme contrasts in hydraulic conductivity is proposed. The validity of the model is established in numerous test problems by comparing the numerical results with the analytical solutions. The results show that the proposed method can improve convergence and numerical stability for solving groundwater flow in unsaturated layered soil with extreme contrasts in hydraulic conductivity.
Acknowledgements
This study was partially supported by the Central Geological Survey of the Ministry of Economic Affairs, R.O.C. The author thanks the Central Geological Survey for their generous financial support.
References
[1] Y.-S. Fok, One-dimensional infiltration into layered soils, J. Irrig. Drain. Div. 96 (1970), 121–129.10.1061/JRCEA4.0000713Search in Google Scholar
[2] D. E. Aylor and J.-Y. Parlange, Vertical infiltration into a layered soil, Soil Sci. Soc. Am. J. 37 (1973), 673–676.10.2136/sssaj1973.03615995003700050015xSearch in Google Scholar
[3] A. Y. Hachum and J. F. Alfaro, Rain infiltration into layered soils: Prediction, J. Irrig. Drain. Div. Am. Soc. Civil Eng. 106 (1980), 311–319.10.1061/JRCEA4.0001321Search in Google Scholar
[4] Z. Samani, A. Cheraghi, and L. Willardson, Water movement in horizontally layered soils, J. Irrig. Drain. Eng. 115 (1989), 449–456.10.1061/(ASCE)0733-9437(1989)115:3(449)Search in Google Scholar
[5] C.-Y. Ku and Y.-H. Tsai, Solving nonlinear problems with singular initial conditions using a perturbed scalar homotopy method, Int. J. Nonlinear Sci. Numer. Simul. 14 (2013), 367–375.10.1515/ijnsns-2013-0029Search in Google Scholar
[6] R. Hanks and S. Bowers, Numerical solution of the moisture flow equation for infiltration into layered soils, Soil Sci. Soc. Am. J. 26 (1962), 530–534.10.2136/sssaj1962.03615995002600060007xSearch in Google Scholar
[7] F. Whisler and A. Klute, Analysis of infiltration into stratified soil columns Proc Wageningen Symp, Proc., IAHS Symposium on Water in the Unsaturated Flow, Wageningen, The Netherlands, 451–470, 1966.Search in Google Scholar
[8] P. Moldrup, D. E. Rolston, and L. A. Hansen, Rapid and numerically stable simulation of one dimensional, transient water flow in unsaturated, layered soils, Soil Sci. 1483 (1989), 219–226.10.1097/00010694-198909000-00009Search in Google Scholar
[9] R. Srivastava and J. T.-C. Yeh, Analytical solutions for one-dimensional, transient infiltration toward the water table in homogeneous and layered soils, Water Resour. Res. 27 (1991), 753–762.10.1029/90WR02772Search in Google Scholar
[10] N. Romano, B. Brunoneb, and A. Santini, Numerical analysis of one-dimensional unsaturated flow in layered soils, Adv. Water Resour. 21 (1998), 315–324.10.1016/S0309-1708(96)00059-0Search in Google Scholar
[11] C. Corradini, F. Melone, and R. E. Smith, Modeling local infiltration for a two-layered soil under complex rainfall patterns, J. Hydrol. 237 (2000), 58–73.10.1016/S0022-1694(00)00298-5Search in Google Scholar
[12] R. Leconte and F. P. Brissette, Soil moisture profile model for two-layered soil based on sharp wetting front approach, J. Hydrologic Eng. 62 (2001), 141–149.10.1061/(ASCE)1084-0699(2001)6:2(141)Search in Google Scholar
[13] L. A. Richards, Capillary conduction of liquids through porous mediums, J. Appl. Phys. 1 (1931), 318–333.10.1063/1.1745010Search in Google Scholar
[14] B. M. Das, Principles of geotechnical engineering, 7th edn, Cengage Learning, 200 First Stamford Place, Suite 400, Stamford, CT 06902, USA, 2010.Search in Google Scholar
[15] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical recipes, the art of scientific computing, Cambridge University Press, Cambridge, 2007.Search in Google Scholar
[16] C.-S. Liu and C.-W. Chang, Novel methods for ill-conditioned linear equations, J. Mar. Sci. and Technol. 17 (2009), 216–227.10.51400/2709-6998.1958Search in Google Scholar
[17] C. Vuik, A. Segal, and J. A. Meijerinky, An efficient preconditioned CG method for the solution of a class of layered problems with extreme contrasts in the coefficients, J. Computat. Phys. 152 (1999), 385–403.10.1006/jcph.1999.6255Search in Google Scholar
[18] C.-S. Liu, Modifications of steepest descent method and conjugate gradient method against noise for ill-posed linear systems, Commun. Numer. Anal. 2012 (2012), 1–24.10.5899/2012/cna-00115Search in Google Scholar
[19] C.-S. Liu, H.-Ki. Hong, and S. N. Atluri, Novel algorithms based on the conjugate gradient method for inverting ill-conditioned matrices, and a new regularization method to solve ill-posed linear systems, CMES: Comput. Model. Eng. Sci. 60 (2010), 279–308.Search in Google Scholar
[20] C.-S. Liu, An optimal multi-vector iterative algorithm in a Krylov subspace for solving the ill-posed linear inverse problems, CMC: Comput. Mater. Continua. 33 (2013a), 175–198.Search in Google Scholar
[21] C.-S. Liu, A two-side equilibration method to reduce the condition number of an ill-posed linear system, CMES: Comput. Model. Eng. Sci. 91 (2013b), 17–42.Search in Google Scholar
[22] C.-Y. Ku, W. Yeih, and C.-S. Liu, Dynamical Newton-like methods for solving ill- conditioned systems of nonlinear equations with applications to boundary value problems, CMES: Comput. Model. Eng. Sci. 76 (2011), 83–108.Search in Google Scholar
[23] M. T. Van Genuchten, A closed-form equation for predicting the hydraulic conductivity of unsaturated soils, Soil Sci. Soc. Am. J. 44 (1980), 892–898.10.2136/sssaj1980.03615995004400050002xSearch in Google Scholar
[24] W. Gardner, Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table, Soil Sci. 85 (1958), 228–232.10.1097/00010694-195804000-00006Search in Google Scholar
[25] A. Warrick, Soil water dynamics, Oxford University Press, New York, 2003.10.1093/oso/9780195126051.001.0001Search in Google Scholar
[26] C.-Y. Ku, A novel method for solving ill-conditioned systems of linear equations with extreme physical property contrasts, CMES: Comput. Model. Eng. Sci. 96 (2013), 409–434.Search in Google Scholar
[27] W. H. Green and G. A. Ampt, Studies on soil physics I. The flow of air and water through soils, J. Agric. Sci. 4 (1911), 1–24.10.1017/S0021859600001441Search in Google Scholar
[28] F. T. Tracy, Analytical and numerical solutions of Richards’ equation with discussions on relative hydraulic conductivity, INTECH Open Access Publisher, Janeza Trdine 9, 51000 Rijeka, Croatia, 2011.10.5772/18502Search in Google Scholar
©2015 by De Gruyter
