高等流体力学-第5讲

合集下载

高等流体力学课件

高等流体力学课件
静止流体满足力的平衡条件,即合力为零。
流体静力学的基本概念
流体静力学是研究流体平衡和压力分布的学 科。
压力分布
静止流体的压力分布与重力场和其他外力场 有关,可以通过静力学方程求解。
流体动力学
总结词
流体动力学的基本概念、一维流动、层流与湍流
一维流动
一维流动是指流体沿着一条线的流动,可以用于 描述长距离管道内的流动或某些对称的流动。
水利工程
机械工程
流体动力学在水力发电、水利枢纽设计、 灌溉系统优化等方面具有广泛应用,为水 利工程提供了重要的技术支持。
流体动力学在机械工程领域的应用也十分 广泛,如内燃机、通风 system等的设计和 优化。
流体在自然界中的应用
气候变化
流体动力学在气候变化研究中发挥着重要作用,如风场、洋流等 对气候的影响研究。
详细描述
连续性方程是流体动力学的基本方程之一,它表达了单位时间内流经某一封闭 曲面微元体的流体质量的增加等于该微元体所受质量源的净增量,用于描述流 体运动的连续性。
动量方程
总结词
描述流体动量守恒的方程
详细描述
动量方程是流体动力学的基本方程之一,它表达了流体动量的变化率等于作用在 流体上的外力之和,包括重力、压力、摩擦力等。
方法
02
常用的线性稳定性分析方法包括特征值分析、傅里叶分析和庞
加莱截面法等。
应用
03
线性稳定性分析在气象、海洋、航空航天等领域有广泛应用,
用于预测和控制流体运动的稳定性。
非线性稳定性分析
定义
非线性稳定性分析是研究流体运动在较大扰 动下的响应,需要考虑非线性效应对流体运 动的影响。
方法
非线性稳定性分析需要求解非线性偏微分方程,常 用的方法包括数值模拟和近似解析法。

高等流体力学-第五讲.

高等流体力学-第五讲.
Mn x nC( x, t)dx
2)浓度分布各阶矩与统计量间的关系
浓浓度度分分布布方质差量:中心2坐标—1 —均(x值:)
μ
2 C(
M1 / M0
x,t)dx
M 0
取: M0 1; M1 0, 即 0 可证明
M2 /
d
dt
M0 2
2
2Dm
q — 在xi方向的单位面积的扩散质量通量;量纲:ML2T 1
Dm — 分子扩散系数;量纲: L2T 1
对三维情况,以矢量表示: q DmC
费克定律说明:在扩散溶液浓度场中的时空点上,单位时间内通过单位
面积的扩散质的质量与该点处扩散溶液浓度的梯度成正比,比例系数为该种
扩散溶液的分子扩散系数;方向与浓度梯度方向相反。
高等流体(水)力学讲稿
11
第五讲 扩散理论
(2)分子扩散的随机游走分析
确定分子游走的概率值:分子的自由程l,经过N次运动后行走X距离 的概率。设在N次行走中,有p次沿x的正向,q次沿x的反向行走;每次行走 是相互独立的。
有:p + q = N;令: p – q = S ,故有:X = Sl
沿x的正向行走X距离的概率P为:
代入移流扩散方程后取时均值,有:
C t
ui
C xi

xi
Dm
C xi
C ui
其中: Cui 是由紊动脉动量产生的沿i方向的质量扩散通量。
仿照分子扩散系数的表示形式,引入紊动扩散系数Dij,
令:
C ui

Dij
C x j
紊动扩散方程可表示为(考虑源、汇项后):
C
C u1 D11

流体力学热能第5章 不可压缩流体动力学基础讲解

