Handling the water content discontinuity at the interface between layered soils within a numerical scheme

View/ Open
Author(s)
Matthews, Christopher
Cook, F.
Knight, John
Braddock, Roger
Year published
2004
Metadata
Show full item recordAbstract
In general, the water content (?) form of Richards' Equation is not used when modelling water flow through layered soil since ? is discontinuous at soil layers. Within the literature, there have been examples of models developed for layered soils using the ?-form of Richards' Equation. However, these models usually rely on an approximation of the discontinuity at the soil layer interface. For the first time, this paper will develop an iterative solution based on Newton's Method to explicitly solve for ? at the interface between two soils within a numerical scheme. The numerical scheme used within this paper is the Method of ...
View more >In general, the water content (?) form of Richards' Equation is not used when modelling water flow through layered soil since ? is discontinuous at soil layers. Within the literature, there have been examples of models developed for layered soils using the ?-form of Richards' Equation. However, these models usually rely on an approximation of the discontinuity at the soil layer interface. For the first time, this paper will develop an iterative solution based on Newton's Method to explicitly solve for ? at the interface between two soils within a numerical scheme. The numerical scheme used within this paper is the Method of Lines, however, the principles of the iterative solution could be used in other numerical techniques. The paper will show that the iterative scheme is highly efficient converging mostly within 2 iterations. The iterative solution will be tested on two main test cases 1) fine over coarse soil, 2) a coarse over fine soil to ensure that the convergence behaviour holds. A test case from the literature is also considered to ensure the accuracy of the solution.
View less >
View more >In general, the water content (?) form of Richards' Equation is not used when modelling water flow through layered soil since ? is discontinuous at soil layers. Within the literature, there have been examples of models developed for layered soils using the ?-form of Richards' Equation. However, these models usually rely on an approximation of the discontinuity at the soil layer interface. For the first time, this paper will develop an iterative solution based on Newton's Method to explicitly solve for ? at the interface between two soils within a numerical scheme. The numerical scheme used within this paper is the Method of Lines, however, the principles of the iterative solution could be used in other numerical techniques. The paper will show that the iterative scheme is highly efficient converging mostly within 2 iterations. The iterative solution will be tested on two main test cases 1) fine over coarse soil, 2) a coarse over fine soil to ensure that the convergence behaviour holds. A test case from the literature is also considered to ensure the accuracy of the solution.
View less >
Conference Title
SuperSoil 2004: Proceedings of the 3rd Australian New Zealand Soils Conference
Publisher URI
Copyright Statement
© The Author(s) 2004 Griffith University. This is the author-manuscript version of the paper. It is posted here with permission of the copyright owners for your personal use only. No further distributions permitted. The Regional Institute Ltd retains the non-exclusive right to publish the proceedings from the SuperSoil 2004, the 3rd Australian New Zealand Soils Conference in electronic format in perpetuity. Reproduction and the making available of this material for personal or non-commercial purposes is authorised, on condition that: attribution is given to the copyright owner; no official connection is claimed; the material is made available without charge or at cost; and the material is not subject to inaccurate, misleading or derogatory treatment.