Numerical Tests of the Lookup Table Method in Solving Richards’ Equation for Infiltration and Drainage in Heterogeneous Soils

The lookup table option, as an alternative to analytical calculation for evaluating the nonlinear heterogeneous soil characteristics, is introduced and compared for both the Picard and Newton iterative schemes in the numerical solution of Richards’ equation. The lookup table method can be a cost-effective alternative to analytical evaluation in the case of heterogeneous soils, but it has not been examined in detail in the hydrological modeling literature. Three layered soil test problems are considered, and the robustness and accuracy of the lookup table approach are assessed for uniform and non-uniform distributions of lookup points in the soil moisture retention curves. Results from the three one-dimensional test simulations show that the uniform distributed option gives improved convergence and robustness for the drainage problem compared to the non-uniform strategy. On the other hand, the non-uniform technique can be chosen for test problems involving flow into initially dry layered soils.