COMSOL混凝土碳化模型

COMSOL 混凝土碳化模型

混凝土表面那层白霜般的碳化层，总让我想起实验室里放了三个月的苏打饼干。这层碳酸钙的生成过程，在COMSOL里建模就像在模拟一场微观世界的化学舞会——二氧化碳分子突破混凝土保护层，与氢氧化钙在孔隙溶液里跳起离子交换的华尔兹。

打开COMSOL新建模型时，先得在全局参数里塞进几个关键数值：

% 材料参数 CaOH2_initial = 0.85; // 初始氢氧化钙浓度 mol/m³ CO2_ambient = 0.03; // 环境CO2体积分数 diffusivity = 1.2e-7; // CO2有效扩散系数 m²/s

这几个数字就像调酒师的秘方，稍微改动就会让整个反应过程变味。扩散系数尤其敏感，试过把1.2e-7改成1.5e-7，碳化前沿推进速度直接快了20%，就像往咖啡里多加了冰块。

化学反应模块的设置是个精细活：

// 碳化反应动力学 rate = k_forward*CO2*CaOH2 - k_backward*CaCO3*H2O;

这里的正向速率常数k_forward需要配合实验数据反复校准。有次偷懒直接用了文献值，结果模拟的碳化深度比实测值多出5mm，被导师笑称建了个"科幻混凝土模型"。

COMSOL 混凝土碳化模型

多物理场耦合才是重头戏。扩散方程和化学反应就像在玩跷跷板：

// CO2传输方程 transport = diffusiveFlux(-D*grad(CO2)) + sourceTerm(-rate); // 氢氧化钙守恒 dCaOH2_dt = -rate;

当扩散速度跟不上反应消耗时，模型会像卡住的齿轮突然停止计算。这时候得去求解器设置里把残差容差从1e-4调到1e-3，相当于给数值计算松了松裤腰带。

边界条件的设定藏着魔鬼细节。暴露面的对流系数如果没考虑风速变化，模拟结果会比实际碳化层厚得像千层饼。有次把对流系数设成随时间正弦波动，意外发现间歇性通风反而加剧碳化，就像用吹风机烘干墙壁反而让水汽渗透更深。

后处理时用切割线数据画出碳化前沿，看着浓度梯度曲线从陡峭变平缓，就像目睹一场持续数十年的风化快放。导出成动画后发现，氢氧化钙的消耗前沿总比CO2渗透层滞后2-3mm，这微妙的相位差解释了为什么有些混凝土表面完好内部却已碳化。

调试模型最崩溃的不是发散，而是那种似是而非的收敛。有次迭代200步后残差曲线美如预期，导出数据却发现氢氧化钙浓度出现负值——原来忘记在变量设置里勾选"物理量非负"选项，导致模型自顾自地造出了反物质混凝土。