Publication Details

Category Text Publication
Reference Category Journals
DOI 10.1002/2015WR017164
Title (Primary) Numerical stability analysis of two-dimensional solute transport along a discrete fracture in a porous rock matrix
Author Watanabe, N.; Kolditz, O. ORCID logo
Source Titel Water Resources Research
Year 2015
Department ENVINF
Volume 51
Issue 7
Page From 5855
Page To 5868
Language englisch
UFZ wide themes RU5;
Abstract

This work reports numerical stability conditions in two-dimensional solute transport simulations including discrete fractures surrounded by an impermeable rock matrix. We use an advective-dispersive problem described in Tang et al. (1981) and examine the stability of the Crank-Nicolson Galerkin finite element method (CN-GFEM). The stability conditions are analyzed in terms of the spatial discretization length perpendicular to the fracture, the flow velocity, the diffusion coefficient, the matrix porosity, the fracture aperture, and the fracture longitudinal dispersivity. In addition, we verify applicability of the recently developed finite element method-flux corrected transport (FEM-FCT) method by Kuzmin (2009) to suppress oscillations in the hybrid system, with a comparison to the commonly utilized Streamline Upwinding/Petrov-Galerkin (SUPG) method. Major findings of this study are (1) the mesh von Neumann number (Fo) ≥ 0.373 must be satisfied to avoid undershooting in the matrix, (2) in addition to an upper bound, the Courant number also has a lower bound in the fracture in cases of low dispersivity, and (3) the FEM-FCT method can effectively suppress the oscillations in both the fracture and the matrix. The results imply that, in cases of low dispersivity, prerefinement of a numerical mesh is not sufficient to avoid the instability in the hybrid system if a problem involves evolutionary flow fields and dynamic material parameters. Applying the FEM-FCT method to such problems is recommended if negative concentrations cannot be tolerated and computing time is not a strong issue.

Persistent UFZ Identifier https://www.ufz.de/index.php?en=20939&ufzPublicationIdentifier=16516
Watanabe, N., Kolditz, O. (2015):
Numerical stability analysis of two-dimensional solute transport along a discrete fracture in a porous rock matrix
Water Resour. Res. 51 (7), 5855 - 5868 10.1002/2015WR017164