POLYMIN - University of Waterloo

POLYMIN 2005
POLYMIN has been developed with some of the most efficient and powerful numerical algorithms
available. Finite elements are employed for high accuracy and to allow deformable domain
geometry, and the Leismann time-weighting scheme (Leismann and Frind, 1989) has been
incorporated to maintain matrix symmetry which saves memory and execution time. The model
includes one of the most efficient preconditioned conjugate gradient solvers available to solve the
matrix equations.
Custom versions, and executables are available to suit individual needs. Please contact the authors
for available updates. A more recent model, BIONAPL/3D, can also simulate density flow along
with multiple-component NAPL dissolution and reactive transport (Molson, 2001). The POLYMIN
model has most recently been applied by Molson et al. (2005) to waste rock piles, based on flow
systems developed by Fala et al. (2005).
1.2 Capabilities, Assumptions and Limitations
The capabilities of POLYMIN, and its major assumptions and limitations are summarized in Table
1.
Table 1. Summary of the capabilities, assumptions and limitations of POLYMIN.
Capabilities
- 2D or 1D domains.
- fully coupled oxygen diffusion, geochemical
speciation and advective-dispersive mass
transport.
- domain can be heterogeneous and
anisotropic.
- deformable elements can conform to
complex geometry.
- versatile boundary condition options.
- computes concentration breakthrough data
(conc. vs time) at selected points.
Assumptions and Limitations
- non-deforming, isothermal aquifer
- fluid is incompressible.
- chemical reactions are assumed at equilibrium
- flow system read separately
- non-fractured or equivalent porous media (2)
(2)
Research versions also available with 1D and 2D fractures (Yang et al, 1996a,b).
6