The analytic element method is a numerical method used for the solution of partial differential equations. It was initially developed by O.D.L. Strack at the University of Minnesota. It is similar in nature to the boundary element method, as it does not rely upon discretization of volumes or areas in the modeled system; only internal and external boundaries are discretized. One of the primary distinctions between AEM and BEMs is that the boundary integrals are calculated analytically.
Mathematical basis
The basic premise of the analytic element method is that, for linear differential equations, elementary solutions may be superimposed to obtain more complex solutions. A suite of 2D and 3D analytic solutions are available for different governing equations. These elements typically correspond to a discontinuity in the dependent variable or its gradient along a geometric boundary. This discontinuity has a specific functional form and may be manipulated to satisfy Dirichlet, Neumann, or Robin boundary conditions. Each analytic solution is infinite in space and/or time. Commonly each analytic solution contains degrees of freedom that may be calculated to meet prescribed boundary conditions along the element's border. To obtain a global solution, a system of equations is solved such that the boundary conditions are satisfied along all of the elements. Notably, the global solution provides a spatially continuous description of the dependent variable everywhere in the infinite domain, and the governing equation is satisfied everywhere exactly except along the border of the element, where the governing equation is not strictly applicable due to the discontinuity. The ability to superpose numerous elements in a single solution means that analytical solutions can be realized for arbitrarily complex boundary conditions. That is, models that have complex geometries, straight or curved boundaries, multiple boundaries, transient boundary conditions, multiple aquifer layers, piecewise varying properties and continuously varying properties can be solved. Elements can be implemented using far-field expansions such that model containing many thousands of elements can be solved efficiently to high precision. The analytic element method has been applied to problems of groundwater flow governed by a variety of linear partial differential equations including the Laplace, the Poisson equation, the modified Helmholtz equation, the heat equation, and the biharmonic equations. Often theses equations are solved using complex variables which enables using mathematical techniques available in complex variable theory. A useful technique to solve complex problems is using conformal mapping which maps the boundary of a geometry, e.g. an ellipse, onto the boundary of the unit circle where the solution is known. In the analytic element method the discharge potential and stream function, or combined the complex potential, is used. This potential links the physical properties of the groundwater system, the hydraulic head or flow boundaries, to a mathematical representation in a potential. This mathematical representation can be used to calculate the potential in terms of position and thus also solving groundwater flow problems. Elements are developed by solving the boundary conditions for either of these two properties, hydraulic head or flow boundary, which results in analytical solutions capable of dealing with numerous boundary conditions. A contemporary student of Strack's who is a proponent of the Analytic Element Method in groundwater modeling applications is Dr. David Steward of Kansas State University.
Comparison to other methods
As mentioned the analytic element method thus not rely on discretization of volume or area in the model, as in the finite elements or finite different methods. Thus, it can model complex problem with an error on the order ofmachine precision. This is illustrated in a study which modeled a highly heterogeneous, isotropic aquifer by including 100,000 spherical heterogeneity with a random conductivity and tracing 40,000 particles. The analytical element method can efficiently be used as verification or as a screening tool in larger projects as it may fast and accurately calculate the groundwater flow for many complex problems. In contrast to other commonly used groundwater modeling methods, e.g. the finite elements or finite different method, the AEM does not discrete the model domain into cells. This gives the advantage that the model is valid for any given point in the model domain. However, it also imposes that the domain is not as easily divided into regions of e.g. different hydraulic conductivity, as when modeling with a cell grid. While there are some solutions that deals with this, e.g. there exist solutions for implementing a vertically varying properties or structures in an aquifer in an AEM model.
