Abstract

A new second-order stabilised numerical method for the Richards' equation with source term is proposed to solve the water infiltration in porous soil. This new method to solve Richard's equation is based on a modification of the head pressure representation with respect to the water content together with the hydraulic conductivity. We show that the modified relations provide stability while preserving the order and accuracy of the numerical solution. Furthermore, a new second-order positivity preserving scheme in time is implemented in substitution of the popular Crank-Nicolson technique. Applications to root suction in critical situations such as water saturation or rarefaction are carried out to highlight the benefits of the method. We also present a strategy to deal with multi-layers problem, by introducing a fixed point strategy to determine the complete solution on the whole domain.

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 call

Disclaimer: 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.