On the behavior of approaches to simulate reactive transport

被引:82
作者
Saaltink, MW
Carrera, J
Ayora, C
机构
[1] Univ Politecn Cataluna, ETSECCPB, Dept Engn Terreny & Cartog, ES-08034 Barcelona, Spain
[2] CSIC, Inst Ciencies Terra Jaume Almera, E-08028 Barcelona, Spain
关键词
hydrogeochemistry; geochemistry; modeling; numerical methods;
D O I
10.1016/S0169-7722(00)00172-8
中图分类号
X [环境科学、安全科学];
学科分类号
08 ; 0830 ;
摘要
Two families of approaches exist to simulate reactive transport in groundwater: The Direct Substitution Approach (DSA). based on Newton-Raphson and the Picard or Sequential Iteration Approach (SIA). We applied basic versions of both methods to several test casts and compared both computational demands and quality of the solution fur varying grid size. Results showed that the behavior of the two approaches is sensitive to both grid size and chemistry. As a general rule, the DSA is more robust than the SIA. in the sense that its convergence is less sensitive to time step size (any approach will converge given a sufficiently small time step). Moreover. the DSA leads to a better simulation of sharp fronts, which can only he reproduced with fine grids after many iterations when the SIA is used. As a consequence, the DSA runs faster than SIA in chemically difficult cases (i.e.. highly non-linear and/or very retarded), because the SIA may require very small time steps to converge. On the other hand, the size of the system of equations is. much larger for the DSA than for the SIA. so that its CPU time and memory requirements tend to be less favorable with increasing grid size. As a result, the SIA may become faster than the DSA for very large. chemically simple problems. The use of an iterative linear solver for the DSA makes its CPU time less sensitive to grid size. (C) 2001 Elsevier Science B.V. All rights reserved.
引用
收藏
页码:213 / 235
页数:23
相关论文
共 32 条
[1]  
Appelo C. A. J., 1993, Geochemistry
[2]   GEOCHEMICAL CALCULATIONS AND OBSERVATIONS ON SALT-WATER INTRUSIONS .1. A COMBINED GEOCHEMICAL MIXING CELL MODEL [J].
APPELO, CAJ ;
WILLEMSEN, A .
JOURNAL OF HYDROLOGY, 1987, 94 (3-4) :313-330
[3]   The genesis of dedolomites: a discussion based on reactive transport modeling [J].
Ayora, C ;
Taberner, C ;
Saaltink, MW ;
Carrera, J .
JOURNAL OF HYDROLOGY, 1998, 209 (1-4) :346-365
[4]   Alternative split-operator approach for solving chemical reaction groundwater transport models [J].
Barry, DA ;
Bajracharya, K ;
Miller, CT .
ADVANCES IN WATER RESOURCES, 1996, 19 (05) :261-275
[5]   CHEMFRONTS:: a coupled geochemical and transport simulation tool [J].
Bäverman, C ;
Strömberg, B ;
Moreno, L ;
Neretnieks, I .
JOURNAL OF CONTAMINANT HYDROLOGY, 1999, 36 (3-4) :333-351
[6]  
DUFF IS, 1977, MA28 SET FORTRAN SUB
[7]   Contaminant transport at a waste residue deposit .2. Geochemical transport modeling [J].
Engesgaard, P ;
Traberg, R .
WATER RESOURCES RESEARCH, 1996, 32 (04) :939-951
[8]   A GEOCHEMICAL TRANSPORT MODEL FOR REDOX-CONTROLLED MOVEMENT OF MINERAL FRONTS IN GROUNDWATER-FLOW SYSTEMS - A CASE OF NITRATE REMOVAL BY OXIDATION OF PYRITE [J].
ENGESGAARD, P ;
KIPP, KL .
WATER RESOURCES RESEARCH, 1992, 28 (10) :2829-2843
[9]   Reactive chemical transport within engineered barriers [J].
Grindrod, P ;
Takase, H .
JOURNAL OF CONTAMINANT HYDROLOGY, 1996, 21 (1-4) :283-296
[10]   Coupled reaction-flow modeling of diagenetic changes in reservoir permeability, porosity and mineral compositions [J].
Le Gallo, Y ;
Bildstein, O ;
Brosse, E .
JOURNAL OF HYDROLOGY, 1998, 209 (1-4) :366-388