渗流PFC_CFD耦合细观数值模拟.pdf

　　　　　　 水文地质工程地质 2008 年第 2 期 　 ·66 · 渗流的 PFCCFD 耦合细观数值模拟 刘 　洋 ,李 　飞 ,柴小庆 , 杜晓伟 , 张 　华 (北京科技大学土木与环境工程学院土木系 ,北京 　100083) 摘要 : 本文针对饱和多孔介质中渗流特点 ,在分析渗流连续介质力学模型的基础上 ,采用描述多孔介质中固相颗粒和液 ( ) 相流体耦合的细观力学模型 ,其中固相颗粒采用离散元的颗粒流理论 PFC 模拟 ,液相流体通过求解平均 NavierStokes 方 ( ) 程的计算流体动力学 CFD 技术计算 ,并利用建议的模型对二维渗流问题进行了模拟验证。数值模拟结果表明 ,采用 PFCCFD 耦合细观力学模型能够描述从低雷诺流到高雷诺流范围很大的流体运动 ,可以用于多孔介质中的渗流问题分 析。 关键词 : 渗流 ; 颗粒流 ; 计算流体动力学 ; 流固耦合模型 + ( ) 中图分类号: p64211 3 ;p6412 　　　　文献标识码 : A 　　　　文章编号 : 2008 　　饱和多孔介质是由固相颗粒集合体以及充填其中 流体的流动以及它们之间的相互作用。 的流体组成的混合物 ,基于颗粒之间以及颗粒与流体 1 　多 - 孔介质的微观模拟 之间的相互作用 ,饱和多孔介质呈现出高度复杂的非 线性力学特性。传统多孔介质的连续耦合模型通常是 11 　固相颗粒的微观模拟 ( ) [1 ] 土体( ) 基于唯象的描述 如混合物理论或Biot 理论 或微观 尤其是砂土 是由不连续的颗粒组成的集合 力学运动方程的均一化来建立的。 体 ,可以采用离散颗粒的集合来有效地模拟。在离散 用于模拟多孔介质的连续模型都是基于固液两相 颗粒模拟中 ,单个颗粒的运动直接采用牛顿第二运动 的唯象描述 ,在场方程和本构关系中有很多假设 ,如初 定律计算 , 粒间接触力通过颗粒之间的接触关系得 始局部孔隙率通常认为与总体平均孔隙率相等 ,并假 到[4 ] 。在岩土力学中最常用的方法就是离散单元法 定其不随时间和空间变化。此外 ,各相之间的动量交 (DEM) 。 换遵循 Darcy 定律 ,而渗透系数却是独立应变的。实 12 　液相流体的微观模拟 际上 ,Darcy 定律是低阶的NavierStokes 方程 ,仅当孔隙 这方面早期的研究主要是解慢速流绕球或圆柱体 流是薄层流、惯性力可以忽略的情况下才有效[2 ] 。在 的 Stokes 方程 ,得到流体作用于单个球体或圆柱体的 较大的水力梯度下产生非薄层流 ,导致大的孔隙流体 拽曳力。在此基础上 ,采用不同的数值方法研究孔隙 速率。随着砂土骨架大的体积变形产生的孔隙率变化 尺度上的流体流动。 [5 ] 也可能导致偏离 Darcy 定律。为了得到对饱和砂土液