Abstract
A computationally simple, numerical algorithm capable of solving a wide variety of two-dimensional, variably saturated flow problems is developed. Recent advances in modeling variably saturated flow are incorporated into the algorithm. A physically based form of the general, variably saturated flow equation is solved using finite differences (centered in space, fully implicit in time) employing the modified Picard iteration scheme to determine the temporal derivative of the water content. The algorithm avoids mass-balance errors in unsaturated regions and is numerically stable. The resulting system of linear equations is solved by a preconditioned conjugate gradient method, which is known to be computationally efficient for the type of equation set obtained. The algorithm is presented in sufficient detail to allow others to implement it easily, and is verified using four published, illustrative sets of experimental data.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
