期刊
PHYSICS AND CHEMISTRY OF THE EARTH
卷 70-71, 期 -, 页码 138-149出版社
PERGAMON-ELSEVIER SCIENCE LTD
DOI: 10.1016/j.pce.2013.11.008
关键词
Reactive transport; Finite elements; Flux correction scheme; Parallelization; Porosity change
资金
- Swiss National Cooperative for the Disposal of Radioactive Waste (NAGRA)
- German Federal Ministry of Education and Research (BMBF) [03E10588 (A-DuR), 03G0797D (CO2BENCH)]
- Helmholtz Association within the POF Renewable Energy program [CO2BENCH]
- Swiss National Cooperative for the Disposal of Radioactive Waste (NAGRA)
- German Federal Ministry of Education and Research (BMBF) [03E10588 (A-DuR), 03G0797D (CO2BENCH)]
- Helmholtz Association within the POF Renewable Energy program [CO2BENCH]
Reactive transport codes that use a Gibbs Energy Minimization (GEM) to solve chemical equilibria are uncommon. We present a new coupling of the Richards flow module of the Finite Element (FE) based OpenGeoSys code with the GEM based chemical solver GEMS3K. The coupled code is highly parallelized using an overlapping domain decomposition approach in combination with execution of multiple threads that solve chemical equilibria in parallel. FE reactive transport schemes are often affected by spurious concentration oscillations. We effectively suppress these oscillations with a linearized algebraic flux corrected transport (FCT) algorithm. An application example is presented which investigates the evolution of material interfaces in a deep geological repository for nuclear waste. The example uses all features of the new coupled code: flow and multi-component transport in variably saturated media, and a very complex chemical setup which makes extensive use of (non-linear) solid solution formulations for mineral phases. (C) 2013 Elsevier Ltd. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据