Premium
Application of the Lanczos Algorithm to the solution of the groundwater flow equation
Author(s) -
Dunbar W. Scott,
Woodbury Allan D.
Publication year - 1989
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/wr025i003p00551
Subject(s) - tridiagonal matrix , lanczos resampling , discretization , finite element method , lanczos algorithm , mathematics , groundwater flow equation , matrix (chemical analysis) , groundwater flow , flow (mathematics) , finite difference method , partial differential equation , finite difference , algorithm , mathematical optimization , mathematical analysis , aquifer , groundwater , geometry , eigenvalues and eigenvectors , materials science , physics , geotechnical engineering , quantum mechanics , composite material , thermodynamics , engineering
The solution to the finite element matrix‐differential equations resulting from the discretization of the groundwater flow equation is normally carried out by a finite difference approximation to the time derivative. The total computational effort in solving a fluid flow problem is then directly related to the number of unknowns and the number of time steps required to obtain accurate and stable solutions. The Lanczos algorithm uses an orthogonal matrix transformation to reduce the finite element equations to a much smaller tridiagonal system of first‐order differential equations. This new system can be solved by a standard tridiagonal solution algorithm with very little computational effort. A matrix‐vector multiplication is then used to obtain the original solution at desired time steps. The algorithm is used to accurately simulate the drawdown of synthetic two‐dimensional aquifers, including ones with substantial hydraulic conductivity contrasts. The method affords an efficient means of solving large problems, particularly when time durations are long.