![地下水流数值模拟基础](https://wfqqreader-1252317822.image.myqcloud.com/cover/960/40936960/b_40936960.jpg)
上QQ阅读APP看本书,新人免费读10天
设备和账号都新为新人
1.1 达西定律
![](https://epubservercos.yuewen.com/83C7AB/21277072501878406/epubprivate/OEBPS/Images/txt002_1.jpg?sign=1739095064-ff1r8g7zSOM0rHWrOnv7irf6QaiREuac-0-9d0d654ad92cfb2f46ff1fca39254917)
用矢量表示渗流速度,形式如下:
![](https://epubservercos.yuewen.com/83C7AB/21277072501878406/epubprivate/OEBPS/Images/txt002_2.jpg?sign=1739095064-Q0gz3vEoDyOKCA6x6VGHQuA8Ydqw1vEs-0-72d894b645dfa2034584115292d9bcb8)
取负号的意思:随着渗流途径的增加,水头逐渐降低。即水头沿着d s的方向d H增量永远是负的,而水力坡度是正的,所以要加一个负号。
Kx、Ky、Kz表示渗透系数在x,y,z 3个方向上的分量。对于给定的多孔介质,Kx、Ky、Kz不等时说明介质是各向异性的;而Kx、Ky、Kz相等时则说明介质是各向同性的。
在三维坐标体系下,渗透系数K是一个具有9个分量的张量:
![](https://epubservercos.yuewen.com/83C7AB/21277072501878406/epubprivate/OEBPS/Images/txt002_3.jpg?sign=1739095064-skJlGxd9rAZIGEiiyCL3wf0iIZ0KB1S9-0-4dc910e4b1676a6853c092658ab30944)
当x,y,z3个坐标轴方向与渗透系数张量主轴方向一致时,K则可简化为主对角线张量:
![](https://epubservercos.yuewen.com/83C7AB/21277072501878406/epubprivate/OEBPS/Images/txt002_4.jpg?sign=1739095064-2Lv8p7xw897zKu7RUO4gLpiuZln67iim-0-ae3c289035dabc0f0049a6200326f735)
适用条件:过去认为达西定律适用于所有做层流运动的地下水,但19世纪40年代以来的多次实验表明,只有雷诺系数小于1~10之间某一数值的层流运动才服从达西定律(图1.1.1)。
![](https://epubservercos.yuewen.com/83C7AB/21277072501878406/epubprivate/OEBPS/Images/txt002_5.jpg?sign=1739095064-5U6lXgCqsyfZ0r6OIimnZYnIIYtUtxRo-0-6edbeb016afa4a9e3dbbe3fe349c4822)
图1.1.1 雷诺系数与达西定律的适用范围(雅各布·贝尔,1982)
在这里,提出一个问题请大家思考:在水力学上,层流与紊流的临界雷诺数约为2100,而在多孔介质中最大也只有100(图1.1.1),为什么两者相差如此之大?
对于这个问题的分析,可以从以下两个要点着手。
(1)在水力学中,雷诺系数2100是通过水流在玻璃圆管中运动的实验得到的,而在多孔介质中,水流的运动空间是曲折和粗糙的,层流与紊流的临界雷诺数自然会比2100小得多。
(2)多孔介质中水流的雷诺数,是参照水力学的定义进行类比得出的:
![](https://epubservercos.yuewen.com/83C7AB/21277072501878406/epubprivate/OEBPS/Images/txt002_6.jpg?sign=1739095064-t2jcQoV3qiu5xfPauyOEnVojeVuHitbr-0-39530bc5886df9d4e6ceff58de9a8bd2)
式中d是多孔介质的某种长度尺寸,在管流试验中,d是圆管直径,在多孔介质中则应是水流通道的直径,但由于无法方便、准确地测定这个参数,在实际计算过程中则是某种代表性尺寸,例如平均粒径、d10等。根据雅可比在《多孔介质流体动力学》中的表述:就一切实际情况而论,只要根据平均粒径计算的雷诺系数不超过1~10之间的某个值,达西定律就是适用的。这也是一般教科书所采用的说法来源。