流体力学:第5章 势流理论-上
流体力学课件(全)
Y 1 p 0 y
欧拉平衡方程
Z 1 p 0 z
p p( , T )
t
1 V V T p
1 V V p T
p p(V , T )
1 t T p
p
p
1 p T
V
p y = pn pz = pn
px = p y = pz = pn = p
28/34
第二章
流体静力学
§1 静压强及其特性 §2 流体静力学平衡方程 §3 压力测量 §4 作用在平面上的静压力 §5 作用在曲面上的静压力 §6 物体在流体中的潜浮原理
29/34
§2流体静力学平衡方程
通过分析静止流体中流体微团的受力,可以建立 起平衡微分方程式,然后通过积分便可得到各种不同 情况下流体静压力的分布规律。 why 因此,首先要建立起流体平衡微分方程式。 现在讨论在平衡状态下作用在流体上的力应满足 的关系,建立平衡条件下的流体平衡微分方程式。
《流体力学》
汪志明教授
5/24
第一章 流体的流动性质
§1 流体力学的基本概念
§2 流体的连续介质假设 §3 状态方程 §4 传导系数 §5 表面张力与毛细现象
《流体力学》
汪志明教授
6/24
§2 流体的连续介质假设
虽然流体的真实结构是由分子构成,分子间有一定的孔隙,但流 体力学研究的并不是个别分子微观的运动,而是研究大量分子组成的 宏观流体在外力的作用下所引起的机械运动。 因此在流体力学中引入连续介质假设:即认为流体质点是微观上 充分大,宏观上充分小的流体微团,它完全充满所占空间,没有孔隙 存在。这就摆脱了复杂的分子运动,而着眼于宏观机械运动。
流体力学第五章 量纲分析和相似理论
第五章 量纲分析与相似原理
5.2 量纲分析与П定理
2. П定理
提议用量纲分析的是瑞利(L.Reyleigh,1877),奠定理论基础的是美国物理
学家布金汉(E.Buckingham,1914):
Π定理
若某一物理过程包含 n 个物理量,即:
f(q1 , q 2,q 3, ……, q n )=0
其中有 m 个基本量(量纲独立,不能相互导出的物理 量),则该物理过程可由 n个物理量构成的 n-m 个无 量纲的关系表达式来描述。即:
5.1 量纲与物理方程的量纲齐次性
1. 物理量的量纲(因次):物理量的本质属性。
2. 物理量的单位:物理量的度量标准。
基本量纲和导出量纲:根据物理量之间的关系把无 任何联系且相互独立的量纲作为基本量纲,可由基本量 导出的量纲为导出量纲。
SI制中的基本量纲:
dim m = M , dim l = L , dim t = T ,dim θ=Θ
第五章 量纲分析与相似原理
5.1 量致性原则,也叫量纲齐次性原理(量纲和谐原理)
物理方程可以是单项式或多项式,甚至是微分方程等,同 一方程中各项的量纲必须相同。
用基本量纲的幂次式表示时,每个基本量纲的幂次应相等,
这就是物理方程的量纲一致性原则,也叫量纲齐次原则或量纲
1. 客观性 2. 不受运动规模的影响 3. 可以进行超越函数运算
整理课件
第五章 量纲分析与相似原理
5.1 量纲与物理方程的量纲齐次性
2. 量纲一的量(无量纲量)
基本量独立性判别条件:
设A、B、C为三个基本量,他们成立的条件是:指数行列式 不等于零。
diB m M 2L 2T 2 diA m M 1L 1T1 diC m M 3L 3T 3
流体力学第五章
5.2 边界层流动
5.2 边界层流动
*
0
u 1 u e e
dy
5.2 边界层流动
**
0
u eue
u 1 u dy e
5.2 边界层流动
平面边界层流动方程
边界层近似假定 1. 纵向偏导数远小于横向偏导数
5.2 边界层流动
边界层分离
理想流体能量转换过程 边界层内粘性对机械能的耗散使得流体微团在逆 压区 MF 段间的某个点处 V 降为零,后来的质点 将改道进入主流区,使来流边界层与物面分离; 在分离点下游区域,受逆压作用而发生倒流。
5.2 边界层流动
边界层分离
分离点:紧邻壁面顺流区与倒流区分界点。 边界层分离的必要条件:粘性、逆压梯度。
湍流边界层摩阻系数大
0.664 C fL Re x
C fT
0.0576 /5 Re 1 x
5.2 边界层流动
边界层分离
边界层流动:流体质点受惯性力、粘性力和压力 作用;粘性力阻滞流体质点运动,使流体质点减 速和失去动能;压力的作用取决于绕流物体形状; 顺压梯度有助于流体加速前进,而逆压梯度阻碍 流体运动。
研究方法:实验、数值(RANS、LES、DNS)
5.1 粘流的基本特性
层流、紊流速度型 紊流粘性应力比层流大
5.2 边界层流动
边界层概念的提出
高 Re流动,惯性力远大于粘性力,研究忽略粘 性的流动有实际意义。 阻力、分离、涡扩散等问题,无粘解与实际相 差甚远。 研究表明:虽然 Re很大,但在靠近物面的薄层 流体内,沿物面法向存在很大的速度梯度,粘 性力与惯性力相当而不可忽略。 Prandtl把物面附近粘性力起重要作用的薄层称 为边界层。
流体力学:第5章势流理论-上
c1
c2
5.2.1 复势的可叠加性 解析函数 W1(z) 1 i1 W2 (的z) 线性2 组i合2 ,
W (z) W1(z) W2 (z)
仍然是解析函数,仍然代表某一种流动的复势。简单 流动组合成复杂流动——叠加法
5.3 平面势流的基本解
目的:求解最简单的流动,为解决复杂势流奠定基础。 内容:均匀流、点源、点涡、偶极。
v 0 (R )
5.1.3 初始条件(initial condition)
初始时刻 t0速度势 (或 )在流 体域内
或边界上满足的条件。
例5-1 半径为R 的固定大球壳中充满不可压缩理想流体,半径为a
的小球以速度V(t) 在其中运动。试建立速度势定解问题。
解 : 取静坐标系o - xyz
z
2 0 (在流体中)
势流问题的数学描述—— Mathematical Model
5.1.1 基本方程——Laplace Equition
v 0
v
0
v
2 0 (in fluid)
Laplace方程是线性方程。要使 解唯一,需给出边界条件、初
v
p(x, y, z,t)
始条件。
R( M )
5.1.2 边界条件(Boundary Condition)
借助复变函数数学工具解平面势流问题。
平面势流:φ和ψ都是调和函数, 2 0, ,且2满足0
x y
y x
(C-R 条件)
5.2.1 复势与复速度(复平面)
1)复势函数:W (z) (x, y) i (x, y)
解析函数
平面势流
2)复速度(导数)与流体速度的关系:
z x iy
dW W W i i u iv Vei
流体力学实验_第五章
§5.4 流动显示的光学方法
1. 适用范围 光学显示方法:利用流场的光学性质,如流体的密 度变化会造成光学折射率或传播速度的变化,通过 适当的光学装置可以显示流体的流动特性。
流场的温度、压力、浓度和马赫数等状态参数与密度 有确定的函数关系,而流体的光学折射率是其密度的 函数,因此下列流动可以采用光学流动显示的方法:
分光镜 补偿片
单色 点光 源
全反镜
风洞实验段
屏幕
40
密度均匀:干涉条纹彼此平行 密度不均匀:干涉条纹发生移动或变形,干涉条纹的改变与
流体密度的变化有关
干涉条纹 41
§5.5 流动显示技术的新发展——定量的流 动显示和测量技术
1. 激光诱导荧光(LIF)技术
激光诱导荧光技术:是一种20世纪80年代发展起来的光 致发光流动显示与测量技术,把某些物质(如碘、钠或 荧光染料等)溶解或混合于流体中,这些物质的分子在 特定波长的激光照射下能激发荧光。
照明光源:高亮度的白光碘钨灯
25
26
27
3. 荧光微丝法
采用直径为0.01 ~0.02mm的合成 纤维丝,经柔化 和抗静电处理, 使微丝染上荧光 物质,粘贴于模 型表面。
光源:采用连续 紫外光源
照相:选用合适 的滤光片
Flourescent minitufts on aircraft wing
在定常流动中,流线、迹线和染色线相同。
但在非定常流动中,是互不相同的。
4
3. 流动显示方法的分类
(1)示踪粒子流动显示:在透明无色的气流或水流中加
入一些可见的粒子,通过可见的外加粒子跟随流体微团的运 动来使各种流动现象显示出来。 固态示踪粒子:
水流(铝粉、有机玻璃粉末或聚苯乙烯小球等) 气流(烟颗粒) 液态示踪粒子:水流(牛奶、染料溶液) 气态示踪粒子:水流(氢气泡、空气泡)
高等流体力学
概念第一章绪论连续介质:但流体力学研究的是流体的宏观运动,不以分子作为流动的基本单元,而是以流体质点为基本单元,把流场看做是由无数流体质点组成的连续体。
流体质点:流场中一个体积很小并可以忽略其几何尺寸,但与分子相比,这个体积可容纳足够多的分子数目的流体元,有一个稳定的平均特性,即满足大数定律理想流体:忽略流体黏性的流体,即μ=0.可压缩流体与不可压缩流体:简单地讲,密度为常数的流体为不可压缩流体,如水、石油及低速流动的气体。
反之,密度不为常数的流体为可压缩流体。
牛顿流体与非牛顿流体:根据流体流动时切应力与流速梯度之间的关系,即牛顿内摩擦定律。
凡是符合牛顿内摩擦定律的成为牛顿流体,如水、空气、石油等。
否则为非牛顿流体,如污泥、泥石流、生物流体、高分子溶液等动力粘度与运动粘度:动力粘度又成为动力黏度系数,动力黏度是流体固有的属性。
运动粘度又称为运动粘性系数,运动黏性系数则取决于流体的运动状态体积力与表面力:体积力亦称质量力,是一种非接触力,即外立场对流体的作用,且外立场作用于流体每一质点上,如重力、惯性力、离心力。
表面力是一种表面接触力,指流体与流体之间或流体与物体之间的相互作用,主要指压力、切应力、阻力等定常流与非定常流:又称恒定流与非恒定流。
若流场中流体质点的所有运动要素均不随时间变化,则这种流动称为定常流;反之只要有一个运动要素随时间变化则为非定常流大气层分为5层:对流层、同温层、中间层、电离层及外逸层第二章流体运动学描述流体质点的位置、速度及加速度的两种方法,即拉格朗日法和欧拉法质点导数:亦称随体导数,表示流体质点的物理量对时间的变化率,亦即跟随流体质点求导数那布拉P9流体质点的运动轨迹称为迹线流线:此曲线上任一点的切线方向就是该点流速方向依照一定次序经过流场中某一固定点的各个质点连线称为脉线,也叫序线。
流体线:在流场中任意指定的一段线,该段线在运动过程中始终保持由原来那些规定的质点所组成。
流体力学第5章 平面势流理论
2π r 2π r
M 1 c o s i s i n M 1 ( c o s i s i n ) ( c o s i s i n )
2 π r
2 π r c o s i s i n
M 1 2π z
工程流体力学
若偶极子放置在 z z0 处,且偶极子中源到汇的方向 同 x 轴,则复势
当均流叠加偶极子组合,会有圆柱流线形成。它们 组合流场的复势为
工程流体力学
W (z)
W1 (z) W2 (z) U 0 z
M 2p
1 z
(M
0)
对于这个组合流场,只要选择适当的偶极子强度 M
和均流速度 U
的大小,使一条零流线与圆柱表面 (r
0
a)
正好重合即可。
首先引入 z rei,得
式中 z xiy, i 1
解析复变函数称为流动的复势。平面势流必然对 应一个确定的复势W(z),而一个复势也代表一种平面 势流。
工程流体力学
5.1.2 几种简单的平面势流复势
1.均匀直线流动(均流)
当流动速度为U 0 ,方向同x轴方向一致时,复势
W ( z ) U 0 x iU 0 y Biblioteka 0 ( x iy ) U 0 z
流线族
U0(1ar22 )rsin
y
U0(1ar22)rsinC
U0
x
如图5.8所示。
图5.8 均流叠加偶极流场
工程流体力学
W(z)
U0z
a2 z
(1)流场的速度分布:
vr r U0(1ar22)cos v rU0(1a r2 2)sin
第5章 理想流体运动
第五章理想流体流动•欧拉运动方程•伯努利方程及其应用•开尔文涡线定理•能量守恒定律•速度势函数与流函数什么是理想流体?为什么要研究理想流体?第一节理想流体的欧拉运动方程式完整的求解一个流动问题有几个未知数?:p压力u:r速度zy x u u :u ,,速度完整的描述此流动问题需要有几个方程?:=∂∂+∂∂+∂∂z u y u x u zy x 质量守恒方程动量方程个分量有矢量方程3,欧拉运动方程柯西方程()()()()T div g v v t v dt v d ρ1+=∇⋅+∂∂=v v v vv ⎟⎟⎠⎞⎜⎜⎝⎛∂∂+∂∂+∂∂+=∂∂+∂∂+∂∂+∂∂z y x f z u u y u u x u u tu zx yx xx x x z x y x x xτττρ1⎟⎟⎠⎞⎜⎜⎝⎛∂∂+∂∂+∂∂+=∂∂+∂∂+∂∂+∂∂z y xf z u u y u u x u u t u zy yy xy y yz yy yx yτττρ1⎟⎟⎠⎞⎜⎜⎝⎛∂∂+∂∂+∂∂+=∂∂+∂∂+∂∂+∂∂z y xf z u u y u u x u u t u zz yz xz z z z z y z x z τττρ1矢量形式()()()p grad g v v tv ρ1−=∇⋅+∂∂v v v v⎟⎠⎞⎜⎝⎛∂∂−=∂∂+∂∂+∂∂+∂∂x p f z u u y u u x u u t u x x z x y x x x ρ1⎟⎟⎠⎞⎜⎜⎝⎛∂∂−=∂∂+∂∂+∂∂+∂∂y p f zu u y u u x u u t u y yz y y y x yρ1⎟⎠⎞⎜⎝⎛∂∂−=∂∂+∂∂+∂∂+∂∂z p f z u u y u u x u u t u z z z z y z xz ρ1矢量形式剪应力全部=0压应力=压强即正应力=-p根据牛顿第二定律得x 方向的运动方程式为()dt du dxdydzdydz x p p dydz p dxdydz X x ρρ=⎟⎠⎞⎜⎝⎛∂∂+−+上式简化后得同理zoyx微元六面体A A1A2dx xPp ∂∂−21dxxP p ∂∂+21pdtdu x p X x=∂∂−ρ1dtdu z p Z dt du y p Y zy =∂∂−=∂∂−ρρ11111xy z du p X x dt du p Y y dt du p Z z dtρρρ∂−=∂∂−=∂∂−=∂对静止流体的欧拉平衡方程式和理想流体的欧拉运动方程式进行对比101010p X x p Y y p Z zρρρ∂−=∂∂−=∂∂−=∂把上式的三个方程依次乘以i、j、k后相加可得理想流体运动方程的矢量形式,即:1d p dt ρ=uf -∇(,,)d dx dy dz dt dt dt dt==r u dz dtdu dy dt du dx dt du dz zpdy y p dx x p Zdz Ydy Xdx z y x++=∂∂+∂∂+∂∂−++)(1)(ρ由于稳定流时流线与迹线重合,质点沿流线运动,由流线上微元矢量(dx,dy,dz)与时间间隔dt所构成的导数便是流体质点的速度,即将欧拉拉运动微分方程式中各式分别乘以dzdy dx ,,相加得(4-4)伯努利方程的推导——分量方法式(4-4)等号右端可变为222211()()22y x z x x y y z z x y z du du du dx dy dz u du u du u du d u u u d u dt dt dt++=++=++=因此)(21)()(1)(2u d dp Zdz Ydy Xdx dz z p dy y p dx x pZdz Ydy Xdx =−++=∂∂+∂∂+∂∂−++ρρ1()()y x z du du du p p pXdx Ydy Zdz dx dy dz dx dy dzx y z dt dt dt ρ∂∂∂++−++=++∂∂∂•思考一下什么情况下左端的项可以消去?–静止流体–稳定流,且沿流线积分–稳定流,且沿涡线积分–稳定流,且为无旋流动•右端三项分别为:重力势能,动能和压力能•可以写成水头的形式,即单位重量流体的能量•利用伯努利方程,如何通过测压力来测量流速?CvpU E =++=22ρ伯努利方程的适用条件第三节开尔文涡线定理•开尔文涡线定理的表述–理想正压流体在有势力场中运动时,连续流场内沿封闭流体线的速度环量不随时间变化–如果理想流体初始状态静止或绕任意封闭流体线的速度环量为0,则流体运动必然是无旋运动–如果理想正压流体在势力场中运动时,如某一时刻无旋,则流场始终无旋。
流体力学-第5章
六. 伯努利方程 的应用举例
%%%%%%%%%%%%
恒定总流伯努利方程表明三种机械能相互 转化和总机械能守恒的规律,由此可根据具 体流动的边界条件求解实际总流问题。
1
%%%%%%%%%%%%
先看一个跌水的例子。取 顶上水深处为 1-1 断面,平 均流速为 v1,取水流跌落高 度处为断面 2-2 ,平均流速 为 v2,认为该两断面均取在 渐变流段中。基准面通过断 面 2-2 的中心点。
Gz dQdt( z2 z1 )
2 2 1 1 u u 2 2 m2u2 m1u1 ( 2 1 ) dQdt 2 2 2 2
外力对系统做功=系统机械能量的增加
2 2 u2 u1 ( p1 p2 )dQdt dQdt( z2 z1 ) ( ) dQdt 2 2
实际流体恒定总流 的伯努利方程
断面 A1 是上游断面,断面 A2 是 下游断面,hl 1-2 为总流在断面 A1 和 A2 之间平均每单位重量流体所损耗 的机械能,称为水头损失。水头损 失如何确定,将在后面叙述。
分析流体力学问 题最常用也是最 重要的方程式
二、恒定总流伯努利方程的几何表示——水头线
u p2 u z1 z2 2g 2g
p1
2 1
2 2
(P57 3-39)
单位重量理想 流体沿元流的 能量方程式
能量方程
•能量方程的
物理意义
z
u2 z Cl 2g p
伯努利方程表示能 量的平衡关系。
单位重量流体所具有的位置 势能(简称单位位置势能) **************** p 单位重量流体所具有的压强 势能(简称单位压强势能) **************** 单位重量流体所具 p z 有的总势能(简称 单位总势能)
流体力学第五章:旋涡理论
Bˊ Aˊ BA
L
E
C
即 Γc=ΓL (与积分路径方向一致时)
§5-2 汤姆逊定理
假设:
(1)理想流体;
(2)质量力有势; (3)正压流体(流体密度仅为压力的数)
汤姆逊定理: 沿流体质点组成的任一封闭流体 周线的速度环量不随时间而变.
即 d 0
dt
(5-14)
证明:
C上微分长 ds 经dt时间后移到C′,移动速度 v '
θ: 是ds与r的夹角
dH的方向: 垂直于ds和r所在的平面,按右手法则确定。
流体力学中毕奥——沙伐尔公式的形式 旋涡强度为J(环量Γ=2J)的ds段涡丝
对于P点所产生的诱导速度:
流场中单一有限长涡丝在P点的诱导速度沿 整个涡丝积分:
该式可算出任意单一涡丝所引起的诱导速度场
流场中多条涡丝可组成一涡面, 每条 涡丝的诱导速度求得后,沿涡面积分就可 求得整个涡面上的诱导速度。流体力学中 速度场可以看成是涡丝诱导出来的。
旋涡场的几个基本概念:
一、涡线,涡管,旋涡强度 涡线(vortex line):
3 2
涡线上所有流体质点在
同瞬时的旋转角速度矢量
与此线相切。
涡线微分方程:
1 ds
取涡线上一段微弧长 ds dxi dyj dzk
该处的旋转角速度 xi y j zk
由涡线的定义(涡矢量与涡线相切),得 涡线微分方程式:
c
dt
c
c2
而积分式
d d dt dt
c
vds
c
dv dt
ds
c
v
d dt
ds
由欧拉方程
第一项积分可写成
c
流体力学-05 微分形式基本方程
控制面净的质量流率
化工流体力学 大连理工大学流体与粉体工程研究设计所 刘志军 2010年9月~10月
控制体内的质量变化率为
质量守恒的微分表达式
化工流体力学 大连理工大学流体与粉体工程研究设计所 刘志军 2010年9月~10月
质量守恒定律表示为
化工流体力学 大连理工大学流体与粉体工程研究设计所 刘志军 2010年9月~10月
为常数的线在该点处斜率的负倒数
Ψ 为常数的线与 Φ 为常数的线是正交的。 为常数的线是正交的
化工流体力学 大连理工大学流体与粉体工程研究设计所 刘志军 2010年9月~10月
5-6.3 6 3 无旋流动和粘度
速度势只在无旋流动中存在,流函数满足连 续性方程,流函数不限定在无旋流动中。 在什么条件下能够形成无旋流动? 开始不旋转的质点在没有角变形时将不会 发展成旋转运动。
化工流体力学 大连理工大学流体与粉体工程研究设计所 刘志军 2010年9月~10月
Surface Left ( x) (Right (+x) = = = = = = = =
化工流体力学 大连理工大学流体与粉体工程研究设计所 刘志军 2010年9月~10月
Bottom = ( y) (Top (+y) Back ( z) (Front (+z) = = =
每一个面的性能参数通过相对于O点的 泰勒级数展开
化工流体力学 大连理工大学流体与粉体工程研究设计所 刘志军 2010年9月~10月
Surface Inside = (-r) ( Right = (+r) Front = ( θ) (Back = (+θ ) Bottom = ( z) (Top = (+z)
上海交通大学精品课程流体力学课件精品文档
k dV 1 V dp
( m2/N )
式中:dV —— 流体体积相对于V 的增量;
V —— 压强变化前(为 p 时)的流体体积;
dp —— 压强相对于p 的增量。
体积(弹性)模量:
K 1 Vdp k dV
K 不易压缩。
( N/m2 )
一般认为:液体是不可压缩的(在 p、T、v 变 化不大的“静态”情况下)。
五、质量力的势函数
对于不可压缩流体, =常数。 令p/ = w,因 p = p ( x, y, z ),则: w = w ( x, y, z )
由综合式有:
d (p/) = fxdx + fydy + fzdz = dw = (w/x)dx + (w/y)dy + (w/z)dz
得: fxd xfyd yfzdz x pd x p yd y p zdz
静压强的全微分
则: d p f x d f x y d f y z dz
—— 欧拉平衡微分方程的综合表达式
此式便于积分。对于各种不同质量力作用下流体 内的压强分布规律,均可由它积分得到。
3、流体的压缩性,一般可用体积压缩系数 k 和体积模量 K 来描述。 在压强变化不大时,液体可视为不可压缩流体。
4、粘性是流体最重要的物理性质。它是流体运动时产生内摩擦力,
抵抗剪切变形的一种性质。不同流体粘性的大小用动力粘度 或 运动粘度 来反映。温度是影响粘度的主要因素,随着温度升高, 液体的粘度下降。理想流体是忽略粘性的假想流体。
§1-1 流体力学研究的内容和方法
流体力学研究的主要内容:
1、建立描述流体平衡和运动规律的基本方程; 2、确定流体流经各种通道时速度、压强的分布 规律; 3、探求流体运动中的能量转换及各种能量损失 的计算方法; 4、解决流体与限制其流动的固体壁面间的相互 作用力。
流体力学的路线图
流体力学的路线图(之一)流体力学基础理论的学习历来被初学者视为畏途,每到学习结束要进入期末考试的时候,老师和学生一样心中难免忐忑,在流体力学这门课上挂科已经成为某种常态。
即使是学习多年的老手也会在具体问题面前感到基础尚不完备,还不够扎实。
这个问题的起源当然与流体运动规律本身的复杂性有关,这个复杂性导致流体力学与大家印象中的“学科”概念有一定的出入。
比如我们在学习高等数学时,很容易发现,数学是一门“咬文嚼字”的学科,里面充满严格定义的概念,不论学习线性代数还是微积分,都是从一些基本公理出发,循着一条严格的逻辑路线,架构起整门课程。
因为数学有这样逻辑严密的特点,所以虽然学起来也不容易,但大家一致认为数学是美的,而且不论谁写的数学书,比如微积分的书,内容都只有程度深浅的差异,而绝没有内容上的巨大差异。
流体力学则有所不同,流体的流动本身是一种连续不断的变形过程,经典的流体力学理论以连续介质假设为基础,将整个流体看作连续介质,同时将其运动看作连续运动。
但是由于流体是复杂的,实际上至今还没有完全掌握其全貌,因此流体力学在建立了基本控制方程后,就开始转而从一些特殊的流动出发,采用根据流动特点进行简化的方式,先建立物理模型,再得到数学模型,进而得到我们在书中经常看到的很多“理论”,比如不可压无旋流、旋涡动力学、水波动力学、气体动力学等等,甚至理论中还包括理论,比如不可压无旋流中还有自由流线理论,等等。
形成一个类似于俄罗斯套娃的学科结构,这种结构容易给人一种支离破碎的印象。
特别是在各个理论之间联系比较薄弱的时候,更容易给人这种印象。
似乎一门课中又包含了很多门“小课”,每门“小课”使用的数学工具也完全不同,甚至很多同行还进一步把自己分成是学气的(比如空气动力学),或者是学水的(比如学船舶的)等等。
就象旅行者要有一张地图才能更高效率地到达目的地一样,如果能有一张流体力学的地图,或者叫路线图(roadmap),应该对初学者有很大帮助。
流体力学-流体力学绪论
外力作用时,变形的方式不同而已。 流体和固体它们有着完全确定而又互不相同的力学特性
流体不是粘弹体、粘塑体
实际上世界上还存在这样的物体,既不能用 流体,也不能用固体来表述它们的力学性质。
如,有些树胶和油漆,长期静置后,会呈现 固体的力学性态;如果加以摇晃或搅拌,就又显 示出流动性。
实验流体力学
采用试验的方法研究解决工 程实际问题,验证理论推导
§1 流体力学的任务
5)流体力学涉及的基础知识面
流体力学涉及的基础知识有高等数学、理论力学、工程热力学、传 热学等。
6)流体力学的学习方法
流体力学作为一门技术基础课,学习时必须着重于掌握基本概念、 基本理论、基本计算方法和实验技能,通过作业了解流体力学的基本理 论在工程实际中的应用。在学习中培养自己分析问题和解决问题的能力。
信息与控制
生物领域
微流体力学是一门新兴学科在生产小到纳米级的微电源,微机械,微传感器,要研究 在这些微系统中的流速、压力和温度的分布,它们是与宏观流体力学有区别的,压考 虑稀薄效应、尺度效应、表面张力的作用。电流变液是一种两相悬浮液,在电场中会 固化,离开电场又会恢复流体特性。
在细胞层次上进行研究人体的生理流动,与工程紧密结合是生物流体力学今后发展趋 势上研究,流体的非牛顿性, 脏器的多孔性,微循环细胞膜的传质,流动的尺度现象, 动物昆虫的仿生流体力学的研究也是发展的方向
§2 流体力学发展简史
流体力学是力学的重要分支之一,它与其它学科一样与人类 的生产活动有着紧密地关系。迄今为止流体力学流体力学经历了 四个重要的发展时期:
1) 古代流体力学发展时期 2) 经典流体力学发展时期 3) 近代流体力学发展时期 4) 现代流体力学发展时期
流体力学第五章(理想不可压缩流体的平面势流)
流体力学——理想不可压缩流体的平面势流内容¾基本方程组,初始条件及边界条件¾速度势函数及无旋运动的性质¾平面流动及其流函¾不可压缩流体平面无旋流动的复变函数表示¾基本的平面有势流动¾有势流动叠加P=Pa , Pa为大气压强。
在直角坐标系中有一个线性的二阶偏微分方程(拉普拉斯方程线性方程的一个优点是解的可叠加性对于定常流:则由伯努利方程得到理想不可压缩无旋流的基本方程为:边界条件静止固壁上自由面上:P = Pa 无穷远处:速度势函数及无旋运动的性质在无旋流中有若已知函数,则可求出若已知速度矢量V,则可由积分求出势函数上式中为任意常数,因此的值相对于不同的Mo点可以差一个,为某一常数,但并不影响流动的实质,因为当求流动的特征量ui, P时,常数的差别便消失不见了,所谓的结果完全一样φ涉及到单值和多值问题在单连通区域 与积分路线无关,而只与起点M0及终点M的位置 有关。
因而势函数为单值函数。
在多连通区域 , 是封闭曲线L绕某一点的圈数, 称为环量 势函数 为多值函数。
速度势函数及无旋运动的性质(已作介绍)内容 ¾ 基本方程组,初始条件及边界条件 ¾ 速度势函数及无旋运动的性质¾ ¾平面流动及其流函数 不可压缩流体平面无旋流动的复变函数表示 基本的平面有势流动 有势流动叠加¾ ¾平面流动及其流函数 平面问题是指 流动在平面内进行,即 u z = 0 ; 垂直平面的垂线上个物理量相 等即适用范围 无限长柱体,它的一个方向的尺寸比其它两个方向的尺寸大得 多,在长方向的速度分量很小,其它物理量的变化也很小。
如:低速机翼表面的压力分布问题的理论计算等,无限长的柱 体平板的绕流等研究平面无旋运动,在平面运动中,涡旋矢量Ω的三个分量为只有 而无旋,可推出存在着速度势函数 使得:速度势函数的性质我们已经讨论过了流函数的意义 如果能够找到某一函数Ψ,满足流动的可能判据 —— 连续性 方程,则称这一函数Ψ为流函数 在平面运动时,不可压缩流体的连续性方程为:若有一函数Ψ(x,y,t)并令 则连续性方程为称为流函数知道了流函数 •若与流速ux ,uy 之间的关系之后 求出流速场已知,可由• 若 ux ,uy 已知,可用积分速度势与流函数 平面流动垂直与z轴的每个平面流动 都相同,称平面流动速度势函数 速度势函数存在的条件∂w ∂v − = 0 ∂y ∂z ∂u ∂w − = 0 ∂z ∂x ∂v ∂u − = 0 ∂x ∂y此条件称 柯西—黎曼条件由高数知识可知,柯西—黎曼条件是使udx + vdy + wdz全微分的充要条件,即成为某一个函数ϕ(x ,y ,z ,t )d ϕ = udx + vdy + wdz而当 t 为参变量, ϕ(x ,y ,z ) 的全微分为∂ϕ ∂ϕ ∂ϕ dϕ = dx + dy + dz ∂x ∂y ∂z比较两式有∂ϕ u = ∂x ∂ϕ v = ∂y ∂ϕ w = ∂z∂ϕ 柱坐标 V r = ∂r 1 ∂ϕ Vθ = r ∂θ ∂ϕ Vz = ∂z把ϕ(x ,y ,z ) 称为速度势函数简称势函数无论流体是否可压缩,是否定常流只要满足无旋条件 ,总有 势函数存在。
工程流体力学-第五章 粘流和边界层流动_完整版
§5.1 粘流的基本特性:(五)流动分离 • 现在,来考虑大雷诺数情况下真实流体绕二维
翼型的流动。 • 如果来流攻角(无穷远处速度与翼弦的夹角)
不大(比如小于 10)0 ,流体平滑地绕翼型流 动而不发生明显的边界层分离。这时,真实流 体效应(粘性)只在紧靠翼型流动而不发生明 显的边界层的。
26 14 作业13 作业12 作业11 作业10
§5.1 粘流的基本特性:(五)流动分离
• 由于在上述情况下,边界层和尾迹都是厚度极 小的薄层,绕翼型的流场(在边界层和尾迹之 外)基本上与理想流体绕同一翼型的流动相 同;
• 翼型表面压力分布和翼型升力系数的实测值与 理想流体位势理论所得结果非常接近。在这种 情况下,翼型所受的阻力主要是摩擦阻力。阻 力的实测值虽不为零,但阻力与升力的比值颇 小。
§5.2.1 附面层的概念
35 14 作业13 作业12 作业11 作业10
§5.2.2 附面层的厚度 ➢ 附面层厚度δ
平板边界层流动示意图
36 14 作业13 作业12 作业11 作业10
§5.2.2 附面层的厚度 • 物面处流体速度 u ,0 物面上方 沿u 方向y 递
18 14 作业13 作业12 作业11 作业10
§5.1 粘流的基本特性:(五)流动分离
v
层紊流1流分 度30分离左度离发右左发生。右生在。在85
19 14 作业13 作业12 作业11 作业10
§5.1 粘流的基本特性:(五)流动分离 • 下面,考虑真实流体绕圆柱的流动。
✓ 根据实验观察,可以发现,在不同的雷诺数范畴,有 完全不同的流动形态。
• 最初是由普朗特在1904年提出来的。当流体流过物体 时,由于物体表面和流体之间的摩擦力对流体产生迟 滞作用,物面上流动速度为0。
流体力学-有旋流动和有势流动
5
涡管 在流场中,取一条
不与涡线重合的封闭曲 线 L,在同一时刻过 L上 每一点作涡线,由这些 涡线围成的管状曲面称 为涡管。
涡管 涡线 Ω
与涡线一样,涡 管是瞬时概念
6
涡通量 通过流场中某曲面 A 的涡量通量 Ω n d A
称为涡通量。
A
n
涡管强度 通过涡管任一截面 A 的涡 通量又可称为涡管强度
这个分类是 很重要的
无旋流动
有旋流动
判别的唯一标准是看流速场的旋度是否为零
3
• 涡量、涡线、涡管和涡通量
涡量 对于有旋流动,将流速场的旋度 称为涡量,它是流体微团旋转角速 度矢量的两倍。涡量场是矢量场。
Ω u 2ω
涡线 涡线是涡量场的矢量线,是某瞬时对应的流场中的曲线,
该瞬时位于涡线上各点对应的涡量都沿着涡线的切向。与流线 一样,涡线是与欧拉观点相对应的概念。
A
dA
Ω
I Ωn d A ( u) n d A 2ωn d A
A
A
A
留下一个问题:为 什么可取任一截面
计算涡管强度
AΩ
7
• 速度环量、斯托克斯定理
速度环量 定义流速矢量 u 沿有向曲线 L 的线积分为速度环量
Γ udl
L
斯托克斯定理 n
Ωn d A u d l
A
L
dA
Ω
封闭曲线 L 是 A 的周界,
4
涡线微分方程 根据定义,涡线的微分方程为
Ω d l 0 其中 dl d xi d yj d zk
i jk dx dy dz 0 x y z
dx dy dz x (x, y, z,t) y (x, y, z,t) z (x, y, z,t)
势流理论
2.经过B,D后又逐渐减小,在C点汇合时速度
又降至零。离开C点后,又逐渐加速,流向后方
的无限远处时再恢复为v0。
柱面上的压力分布:
定常,不计质量力的拉格朗日积分式为:
p
V 2
2
p0
V02
2
无穷远均匀流中压力
将(6-22)式代入即得圆柱表面上压力分布:
p
p0
V02
2
(1 4sin2 )
这些圆与ψ =const正交 注意:
偶极子的轴线和方向
轴线:源和汇所在的直线
方向:由汇指向源的方向
图6-8(b)
偶极子的方向
为x轴负向
四、点涡(环流)
点涡:无界流场中坐标原点处一无穷长直线涡,
方向垂直于x0y平面,与xoy平面的交点 诱导速度沿点涡为中心的圆周切线方向,大小
与半径成反比:
vs 2 r
Q
2 r
dr
d
(
r
dr
d )
VsdrrVrd
Q
2
d
Q 2
ln
r
Q 2
(6-8)
流线为θ =const,为原点 引出的 一组射线
等势线为r=const,流
线为同心圆,相互正交。
图6-6
当Q>0,则 Vr>0为点源,反之为点汇。 对于扩大(收缩)流道中理想流体的流动,
利用泰劳展开: ln(1 z) z z2 z3
23
令 z x cos1
r2
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
W ( z)
M 2z
位于(0,0)偶极:
M x M cos 2 x 2 y 2 2 r
M y M sin 2 x 2 y 2 2 r
W ( z)
M 2z
5.3.4 平面点涡 (vortex)
位于(0,0)点涡:
2 ln r 2
5.4.3
沿轴正向均匀流与偶极的叠加
偶极位于(0,0),方向沿 - 轴:
M cos 2 r M sin V0 r sin 2 r
V0 r cos
W ( z ) V0 z
M 2z
速度分布: vr V0 cos (1
仍然是解析函数,仍然代表某一种流动的复势。简单 流动组合成复杂流动——叠加法
5.3 平面势流的基本解
目的:求解最简单的流动,为解决复杂势流奠定基础。 内容:均匀流、点源、点涡、偶极。 方法:利用已知流动的特征,“凑”。
V0
m
M
5.3.1 均匀流 (uniform stream)
≠ 0 时:
m xb x b u V0 2 ( x b) 2 y 2 ( x b) 2 y 2 m y m y v 2 2 2 2 2 ( x b) y 2 ( x b) y
5.4.2 均匀流和一对等强度源汇的叠加
驻点位置:
v d r vr r d
W ( z)
y
m 2
r
x
m m m ln r i ln z 2 2 2
ψ=const φ=const
5.3.2 平面点源、点汇 (source and sink)
点源位于(x0,y0):
m ln ( x x0 ) 2 ( y y 0 ) 2 2
x0 0
M 1 y y0 M W ( z) 2 2 2 z z 0 2 ( x x0 ) ( y y0 )
5.3.3 平面偶极 (dipole)
位于(0,0)偶极:
M x M cos 2 x 2 y 2 2 r
M y M sin 2 x 2 y 2 2 r
z
解 : 取静坐标系o - xyz
2 0 (在流体中)
内边界条件:小球表面方程为
o
y
x
V(t) x0
F ( x x0 ) 2 y 2 z 2 a 2
( x x0 )V (t ) ( x x0 ) y z 0 x y z
外边界条件:大球面方程为 F x 2 y 2 z 2 R 2 ,得
F F 0 t
若物面静止不动: v b 0 ,则物面边界条件简化为
0 n
F F F 0 ( on S ) x x y y z z
2) 无穷远边界条件
(1)大地坐标系:
v0
( R )
0
(2)随体坐标系:若物体以V0 运动,则问题转化为物体不动,
势流问题的数学描述—— Mathematical Model 5.1.1 基本方程——Laplace Equition
v 0 v 0 v
2 0 (in fluid)
Laplace方程是线性方程。要使 v 解唯一,需给出边界条件、初 始条件。
p( x, y, z, t )
2 2 2
vb
— 物面运动速度
v
— 流体质点的速度
n — 物面的单位外法向量
v n vb n
(on S)
n
F F
vb n n
5.1.1 基本方程——Laplace Equition
v n vb n
F v b F (on S)
W i ( i ln r ) ln z 2 2
v
Γ顺时针方向,若逆时针,上式加负号。
不难验证,上述基本解满足 Laplace 方程和相应的无穷远条件的。 另外,在源、涡和偶极的位置上存在奇异性(奇点)。可见,点源 (汇)、偶极以及点涡都是奇点,均匀流是一个特殊的奇点。
作业
5-4
W ( z ) (1 i ) z
补充题:已知复势为:
z 1) w ( z ) (1 i ) ln z4
试分析以上流动的组成,绘出流线图.
5.4 平面势流基本解的叠加
1. 由基本解构造复杂流动的解 —— 基本解(奇点)叠加法。 2. 基本解叠加 代表何种物理流动?
y
V V0e
i
u V0 cos , v V0 sin
o
平板
V0
u d x v d y V0 x cos V0 y sin
v d x u d y V0 x sin V0 y cos
x
W ( z) V0 z cos iV0 z sin V0 ze i
y x
V0
m
2 a
均匀流和源叠加可模拟绕弹形物体的流动。调整源强m和速度V0, 改变流线的形状。
5.4.1 均匀流和点源的叠加 流场中压力分布
p ( p0 1 2 1 v0 ) V 2 2 2
压力系数
V 2 p pБайду номын сангаас cp 1 v2 1 2 0 v0 2
位于(x0,y0),沿 -x 轴方向:点源
( x0 , y0 ) ,点汇 ( x0 x0 , y0 )
m m 2 2 ln{( x x0 ) ( y y0 ) } ln{[ x ( x0 x0 )]2 ( y y0 ) 2} 4 4
x x0 M 2 ( x x0 ) 2 ( y y0 ) 2
y y0 m arcg 2 x x0
W ( z)
m ln( z z 0 ) 2
5.3.3 平面偶极 (dipole)
偶极强度:设强度为m 的源和汇相距 x0
x0 0
lim m x0 M
x
+m
-m
这对源汇构成一新的奇点为偶极,方向由汇指向源。 偶极既有大小,又有方向。
若物面运动:对 F ( x, y, z, t ) 0 求全(物质)导数
dF F F dx F dy F dz F v b F 0 dt t x dt y dt z dt t
F F F F 0 t x x y y z z
( R )
v
p
F
奇点叠加法;保角变换法(平面流)。 数值解:复杂边界问题。
CFD — Computational Fluid Dynamics
5.2 复势(complex potential )
借助复变函数数学工具解平面势流问题。
平面势流:φ和ψ都是调和函数, 2 0, 2 0,且满足
驻点位置
m x a, y 0 或 2V0
r a
5.4.1 均匀流和点源的叠加 过驻点(a, 0)流线方程:
v0 y
m C v 0 a 2
y a( )
y a
。
r 时, 0 ,流线Ⅰ在无穷远处的半宽为
sin 2 sin surface
2
5.4.2 均匀流和一对等强度源汇的叠加
x方向均匀流
+
等强度源汇:源(-b,0)、汇(b,0)
2b
+
V0
m
m
+
V0 x
m m ln ( x b) 2 y 2 ln ( x b) 2 y 2 2 2
流体力学 第5章 势流理论 (Chapter 5. Potential Flow Theory)
本章内容: 研究不可压理想流体无旋运动流场的 速度分布、压力分布及作用于物体上的力。
Background:
Aviation, ship & ocean eng. water waves.
5.1 势流问题的基本方程和边界条件
dW u iv dz
z x iy
dW u2 v2 v V dz
5.2.1 复势与复速度(复平面)
3)复速度的环路积分与速度环量和流量的关系:
l
dw dz dw d id l iQl l l dz
dw l Re l dz
x m 2b b2 , 2 v0 y0
过驻点流线: m 1 2by vo y tg ( 2 )0 2 2 2 x y b
V0
o
+
y
x
2v0 y 2by tg ( ) 2 m x y 2 b2
点源推开流线,点汇收回流线。
将流线替换成物面,该解模拟流体绕卵形体的外部流动。
5.4.1 均匀流和点源的叠加
m V0 r cos ln r 2 W ( z ) V0 z m ln z 2
+
dW m x m y V0 i dz 2 x 2 y 2 2 x 2 y 2
5.4.1 均匀流和点源的叠加 流体速度:
u V0
m x m cos V 0 2 x 2 y 2 2 r m y m sin v 2 2 2 x y 2 r
而流体从无穷远处以-V0 流来 —— 绕流问题。
v 0
( R )
5.1.3 初始条件(initial condition)
初始时刻 t 0 速度势 (或 )在流 体域内或边界上满足的条件。