Numerical Solution of Water Flow in Unsaturated Soil with Evapotraspiration

Flow movement in unsaturated soil can be expressed by Richards equation. This equation can be obtained by applying the mass conversation law and the Darcy law. In this work, we solve one-dimensional Kirchhoff transformed Richards equation with loss of water due to the evaporation of unsaturated porous media (soils) and transpiration of plants numerically using Crank-Nicolson method. The result has compared with evapotranspiration function and without it in the governing equation. It has found that an additional work in time and flow movement is needs to reach the given boundary condition for the model without evapotranspiration.