Semi-discrete finite difference multiscale scheme for a concrete corrosion model: approximation estimates and convergence

Vladimír Chalupecký,Adrian Muntean
DOI: https://doi.org/10.1007/s13160-012-0060-6
2011-06-24
Abstract:We propose a semi-discrete finite difference multiscale scheme for a concrete corrosion model consisting of a system of two-scale reaction-diffusion equations coupled with an ode. We prove energy and regularity estimates and use them to get the necessary compactness of the approximation estimates. Finally, we illustrate numerically the behavior of the two-scale finite difference approximation of the weak solution.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the numerical approximation and convergence of the concrete corrosion model. Specifically, the author proposes a semi - discrete finite - difference multiscale scheme for simulating a concrete corrosion model coupled by a two - scale reaction - diffusion equation and an ordinary differential equation (ODE). This model aims to study the influence of microstructure effects on macroscopic observables and verify whether the proposed two - scale finite - difference scheme converges, that is, whether it can approximate the weak solution of the two - scale system. ### Background of the Paper and Problem Statement 1. **Problem Background**: - Concrete corrosion is a complex biochemical process. Especially in the wastewater environment, bacteria - mediated hydrogen sulfide gas is formed and converted into sulfuric acid, which then attacks concrete and steel bars. - Sulfuric acid has a devastating impact on the surfaces of concrete and steel exposed to the harsh wastewater environment, resulting in high maintenance costs. 2. **Mathematical Modeling**: - The model describes the process of sulfate attacking concrete, involving a one - dimensional macroscopic region $\Omega=(0, L)$ and a typical microstructure $Y = (0,\ell)$. - The model consists of reaction - diffusion equations at two scales and an ordinary differential equation, which describe the changes in the concentrations of different chemical substances over time and space. 3. **Main Problems**: - Verify whether the proposed two - scale finite - difference scheme converges, that is, whether it can approximate the weak solution of the two - scale system. - Specifically, the author focuses on how to approximate the weak solution in a simple and efficient way while maintaining the regularity of the model structure, data, and parameters. ### Mathematical Formula Representation The main equations involved in the model are as follows: - Macroscopic - scale equation (change in hydrogen sulfide concentration $u_1$): \[ \frac{\partial u_1}{\partial t}-d_1\frac{\partial^2 u_1}{\partial x^2}=d_2\left.\frac{\partial u_2}{\partial y}\right|_{y = 0},\quad\text{in }\Omega \] - Microscopic - scale equations (changes in dissolved hydrogen sulfide concentration $u_2$ and sulfuric acid concentration $u_3$): \[ \frac{\partial u_2}{\partial t}-d_2\frac{\partial^2 u_2}{\partial y^2}=-\zeta(u_2, u_3),\quad\text{in }\Omega\times Y \] \[ \frac{\partial u_3}{\partial t}-d_3\frac{\partial^2 u_3}{\partial y^2}=\zeta(u_2, u_3),\quad\text{in }\Omega\times Y \] - Gypsum generation rate (change in gypsum concentration $u_4$): \[ \frac{\partial u_4}{\partial t}=\eta(u_3|_{y=\ell}, u_4),\quad\text{in }\Omega \] Boundary conditions and initial conditions also describe the behavior of these equations in detail. ### Summary The core problem of the paper is to verify whether the proposed two - scale finite - difference scheme can converge to the weak solution, and the effectiveness of this method is demonstrated through numerical experiments. This not only helps to understand the microscopic mechanism of concrete corrosion but also provides theoretical support for practical engineering applications.