COMSOL多物理场建模在混凝土碳化预测中的应用

1. 混凝土碳化现象与工程意义

混凝土碳化是钢筋混凝土结构耐久性研究中的经典课题。当环境中的二氧化碳渗透进混凝土孔隙，与水泥水化产物氢氧化钙发生化学反应，生成碳酸钙和水，这个过程会降低混凝土的碱性。当碳化前沿到达钢筋表面时，钢筋表面的钝化膜遭到破坏，在氧气和水分共同作用下就会引发钢筋锈蚀。

我在参与某跨海大桥耐久性评估项目时，曾亲眼见证过碳化深度超标导致的保护层剥落案例。当时在浪溅区取芯检测发现，设计寿命50年的结构在服役12年后碳化深度已达38mm，远超25mm的保护层厚度。这个案例让我深刻意识到碳化预测对工程维护的重要性。

传统碳化深度预测多采用经验公式法，如欧洲规范推荐的fib Model Code公式。但这类方法无法反映环境湿度、温度、CO2浓度等多场耦合效应。而基于COMSOL的多物理场建模，可以构建更精确的碳化反应-扩散耦合模型，这对重要基础设施的寿命预测具有突破性意义。

2. 模型理论基础与方程推导

2.1 碳化反应动力学机制

混凝土碳化的核心是气-液-固三相反应：

1. 气相CO2通过孔隙扩散进入混凝土
2. CO2溶解于孔隙液中形成H2CO3
3. H2CO3与Ca(OH)2发生中和反应

其反应速率遵循二级动力学方程：
$$
r = k \cdot [CO_2] \cdot [OH^-]
$$
其中速率常数k与温度关系符合阿伦尼乌斯公式：
$$
k = A e^{-E_a/RT}
$$

在COMSOL中，我们需要将这个非线性反应源项耦合到物质传递方程中。实测数据显示，普通硅酸盐水泥的k值在20℃时约为3.5×10^-6 m/s，而矿渣水泥的k值会低1-2个数量级。

2.2 控制方程建立

完整的碳化模型需要求解以下耦合方程组：

1. 物质守恒方程（CO2传输）：
   $$
   \frac{\partial C_{CO_2}}{\partial t} = \nabla \cdot (D_{eff} \nabla C_{CO_2}) - r
   $$
   其中有效扩散系数D_eff需考虑孔隙率ε和曲折度τ：
   $$
   D_{eff} = D_0 \cdot ε^{1.8} / τ
   $$

2. 反应进度变量方程：
   引入碳化程度η（0≤η≤1）：
   $$
   \frac{\partial η}{\partial t}