流体力学热能第5章 不可压缩流体动力学基础讲解
第七章 不可压缩流体动力学基础
本章讨论三元流动,主要内容是有关流体运动的基本概念和基本原理, 以及不可压缩流体流动的基本方程。 积分形式的基本方程用于解决控制面上的流动参数问题。 微分方程可用于解决流 动参数在流场中的分布问题。
一、运动形式
§7-1 流体微团运动的分析
1、流体微团:指体积微小,随流体一起运动的一团流体物质。与流体质点不 同,虽体积微小,但包含无数个流体质点。各质点间存在着相对位置的变化。
?x
?
?ux ?x
?y
?
?uy ?y
差值为正,发生伸长变形。
?z
?
?uz ?z
3、旋转角速度
逆时针为正
对角线EMF 的旋转角速度定义为
A
整个流体微团在oxy 平面上的旋转角速度。 E
? ?
?
?
z
?
1
?u (
y
2 ?x
?
?ux ) ?y
?
??
?
y
?
1 (?ux 2 ?z
?
?uz ) ?x
?
??
?
x
?
2、基本运动形式
平移运动
旋转运动
线变形、
变形运动
角变形
BF
二、运动分析
以二元流动的情况为例,研究几种
A
uy
ux C
dy
基本运动形式的速度表达式。
M
E
如图,方形流动微团
D
dx
各侧边中点A、B、C、D的流速分量分别为
M
A
ux
ux
?
?ux ?x
? dx 2
B
ux
?
?ux ?y

高等流体力学第五章(1)

高等流体力学第五章(1)

1 dl 2 Pl ( x) l x 1 l 2 l! dx
其前3项分别是,


P0 ( x) 1
P ( x) x 1
P2 ( x) 1 3x 2 1 2


5.4
势函数
均匀流
沿 x 方向均匀流,速度为 U,P点的势函数 ,
Ux
x r cos

P
U
r

Ur cos
称偶极子的强度。请注意在求上述偶极子势函数过程中,点汇在 x
轴正方向放置,点源在 x 轴负方向放置,相互无限靠近。
x
o


x

Q
5.6
偶极子流动
流函数
依据流函数与势函数之间的关系式求偶极子流的流函数,
1 1 2 cos 2 3 r r sin 2r r sin
er 1 r 2 sin r ur
u 0
re ru
r sin e e ru u r r r 0
e r 1 2 0 r r r sin r r sin

流函数也可利用势函数与流函数关系式求得.
5.6
势函数
偶极子流动
P
求一对相等强度的点源和点汇 在 P 点的势函数,

Q Q 4 r 4 (r r )
r
Q
r
r r
当 r / r 1 时
2 2 Q r r Q r r O O 1 1 4 r r r 4 r r r

【计算流体力学】第5讲-差分方法3

【计算流体力学】第5讲-差分方法3

通量差分分Байду номын сангаас (FDS): 耗散低、分辨率高
Step 1: 运用差分格式,计算
U ,U L j 1/ 2
R j 1/ 2
Step 2: 运用Riemann解, 计算
F j 1/ 2
F
(U
L j 1/
2
,U
) R
j 1/ 2
Step 3: F Fj1/2 Fj1/2
x
x
U ,U L j 1/ 2
R j 1/ 2
f
j 1/
2
x
x
f f + f =+
x x x
13
3. 特征重构方法
常系数方程组:
U t
A U x
0
U t
S1ΛS U x
0
V t
Λ V x
0
vk t
k
vk x
0
变系数情况—— 局部冻结系数
完全 解耦
U f(U) 0 U A U 0
t x
t x
在基架点上系数 A j 不变
U t
u j u j1 u j1 u j
二阶精度区
TVD区
二阶精度TVD区(二 者交集)
1
通量分裂技术: 模型方程 NS/ Euler 方程
Step 1 针对模型方程构造差分格式
u a u 0 t x
u uˆ j1/2 uˆ j1/2
x
x
a0
uˆ j1/2 =......
格式1
a0
uˆ j1/2 =......
6
➢ Steger-Warming 具体步骤 (以一维为例)
u a u 0
t x

高等流体力学习题

高等流体力学习题

第一讲绪论习题:1.综述流体力学研究方法及其优缺点。

2.试证明下列各式:(1)grad(φ±ψ)=grad(φ)±grad(ψ)(2) grad(φψ)=ψgrad(φ)+φgrad(ψ)(3)设r= x i+y j+ z k,则=(4) 设r= x i+y j+ z k,求div(r)=?(5) 设r= x i+y j+ z k,则div(r4r)= ?3.给定平面标量场f及M点处上已知两个方向上的方向导数和,求该点处的grad f 第二讲应力张量及应变张量例2-1试分析下板不动上板做匀速运动的两个无限大平板间的简单剪切流动,,式中k为常数,且k=u0/b。

解:由速度分布和式(2-14、16和17)可得再由式(2-18)可得所以II=k=u0/b。

