COMSOL瓦斯抽采二氧化碳封存躯替多物理场耦合方面 论文复现下午三点实验室的咖啡机又罢工了我对着屏幕上扭曲的流速云图发呆。做煤层气抽采模拟最头疼的就是这个——达西流和固体变形较着劲渗透率参数像过山车似的忽上忽下。突然想起上个月在COMSOL年会上偷师的那个骚操作抄起建模区开始搞事情。在煤层模块里先摆个达西流物理场参数栏直接敲k0 2e-13; //初始渗透率[m²] sigma_eff solid.sxx solid.syy solid.szz - 3*p; //Terzaghi有效应力 k k0 * exp(-0.05*sigma_eff); //指数型渗透率衰减这段代码把固体力学模块的应力场和流体压力场勾搭上了。跑个瞬态分析看着渗透率云图像被揉皱的锡纸孔隙压力下降时煤岩骨架果然开始喘气——应力重分布导致的新渗流通道逐渐显现。COMSOL瓦斯抽采二氧化碳封存躯替多物理场耦合方面 论文复现做二氧化碳封存更刺激。在地下水库模块里把传统的达西定律改成mu_CO2 (1.0 0.035*(T-313)) * 1e-3; //温度依赖粘度[Pa·s] rho_CO2 1.98 * p / (1e5); //可压缩流体密度这还没完接着在化学物质传输接口插入反应速率 1e-6 * c_CO2^2 * (1 - exp(-0.1*t)); //碳酸盐岩溶解动力学结果在2号注入井周围矿物饱和度像病毒扩散似的蔓延开。特别要注意孔隙率变化引发的反馈phi phi0 0.02*(1 - exp(-5*反应速率)); //孔隙率动态更新这时候多物理场耦合的蝴蝶效应就来了——孔隙结构改变引发流场畸变流场变化又影响反应物传输整个模型像被丢进滚筒洗衣机的乐高积木。调试这种模型要善用探针功能。某次在井筒监测点抓取到诡异的压力震荡最后发现是固体力学求解器的阻尼系数设太大导致应力波传播失真。改完参数后顺手加个with(模型.结果().导出().动画()), framerate(24), quality(高清), filename(打工人の救赎.mp4) );毕竟老板要看的是会动的云图不是静态的论文配图。这种数值模拟就像在虚拟地层里玩华容道既要保证物理规律严丝合缝又得让计算结果讲人话。有时半夜盯着收敛曲线感觉自己像在跟偏微分方程谈恋爱——它虐我千百遍我待它如初见。