OpenGeoSys-Gem: A numerical tool for calculating geochemical and porosity changes in saturated and partially saturated media |
| |
| Institution: | 1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China;2. The key laboratory of safety for geotechnical and structural engineering of Hubei province, School of Civil Engineering, Wuhan University, Wuhan 430072, China |
| |
| Abstract: | 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. |
| |
| Keywords: | Reactive transport Finite elements Flux correction scheme Parallelization Porosity change |
| 本文献已被 ScienceDirect 等数据库收录！ |
|
