Ph.D. projects


Improved prediction of fluid-solid reaction kinetics using multiple numerical approaches

Jonas Schabernack
Eu atoms (green) at an etch pit surface on the
base of muscovite (kMC simulation)
Ph.D. student:

Jonas Schabernack

Supervisor:

PD Dr. Cornelius Fischer (HZDR)

Division:

Reactive Transport

Duration:

02/2020–01/2023

Recent numerical investigations revealed that the heterogeneity of the dissolution rate observed in numerous experiments cannot be explained by fluid transport effects. This heterogeneity is attributed to intrinsic surface reactivity. Therefore, reactive transport models (RTM) require parameterization of the surface reactivity for accurate predictions. For this purpose, a nanotopographic parametrization based on surface slope has been recently suggested. In this study, we utilize and improve this parametrization for RTMs of pore-scale systems, from the crystal surface to the single crystal geometry, going beyond the previous reactivity parametrization. 2D and 3D RTMs were developed using COMSOL Multiphysics for calcite systems based on experimental measurements. We compared the results between classically parameterized RTMs, RTMs with new slope parameterization, and experimental data. The effect of flow on dissolution under conditions far-from-equilibrium is found to be negligible, highlighting the importance of surface reactivity in the dissolution reaction. For the first time, the new slope factor was able to accurately reproduce the experimental results on a crystal surface with large field-of-view, large height variability of the topography, and over a long-term reaction period. The new parameterization had greatly improved sensitivity for intermediate reactivity ranges compared to the previous parameterization. A 3D model is used to present the general applicability of the parameterization for use in realistic geometric data sets. Thus, we also show that neglecting surface reactivity in an RTM leads to incorrect predictions regarding the porosity, pore geometry, and surface topography of the system. Our new slope factor can successfully serve as a first-order proxy for the distribution of surface reactivity in 3D pore-scale rock systems. The description of surface reactivity is crucial for accurate long-term modeling of natural rock systems. For details see: https://doi.org/10.1016/j.gca.2022.08.003

In deep geological repositories for nuclear waste, the surrounding rock formation serves as an important barrier against radionuclide migration. Multiple potential host rocks contain phyllosilicates, which have shown high efficiency in radionuclide sorption. Recent experimental studies report a heterogeneous distribution of adsorbed radionuclides on nanotopographic mineral surfaces. In this study, the energetic differences of surface sorption sites available at nanotopographic structures such as steps, pits, and terraces are investigated. Eleven important surface sites are selected and the energies of ad- and desorption reactions are obtained from density functional theory calculations. The adsorption energies are then used for the parameterization of a kinetic Monte Carlo model simulating the distribution of adsorbed europium on a typical nanotopographic muscovite surface. On muscovite, silicon step sites are favorable for europium sorption and lead to an increased adsorption in regions with high step concentrations. Under identical chemical conditions, sorption on typical nanotopographic surfaces is increased by a factor of three compared to atomically flat surfaces. Desorption occurs preferentially at terrace sites, leading to an overall 2.5 times increased retention at nanotopographic structures. This study provides a mechanistic explanation for heterogeneous sorption on nanotopographic mineral surfaces due to the availability of energetically favorable sorption sites.

For details see: https://doi.org/10.1002/adts.202300406