Numerical model for saturated-unsaturated flow in deformable porous media: 2. The algorithm

• Published in: Water resources research 1978-04, Vol.14 (2), p.255-261
• Main Authors: Narasimhan, T. N., Witherspoon, P. A., Edwards, A. L.
• Format: Article
• Language: English
• ISSN: 0043-1397; 1944-7973
• DOI: 10.1029/WR014i002p00255

An integrated finite difference algorithm is presented for numerically solving the governing equation of saturated‐unsaturated flow in deformable porous media. In recognition that stability of the explicit equation is a local phenomenon a mixed explicit‐implicit procedure is used for marching in the time domain. In this scheme the explicit changes in potential are first computed for all elements in the system, after which implicit corrections are made only for those elements for which the stable time step is less than the time step being used. Time step sizes are automatically controlled in order to optimize the number of iterations, to control maximum change in potential during a time step, and to obtain desired outputs. Time derivatives, estimated on the basis of system behavior during two previous time steps, are used to start the iteration process and to evaluate nonlinear coefficients. Boundary conditions and sources can vary with time or with the dependent variable. Input data are organized into convenient blocks. Accuracy of solutions can be affected by modeling errors, different types of truncation errors, and convergence errors. The algorithm constitutes an efficient tool for analyzing linear and nonlinear fluid flow problems in multidimensional heterogeneous porous media with complex geometry. An important limitation is that the model cannot conveniently handle arbitrary anisotropy and other general tensorial quantities.