地下水运动方程理论分析及其在cfd中的应用 theoretical analysis of groundwater movement equation and its application in cfd.pdfVIP

第28卷第1期 北方工业大学学报 V01．28NO．1 2016年3月 J．NORTHCHINAUNIV．OFTECH． Mar．2016 地下水运动方程理论分析及其在CFD中的应用* 宋 伟 代 云 郁文红 (北方工业大学土木工程学院，100144，北京) 摘要在单井循环地下换热器仿真研究中，对于多孔介质模型参数设置常根据经验确 定，相关理论研究较为缺乏，导致仿真研究结果与实验结果偏差较大．因此，本文对地下水运动 方程进行理论分析，利用纳维一斯托克斯方程与米斯马赫方程推导得出不可压缩流体在不可变 形多孔介质中的运动方程；利用伯努利方程和米斯马赫方程，将能量方程方程中的水头损失与 流体在多孔介质中流动的阻力损失联系起来，使伯努利方程在多孑L介质中的得到应用；根据实 验结果与仿真数据确定了Fluent中多孔介质模型的2个经验系数，C。和C。分别为2．5×106 和1．2；最后，应用该经验系数，仿真得到3种单井循环地下换热器影响下的地下水速度流线． 关键词 地下水运动方程；单井循环系统；多孔介质；米斯马赫方程；CFD 分类号 TK529 含水层作为地下水运动的载体，是一个或一 Darcy)通过研究水在直立均匀砂柱中的渗流实 组饱和的具有储存和产出可用淡水水量的地质 验]，得到了著名的达西定律，建立了流量(或流 层n]，导水、储水和给水度是含水层的主要性 速)与水头之间的关系，反映了多孔介质中水流 质口]．作为多孔介质的含水层[3]，流体的运动发 的宏观运动规律．由于渗流速度与水力坡度成正 生在骨架的孔隙和缝隙中，流体在以孔隙或缝隙 比，也称它为线性渗透定律，如式(1)和(2)所示． 壁面为边界的小通道中运动¨]．按照其是否含有 多年来，研究者们做了大量的实验，证实了线性 潜水面而分为潜水含水层和承压含水层[5]．单井 渗透定律具有一定的适用范围[8’9]．雅克比·贝 循环地下换热器多应用于承压含水层中，本文以 尔把地下水渗流按照渗流速度从低到高分为3 假想的连续系统来代替真实的多孔介质系统[6]， 种情、况[10]，如表1所示．实验证明，只有当雷诺 对涉及该系统的流动特性的相关数学模型展开 数R。10的条件下，经过多孔介质的流体作层 研究．分别从动量方程(纳维一斯托克斯方程)和 流运动，渗流才满足达西定律． 能量方程(伯努利方程)出发，对含水层中地下水 (1) Q7一KA旦尝：KAJ ．) 井流运动方程进行理论推导，深人剖析其与达西 n 定律和米斯马赫方程之间内在联系，从而指导 U一旱一KJ (2) ^ CFD仿真实践． 式中：Q’～渗流量(m3／s)；K一渗透系数 (m／s)；A一砂柱的横截面面积(m2)；H。一H2一 地下水运动方程理论分析 砂柱顶底断面上的测压水头差(mHzO)；S一渗 1．1 达西定律