流动的旋转张量R的分量不全为零说明流动是有旋流动,I=tr A=0表明流动为不可压缩流动,II=k表明了流场的剪切速率为常数。

第三讲流体的微分方程习题:试由纯粘流体的本构方程和柯西方程推导纳维尔-斯托克斯方程(N-S方程)。

第四讲流动的积分方程【例3-1】在均匀来流速度为V的流场中放置一个垂直于来流的圆柱体,经过若干距离后测得的速度分布如图所示,假设图示的控制体边界上的压力是均匀的,设流体为不可压缩的,其密度为ρ,试求:(1)流线1-2的偏移量C的表达式;(2)单位长度圆柱体的受力F的表达式。

解:(1)无圆柱体时流管进出口一样大(即流线都是直线,无偏移),进出口的流速分布也是相同的,而放入圆柱体之后出口处的流速分布变成图示的那样,即靠近中心线部分的流速变小,由于已经假定流体是不可压缩的流体,若想满足进出口流量相同——连续性方程,必然会导致流管边界会向外偏移,也就是说出口处流管的截面会增大。

因此,求解时可由进出口流量相等入手,设入口处平均流速为V,取宽度为L,所得的连续性方程应为:求得C=a/2(2)在流管的进出口截面1-1与2-2之间使用动量方程,即圆柱体的阻力应等于单位时间内流出2-2面的流体的动量与流入1-1面的流体的动量差,列x方向的动量方程可表示为则,F=-R【例3-2】试求如图所示的射流对曲面的作用力。

高等流体力学第5讲

高等流体力学第5讲

f (z) (x, y) i (x, y)
z(x,y) r
式中 i 1 ,

O
x
z x iy rei 为复变量。
x, y是平面直角坐标, r, 是平面极坐标的向径和幅角。
设函数 f (z) 在点 z0的某邻域内处处可导,则称函数 f (z) 在 点 z0 处解析;又若 f(z) 在区域 D内的每一点解析,则称 f (z) 在 区域 D 内是 解析函数。
x y 可得 f f
x (iy)


y



x
三、 复速度
用复势描述势流运动,需要建立复势f (z)和速度矢量 v 的关系。
复势f (z)的导数为
df i u iv
dz x x
称复势的导数为复速度,其实数为x向的分速度, 其虚数为 y 向分速度的负值。
其中记号


2



2 x2
2 2 y2 z2
称为拉普拉斯算子(Laplace operator)
笛卡尔坐标系下 2 2 2
x2 y2 z2 0
对于不可压缩均质流体 = const,有 1 p p


对于外力有势且质量力只有重力的情况,有 f gk (gz)
代入理想流体 Lamb 形式的运动方程
v t


v v
2


v

f

1

p
积分可得

v v gz p f (t)
t 2

Lagrange积分方程
对于理想不可压缩均质流体的势流运动,其控制方程为

中科院计算流体力学最新讲义CFD2011-第5讲-差分方法3

中科院计算流体力学最新讲义CFD2011-第5讲-差分方法3

(1)
2
保证“系数非负”
1 un j
a u
k k
n j k
ak 0

Cn 0, j 1
Dn 0, j 1
2
Cn Dn 1 j 1 j 1
则格式(1)是TVD格式
1 n n n n n n un un j j C j 1 (u j 1 u j ) D j 1 (u j u j 1 )
t
a
n un j 1 u j 1
2x

a 2 t n n (u j 1 2u n j u j 1 ) 2 2x
1 n n n n n n un un j j C j 1 (u j 1 u j ) D j 1 (u j u j 1 )
2 2
6
Copyright by Li Xinliang
2) 重要概念: 网格 Reynolds数 以网格尺度度量的Reynolds数
Re x Rex
j+1
含义: 数值振荡—— 流动尺度为网格尺度 网格 Reynolds数小,该尺度的能量 被耗散掉—— 不发生振荡
网格足够小:不会发生振荡; 网格小于激波的实际厚度,则不会振荡
有效网格点数: 一个波长里 面的网格点数 (PPW: Point per Wavelength)
PPW

x

2 2 k x
GVC2 格式
u j 1 / 2 (3u j u j 1 ) / 2 (u j u j 1 ) / 2 when u j u j 1 u j 1 u j when u j u j 1 u j 1 u j
(b1x) 2 I
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

