2非饱和水流运动基本方程精选文档PPT课件
非饱和水流运动基本方程
Darcy’s Law of Soil Water Flow in Unsaturated Zone
非饱和土壤水分运动和饱和土壤水分运动一样,水分从水势高 处向水势低处运动。一般认为,适用于饱和水流动的达西定律 在很多情况下也同样适用于非饱和土壤水分流动。 1931年,Richards最早将达西定律引入非饱和土壤水流动。非 饱和土壤水分流动的达西定律:
则 (x+x)面流速为
(vx
vx x
x);
y
vy vx
vz
Δx
Δz x
同理(y+y) 、(z+z) 面流速为:
(vy
vy y
y)
(vz
vz z
z)
z
mxin vxyzt
vy
mxout
(vx
vx x
x)yzt
mx (in
out)
vx x
xyzt
vx vz
Δx
Δz
x
同理:
y
my (in
out)
vy y
vx
kx ( )
H x
vy
k y
(
)
H y
vz
kz ( )
H z
将上式代入
t
vx x
vy y
vz z
t
x
k
x
(
)
H x
y
k
y
(
)
H y
z
k
z
(
)
H z
——非饱和土壤水运动基本方程,可简写为:
通量q h=-1h0c=0m,k=ks
k (h)
非饱和土力学
根据Ja的定义可知,通过单位面积土的空气质
Va J a a a v a t Va为通过的空气体积;va为通过的空气体积流速。 上两式相等得
量流量可用下式表示
v a k a i ay
k a D* g a
ka称为空气在土中的渗透系数
9.2.2 水流动 -广义达西定律
饱和土的达西定律
参数的测定
f c ( u a ) tg stg
' 可由饱和土的常规CU试验测定。为了 c’和
测定 '' ,应取若干相同初始孔隙比的试样进 行常含水量剪切试验。 试验中施加不同周围压 σ3,并调整ua值,以保 持所有试样的(σ3- ua)值为某一选定 的常数。 在施加附加轴向压力时,仍应随时调整ua值, 始终保持(σ3- ua)值不变,同时测读孔隙水压 力uw,直至试样剪破。于是可得一套具有相 同(σ3- ua)值、不同s=ua-uw值的极限应力圆, 如下图。
毛细粘聚力
粒子间的结合力,是影响土的抗剪强度
的重要因素之一,特别是粘性土。 然而,随着饱和土中弯液面的消失,该 力也随之消失,所以由水的表面张力产 生的粘聚力有时也称为毛细粘聚力。 大家都可能有这样的经验,在砂滩上堆 起的砂堆中挖隧道,当砂处于饱和和完 全干燥的状态时都是不可能的,只有在 适当湿的砂堆中才能容易完成。这是因 为水的表面张力即吸力产生的毛细粘聚 力在起作用。
u u a (u a uw )
Bishop (1959)的有效应力与强度
为了考虑ua 和uw 对非饱和土变形和强度特 性的影响,Bishop引进等效孔隙压力概 念,试图把适用于饱和土的有效应力原 理直接引伸到非饱和土,即
u u a (u a u w ) u a s
第三章 流体流动的基本方程 副本PPT课件
海南大学机电学院
工程流体力学
1.1 描述流体运动的方法
x=x (a,b,c,t) 1. 流体质点的位置坐标: y=y (a,b,c,t) 流体质点的运动方程
z=z (a,b,c,t)
ux(a,b,c,t)u(a,b,c,t) t
2. 速度:
vy(a,b,c,t)v(a,b,c,t) t
流体的出流
工程流体力学
1.2 描述流体运动的基本概念
一、流动的分类
•一维流动、二维流动和三维流动
1. 定义
流动参量是几个坐标变量的函数,即为几维流动。
一维流动
vv(x)
二维流动
vv(x,y)
三维流动
vv(x,y,z)
2 .实际流体力学问题均为三元流动。工程中一般根据具体情况加以 简化。
海南大学机电学院
不可能出现折点。 (3)定常流动时流线形状不变,
非定常流动时流线形状发生变化。
(4)流线簇的疏密反映了速度的大小;
强调的是空间连续质点而不是某单个质点 形成是在某一瞬间而不是一段连续时间内 表示的是质点的速度方向而不是空间位置连线
v1 交点 v2
s1
s2
v1 折点 v2
s
海南大学机电学院
工程流体力学
质量交换
动量和能 量交换
系统
无 有 单纯界面相互作用
海南大学机电学院
控制体
有
有 界面相互作用 流体流入流出作用
工程流体力学
海南大学机电学院
工程流体力学
2 连续性方程
流体连续地充满所占据的空间,当流体流动时在其内部不形成空隙, 这就是流体运动的连续性假设。
质量守恒定律(conservation of mass) :
水 文 学 原 理(五土壤水1)
z x y
(θρ w ) ( ρ wυ x ) ( ρ wυ y ) ( ρ wυ z ) = + + t x y z
单位时间内,流入控 单位时间内, 制体的水量 - 流出控 制体的水量= 制体的水量= 控制体 内土壤水的改变量
r θ = υ t
§4 土壤水运动的控制方程
3 非饱和水流运动的基本微分方程
固体密度: 固体密度:
ρ s = M s Vs
ρs = 2.6 ~ 2.7g / cm3 一般土壤
干 容 重:
Ms Ms ρb = = Vt Vs + Vw + Va
§1 土壤的质地结构及"三相"关系 土壤的质地结构及"三相"
b 表示土壤中液相比例的物理量: 表示土壤中液相比例的物理量:
质量含水率: 质量含水率: ω
§3 土壤水的能量状态
1 土水势
——土壤水所具有的势能, ——土壤水所具有的势能,称为土水势 土壤水所具有的势能
标准参照状态: 在大气压下,与土壤同温度, 标准参照状态: 在大气压下,与土壤同温度,具有固定高度的一个
假想纯自由水面的储水池. 假想纯自由水面的储水池.
土水势分析的关键点: 1 标准参照面的确定 土水势分析的关键点:
孔隙度: 孔隙度:
Va + Vw f = = Vt Va + Vw + Vs Vf
e= Vf Vs = Vf Vt V f
孔隙比: 孔隙比:
Va Va = 充气孔隙度: 充气孔隙度: f a = Vt Vs + Vw + Va
§2 土壤水的存在形态
1 土壤水作用力
分子力: 分子力: 土壤颗粒表面的分子对水分子的吸引力称为分子力 毛管力: 由于浸润性液体表面张力作用引起的水分在土壤毛 毛管力:
第2节 流体流动的基本方程PPT课件
宾汉塑性流体剪应力与速度梯度的关系
四、连续性方程
在稳态流动系统中,对直径不同的管段做物料衡算
衡算范围:取管内壁截面1-1’与截面2-2’间的管段。对于
稳态流动:
ms1 ms2
m sVsuA
u1A 1 1u2A 2 2
如果把这一关系推广到管路系统的任一截面,有:
m S u 1 A 1 1 u 2 A 2 2 u A 常 数
AA
对于圆形管道, A d 2
4
u VS d2
4
d 4VS
u
——管道直径的计算式
二、稳态流动与非稳态流动
稳态流动:运动流体的流速、压强、密度等有关物理量 仅随位置而改变,而不随时间而改变
非稳态流动:上述物理量不仅随位置而且随时间变化的流 动。
三、牛顿粘性定律与流体的粘度
1. 牛顿粘性定律
粘流指数: n>1
涨塑性流体包括玉米粉、糖溶
液、含细粉浓度很高的水浆等
0
d u /d y
胀塑性流体剪应力与速度梯度的关系
3. 宾汉塑性流体
流体的应力与应变成线性关系,但存在一屈服应力 表观粘度值为一常数
τ
0
K
du dy
粘流指数:n=1
常见的宾汉塑性流体如牙 膏、肥皂、纸浆等。
0
d u /d y
③ 一般气体的粘度值远小于液体的粘度值
④ 流体的粘度是温度T的函数
气体:T↑,粘度↑ 液体:T↑,粘度↓
?
⑤ 流体的粘度值一般不随压力而变化
流体的分类:
非饱和土力学ppt课件.ppt
• 非饱和土基本特性的学习/2、非饱和土的吸力特性
土-水特征曲线形态的重要参数
由于土中的水分可以有 结晶水、吸着水、结合水(薄膜水)和自由水等
具有不同属性的不同类型。 含水量变化时,土中水有不同的类型,气有不同的连通,
孔隙水压力和孔隙气压力分别在土的孔隙水体 和孔隙气体中是各向等压的静水压力型应力
孔隙水压力和孔隙气压力 各自作用在其与土颗粒接触部分的表面上, 其差值对土骨架的作用不会是各处相等的。
当孔隙水为弯液面环状水时,吸力只在接触点的 法向上作用;当孔隙水为有弯液面的体积水时, 所产生的吸力必然有法向和切向两个方向上分力 的作用。国内也出现了湿吸力与牵引力的提法(汤连生)。
单一有效应力型的应力状态变量
人们在寻求非饱和土的应力状态变量时,首先想到了 单一有效应力型的应力状态变量
它不是一般的纯力学量,而是一个材料有关的力学量,与材料 的本构关系有着密切的联系(如饱和土力学中的有效应力)。 研究提出具有真实合理性的有效应力表达式是当前的主要任务。
对已经提出的各种表达式还需要作出认真的选择与检验。
导致了非饱和土十分复杂的力学性质。
• 非饱和土基本特性的学习/2、非饱和土的吸力特性
2、非饱和土的吸力特性
非饱和土的土水势一般包括 温度势、压力势、重力势、基质势和溶质势
在等温、等压、等高(不计重力)的情况下, 土中水的温度势、压力势、重力势保持不变,
自由能的变化只有基质势和溶质势的变化。
如将它们分别称之为基质吸力和溶质吸力, 它们之和,即此时的自由能,称为总吸力,则有
应该取决于各自的相对压缩性。
在孔隙流体不能排出的条件下,土受力后的孔隙水压力 和孔隙气压力的增量是一种超孔隙压力
非饱和水流基本方程2-方程及通量法
第二章非饱和水流基本方程直角坐标非饱和水分基本方程柱坐标系非饱和水分基本方程水分运动基本方程的定解条件 2.1.1基本方程推导z理论基础:)达西定律;)()(∂∂−=∂∂−=xhk x H k v x θθz h H nHk v +=∂∂−=;)(θUnsaturated Flow32.1直角坐标非饱和水分基本方程2.1.1基本方程推导z 取微分单元体,体积:Δx Δy Δzz设(x,y,z)流速为:v x ,v y ,v z ,z右面(x+Δx)流速为;;)( )(z z vv y y v v z z y y Δ∂∂+Δ∂∂+;)(x xv v x x Δ∂∂+Unsaturated Flow4tz y v x x ΔΔΔ=ρ方向流入质量t z y x zv y v x v zy x ΔΔΔΔ⎟⎟⎠⎞⎜⎜⎝⎛∂∂+∂∂+∂∂−=−ρout in outin m -m m m 流出质量流入质量t z y x x v v x xx ΔΔΔΔ∂∂+=)(ρ方向流出质量tz y x x vx x ΔΔΔΔ∂∂−=−ρ质量流出流入方向)(tz y x z vz t z y x y v y z yΔΔΔΔ∂∂−=−ΔΔΔΔ∂∂−=−ρρ质量流出流入方向质量流出流入方向)()(水质量为:根据质量守恒原理有t z y x tm ΔΔΔΔ∂∂=Δθρtz y x z v y v x v z y x ΔΔΔΔ⎟⎟⎠⎞⎜⎜⎝⎛∂∂+∂∂+∂∂−=ρout in m -m []H k tz H k z y H k y x H k x t ∇∇=∂∂⎥⎦⎤⎢⎣⎡∂∂∂∂+⎥⎦⎤⎢⎣⎡∂∂∂∂+⎥⎦⎤⎢⎣⎡∂∂∂∂=∂∂)()()()(θθθθθθ可简写为:代入上式得:z k z h x t yy ∂∂±⎥⎦⎤∂∂⎢⎣∂∂∂=∂)()θθ1。
以含水量为变量的垂向一维流基本方程)()()(向上为正z k h k ∂+⎟⎞⎜⎛∂∂=∂θθθzk z h k z y h k y x h k x t ∂∂±⎥⎦⎤⎢⎣⎡∂∂∂∂+⎥⎦⎤⎢⎣⎡∂∂∂∂+⎥⎦⎤⎢⎣⎡∂∂∂∂=∂∂)()()()(θθθθθz∂+⎟⎠称为扩散度。
流体运动的基本概念与基本方程课件
人类如何与动物和自然环境共存?
随着科技的进步和人类文明的发展,人口数量不断增加,对自然资源和动物栖息地的需求也越来越大,造成了严重的环境问题。
因此,如何与动物和自然环境共存成为摆在我们面前的一个大问题。
本文将从以下几个方面展开讨论。
一、保护自然环境
自然环境是人类赖以生存的基础,保护自然环境是维护人类生存的关键。
人类应该减少对自然环境的破坏,降低碳排放量,减少污染,垃圾分类,开展环保宣传等措施来保护自然环境。
二、保护动物栖息地
许多动物物种的栖息地遭到破坏,导致它们的生存面临威胁。
如果我们不保护动物的栖息地,那么这些动物的数量将会逐渐减少,最终可能会走向灭绝。
保护动物栖息地是非常必要的,这不仅仅是对动物本身的保护,也会保护我们的生态环境。
三、合理利用自然资源,是指在保护自然环境的基础上,科学地开发和利用自然资源。
尽可能地减少浪费,提高资源的利用效率,这有利于稳步地推进人类文明的发展,同时也有利于保护自然环境。
四、环保教育的普及
环境保护是人类的共同责任。
在进行环保教育时,我们应该尝试将其融入各个领域,如在学校、工厂、商场等,让更多的人们了解环境保护的重要性,并树立起环境保护意识。
五、遵守相关法规
政府对环境保护和动物保护有一定的法规和规定。
对于不遵守规定的人员,要进行相应的惩罚以保障环境和动物的权益。
因此,在正常的生产和生活中,人们也要遵守相关法规,并自觉地履行自己的环保义务。
总之,人类与动物和自然环境的共生是非常重要的,我们应该从各个方面出发,不断探索更好的共生方式。
只有这样,我们才能实现美好的环境和社会。
描述液体运动的两种方法及液体运动的基本概念ppt课件
dt
式中, (x, y, z, t ) : 欧拉变数 (ux uy uz) : 通过固定点的流速分量
.
拉格朗日法
(a, b, c) :
t
:
(x, y, z) :
(a, b, c , t ) :
质点起始坐标 任意时刻 质点运动的位置坐标 拉格朗日变数
欧拉法
(x, y, z) : t:
(x, y, z , t ) :
u
y
dy ( a , b , c , t ) dt
u
z
dz
(a , b , c , t ) dt
.
u
x
dx
(a , b , c , t ) dt
u
y
dy
(a , b , c , t ) dt
u
z
dz
(a , b , c , t ) dt
d 2 x (a,b,c,t)
ax
dt 2
a y
.
上式对t 求导,得到液体质点的速度
x x(a,b,c,t)
ux
d
x(a, b, c, t ) dt
d dt
ቤተ መጻሕፍቲ ባይዱ
y
y(a, b, c, t )
z z(a,b,c,t)
uy
d y(a,b,c,t) dt
d z(a,b,c,t)
uz
dt
.
x x (a ,b,c,t )
u
x
d
x (a , b , c , t ) dt
.
z 空间固定点
O
x
考察不同固定点上、不同液体质点通过时的运动情 况。这句话包含两层意思:
.
z 空间固定点
2-非饱和水流运动基本方程PPT课件
呈直线关系,但其斜率(即水力传导度)随平均负压而变。
通量q h=-1h0c=0m,k=ks
h=-30cm h=-50cm
负压梯度△h/△x
k (h)
.
h
11
容水度(或比水容量)
单位基膜势(负压值)变化所引起土壤含水率的变化,一 般称为容水度或比水容量(C),可以下式表示:
Ch d
dh
表示在单位压力水头降低时自单位体积土壤中所释放 出来的水的体积,它与饱和土壤的给水度相似。 用测水分特征曲线的方法来测定
y)
(vz
vz z
z)
.
16
设六面体土壤含水量为,则t内六面体内土壤水质量变 化量为:
m xyzt
t
根据质量守恒原理有 p
min
- mout
vx x
v y y
vz z
xyzt
xyzt
t
即
t
vx x
vy y
vz z
.
17
根据达西定律有:
vx
kx ( )
H x
vy
ky ( )
H y
vz
壤水通量q可由l端补给量或2端
溢出量测得,两者相等时,水流 处于稳定状态。
➢ 非饱和土壤水力传导度可由达西 定律求得。
k(h) q H
q H
qX H
X
△x 多孔板
2
1
H1
平水 装置
H2
△H
计算的k是平均的k
试样中各点的是不相同的
.
10
k(h) q H
q H
qX H
X
在不同的平均负压(吸力)值下,通量与负压梯度成正比,两者
[理学]3第三章-水动力学基础PPT课件
《水力学》
第三章 水动力学基础 §1 描述液体运动的两种方法
1 描述液体运动的两种方法
拉格朗日法
z
以研究单个液体质点的运动过程作为基 础,通过对每个液体质点运动规律的研 究来获得整个液体运动的规律性。
t
(x,y,z)
(t0)
O M (a,b,c)
x
x x(a,b,c,t) y y(a,b,c,t) z z(a,b,c,t)
第三章 水动力学基础
§4 恒定总流的能量方程
4 恒定总流的能量方程
方程应用的注意点:
z1p 12 1V g12z2p 22 2V g22hw
选取高程基准面; 选取两过水断面;
所选断面上水流应符合渐变流的条件,但两个断面之
间,水流可以不是渐变流。
选取计算代表点;
选取压强基准面;
动能修正系数一般取值为1.0。
Yangzhou Univ
-
19
《水力学》
第三章 水动力学基础
§4 恒定总流的能量方程
例题:采用直径d=350mm的虹吸管将河水引入堤外供给灌溉。已知H=4m, h=3m,若不计损失,试确定该虹吸管的输水量,并计算图中1、2、3、4各 点的位置水头z、压强水头p/γ、流速水头V2/2g及总水头。
-
27
《水力学》 能量方程(伯努利方程)应用的拓展
恒定总流的能量方程(又称为伯努利方程)是 应用最广的流体动力学基本方程。它不仅可以 用于液体运动,还可以应用于气体运动。
特别提醒
伯努利方程式是对不可压缩
液体导出的,而气体是可压缩流体,因此,将
伯努利方程式应用于气流时,要求流速不很大
(V<50m/s),压强变化不很大,密度的变
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
同理:
y
my(i
nou )tvyxyzt
y
m z(in o)u t v zz x y z t
m in-m out v xx v yy v zz x y z t
设六面体土壤含水量为,则t内六面体内土壤水质量变 化量为:
mxyzt
t
根据质量守恒原理有
m in-m ou t v x x v y y v zz x y z t
第二章 非饱和土壤水流的基本方程
Chapter.2 Basic Equations of Soil Water Dynamics
§1 非饱和土壤水流的达西定律 §2 土壤水分运动参数及其测定方法 §3 非饱和土壤水运动的基本方程 §4土壤水分运动基本方程的定解条件 §5 土壤水分通量法
§1 非饱和土壤水流的达西定律
ks ≥k(θ)
θ ↗⇒k(θ ) ↗
§2 土壤水分运动参数及其测定方法
Parameters of Soil Water Movement and its measurement methods
3.1 非饱和/饱和水力传导度 3.2 容水度 3.3 土壤水分扩散度
3.1 非饱和/饱和水力传导度
非饱和水力传导度及其测定 饱和水力传导度及其测定
Guelph土壤入渗仪
3.2 容水度(或比水容量)
单位基膜势(负压值)变化所引起土壤含水率的变化,一 般称为容水度或比水容量(C),可以下式表示:
Ch d
dh
表示在单位压力水头降低时自单位体积土壤中所释放 出来的水的体积,它与饱和土壤的给水度相似。 用测水分特征曲线的方法来测定
3.3 土壤水分扩散度D
则 (x+x)面流速为
(vx
vx x
x);
y
vy vx
vz
Δx
Δz x
同理(y+y) 、(z+z) 面流速为:
(vy
vy y
y)
(vz
vz z
z)
z
m xi n vxyzt
m xou t(vx v xx x) y z t m x(in o)u t v xx x y z t
vy vx
vz
Δx
Δz
Unsaturated Zone
理论基础:
➢ 达西定律 qk(h)H
➢ 质量守恒定律(水流连续原理)
3.1 直角坐标非饱和水分基本方程 3.2 柱坐标系非饱和水分基本方程(自学)
3.1直角坐标非饱和水分基本方程
基本方程推导
z
取微分单元体,体积:
xyz
设沿x、y和z方向流速分别为:
vx,vy,vz,
Darcy’s Law of Soil Water Flow in Unsaturated Zone
非饱和土壤水分运动和饱和土壤水分运动一样,水分从水势高 处向水势低处运动。一般认为,适用于饱和水流动的达西定律 在很多情况下也同样适用于非饱和土壤水分流动。 1931年,Richards最早将达西定律引入非饱和土壤水流动。非 饱和土壤水分流动的达西定律:
q K ( m ) 或 q K ( )
饱和土壤水分流动的达西定律:
qKsH
水势组成: 流动准则:
ψg :
饱和流
非饱和流
ψ =ψg +ψp
总水头
高
低
ψ =ψg +ψm
土水势
高
低
相对参考平面之高度
ψp:
至地下水面的高度
ψp= 0
ψm :
ψm = 0
ψm 取决于土壤的干湿程度
K:
ks,渗透系数=const. k(θ ) 为土壤含水率的函数
扩散率D
土壤水分扩散度为单位含水率梯度下,通过单位面积的 土壤水流量,其值为土壤含水率的函数,即
DkhC k
土壤水分扩散度与土壤的关 系可用以下经验公式表示
Dabe
土壤体积含水量θ
§3 非饱和土壤水运动的基本方程 Basic Equations of Soil Water Movement in
非饱和水力传导度的测定
非饱和达西实验
➢ 在水平土柱两端有多孔板,分别 由平水箱保持一定水位,使其负 压为h1和h2,在梯度作用下,土 柱中土壤水从l端向2端运移。土
壤水通量q可由l端补给量或2端
溢出量测得,两者相等时,水流 处于稳定状态。 ➢ 非饱和土壤水力传导度可由达西 定律求得。
k(h)qHqHqH X X
水率降低而减小。
不同质地土壤的非饱和 水力传导度
➢ 具有大孔隙粗质土壤,在吸力 作用下孔隙中水分很快排除, 导水率迅速下降;
➢ 而粘质细颗粒土壤,在较高吸 力下,许多小孔隙仍充满水, 仍具有一定的导水性能,导水 率下降较缓慢。
➢ 所以,同一吸力条件下,粘性 土的导水率可以大于砂性土的 导水率。
土壤非饱和导水率K
△x 多孔板
2
1
H1
H2
△H
平水 装置
计算的k是哪点的k? 试样中各点的是否相同?
通量 q h=-10hc=m0,k=k s
k(h)qHqHqH X X
在不同的平均负压(吸力)值下,通量与负压梯度成正比,两者 呈直线关系,但其斜率(即水力传导度)随平均负压而变。
k (h)
h=-30cm h=-50cm
3.1.1 非饱和水力传导度及其测定
非饱和水力传导度的概念及特征
是反映土壤水分在压力水头差作用下流动的性能。一般在 饱和土壤中导水率称为渗透系数,为常量。
定义:在单位水头差作用下,单位断面面积上流过的水流
通量。
qk(h)H
k(h) q H
qk()H
k() q
H
在非饱和土壤中,导水率是负压或含水率的函数,随着含
砂质土 粘质土
土壤吸力S
Parametric Models for K(θ) \K(h)
K(h)
k(h)ksech
k(h) ahn a
k(h) h n b
van Genuchten-Mualem model:
K(θ)
k()ksecs
k()kBrooks and Corey (BC) model:
h
负压梯度△h/△x
3.1.2 饱和水力传 导度及其测定
双环入渗仪
The assumption is that the soil layer immediately below the ponded area is fully saturated and thus the matric potential is essentially zero. Common Steady Flow Analysis (Unit Gradient): Accounts only for the flow component due to gravity.