COMSOL粗糙交叉裂隙注浆模拟 [1]模型简介：本案例使用数值模拟软件COMSOL，对水平—竖直交叉裂隙采用粗糙裂隙方法进行建模，然后采用层流水平集方法追踪了浆液和水运动的两相流过程 结果表明:以粗糙裂隙方法构建的水平及竖直裂隙能有效捕捉裂隙表面的不规则和复杂性特征 此外粗糙表面的不规则性显著影响了流体的运动模式，导致流体在交叉点附近出现了“偏流”效应 [2]案例内容：COMSOL数值模型一份 [3]模型特色：（1）掌握在COMSOL中粗糙裂隙的建模方法；（2）掌握水平集追踪两相流过程 COMSOL6.2版本

先搞点粗糙的

在岩土工程现场，裂隙从来都不是规整的几何线条。最近在COMSOL6.2上捣鼓了个粗糙裂隙注浆模型，发现用参数化曲面搞粗糙面特别有意思。比如在交叉裂隙的建模中，我往壁面坐标里塞了个高斯随机场：

surface_z = 0.1*randn(x,y) + 0.05*sin(3*x)*cos(2*y)

这种操作让原本光滑的裂隙表面瞬间长满"青春痘"。运行后发现，当浆液流经这些凸起时，速度场会出现明显的蛇形走位——就像早高峰地铁换乘通道里的人流，遇到障碍物自动分岔。

水平集的黑魔法

两相流界面追踪用层流+水平集组合拳效果拔群。设置相场参数时有个坑：界面厚度系数得和网格尺寸匹配。有次手滑设成0.5mm，结果浆液前锋直接表演量子隧穿，把隔壁网格都染红了。后来按这个公式调整才正常：

epsilon = 2*max(mesh_size)

物理场设置里最骚的操作是在交叉点区域加了个瞬态脉冲扰动。代码里写了个时间触发器，当浆液到达交叉点0.5mm范围内时，自动激活局部湍流模型。这招让模拟结果出现了文献里说的"浆液抢道"现象。

COMSOL粗糙交叉裂隙注浆模拟 [1]模型简介：本案例使用数值模拟软件COMSOL，对水平—竖直交叉裂隙采用粗糙裂隙方法进行建模，然后采用层流水平集方法追踪了浆液和水运动的两相流过程 结果表明:以粗糙裂隙方法构建的水平及竖直裂隙能有效捕捉裂隙表面的不规则和复杂性特征 此外粗糙表面的不规则性显著影响了流体的运动模式，导致流体在交叉点附近出现了“偏流”效应 [2]案例内容：COMSOL数值模型一份 [3]模型特色：（1）掌握在COMSOL中粗糙裂隙的建模方法；（2）掌握水平集追踪两相流过程 COMSOL6.2版本

求解器也得跟着糙

粗糙表面的存在直接干碎了常规求解策略。原本用稳态求解器跑10分钟就能出结果，现在必须切到瞬态模式，时间步长还得动态调整。这里有个小技巧：在交叉点区域设置局部时间步长缩放因子：

dt_local = dt_global * (1 + 0.3*rand())

随机扰动项模仿了真实工况的不确定性。后处理时发现，当浆液黏度超过50Pa·s时，粗糙面会形成微型"水坝"，导致注浆压力出现毛刺状波动——这现象在光滑模型里根本看不到。

偏流效应真人秀

最带劲的结果出现在流速矢量图上。原本预期浆液会在交叉点均匀分流，实际模拟显示约63%的浆液选择水平裂隙优先通过。用截面探针抓取数据后发现，竖直裂隙入口处出现了负压漩涡，这可能是导致偏流的元凶。有个反直觉的现象：当表面粗糙度超过临界值0.8mm后，偏流效应反而减弱——可能因为流体开始摆烂，直接走之字形路线了。

模型文件里藏着个彩蛋：在裂隙末端设置了自动报警变量，当注浆压力超过岩石抗拉强度时，会触发Python回调函数发送邮件提醒。这个设计后来真帮用户避免了两次模拟中的"虚拟井喷"事故。