Shooting the Numerical Solution of Moisture Flow Equation with Root Water Uptake Models: A Python Tool

Fabio V. Difonzo, Costantino Masciopinto, Michele Vurro and Marco Berardi ()
Additional contact information
Fabio V. Difonzo: Czech Technical University in Prague
Costantino Masciopinto: Consiglio Nazionale delle Ricerche
Michele Vurro: Consiglio Nazionale delle Ricerche
Marco Berardi: Consiglio Nazionale delle Ricerche

Water Resources Management: An International Journal, Published for the European Water Resources Association (EWRA), 2021, vol. 35, issue 8, No 16, 2553-2567

Abstract: Abstract Modeling the water uptake by plant roots is a key issue in studying soil processes, which are governed by water dynamics: a comprehensive understanding and forecast of such dynamics is a relevant issue in managing water resources. Typically, movement of water in soils and uptake by roots are described by the Richards’ equation with a sink term, and numerical treatment of this problem is still a challenge, together with its practical implementations in user-friendly softwares. In order to tackle this problem, in the present paper we propose a simple and computationally fast algorithm developed as a Python code, implementing a numerical approach based on the shooting method, a classical tool for handling boundary value problems (BVPs) arising here from a discretization recently introduced for Richards’ equation: such a method is applied to the linearized Richards’ equation with Gardner’s hydraulic functions. This method is implemented also in MATLAB, in order to accomplish comparisons with built-in MATLAB solver for parabolic partial differential equations. The Python code is made available to readers, and is intended to be an easy tool for handling this problem in the framework of Gardner’s constitutive relations, filling the gap of other commercial codes, which do not provide choice of Gardner functions. Many numerical simulations are performed: the results are promising, since the proposed method behaves efficiently and in some cases it is able to converge even when the MATLAB solver fails; mass balance properties and order of accuracy issues are also investigated.

Keywords: Richards’ equation; Root water uptake models; Shooting method; MATLAB software; Python (search for similar items in EconPapers)
Date: 2021
References: View references in EconPapers View complete reference list from CitEc
Citations: View citations in EconPapers (3)

Downloads: (external link)
http://link.springer.com/10.1007/s11269-021-02850-2 Abstract (text/html)
Access to the full text of the articles in this series is restricted.

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:spr:waterr:v:35:y:2021:i:8:d:10.1007_s11269-021-02850-2

Ordering information: This journal article can be ordered from
http://www.springer.com/economics/journal/11269

DOI: 10.1007/s11269-021-02850-2

Access Statistics for this article

Water Resources Management: An International Journal, Published for the European Water Resources Association (EWRA) is currently edited by G. Tsakiris

More articles in Water Resources Management: An International Journal, Published for the European Water Resources Association (EWRA) from Springer, European Water Resources Association (EWRA)
Bibliographic data for series maintained by Sonal Shukla () and Springer Nature Abstracting and Indexing ().