高等流体力学_第一讲.共40页
高等流体力学_第一讲.
曲面所围体积之比的极限值；
div
a
lim
a
S
nds
V 0 V
封闭曲旋线度所（张cu的rl面or积r比ota值tio的n极）限：；向量场中围绕一点的封闭曲线的环量与该
a dl
rot a lim
S S 0
北京工业大学市政学科部——马长明
高等流体(水）力学讲稿
课程简介
一、课程名称:
高等流体力学——水利水电工程 高等水力学——给排水工程（土木工程）
——学什么？
二、教材：
1、高等流体力学？天津大学——新世纪研究生适用教材
相对于本科“水力学”或“流体力学”，在相关问题上进行更深入的理论分析 和论述，以满足现代水力工程对流体力学的要求，有助于提高理论修养，深入理解现代 流体力学的内容。是水力工程以及学科各硕士专业的学位课。
8、地下水中的弥散
北京工业大学市政学科部——马长明
高等流体(水）力学讲稿
2
课程简介
三、内容
环境流体力学——董志勇（共8章）
1、绪论； 2、迁移扩散理论； 3、剪切流离散； 4、射流、羽流和浮射流； 5、水质模型； 6、地下水污染模型； 7、分层流； 8、生态水力学引论。
北京工业大学市政学科部——马长明
五、教、学与评价探讨
课程特点： 1）要求数学知识多；方程、公式多，推演论证繁琐；解题 难度大。 2）学时少（32），所留自学时间也少，而教学内容多。
北京工业大学市政学科部——马长明
高等流体(水）力学讲稿
5
数学基础知识
一、正交曲线坐标系
1、直角坐标系、柱坐标系与球坐标系 1）坐标线与坐标面 2）坐标系间的转化
高等流体力学课件 高等流体力学(1)
ti ijn j
aij bikckj d eij kk
在方程同一项中重复出现的指标称为哑指标，哑指标在作 求和运算后就消失了，因此改变哑指标的字母不改变表达 式的内容。
在方程同一项中只出现一次的指标称自由指标，在同一 方程的所有项中出现的自由指标必须相同。
为避免混淆，同一项中相同指标出现的次数不能多于2。
散度的性质： diva 0 为无源场，无源场的性质见P14。自己看书。
4.旋度
矢量 a 沿任意曲线L的线积分，即 a 沿任意曲线L的环量：
a • dr axdx aydy azdz
L
L
a • dr
设张于L上的曲面为S，则定义 lim L
为矢量 a
S0 S
的旋度矢量 rota 在法线 n 上的投影，
12月31日,1月8日：粘性不可压缩流体的流动
12月31日（第13周）
N-S方程的精确解
1月8日（第14周）
小雷诺数流动的近似解
1月15日（第15周）
大雷诺数下的边界层理论
绪论
人类生活在流体环境中，人们对一些流体运动 现象却缺乏认识，比如：
1. 高尔夫球 ：表面光滑还是粗糙？ 2. 汽车阻力 ：来自前部还是后部？ 3. 机翼升力 ：来自下部还是上部？
(1) v w； (2) v w； (3) v v；
(4) e1 v； (5) e2 v； (6) r v, r xi yj zk 是位置矢量。
解:
v w vi wi
v1w1 v2 w2
v3w3
1 3 2 1 5 1 4
i jk v w 1 2 5 i (2 1 1 5) j(3 5 11) k (11 2 3)
高等流体力学
高等流体力学课件 第一章 流体力学的基本概念
J 0
x y z x0 x0 x0 J x y z 0 y0 y0 y0 x y z z0 z0 z0
有限大的正数
rr r0 , r
互为反函数。
§1.1 拉格朗日参考系和欧拉参考系
三、两个参考系间的相互转换
2.两个参考系间的相互转换
若已知流线经过点 (x0,y0,z0) ,则参数方程的初始条件可定为，
《高等流体力学》电子课件
§1.1 拉格朗日参考系和欧拉参考系
一、拉格朗日参考系
1.流动的描述
着眼于流体质点。 描述每个流体质点自始至终的运动，即位置随时间的变化。
r r r r(x 0,y0,z0,t) 式中x0 , y0 , z0 是t =t0 时刻流体质点的空间坐标，用来区分不同的流体质点。
二、流线
1.定义
某时刻，流场中的一条曲线，曲线上各点的速度矢量方向和曲线在 该点的切线方向相同。
2.流线方程的微分方程
d r d i x d j y d k z u u i v j w k
i dru dx u
j dy v
k dz0 w
2.流动物理量随时间的变化
加速度：
ai
ui t
uj
ui xj
其他物理量：
d dt
t uj
xj
dp dt
p t
uj
p xj
§1.1 拉格朗日参考系和欧拉参考系
三、两个参考系间的相互转换
1.两个参考系间相互联系——雅可比行列式
0 初始时刻流体微团体积 T时刻变形后流体微团体积
1.流动的描述
着眼于空间点。 描述流过每个空间点上的流体质点的运动。
高等流体力学讲义
高等流体力学授课提纲第一章概论§1.1 流体力学的研究对象§1.2 流体力学发展简史§1.3 流体力学的研究方法§1.3.1 一般处理途径§1.3.2 应用数学过程§1.3.3 流体力学方法论：一般方法§1.3.4 流体力学方法论：特殊方法●Lagrange描述和Euler描述●无量纲化●线性化●分离变量法●积分变换法●保角映射法●奇点法（孤立奇点法、分布奇点法、Green函数法）●控制体积法●微元法第一章概论§1.1 流体力学的研究对象（1）物质四态：●四态：固态—液态—气态—等离子态；等离子体=电离气体●界限：彼此无明确界限（高温下的沥青；冰川），取决于时间尺度；●流体力学的具体研究对象：液体、气体、等离子体（电磁流体力学、等离子体物理学）；●液体与气体的差别：液体—有固定容积、有自由面、不易压缩、有表面张力；气体—无固定容积、无自由面、易压缩、无表面张力。
（2）流体的基本性质：易流动性：静止流体无剪切抗力；压缩性（膨胀性）：压差、温差引起的体积改变，判据：马赫数；粘性：运动流体对剪切的抗力，判据：雷诺数；热传导性：温差引起的热量传递，普朗特数。
（3）流体的分类：i）按有无粘性、热传导性分：真实流体（有粘性、有热传导、与固体有粘附性无温差）；理想流体（无粘性、无热传导、与固体无粘附性有温差）；ii）按压缩性分：不可压缩流体，可压缩流体；iii）按本构关系分：牛顿流体（牛顿粘性定律成立），非牛顿流体（牛顿粘性定律不成立），下分纯粘性流体（拟塑性流体，涨塑性流体）；粘塑性流体（非宾汉流体、宾汉流体）；时间依存性流体（触变流体、振凝流体）；粘弹性流体拟塑性流体（剪切流动化流体）：剪切应力随剪切速度增加而减小，如淀粉浆糊、玻璃溶液、高分子流体、纤维树脂；涨塑性流体（剪切粘稠化流体）：剪切应力随剪切速度增加而减小，如淀粉中加水、某些水-砂混合物；粘塑性（非宾汉和宾汉流体）：存在屈服应力，小于该应力无流动，如粘土泥浆、沥青、油漆、润滑脂等，所有粘塑性流体为非宾汉流体，宾汉流体为近似；触变流体（摇溶流体）：粘性或剪切应力随时间减小，如加入高分子物质的油、粘土悬浊液；振凝流体：粘性或剪切应力随时间增大，如矿石浆料、膨润土溶胶、五氧化钒溶液等；粘弹性流体：兼有粘性和弹性性质的流体，能量不像弹性体守恒，也不像纯粘性体全部耗散。
高等流体力学第1讲
第一讲绪论一、参考教材1.流体力学，周光炯等编写，高等教育出版社2.流体力学，吴望一编写，北京大学出版社3.流体力学的先期课程：数学（微积分、线性代数、复变函数、数理方程、场论、张量分析、数值分析、偏微分方程数值解法乃至泛函分析等等）、力学（分析力学）基础。
二、流体力学的研究方法实验方法:同物理学等其它的自然科学学科的研究方法一样，非牛顿流体力学的研究方法包括理论方法和实验方法。
理论方法就是根据流动的物理模型和物理定律建立描写流体运动规律的封闭方程组以及相应初始条件和边界条件，运用数学方法准确或近似地求解流场，揭示流动规律；实验方法就是运用模型实验理论设计试验装置和流程，直接观察流动现象，测量流体的流动参数并加以分析和处理，然后从中得到流动规律。
在非牛顿流体力学的发展过程中，实验方法是最先采用的方法，也是最基本的方法。
即使到现在，不使用实验方法，航空航天、大型水利枢纽、聚合物驱油等复杂系统的研究几乎是不可能的。
实验方法主要包括以下几个步骤：○1运用相似理论，针对具体的研究对象确定相似准数和相似准则；○2依据模型律来设计和制造模型，确定测量参数，选择相应的仪器仪表，建立实验装置；○3制定实验方案并进行实验，观察流动现象，测量流动参数；○4运用量纲分析等方法整理和分析实验数据，与其它方法或著作所得的结果进行比较，从中总结出流动规律。
实验研究方法的优点：能够直接解决工程实际中较为复杂的流动问题，能够根据观察到的流动现象，发现新问题和新的原理，所得的结果可以作为检验其他方法的正确性和准确性。
实验研究方法的缺点主要是对于不同的流动需要进行不同的实验，实验结果的普遍性稍差。
解析方法:解析方法是非牛顿流体力学各种研究方法中最为准确的和最为理想的方法。
解析方法主要包括：○1详细分析问题的物理学本质，通过适当的简化建立物理模型；○2运用物理定律建立数学模型，通常是建立起微分方程或微分方程组，确定流动方程边界条件和初始条件；○3运用数学方法求解出流动方程的解析解；○4列举计算实例，然后再与其他方法所得的结果进行比较，以检验物理模型和数学模型的合理性。
高等流体力学PPT课件
2
aij ijkk
uD S r
表示由于流体微团变形而产生的 M 点相对于M点 的速度变化。
uR
1 ur
2
表示由于流体微团绕瞬时轴旋转而产生的 对于M 点的速度变化。
M 点相
26
26
欧拉和拉格朗日参考系中的时间导数
欧拉参考系：
u t x,y,z
u
u(x,
y,
z,
t)
某一空间点上的流体速度随时间的变化，称当地导 数或局部导数。
拉格朗日参考系：u u(x0, y0, z0,t)
u
t
x0 , y0 ,z0
流体质点速度随时间的变化，即加速度。
在欧拉参考系下用 Du 表示流体质点的速度变化。
25
速度分解定理，应变率张量和旋转率张量
速度分解定理
ui
ui x j
xj
1 2
ui x j
u j xi
1 2
ui x j
u j xi
xj
sij x j aij x j S r A r
Sr 1 ur
2
u uD uR
aij x j ijk x jk r
物质导数
以矢量和张量下标形式表示的物质导数
D
Dt
t
uk
xk
D
Dt
t
u
t
u
算符
u
ui vj wk
i
x
j
y
k
z
u v w x y z
13
13
物质导数物理意义
D Dt t uk xk
D 物质导数，质点导数，随体导数；
Dt
欧拉参考系中的时间导数，称局部导数或就地导数，表示空
高等流体力学第一讲.ppt
v v v v a b (a2b3 a3b2 )e1 (a3b1 a1b3 )e2 a1 a2 v (a1b2 a2b1 )e3 b1 b2
3
v e1
v e2ቤተ መጻሕፍቲ ባይዱ
v e3 a3 b3
第一讲，附录部分：数学基础
二、场的概念，梯度及方向导数
v v v ai a1e1 a2e2 a3e3
a11 a ij a 21 a 31
a12 a 22 a 32
a13 a 23 a 33
2.求和约定
①在同一项中如有两个指标相同时，就表示对该指标从1到3求和:
aibi a1b1 a2b2 a3b3
n为自由指标 m为哑指标
北京工业大学市政学科部——马长明
高等流体(水）力学讲稿
10
第一讲，附录部分：数学基础
3.张量的基本运算规则
（1）克罗内克（Kroneker）符号δ
ij
1 i j ij 0 i j
是二阶单位张量。 符号具有以下重要性质：
v v ij ei e j
两矢量的点积可表示为：
ai bj aiei bj e j aibjij aibi a jbj
11
第一讲，附录部分：数学基础
1 i j 符号具有以下重要性质： ij 0 i j
ij jk i11k i 22k i33k
12
第一讲，附录部分：数学基础
（2）里奇（Ricci）置换符号ε
ijk
ijk
1 1 0
偶排列，即：123,231,312； 奇排列，即：213,321,132 有两个或两个以上指标相同。
高等计算流体力学讲义(1)
(8)
∂φ ∂φ ∂ 2φ ∂ 2φ ∂ 2φ ∂ 2φ = ξ xx + η xx + ξ x [ 2 ξ x + ηx ] +ηx[ ξx + 2 ηx ] ∂ξ ∂η ∂ξ ∂ξ∂η ∂ξ∂η ∂η ∂φ ∂φ ∂ 2φ ∂ 2φ ∂ 2φ 2 = ξ xx + η xx + 2 (ξ x ) + 2 ξ xη x + 2 (η x ) 2 ∂ξ ∂η ∂ξ ∂ξ∂η ∂η
2、度量系数及其计算方法
在导数的坐标变换公式中涉及到下列坐标变换系数： ξ x , ξ y ,η x ,η y 。这些系数 称为坐标变换公式(5)对应的度量系数（metrics）。我们看到，为了求解计算平 面中的偏微分方程，如(9)式，必须确定度量系数（有时还包括 ξ xx , ξ xy , ξ yy ,η xx ,η xy ,η yy 等）的离散值。那么，这些度量系数如何计算呢？由于一 般情况下，我们只知道坐标变换关系(5)、(6)的离散表达式，度量系数一般也要 通过有限差分方法近似计算。但是，直接构造 ξ x , ξ y ,η x ,η y 的差分近似是不容易 的。以 ξ x 为例，根据偏导数的意义， ξ x 为 y 保持不变时 ξ 随 x 的变化，如图 2 所示，网格点 P 处的 ξ x 的计算公式应为：
不计质量力的情况下，在直角坐标系中，守恒型 N－S 方程可以写为下列 向量形式： ∂U ∂ ( F − Fv ) ∂ (G − G v ) ∂ ( H − H v ) + + + =0， (1) ∂t ∂x ∂y ∂z 其中
ρu ρv ρw ⎞ ⎛ ⎞ ⎛ ⎞ ⎛ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ 2 ⎜ ρu + p ⎟ ⎜ ρ vu ⎟ ⎜ ρ uw ⎟ F = ⎜ ρ uv ⎟ G = ⎜ ρ v 2 + p ⎟ H = ⎜ ρ vw ⎟ ， ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ 2 ⎜ ρ uw ⎟ ⎜ ρ vw ⎟ ⎜ ρw + p ⎟ ⎜ ( ρ E + p)u ⎟ ⎜ ( ρ E + p )v ⎟ ⎜ ( ρ E + p) w ⎟ ⎝ ⎠ ⎝ ⎠ ⎝ ⎠ 0 ⎛ ⎞ 0 ⎛ ⎞ ⎜ ⎟ ⎜ ⎟ τ xy τ xx ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ τ ⎜ ⎟ τ xy yy G = Fv = ⎜ v ⎜ ⎟， ⎟ τ τ yz xz ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ T ∂ ∂T ⎜ uτ xy + vτ yy + wτ yz + k ⎟ ⎜ uτ xx + vτ xy + wτ xz + k ⎟ ∂y ⎠ ∂x ⎠ ⎝ ⎝ 0 ⎛ ⎞ ⎜ ⎟ τ xz ⎜ ⎟ ⎜ ⎟ τ zy Hv = ⎜ ⎟。 τ zz ⎜ ⎟ ⎜ ∂T ⎟ ⎜ uτ xz + vτ zy + wτ zz + k ⎟ ∂z ⎠ ⎝ 如果忽略 N－S 方程中的粘性和热传导，得到的简化方程为 Euler 方程：