牛顿平板实验 实验结论: 速度呈线性分布 摩擦力与速度梯度呈正比
U τ =μ h
du τ =μ dy
牛顿流体与非牛顿流体
2.2 流体的连续介质假说
2.2.3 流体的性质及分类
¾ 粘性 流体抵抗剪切变形(或相对运动)的一种属性。 粘性系数 ——理想流体和粘性流体 动力粘性系数μ,单位:Pa.s, U −1 ⎤ ⇒ [ μ ] = [ Pa ⋅ s] ⇒ [ Pa] = μ ⎡ s τ =μ ⎣ ⎦ h 运动粘性系数ν, 单位: m2/s Pa ⋅ s Pa ⋅ s N / m2 ⋅ s kg ⋅ m / s2 ⋅ s 2 = = ν =μ/ρ = = 3 3 = m /s 3 kg / m kg / m kg / m kg / m 20℃时,空气ν =0.151×10-4 m2/s, 水ν = 0.1×10-5 m2/s 粘性系数与温度的关系 液体的粘性主要取决于液体分子间的距离和分子间的吸引 力,所以温度升高液体粘性下降。 气体的粘性主要取决与气体分子热运动所产生的动量交 换,所以温度升高气体粘性增大。
1 2
mv2 ≈ ΔE
液体:分子间形成不确定、不稳定的对偶结构
2.2 流体的连续介质假说
2.2.2 连续介质假设—宏观
¾ 连续介质模型 假设组成流体的最小物质实体是流体 质点,流体由无限多的流体质点连绵 不断地组成,质点之间不存在间隙。 ¾ 流体质点、流体微团 相对于宏观很小, 相对于微观,又包含足够多的分子。 ¾ 物理参数 微团内大量分子所表征物理参数 的平均值。 ¾ 极限过程 基于宏观意义上的连续性和可微性。
ρ=
ΔV → 0 ΔV →ΔV ′
lim
Δm ΔV
2.2 流体的连续介质假说
2.2.2 连续介质假设—宏观
— 连续介质假说有三个方面含义:
(1) 物质体介质在某点的物理参数值是指位于该点 — 宏观微团上大量分子相应物理参量的平均值; (2) 物质体介质的连续性是宏观意义上的连续是指 — 表征大量分子团的平均物理参数在空间上的连续性; (3) 所谓极限过程 — 宏观尺度上的无限逼近。
G G r = r ( a , b, c , t )
¾ 不同的t + 相同 的(a, b, c) —空间曲线 : 质点随时间变化在空间中运动而 留下轨迹线 — 迹线 ¾ 不同的t + 不同的拉格朗日变数(a, b, c) : 不同质点的不同迹线。 ¾ 相同的t + 不同的拉格朗日变数(a, b, c) :某一时刻下不同质点的空间位置。
2.3 描述流体运动的基本观点和方法
2.3.1 拉格朗日方法
于是 速度: 加速度:
G G ∂r (a, b, c, t) V (a, b, c, t) = ∂t 2G • G ∂ r (a, b, c, t) V (a, b, c, t ) = ∂t 2
G G 向径 r = r ( a, b, c, t )
对拉格朗日变数 (a, b, c) 不 一定连续可微
·直角坐标系下 质点向径(坐标)
x = x ( a , b, c , t )
y = y (a , b, c, t )
z = z ( a , b , c ,t )
∂x u = ∂t ∂y v= ∂t ∂z ω = ∂t
质点速度
质点加速度
∂2x u= ∂t 2
2.1.1 流体力学研究对象[52]
固体力学 理论力学、材料力学、弹 塑性力学、断裂力学等
流体,是指在任何 剪切力持续作用下 连续变形的物质。
物质的机械运动以及机械 运动和其他运动的关系
力 学
物 质
固体 流体
流体力学
工程流体力学、空气动力 学等
流体静力学
运动学
动力学
静止流 体
流体的机械运动
各种物理参量与运动 状态之间的相应关系
2.2 流体的连续介质假说
2.2.1 物质分子运动论
— 分子平均动能与结合能的比较 (1) 分子平均动能远大于其结合能:
1 2
mv2 >> ΔE
r0
rc
气体:分子间不形成对偶结构 ; (2) 分子平均动能远小于其结合能:
1 2
mv2 << ΔE
固体:分子间形成稳定的对偶结构; (3) 分子平均动能约等于其结合能:
理论求解 物理 解释
离散化
数值解
有限差分法 有限元素法 边界元法 有限基本解
计算流体力学
风洞 水洞 水槽 水池
实验模型
实验参数的测定
数据处理及解释
2.2 流体的连续介质假说
2.2.1 物质分子运动论[54]
— 物质的组成与形态 物质的状态(气、液、固)取决于组成物质的这些分子间的间 距和它们之间的排列形式。 固体:分子间间距小,分子按空间点阵排列(如晶格),以各自的点 阵位置为中心作无规则的振动,并受到相邻分子的约束。 液体:分子间间距较固体大,没有一定空间阵列,分子间有约束,但 受扰动能摆脱这种约束,同时分子的运动速度也比固体大,以 致于无法维持分子间的固定结构。 气体:分子间完全没有约束,间距很大,能在空间中自由运动。
小 结
¾
流体力学研究的对象:
流体
¾
流体力学的三个分支:
理论流体力学、实验流体力学、计算流体力学
¾
连续介质假说
宏观假设,物理参数是大量分子的平均值
¾
流体的性质及分类:
流动性、粘性、可压缩性
2.3 描述流体运动的基本观点和方法
拉格朗日方法 —— 欧拉方法
2.3.1
拉格朗日方法
拉格朗日变数 — 有序实数三 位(a, b, c); —质点坐标; —不同(a, b, c) 对应不同质点。 质点空间定位
2.1 流体力学研究的对象、方法及其应用
2.1.2 流体力学研究的方法及其应用[52]
介质运 动过程中 的力学特 征和物理 定律
理论流体力学 物理 模型 物理模型的量化描述: ( 运动参数、状态参数以及这 些参数间的定量关系的确定 )
初边条件 实验流体力学
量纲分析 相似原理
数学模型 运动微分方程 离散化微分方程
r0
rc
第二章 流体力学基本概念
2.1 流体力学研Байду номын сангаас的对象、方法及其应用
2.1.1 流体力学研究的对象
2.2 流体的连续介质假说 2.3 描述流体运动的基本观点和方法 2.4 流体微团运动学分析 2.5 流体微团的受力分析 2.6 物质元素及物质积分的随体导数
2.1 流体力学研究的对象、方法及其应用

