Premium
Numerical stability of the BEM for advection‐diffusion problems
Author(s) -
Peratta Andrés,
Popov Viktor
Publication year - 2004
Publication title -
numerical methods for partial differential equations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.901
H-Index - 61
eISSN - 1098-2426
pISSN - 0749-159X
DOI - 10.1002/num.20009
Subject(s) - discretization , boundary element method , hydraulic conductivity , advection , mathematics , matrix (chemical analysis) , fracture (geology) , stability (learning theory) , diffusion , boundary (topology) , partial differential equation , flow (mathematics) , porous medium , mathematical analysis , finite element method , boundary value problem , conductivity , porosity , geometry , geotechnical engineering , geology , computer science , materials science , engineering , physics , structural engineering , machine learning , soil science , soil water , thermodynamics , quantum mechanics , composite material
Abstract A boundary element method (BEM) approach has been developed to solve the time‐dependent 1D advection‐diffusion equation. The 1D solution is part of a 3D numerical scheme for solving advection‐diffusion (AD) problems in fractured porous media. The full 3D scheme includes a 3D solution for the porous matrix, which is coupled with a 2D solution for fractures and a 1D solution for fracture intersections. As the hydraulic conductivity of the fracture intersections is usually higher than the hydraulic conductivity of the fractures and by at least one order of magnitude higher than the hydraulic conductivity of the porous matrix, the fastest flow and solute transport occurs in the fracture intersections. Therefore it is important to have an accurate and stable 1D solution of the transient AD problems. This article presents two different 1D BEM formulations for solution of the AD problems. The particular advantage of these formulations is that they provide one of the most straightforward and simplest ways to couple multiple intersecting 2D Boundary Element problems discretized with linear discontinuous elements. Both formulations are tested and compared for accuracy, stability, and consistency. The analysis helps to select the more suitable formulations according to the properties of the problem under consideration. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004