Physical Sciences and Mathematics Commons^™

Finite Element Simulation Of Saturated-Unsaturated Subsurface Flow, George E. Blandford

KWRRI Research Reports

A two-dimensional transient model for flow through saturated-unsaturated porous media is developed, The model numerically solves the pressure head dependent or moisture content dependent form of Richard's equation. The model code uses isoparametric quadratic triangular and/or quadrilateral finite elements for the geometric representation and for the weak Galerkin spacial integrations. An implicit, unconditionally stable single-step numerical time integration scheme with an oscilliatory noise reduction option is utilized for the temporal discretization. The highly efficient symmetric skyline (profile) solution scheme is used to solve the resulting simultaneous equations. The nonlinear subsurface flow parameters are approximated using cubic spline interpolation. The element …