混凝土碳化数值模拟与COMSOL工程应用

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

混凝土碳化是影响建筑结构耐久性的关键因素之一。当环境中的二氧化碳渗透进混凝土孔隙，与水泥水化产物氢氧化钙发生化学反应，生成碳酸钙和水，这个过程我们称之为碳化反应。从微观上看，这个反应会改变混凝土的孔隙结构和碱性环境，直接导致内部钢筋失去钝化膜保护，进而引发锈蚀膨胀和混凝土开裂。

在实际工程中，我们评估一座桥梁或建筑的剩余使用寿命时，碳化深度往往是首要测量指标。以某滨海城市高架桥检测数据为例，浪溅区的桥墩在服役15年后碳化深度达到保护层厚度的60%，这意味着钢筋锈蚀风险已经进入预警阶段。因此，准确预测碳化进程对制定科学的维护策略至关重要。

传统实验室测定碳化深度的方法存在明显局限：

• 需要钻芯取样破坏结构
• 加速碳化试验与实际情况存在差异
• 无法反映复杂环境下的动态变化

这正是我们需要建立数值模型的原因。通过COMSOL Multiphysics这类多物理场仿真平台，可以构建包含气体扩散、化学反应、湿度传导等多重机制的耦合模型，实现从实验室数据到工程预测的跨越。

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

2.1 碳化反应动力学机制

碳化过程本质上是气-液-固三相反应体系。二氧化碳首先通过混凝土孔隙扩散，溶解在孔隙液中形成碳酸，继而与固相氢氧化钙反应。这个过程的反应速率受三个关键因素控制：

1. 扩散控制阶段：CO₂在孔隙网络中的传输

   • 遵循Fick第二定律
   • 有效扩散系数Deff与孔隙率关系：Deff=D0·ε^1.8
     （其中D0为自由空气扩散系数，ε为孔隙率）

2. 溶解平衡阶段：CO₂在孔隙液中的溶解

   • 亨利定律：Caq=KH·Pgas
   • 温度修正系数：KH(T)=KH(298K)·exp[2400(1/T-1/298)]

3. 表面反应阶段：Ca(OH)₂的消耗

   • 反应速率方程：r=k·Cca·Cco2
   • 反应速率常数k与温度关系：k=A·exp(-Ea/RT)

2.2 控制方程建立

在COMSOL中我们需要建立耦合的偏微分方程组：

质量守恒方程（CO₂传输）：
∂(εCco2)/∂t = ∇·(Deff∇Cco2) - R

化学平衡方程（Ca(OH)₂消耗）：
∂Cca/∂t = -R·Mca/ρca

其中反应源项R采用以下