COMSOL中基于浆液黏度时空变化的水平裂隙岩体注浆扩散数值模拟

COMSOL基于浆液黏度时空变化的水平裂隙岩体注浆扩散数值模拟

岩层注浆这事儿挺有意思的，就像给地下裂缝打葡萄糖浆似的。最近用COMSOL折腾了个动态黏度的注浆模型，发现浆液这玩意儿根本不是老实人——它的黏度会在流动过程中随时空变化，把牛顿流体模型按在地上摩擦。

先看裂隙几何设定，在COMSOL里用窄矩形域模拟水平裂隙，宽度控制在0.5mm到3mm之间。这里有个坑：网格得沿着流动方向拉伸，否则计算量直接爆炸。我习惯用自由四面体网格配合边界层，这样既能捕捉近壁面效应，又能保证计算速度：

mesh_params = {
    'element_size': 0.2,  # 基础网格尺寸(mm)
    'boundary_layers': 3,  # 边界层数
    'stretching_ratio': 1.5  # 流动方向拉伸系数
}

黏度动态变化才是重头戏。实测数据表明浆液黏度随时间呈指数衰减，同时受流速影响。在材料属性里搞了个自定义函数：

% COMSOL内置的黏度模型扩展
mu = mu0 * exp(-alpha*t) * (1 + beta*ux^2); 
% mu0初始黏度，alpha时衰系数，beta流速影响因子

这个非线性项让方程组变得调皮捣蛋。求解器设置得用全耦合+自动牛顿法，时间步长建议用自适应，初始步长别超过0.1秒。有次手贱设了固定步长1秒，结果残差曲线比过山车还刺激。

边界条件方面，注浆口用流量入口，出口设零压力。特别注意接触角效应，用壁面滑移条件处理浆液-岩壁相互作用：

// 滑移速度计算
u_slip = k_slip * (tau_wall - tau_yield);
// k_slip滑移系数，tau_yield浆液屈服应力

算完发现个反直觉的现象：低黏度区域反而出现流速减缓，仔细检查才发现是剪切稀化效应和压力梯度耦合搞的鬼。这时候后处理就得祭出粒子追踪模块，看单个浆液微团的运动轨迹，比看速度云图直观多了。

最后说个实用技巧：在结果分析里加个黏度场动画输出，能清楚看到浆液前锋的黏度突降过程。这招在项目汇报时特别好使，甲方爸爸看到动态变化的彩色云图，比看十页公式都开心。