考虑空气阻力影响的降雨入渗过程分析研究.pdfVIP

考虑空气阻力影响的降雨入渗过程分析研究 朱岳明孙冬梅 ， 。 (河海大学水利水电工程学院，江苏南京210098) 摘要：本文根据多相流理论，结合水、气的质量守恒定律和达西定律，建立水一气二相流模型，模型的求解采 用积分有限差分法和Newton-Raphson迭代方法，通过变换主要变量来表达相态的变化，实现了水相、气相边 界条件及降雨入渗边界的精确模拟。利用上述模型对一土柱试验进行模拟，从而验证了模型的正确性，然后研 究一均质土层的降雨入渗过程，得到孔隙水压力、孔隙气压力和毛细压力的变化过程，得到入渗率与地表孔隙 气压力的变化关系，表明孑L隙气压力的增大对人渗水流产生阻滞作用，在求解土壤水流人渗的非稳定渗流问题 时，需考虑空气阻力的影响。 关键同：水一气二相流饱和一非饱和渗流积分有限差分法降雨人渗 一般认为非饱和土是南同体颗粒、孔隙气和孔隙水所组成的三相体【l】，水流在非饱和带中的入渗过程实 质上是水在下渗的过程中驱替空气的过程。在目前的研究中，一般仅限于水相，气相的流动通常被忽略，认 为气相压力与大气压力平衡[21。而在实际的强降雨人渗过程中，孑L隙气的存在对水相的流动将产生不可忽略 的影响，大量试验都证明了这一点【34】。 ． 本文以土壤的孔隙水和孔隙气作为研究对象，建立水一气二相流数学模型。水一气二相流模型的求解采 用积分形式的有限差分法，然后利用该模型分析研究均质土壤在降雨入渗下非稳定渗流问题，分析气相的存 在对水流下渗的作用。’ 1 水一气二相流数学模型 对于常见的岩土工程，可认为渗流是发生在绝热的恒温系统中，水相和气相分别满足各自的质量连续性 方程，假设固体骨架不变形，处于平衡状态的气相(g)和液相(I)在多孔介质孔隙中流动，满足达西定律，并且 考虑毛细作用力，气相可看做是理想气体，满足理想气体定律。在模型中考虑空气与水的相互溶混，液相包 含液态水和溶解其中的空气，通常空气的溶解性很小，符合Henry定律；气相包含干空气和水蒸气，气相压 “ ‘ 力满足Dairon气体分压力定律，并假设水蒸气的分压力等于饱和水蒸气压力【们。 1．1基本控制方程 将任意体积内空气和水所满足的质量守恒方程作为水一气二相流模型的基本控制方程‘5】： ． ．1 ，I ． LF。×疗dF一+￡g。dr．2盖LM”dK (1) 式中：圪为单元体；F。为圪的封闭边界；上标K表示组分(水W或空气口)；为法向量余弦，指向单元体内 相g，≯为孔隙率，-％为∥相的饱和度，乃为∥相的密度，k∥m3；巧为r在夕相中所占的质量分数；F。为 的单位流量，kg／s，流人为正，不考虑空气在水中的扩散，有F。=∑占巧易，’易=野坳为∥相的单位流量， kg／s，“8为∥相的流速，m／s，遵循达西定律：， 厶 ％=一七孚(即声一岛g) (2) p8 式中：后为固有渗透系数张量，m2；t。为液相或气相的相对渗透系数；∥。为液相或气相的动力黏滞系数， 合作项目资助(BZ2003046) 第一作者简介：朱岳明，1958年生，男，教授，博士，主要从事混凝土渗控、温控和稳控和反问题研究。E-mail： hhuzym@126．com