COMSOL粗糙交叉裂隙注浆模拟 [1]模型简介：本案例使用数值模拟软件COMSOL

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回调函数发送邮件提醒。这个设计后来真帮用户避免了两次模拟中的虚拟井喷事故。