岩体单裂隙高速非达西流理论试验研究及数值模拟.docx
- 文档编号:26509943
- 上传时间:2023-06-20
- 格式:DOCX
- 页数:17
- 大小:27.19KB
岩体单裂隙高速非达西流理论试验研究及数值模拟.docx
《岩体单裂隙高速非达西流理论试验研究及数值模拟.docx》由会员分享,可在线阅读,更多相关《岩体单裂隙高速非达西流理论试验研究及数值模拟.docx(17页珍藏版)》请在冰豆网上搜索。
岩体单裂隙高速非达西流理论试验研究及数值模拟
岩体单裂隙高速非达西流理论、试验研究及数值模拟
第1章绪论2
1.1课题研究目的和意义2
1.2研究历史及现状2
1.2.1实验研究3
1.2.2计算模型6
1.2.3数值方法8
1.3本文研究内容9第2章岩体单裂隙渗流理论10
2.1控制方程10
2.1.1连续性方程11
2.1.2动量方程12
2.1.3本构方程13
2.2定解条件14
2.2.1边界条件14
2.2.2初始条件15
2.3裂隙渗流层流规律15
2.3.1立方定律15
2.3.2线性变隙宽渗流规律17
2.4裂隙渗流紊流规律17
2.4.1层流与紊流的判别雷诺数18
2.4.2紊流时均连续性方程和雷诺方程18
2.4.3紊流的半经验理论20
2.4.4裂隙紊流流速分布21
2.5本章小结24第3章高速非达西渗流规律的实验探讨25
3.1实验目的25
3.2实验介绍25
3.2.1实验装置25
3.2.2实验数据测量方法29
3.2.3实验步骤31
3.3实验结果及分析32
3.3.1层流实验结果32
3.3.2高速非达西实验结果33
3.4本章小结40第4章Navier-Stokes方程数值解41
4.1CFD简介41
4.2Fluent简介42
4.2.1Fluent建模求解步骤42
4.2.2Fluent求解方法43
4.3物理模型45
4.3.1网格模型45
4.3.2控制方程46
4.3.3紊流模型46
4.3.4壁面区处理47
4.3.5定义材料48
4.3.6边界条件48
4.4Fluent求解48
4.4.1控制方程的离散48
4.4.2控制方程的求解49
4.4.3松弛因子50
4.4.4收敛标准51
4.5数值计算结果51
4.6本章小结55第5章总结与展望56
5.1全文小结56
5.2存在的问题及研究展望56第1章绪论
1.1课题研究目的和意义岩体中的渗流问题作为岩土工程问题,由于其制约着基础岩土体的承载能力.影响着岩质边坡的稳定性.决定了地下洞室的渗涌水量等,在水利水电.交通.工民建以及国防等工程中越来越受到重视。
锦屏二级水电站引水隧洞长探洞揭示的地质情况和已基本查明的工程区岩溶水文地质条件表明,该工程区涌水具有高水头.大流量的特点。
而未来在地下厂房开挖过程中的稳定涌水量和单点瞬时涌水量均可能比探洞中要大,必然影响工程安全及稳定。
对于类似于锦屏二级水电站区域岩溶裂隙围岩介质,按水动力特征可将其含水空间分为三类:
1.主要起储水.释水作用的遵循达西定律的基质岩块孔隙和微裂隙;
2.主要起导水作用.其水流仍服从达西定律的中宽裂隙空隙介质;
3.起导水作用但其水流的流态已属于紊流的岩溶管道和宽裂隙空隙介质。
目前国内外对于裂隙岩体渗流分析所依据的基本定律还是达西定律,而对于不再符合达西定律的深埋隧洞大流量涌水问题,如果仍采用和依据达西定律来分析,就会导致工程不良的后果,高速非达西渗流问题的研究已显得非常迫切和重要。
单条裂隙作为岩体裂隙网络的基本单元,决定了地下水在岩层中的基本渗透特征,是各种渗透理论模型的基础,因此单裂隙渗流规律的研究是最基础的研究。
1.2研究历史及现状渗流的研究开始于土体。
1856年,法国工程师Darcy根据直立均质沙柱模型渗流实验,总结出线性渗流定律流速与水力梯度呈正比,即著名的达西定律,这标志着经典渗流理论的诞生[1]。
早期,俄国对渗流的研究处于领先地位。
1889年,N.E.儒可夫斯基导出了渗流微分方程,指出在数学形式上渗流微分方程相似于热传导方程。
100多年来,工程实践中的渗流分析都是采用经典渗流理论,渗流力学已经成为流体力学的一个重要分支。
对渗流问题采用经典渗流理论进行分析是有一定条件的,那是因为经典渗流理论以连续介质假定为基础,只有表征体积单元与研究域的尺寸相比足够小,才能使用连续介质渗流理论。
对于土体渗流,表征体积单元相比于研究区域足很小,连续介质渗流理论可以直接应用。
岩体渗流不同于土体渗流。
对同一种岩石,致密的岩块渗透系数一般小于10-7~10-8cm/s,岩体的渗透系数要比岩块渗透系数大104~107倍[2]。
产生这种差异的原因是岩体中一般分布有产状不同.大小不等的裂隙。
由此可见,对岩体裂隙导水性的研究更为重要。
岩体裂隙从其形式看在二维方向上延伸,其两个方向上的尺寸比第三个方向上的尺寸大得多,另外裂隙尺寸范围可以从几厘米到几千米,甚至延伸到整个研究域,这样就不一定存在与研究域的尺寸相比很小的表征体积单元,所以不是所有岩体渗流问题都可以使用连续介质渗流理论进行分析。
但历史上相当长一段时间工程设计中仍把岩体视为与土体类似的多孔介质,用连续介质渗流理论进行渗流分析。
直到法国Malpasset拱坝溃决和意大利Vajiont拱坝失事,人们才意识到岩体渗流不能简单类比于土体渗流,必须对岩体渗流进行研究。
1.2.1实验研究理论分析.数值计算与实验研究(包括现场观测和物理模型试验)是目前研究岩体渗流的主要手段,其中物模实验法是研究裂隙岩体渗流的一个最重要且最直接的途径[3]。
前人对岩体单裂隙渗流规律进行了大量的实验研究,主要成果如下:
在假定裂隙光滑.平直.无限长,裂隙间流体为定常层流及裂隙流体为不可压缩.粘性的条件下,Snow[4]通过Navier-Stocks公式导出立方定律单裂隙渗流规律,其表达式为:
(1-1)式(1-1)中:
为单宽流量;为裂隙开度;为水流运动粘滞系数;为重力加速度;为水力梯度。
之后,许多学者进行了单裂隙渗流实验研究,并证明了立方定理的有效性,如苏联学者Lomize.法国学者Louis等。
Romm[5]通过对微裂隙和极微裂隙的研究,提出只要隙宽大于0.2微米,立方定律总是成立的。
速宝玉[6]亦进行了微裂隙渗流实验,指出此时水流已呈现出类似非牛顿流体特性。
然而天然岩体裂隙均为粗糙裂隙,并不满足平行板裂隙的假定。
于是,许多学者进行了人工粗糙裂隙的渗流实验研究:
Lomize[7]采用两块薄而窄的玻璃板形成裂隙,裂隙长度20厘米,裂隙张开度变化范围0.05~
1.0厘米,实验中裂隙面形状为波浪行.楔形等;Louis[8]用两块长200~280厘米,宽80厘米的混凝土平板形成粗糙裂隙;速宝玉[9]在玻璃板壁面上粘贴粒径0.5~0.63毫米沙子来形成粗糙裂隙。
并且各自得到了立方定律的修正公式:
(1-2)式(1-2)中:
为平均裂隙开度;为修正系数,与裂隙面粗糙度及裂隙开度有关。
Lomize[7]提出,Louis[8]提出,速宝玉[9]则认为。
为裂隙粗糙度;为裂隙开度均方差。
Amadei[10]亦得到式(1-2)相似的结论,其中。
公式(1-2)是基于立方定律的假设对立方定律进行修正而得,许多学者则考虑别的因素对粗糙裂隙流规律进行总结。
Barton[11]通过大量实验,提出(节理粗糙度系数)修正法:
(1-3)式(1-3)中:
为节理粗糙度系数,为等效水力裂隙开度。
耿克勤[12]选用天然裂隙,裂隙起伏差
2.3~
2.8毫米,机械裂隙开度最大3毫米左右,进行渗流实验研究,得出经验公式:
(1-4)式(1-4)中:
,为公式拟合系数,为机械裂隙开度。
对于小开度裂隙层流,;对于中开度,;对于大开度,。
Iwai[13]发现裂隙面粗糙性对裂隙水流规律的影响主要与隙面面积接触率(隙面凸起接触面积与总面积之比)有关:
(1-5)式(1-5)中:
为相应于接触面积率的流量;为的流量;为经验系数。
以上结论公式都是基于裂隙流为层流的实验基础上建立起来的,对于裂隙高速非达西流规律还未充分认识,有待进一步研究。
Sharp通过对花岗岩中的天然裂隙进行渗流实验,裂隙开度范围0.27~
1.54毫米,提出了非达西流的存在。
Lomize[7]和Louis[8]各自提出了紊流范围的计算公式。
Lomize[7]紊流公式:
(1-6)适用上限:
(1-7)Louis[8]紊流公式:
(1-8)适用上限:
(1-9)式(1-6)和式(1-8)中:
为单宽流量;为裂隙开度;为水流运动粘滞系数;为重力加速度;为水力梯度。
式(1-7)和式(1-9)中的均为:
(1-10)速宝玉[6]指出Lomize[7].Louis[8]对裂隙紊流的总结,是移植了水力学管流理论。
裂隙形状不同于圆管,圆管的摩擦阻力经验系数必然不同于圆管,利用圆管摩擦阻力经验系数得到的裂隙紊流公式必然不能准确地描述裂隙水流规律。
速宝玉[6]从Prandtl紊流半经验理论出发,经数学推导,并应用实验结果确定表达式中的待定系数,得到以下半经验公式:
(1-11)式(1-11)中:
为缝隙高度(cm),为缝隙开度(cm),为重力加速度(980cm/s2),为水力梯度,为水流运动粘滞系数(cm2/s),为缝隙的总过流量(cm3/s)。
但速宝玉的实验中裂隙开度在
2.85mm以下,且达到的最高水力梯度较小(裂隙开度2mm时最高水力梯度小于0.5),因此其半经验公式也具局限性,有待进一步研究。
1.2.2计算模型裂隙成因各不相同,极为复杂,导致裂隙的产状.分布密度.张开程度.隙面粗糙度等变化极大。
而裂隙的导水性决定了岩体的导水性,因此,对渗流模型而言,能正确反映裂隙特性是极为重要的。
目前的渗流计算模型根据其是否考虑了岩体中裂隙系统和岩块孔隙系统之间的水交替过程来分,可以分为两类:
裂隙-孔隙双重介质模型和非双重介质模型[14~16]。
(1)裂隙-孔隙双重介质模型裂隙-孔隙双重介质模型认为裂隙岩体是由裂隙系统和岩块孔隙系统共同构成的统一体。
在孔隙系统和裂隙系统构成的统一体中,存在两个渗流场,并且在两个系统之间可实行水交替。
其基本思想比较符合实际。
该模型是通过以下方式来实现渗流计算的:
1.基于达西定律分别对两个系统建立控制方程,
2.利用两类系统之间的水交替方程将两个系统的水流输运方程联系起来。
然而该模型比较复杂,实现比较困难。
根据水交替方程的建立方法,又可将其分为拟稳态流模型和非稳态流模型。
①拟稳态流模型拟稳态流模型中两类系统间的水交替方程是根据裂隙系统与岩块孔隙系统的水交替量和两类系统中的水头差成正比而建立的。
由于水交替量是时间的隐式,所以称为拟稳态流模型。
该模型的主要代表人物有Barenblatt[17].Warren和Rott[18]等。
②非稳态渗流模型非稳态流模型虽然考虑了两类系统的水交替过程,但只考虑岩块孔隙水向裂隙流动而忽略裂隙水向孔隙流动。
该模型根据岩块孔隙中的水流运动规律来建立水交替方程。
由于水交替量是关于时间的显式,因此称之为非稳态流模型。
根据模型中裂隙分布假定不同,可分为平行裂隙非稳态流模型[15~16]和组裂隙非稳态流模型[17]。
(2)非双重介质模型非双重介质模型不考虑岩体中孔隙系统与裂隙系统的水交替过程,忽略孔隙的导水作用,着重研究裂隙的导水作用。
该模型进行了工程允许范围内的简化,是目前岩体渗流计算使用最多的模型。
非双重介质模型主要包括等效连续介质模型.离散裂隙网络模型和等效-离散耦合模型。
①等效连续介质模型等效连续介质模型采用平均的思想将裂隙的导水作用平均到整个岩体裂隙,将岩体模拟为具有对称渗透张量的各向异性连续体,采用经典的渗流理论进行分析。
Snow[18].Long[19].Oda[20].张有天[21]和田开铭[22]等对此模型进行了研究。
该模型的突出优点是可以使用经典渗流理论进行分析,具有雄厚的理论基础和丰富的工程经验。
另外对岩体渗流进行等效平均后,就不需知道每条裂隙的确切位置和水力特性,大大降低了工作量。
然而此模型的适用性.有效性不一定得到保证,即裂隙岩体连续性假设不一定成立。
另外,使用此模型时,在确定裂隙岩体等效渗透张量方面亦存在着难度。
②离散裂隙网络模型离散裂隙网络模型是在确定所有裂隙产状.裂隙开度后,以单裂隙渗流规律为基础,利用任意裂隙交叉点流量平衡来建立方程。
Wittke[23].Louis[24].毛昶熙[25~26].Wilson[27].Long[28].Nordqvist[29].Dershowitz[30].万力[31]等对此模型进行了研究。
离散裂隙网络模型试图模拟岩体裂隙网络中每条裂隙,模拟每条裂隙中水流输运过程,以反映裂隙岩体中的真实渗流场,具拟真性好.精度高等优点。
然而岩体中普遍发育多组裂隙,要全部模拟工作量相当大,特别是三维问题,几乎不能实现。
岩体中的裂隙分布还具有很大的不确定性,要建立真实的裂隙网络也相当困难。
由于离散裂隙网络模型的实现比较困难,还不能广泛运用到实际工程中。
③等效-离散耦合模型等效-离散耦合模型结合了等效连续介质模型和离散裂隙网络模型优点,避免了离散裂隙网络模型模拟每一条裂隙的繁琐,又能保证等效连续介质模型的有效性。
该模型将渗流区域划分为两类,即采用经典渗流理论进行分析的等效连续介质和基于单裂隙渗流规律的裂隙网络系统。
建立等效-离散耦合模型的原则是:
1.等效连续介质和离散介质接触处的水头相等,
2.流进.流出等效连续介质区域和离散介质区域的流量相等。
王媛[32]提出:
根据区域裂隙密度大小来分类,裂隙密度小的区域采用离散裂隙网络模型,裂隙密度大的区域采用等效连续介质模型;或者根据裂隙开度来分类,对数目不多的大中型裂隙采用离散裂隙网络模型,对相对于比较密集的小裂隙区域采用等效连续介质模型。
再根据两类介质接触处的水头相等以及结点流量平衡来建立方程。
1.2.3数值方法对于在求解域内建立的控制方程,理论上是有真解的。
但由于所处理的问题自身的复杂性或者复杂的边界条件,造成很难获得方程的真解,这时可用数值方法来得到近似解。
渗流数值计算方法主要有有限差分法.有限元法.有限体积法.边界元法等计算方法[33]。
(1)有限差分法(FiniteDifferentMethod,简称FDM)有限差方法是数值解法中最经典的方法。
它是将求解域划分为差分网格,用有限个网格节点代替连续的求解域,然后将偏微分方程(控制方程)的导数用差商代替,推导出含有离散点上有限个未知数的差分方程组。
求差分方程组(代数方程组)的解,就是微分方程定解问题的数值近似解,这是一种直接将微分问题变为代数问题的近似数值解法。
这种方法发展较早,比较成熟,适合于求解双曲型和抛物型问题。
用它求解边界条件复杂.尤其是椭圆型问题不如有限元法或有限体积法方便。
(2)有限元法(FiniteElementMethod,简称FEM)有限元法是将一个连续的求解域任意分成适当形状的许多微小单元,并于各小单元分片构造插值函数,然后根据极值原理(变分或加权余量法),将问题的控制方程转化为所有单元上的有限元方程,把总体的极值作为各单元极值之和,即将局部单元总体合成,形成嵌入了指定边界条件的代数方程组,求解该方程组就得到各节点上待求的函数值。
有限元法的基础是极值原理和划分插值,它吸收了有限差方法中离散处理的内核,又采用了变分计算中选择逼进函数并对区域进行积分的合理方法,是这两类方法相互结合.取长补短发展的结果。
它具有很广泛的适应性,特别适用于几何及物理条件比较复杂的问题,而且便于程序的标准化。
对椭圆型方程问题有更好的适用性。
有限元法的不足是求解速度较有限差方法和有限体积法慢。
(3)有限体积法(FiniteVolumeMethod,简称FVM)有限体积法又称为控制体积法。
其基本思路是:
将计算区域划分为网格,并使每个网格点周围有一个互不重复的控制体积。
将待解的微分方程(控制方程)对每一个控制体积积分,从而得出一组离散方程。
其中的未知数是网格点上的因变量。
为了求出控制体积的积分,必须假定因变量在网格点之间的变化规律。
从积分区域的选取方法看来,有限体积法属于加权余量法中的子域法,从未知解的近似方法看,有限体积法属于采用局部近似的离散方法。
简言之,子域法加离散,就是有限体积体积法的基本方法。
有限体积法兼有有限差分法和有限单元法的特点,能够适应结构化网格和非结构化网格。
虽然出发点不同,在结构化网格下,FVM离散结果往往与FDM相同;在非结构网格下,又有与FEM类似之处。
FVM兼收了FDM的效率和FEM的灵活。
(4)边界元法(BoundaryElementMethod,简称BEM)边界元法在求解域边界划分单元,以边界积分方程为控制方程,通过对边界分元插值离散化为代数方程求解。
边界元法所利用的微分算子基本解能自动满足无限远处的条件,因而适用于处理无限域以及半无限域问题。
对于边界变量变化梯度较大的问题,其精度更是高于有限元法。
边界元法的主要缺点是它的应用范围以存在相应微分算子的基本解为前提,对于非均匀介质等问题难以应用。
另外其所建立的系数矩阵是不对称满阵,占据相当大的内存。
1.3本文研究内容岩体渗流问题是影响工程设计.施工以及安全使用的重要因素。
当前,对于在地下工程.建筑.水利水电等领域遇到的高速非达西渗流问题,依然依据采用达西定律来进行分析设计是不够科学的,很有可能造成工程的不良后果,对于高速非达西渗流的研究已经格外迫切。
为此,本文研究的主要内容有以下几个方面:
1.从流体力学中的系统及控制体基本概念出发,依据质量守恒定理.动量守恒定理及流体的本构关系,推导了裂隙渗流的控制方程连续性方程和Navier-Stocks方程,并结合定解条件,得到了单裂隙定隙宽.变隙宽渗流层流规律;
2.介绍了层流与紊流的衡量标准,对紊流瞬时控制方程进行时均运算,得到时均连续性方程和雷诺方程,依据紊流半经验理论Prandtl混合长理论得到了裂隙高速非达西流的流速分布,并在此假设基础上进一步得到了紊流断面内平均流速与水力梯度.裂隙开度的关系;
3.根据岩体单裂隙高速非达西渗流特点,建立光滑岩体单裂隙大开度高水力梯度渗流的室内实验模型,开展高速非达西渗流水力特性的室内实验,通过实验方法寻找单裂隙高速非达西渗流规律;
4.介绍CFD软件Fluent的建模求解步骤以及求解方法,并用它建立对应于实验的数值模型进行数值计算,然后将数值模拟结果与实验结果.Lomize紊流式和速宝玉半经验公式进行比较分析。
第2章岩体单裂隙渗流理论本章从流体力学中的系统及控制体基本概念出发,依据质量守恒定律.牛顿第二定律及牛顿流体本构方程推导单裂隙渗流控制方程,并结合定解条件,寻找裂隙渗流层流以及紊流规律;在紊流半经验理论Prandtl混合长理论基础上寻找裂隙紊流断面平均流速与水力梯度.裂隙开度的关系。
2.1控制方程流体力学中,系统指由确定的流体质点所组成的流体团。
系统的边界是把系统和外界分开的边界面。
控制体是相对于所选定的坐标系来说固定的空间体积。
控制体的边界面称为控制面,控制面总是封闭面。
若代表岩体单裂隙内单位体积流体的质量.动量或其他物理量。
代表系统内流体的总质量.总动量或其他物理总量。
代表控制体内流体的总质量.总动量或其他总物理量。
.可以用积分分别表示为:
(2-1)(2-2)如图
2.1所示,有一相对于坐标系静止不动的任意控制体积,控制面用实线表示。
流体穿越控制面流动,带箭头的实线为流线。
在时刻(图a)流体系统边界与控制面重合。
时刻(图b)系统运动到新的位置(区域Ⅱ+Ⅲ),而控制体积仍留在原处(区域Ⅰ+Ⅱ)。
由图(b)可见,区域Ⅰ内的流体是时间内通过部分控制面流入控制体积的,而区域Ⅲ中的流体则是时间内穿越部分控制面流出控制体积的。
时刻,控制体积由区域Ⅰ和Ⅱ组成。
系统总量对时间的变化率按跟踪系统的拉格朗日概念应为:
(2-3)图
2.1系统及控制体示意图由图
2.1可知:
将它们代入式(2-3),进行整理:
开曲面与合成封闭曲面即控制面,于是可得:
(2-4)
2.1.1连续性方程裂隙中的流体为连续介质且作连续流动,在无汇源补给情况下,流体系统运动时其质量守恒,即。
令式(2-4)中的代表密度,有:
(2-5)由质量守恒得到:
(2-6)式(2-6)中的第二项:
那么式(2-6)可化为:
(2-7)由于控制体积是任取的,因此在任意空间点处均有(2-8)式(2-8)即为岩体单裂隙渗流的连续性方程。
在推导过程中,除考虑为单一成分的流体外,并未对流体性质及运动情况加上任何限制,因而是最一般的形式。
它表达了速度场与密度场之间相互制约的关系。
流体连续介质在运动过程中,只要连续性不破坏,连续性方程就必须满足。
如果裂隙流是恒定流,则不论流体可否压缩,均有,由式(2-8)得:
(2-9)此即为恒定流连续性方程。
2.1.2动量方程对于岩体单裂隙中的流体系统,由牛顿第二定律,可得:
(2-10)式(2-10)中:
为系统的质量力,为系统的表面力。
以代替式(2-4)中的,得:
(2-11)式(2-11)左端:
由流体连续性方程,式(2-11)左端可简化为:
式(2-11)右端:
则式(2-11)经整理得:
(2-12)式(2-12)左右两端为同一控制体积的体积分,该控制体积又是任取的,所以可得:
(2-13)式(2-13)即为岩体单裂隙渗流的动量方程。
如果裂隙流是恒定流,则由式(2-13)得:
(2-14)此即为恒定流动量方程。
2.1.3本构方程Stocks早在1845年就作出的三项假设,至今仍是讨论牛顿流体应力与应变率的关系即流体变形律或本构方程的基础。
其内容为:
(1)流体是连续的,且应力分量通常是应变率分量的线性函数;
(2)流体是各向同性的,即其性质与方向无关,因而流体变形律的表达式与坐标系的选择无关;(3)当应变率全为零,即流体静止时,变形率必须退化为流体静力条件,即(为流体静压强)。
假设
(1)是建立在实验基础上的牛顿内摩擦定律的三维推广。
假设
(2)虽对适用范围作了进一步限制,但却使流体变形律即本构方程得到很大的简化。
实际上,大多数流体是各向同性的,包含长链分子的悬浮液或溶液可能呈现某种方向性,这样的液体不在本文考虑范围内。
假设(3)应当是自然满足的,静止(或平衡)本来就是运动的特殊情况,应当包含在一般运动状态下变形律之内。
基于Stocks的三项假设便可得到牛顿流体的本构方程,式(2-15)。
许多流体力学书籍中都有其推导过程,本文不再赘述。
(2-15)将不可压缩流体的本构方程(2-15)代入动量方程(2-13)即得不可压缩流体的Navier-Stocks方程:
(2-16a)即(2-16b)连续方程(2-8)联立Navier-Stocks方程(2-16)即为裂隙渗流的控制方程。
(2-17)
2.2定解条件要解上述控制方程,还需要知道一些特定条件,即定解条件。
对于稳定渗流只需要边界条件即可,而对于非稳定渗流还需要给出初始条件。
2.2.1边界条件Navier-Stocks方程是三个椭圆型二阶偏微分方程式联立而成的方程组,因此其边界条件应是在一个封闭边界上的狄利克雷(PeterGustarLejeuneDirichlet)条件或诺伊曼条件。
常见的边界条件有两类:
压力边界条件和壁面边界条件。
(1)压力边界条件当渗流区域的某一部分边界()上的压力已知时,边界条件为:
(2-18)
(2)壁面边界条件在固体边界处,如果固体边界()的速度为,则此处流体流速为:
(2-19)
2.2.2初始条件如果研究非稳定渗流,则还需要有初始条件。
通常流场中的压力分布在初始时刻对整个流场起支配作用。
所以进行非稳定渗流计算时,必须先求得开始时刻流场的压力分布作为初始条件。
即(2-20)
2.3裂隙渗流层流规律
2.3.1立方定律如图
2.2所示,假定裂隙光滑.平直.无限长,裂隙间流体为定常层流且裂隙流体不可压缩.粘性。
裂隙开度不变,为。
取运动方向为轴,竖直向上为轴,坐标原点在中线上。
图
2.2裂隙渗流层流示意图
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 岩体单 裂隙 高速 非达西流 理论 试验 研究 数值 模拟