Premium
Solving the steady‐state groundwater flow equation for finite linear aquifers using a generalized Fourier series approach in two‐dimensional domains
Author(s) -
Capilla Jose E.,
PulidoVelazquez David,
Sahuquillo Andrés,
Andreu Joaquín
Publication year - 2009
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.2552
Subject(s) - fourier series , groundwater flow equation , fourier transform , mathematics , mathematical analysis , discrete fourier series , spectral method , flow (mathematics) , finite element method , fourier transform on finite groups , matrix (chemical analysis) , series (stratigraphy) , hydraulic head , fourier analysis , groundwater flow , geometry , fractional fourier transform , aquifer , physics , groundwater , geotechnical engineering , materials science , biology , composite material , thermodynamics , engineering , paleontology
A method to solve steady linear groundwater flow problems using generalized Fourier Series is developed and particularized for multiple Fourier series in two‐dimensional domains. It leads to a linear vector equation whose solution provides a finite number of generalized Fourier coefficients approximating the hydraulic head field. Its implementation is shown and two relevant properties are found for the system matrix. It is always symmetric and, once computed, if additional Fourier terms are needed for a better approximation of the hydraulic head field, previously computed matrix elements remain invariant, i.e. only new rows and columns are added to the system matrix. The method is demonstrated in three simple cases with different geometries and transmissivity fields, where solutions are compared with analytical and finite element method results. Thus, the method is verified as an alternative to other flow solvers. Additionally, it provides a direct way to obtain the spectral form of the flow equation solution, given a spectral representation of transmissivity, and can be easily extended to obtain continuous velocity fields and their approximated spectral expressions. Copyright © 2009 John Wiley & Sons, Ltd.