AN ADAPTIVE DOMAIN DECOMPOSITION METHOD FOR SIMULATION OF TRANSPORT IN POROUS-MEDIA

A Galerkin finite element method is used along with a self-adaptive strategy of domain discretisation to model dispersion in an axisymmetric cylindrical porous medium. A solution strategy is proposed based on the use of a Gear scheme for the time stepping and partial vectorisation of the code. The domain is highly discretised in the area of the sharp transient front, while the remainder is coarsely discretised. The area covered by the fine mesh is determined by the value of the local concentration gradients. Numerical results are presented for the one and two dimensional cases.