The Elder problem with reactive infiltration effects

0580180 - GFÚ 2024 RIV DE eng J - Journal Article
Hurtiš, J. - Guba, P. - Kyselica, Juraj
The Elder problem with reactive infiltration effects.
Transport in Porous Media. Roč. 150, č. 3 (2023), s. 627-652. ISSN 0169-3913. E-ISSN 1573-1634
Institutional support: RVO:67985530
Keywords : Elder problem * buoyancy-driven convection * porous medium * reactive infiltration * oscillatory convection
OECD category: Applied mathematics
Impact factor: 2.7, year: 2022
Method of publishing: Open access
https://link.springer.com/article/10.1007/s11242-023-02026-z

The occurrence of buoyancy-driven flow and reactive solute transport in a fluid-saturated porous medium can be induced by either natural processes or human activities. Typical examples include the groundwater salinization in carbonate-rock aquifers and the acid treatment of oil wells in petroleum drilling industry. In this paper, the classical Elder problem of buoyancy-driven convection in two-dimensional porous media is extended to include the local chemical interactions between the solute in the pore liquid (e.g. salt such as NaCl or acids such as HCl and HCl/HF mixtures) and the solid space of the porous medium (e.g. minerals such as calcite and dolomite). Effects of the geochemical processes on the flow and mass transport are investigated. For the reactive, strongly solute-driven case in the regime dominated by the diffusive mass transport, a decrease in the net solute concentration is found as compared to the non-reactive case. This decrease is pronounced at higher values of the Damkohler number when the solute reaction rate becomes larger than the solute diffusion rate. Furthermore, the flow structure is affected by products generated by the chemical reaction when the Rayleigh number for the products is sufficiently high. In that case, numerical simulations show the formation of diluted fluid tongues exhibiting damped periodic oscillations about a nonzero mean. The results are obtained using the pseudospectral numerical method, verified against the analytical solution for the non-reactive, purely diffusive case.
Permanent Link: https://hdl.handle.net/11104/0348942