做网站一般用什么几号字广西建设网注册中心

做网站一般用什么几号字,广西建设网注册中心,南昌网站设计公司,网站技能培训机构二维钻孔封孔效果模拟案例钻孔封孔效果模拟这事挺有意思的。咱今天拿个简化版的二维模型练手#xff0c;用Python搞个渗流场可视化。先说场景#xff1a;地下50米有个直径0.2米的钻孔#xff0c;现在要往里面注水泥浆#xff0c;得看看封堵效果咋样。先整网格。用numpy搞个…二维钻孔封孔效果模拟案例钻孔封孔效果模拟这事挺有意思的。咱今天拿个简化版的二维模型练手用Python搞个渗流场可视化。先说场景地下50米有个直径0.2米的钻孔现在要往里面注水泥浆得看看封堵效果咋样。先整网格。用numpy搞个非均匀网格靠近钻孔的位置网格加密import numpy as np radius 0.2 # 钻孔半径 domain_size 5.0 # 模拟区域边长 r np.concatenate([np.linspace(radius, 0.5, 30), np.linspace(0.5, domain_size, 70)[1:]]) theta np.linspace(0, 2*np.pi, 120) R, Theta np.meshgrid(r, theta)这里用极坐标系更贴合钻孔形状。靠近钻孔的0.2米到0.5米区间安排了30层网格外围则用较稀疏的70层总网格数控制在合理范围。注意[1:]的用法是为了避免端点重复。接着上物理模型考虑浆液扩散的非牛顿流体特性def calc_viscosity(shear_rate): # 赫歇尔-巴尔克利模型 tau_0 120 # 屈服应力 k 0.5 # 稠度系数 n 0.3 # 流动指数 return tau_0/(shear_rate 1e-6) k*(shear_rate)**(n-1)这个粘度模型是关键避免了牛顿流体假设带来的误差。1e-6的小量防止除零实际工程中常见处理方式。参数取值根据实际水泥浆配比调整这里用典型值示意。边界条件设置要讲究注浆压力随时间变化injection_pressure np.interp(time, [0, 10, 30], [5e6, 8e6, 6e6]) # 分段线性加压模拟真实注浆过程——初始阶段逐步升压中期维持高压后期缓慢降压。用np.interp做线性插值比写循环更高效适合向量化计算。求解器部分用显式差分注意稳定性条件dt 0.1 * (dx**2) * rho_max / (viscosity_max) # CFL条件 for _ in range(steps): new_p p.copy() # 极坐标下的扩散方程离散 new_p[1:-1,1:-1] p[1:-1,1:-1] dt*( ... ) p new_p这里时间步长根据最大粘度和网格尺寸动态调整避免计算发散。扩散项的具体离散需要考虑极坐标的1/r项代码中简写为(...)实际需要展开拉普拉斯算子。最后来个动态可视化plt.figure(figsize(10,6)) contour plt.contourf(X, Y, pressure_field, levels20, cmapjet) plt.plot(wellbore_x, wellbore_y, w-, linewidth2) # 绘制钻孔 plt.colorbar(labelPressure (Pa)) plt.title(fGrout propagation t{time:.1f}s)用极坐标转笛卡尔坐标后的X,Y做显示更直观。颜色映射选jet突出压力梯度白线标出钻孔位置。动态更新时只需要更新contour的数据比反复创建新图快得多。跑完模拟发现个有趣现象在注浆中期浆液在钻孔北侧形成指进现象。检查代码发现是网格各向异性导致的后来改用自适应网格后改善。做工程仿真时这种数值假象要和物理真实现象区分开得反复验证。