∂2y v= ∂t 2 • ∂2z ω= 2 ∂t

1 ⎛ dV / m ⎞ 1 ⎛ dκ 1 ⎛ dV ⎞ =− ⎜ ⎜ ⎟ =− ⎜ ⎟ V ⎝ dp ⎠T V / m dp κ dp
体积弹性模量Eν
dp dp =ρ Eκ = = −κ β dκ dρ
1
根据流体是否可压缩,分为可压缩流体和不可压缩流体 尽管所有流体都具有一定的可压缩性,但一般当作不可压缩流体进行处 理,所引起的误差不大。
需要声明:连续介质假设不适用于稀薄气体的情况, 气体分子的自由程与问题所涉及的物理尺度相比较已是 不可忽略的了。
2.2 流体的连续介质假说
2.2.3 流体的性质及分类
¾ 流动性 流体,只能承受压力,不能承受拉力或剪切力 任何剪切力的持续作用都会使流体产生任意大的连续变形。 ¾ 粘性 流体抵抗剪切变形(或相对运动)的一种属性。
2.2 流体的连续介质假说
2.2.3 流体的性质及分类
¾ 可压缩性 所有流体都具有一定的可压缩性
压缩性系数β
β =−
⎞ ⎟ 1 ⎛ dκ ⎞ 1 ⎛ d ρ ⎞ ⎝ ⎠T ⎝ ⎠T ⇒β =− ⎜ ⎟ = ⎜ ⎟ dκ dρ κ ⎝ dp ⎠T ρ ⎝ dp ⎠T ρκ = 1 ⇒ ρdκ + κ d ρ = 0 ⇒ = − κ ρ
质点跟踪 — 着眼于流体介质中的每个质点,这就需 要对流体介质中的每个质点进行辨识 — 三维标号 — 通常用初始时刻流体质点在空间中所处的不同位置作 为区别流体各质点的标记:
t = t 0 : ( a , b, c )
对于给定的拉格朗日变数(a, b, c),空间中某个质 点在任意时刻t 所在的空间位置可表示为:
相关文档
最新文档