基于四参数随机生长法重构土体的格子玻尔兹曼细观渗流研究

多孔介质模型的重构问题是土体细观渗流机理研究的基础和关键。由四参数随机生长法(QSGS)构建土体模型,采用格子玻尔兹曼方法(LBM),通过MATLAB自编程序研究重构土在不同条件下的细观渗流机理。结果表明:随模型尺寸增大,孔隙连通程度显著提高,300×300格点大小的模型连通孔隙率增长幅度(34.38%)最大,继续扩大模型尺寸发现增加不明显;流体粒子在孔隙连通性好、孔径大的区域,会形成主渗流通道,且存在指进效应,孔道中间流速最大,可达0.0324,越靠近孔壁流速越小;大孔隙率土的流速比小孔隙率土大,而低孔隙率土中的流速相比大孔隙土更稳定;LBM模拟渗透率与经典K-C公式计算结果对比发现,孔隙率越高计算渗透率越准确(n=0.78,误差为10.22%);土颗粒越小,渗流孔道越细窄、分布越密集,对应的速度场分布更为均匀,同时流速也更小。该研究成果能较好地揭示重构土的细观渗流机理,也可为现有细观土体孔隙流研究提供一定借鉴。     

非饱和花岗岩残积土水-气两相驱替过程数值模拟

水流在非饱和土体中的入渗过程实质上是水在下渗的过程中驱替空气的两相流问题。为揭示非饱和花岗岩残积土水-气两相驱替动态渗流机理,选取福州某地原状花岗岩残积土作为研究对象,基于工业CT扫描图像与Level Set方法,研究了原状土样两相驱替的动态特征。结果表明:对于细观尺度水-气两相驱替模拟,Level Set法能很好地捕捉两种不混溶流体间的界面位置;水-气两相驱替过程存在大孔隙优先流特征,且“绕流”现象一般易于出现在孔隙成圆度较高处;两相渗流速度主要受孔道迂回度控制,笔直、较宽孔道,渗流速度相对较高,同时存在明显的“优势通道”,且随渗流时间增大以先急后缓的特征呈正相关变化,最大增速率为 10.77%,最小仅 1.90%;孔道横截面速度大小分布与孔隙结构有关,“回流”和“绕流”现象会使驱替速度骤降,降低幅度可达21.62%;驱替阻力最大出现在孔壁处,孔道越窄,阻力越大;驱替效率与驱替压差成正比关系,且初期加压增速效果显著,可达25.49%,后期仅为1.47%。该研究成果可丰富降雨型滑坡理论基础并预防灾害产生,具有重要的理论价值及工程意义。