Premium
Numerical stability analysis of two‐dimensional solute transport along a discrete fracture in a porous rock matrix
Author(s) -
Watanabe Norihiro,
Kolditz Olaf
Publication year - 2015
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1002/2015wr017164
Subject(s) - finite element method , discretization , matrix (chemical analysis) , fracture (geology) , discontinuous galerkin method , instability , mechanics , stability (learning theory) , flow (mathematics) , galerkin method , mathematics , geotechnical engineering , materials science , mathematical analysis , geology , structural engineering , physics , computer science , engineering , composite material , machine learning
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 ([Kuzmin, D., 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.