Reactive transport modelling has wide applications in geosciences. In the field of hydrogeology, it has been utilised to simulate the biogeochemical processes that disperse and degrade contaminants in the aquifer. For geotechnical applications, such as geological CO2 sequestration, the reaction of CO2 with the ambient saline aquifer determines the final success of storage. In a radioactive waste repository, scientists rely on reactive transport models to predict the mobilisation of hazardous radionuclides within space and time. In this work, the multi-component mass transport code OpenGeoSys, was coupled with two geochemical solvers, the Gibbs Energy Minimization Selektor (GEM) and the Biogeochemical Reaction Network Simulator (BRNS). Both coupled codes were verified against analytical solutions and simulation results from other numerical models. Moreover, the coupling interface was developed for parallel simulation. Test runs showed that the speed-up of reaction part had a very good linearity with number of nodes in the mesh. However, for three dimensional problems with complex geochemical reactions, the model performance was dominated by solving transport equations of mobile chemical components. OpenGeoSys-BRNS was applied to a two dimensional groundwater remediation problem. Its calculated concentration profiles fitted very well with analytical solutions and numerical results from TBC. The model revealed that natural attenuation of groundwater contaminants is mainly controlled by the mixing of carbon source and electron donor. OpenGeoSys-GEM was employed to investigate the retardation mechanism of radionuclides in the near field of a nuclear waste repository. Radium profiles in an idealised bentonite column was modelled with varying clay/water ratios. When clay content is limited, Ba-Sr-Ra sulfate solid solutions have a very strong retardation effect on the aqueous radium. Nevertheless, when clay mineral is abundant, cation exchange sites also attract Sr and Ba, thus dominates the transport of Ra.

Reactive transport modelling has wide applications in geosciences. In the field of hydrogeology, it has been utilised to simulate the biogeochemical processes that disperse and degrade contaminants in the aquifer. For geotechnical applications, such as geological CO2 sequestration, the reaction of CO2 with the ambient saline aquifer determines the final success of storage. In a radioactive waste repository, scientists rely on reactive transport models to predict the mobilisation of hazardous radionuclides within space and time. In this work, the multi-component mass transport code OpenGeoSys, was coupled with two geochemical solvers, the Gibbs Energy Minimization Selektor (GEM) and the Biogeochemical Reaction Network Simulator (BRNS). Both coupled codes were verified against analytical solutions and simulation results from other numerical models. Moreover, the coupling interface was developed for parallel simulation. Test runs showed that the speed-up of reaction part had a very good linearity with number of nodes in the mesh. However, for three dimensional problems with complex geochemical reactions, the model performance was dominated by solving transport equations of mobile chemical components. OpenGeoSys-BRNS was applied to a two dimensional groundwater remediation problem. Its calculated concentration profiles fitted very well with analytical solutions and numerical results from TBC. The model revealed that natural attenuation of groundwater contaminants is mainly controlled by the mixing of carbon source and electron donor. OpenGeoSys-GEM was employed to investigate the retardation mechanism of radionuclides in the near field of a nuclear waste repository. Radium profiles in an idealised bentonite column was modelled with varying clay/water ratios. When clay content is limited, Ba-Sr-Ra sulfate solid solutions have a very strong retardation effect on the aqueous radium. Nevertheless, when clay mineral is abundant, cation exchange sites also attract Sr and Ba, thus dominates the transport of Ra.