scientific article; zbMATH DE number 1065008

From MaRDI portal
Publication:4356307





The authors develop a cell-centred finite difference method as an approximation to an expanded mixed finite element method which is found suitable for groundwater flows and transport problems in general geometry with tensor permeabilities or dispersivities. The method is both locally (element by element) mass conservative and highly accurate. General geometry is handled by a mapping between the computational and physical domain. By using a rectangular mesh, the authors apply special quadrature rules to the expanded finite element method to transform it into a cell-centred finite difference method. The solution of the logically rectangular cell-centred finite difference scheme converges to the true solution at the optimal order as the mesh is refined. Moreover, super-convergence is attained by the velocity at certain discrete points of the domain on \(C^2\)-smooth grids away from the boundary. It is observed that grid geometry strongly affects the approximation error, so it is necessary to define and refine grids in a J\(C^2\)-smooth manner. The authors introduce Lagrangian multiplier pressures or infinitely thin cells along faces where tensor geometry changes discontinuously. This allows the adjusted pressure gradient \(\widetilde u\) to be discontinuous, so that the coefficient of \(\widetilde u\) approximates well the continuous velocity. Mixed methods accurately approximate highly nonlinear two-phase flow. The main difficulty is found to be that the true solutions are typically lacking in regularity. The mixed method presented herein is new and of great interest to these who are working on `subsurface simulation'. For the entire collection see [Zbl 0869.00055].











This page was built for publication:

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4356307)