流体力学-讲义2

合集下载

高等计算流体力学讲义(2)

高等计算流体力学讲义(2)

高等计算流体力学讲义(2)第二章 可压缩流动的数值方法§1. Euler 方程的基本理论 0 概述在计算流体力学中,传统上,针对可压缩Navier -Stokes 方程的无粘部分和粘性部分分别构造数值方法。

其中最为困难和复杂的是无粘部分的离散方法;而粘性项的离散相对简单,一般采用中心差分离散。

所以,本章主要研究无粘的Euler 方程的解法。

在推广到Navier -Stokes 方程时,只需在Euler 方程的基础上,加上粘性项的离散即可。

Euler 方程是一种典型的非线性守恒系统。

下面我们将讨论一般的非线性守恒系统以及Euler 方程的一些数学理论,作为研究数值方法的基础。

1非线性守恒系统和Euler 方程一维一阶非线性守恒系统(守恒律)可写为下列一般形式=∂∂+∂∂xF tU ,0,>∈t R x(1)其中U 称为守恒变量,是有m 个分量的列向量,即T m u u u U ),...,(21=。

T m f f f F ),...,(21=称为通量函数,是U 的充分光滑的函数,且满足归零条件,即:0)(lim=→U F U即通量是对守恒变量的输运,守恒变量为零时,通量也为零。

守恒律的物理意义设U 的初始值为:0(,0)(),U x U x x =∈R 。

如果0()U x 在x ∈R 中有紧支集(即0U 在有限区域以外恒为零),则0(,)()U x t dx U x dx =⎰⎰RR。

即此时虽然(,)U x t 的分布可以随时间变化,但其总量保持守恒。

多维守恒律可以写为)(=++∙∇+∂∂k H j G i F tU(2)守恒律的空间导数项可以写为散度形式。

守恒系统(1)可以展开成所谓拟线性形式)(=∂∂+∂∂xU U A tU (3)A 是m m ⨯矩阵,称为系数矩阵或Jacobi 矩阵,其具体形式为111122221212.........m m m m mm f f f u u u f f f u u u A f f f u u u ∂∂∂⎡⎤⎢⎥∂∂∂⎢⎥∂∂∂⎢⎥⎢⎥∂∂∂=⎢⎥⎢⎥⎢⎥∂∂∂⎢⎥∂∂∂⎢⎥⎣⎦(4),容易验证:F U Axx∂∂=∂∂,通常也记F A U∂=∂。

流体力学2章讲稿

流体力学2章讲稿

第二章 流体运动学只研究流体运动, 不涉及力、质量等与动力学有关的物理量。

§2.1 流体运动的描述 两种研究方法:(1)拉格朗日(Lagrange)法: 以流场中质点或质点系为研究对象, 从而进一步研究整个流体。

理论力学中使用的质点系力学方法,难测量,不适用于实用理论研究。

(2)欧拉(Euler)法: 将流过空间的流体物理参数赋予各空间点(构成流场),以空间各点为研究对象,研究其物理参数随时间t ,位置(x ,y ,z )的变化规律。

易实验研究,流体力学的主要研究方法。

两种研究方法得到的结论形式不同,但结论的物理相同。

可通过一定公式转换。

1. 拉格朗日法有关结论质点: r=r (t ) dt d rV = dtd dt d V r a ==22x=x (t ) dt dxu = 22dtx d a x =y=y (t ) dtdyv = 22dt y d a y =p=p (t ) T=T (t ) .. .. .. .. .. .. .. .. 质点系:x=x (t,a,b,c ) p=p (t,a,b,c ) T=T (t,a,b,c ) .. .. .. .. .. .. .. ..(a, b, c)是质点系各质点在t =t 0时刻的坐标。

(a, b, c)不同值表不同质点2. 欧拉法物理量应是时间t 和空间点坐标x, y,z 的函数u =u(x, y, z, t) p =p(x, y, z, t) T =T(x, y, z, t) 3. 流体质点的随体导数!!流体质点的随体导数:流体质点物理参数对于时间的变化率。

简称为质点导数。

例:质点速度的随体导数(加速度)dt d V 质点分速度的随体导数dtdu质点压力的随体导数dtdp质点温度的随体导数dt dT.. .. .. .. .. .. 质点导数是拉格朗日法范畴的概念。

流体质点随体导数式---随体导数的欧拉表达式dt d V =z wy v x u t t∂∂+∂∂+∂∂+∂∂=∇⋅+∂∂V V V V V V Vdt du =z u w y u v x u u t u u tu∂∂+∂∂+∂∂+∂∂=∇⋅+∂∂Vdt dT =z T w y T v x T u t T T tT∂∂+∂∂+∂∂+∂∂=∇⋅+∂∂V普遍形式: dt dF =z F w y F v x F u t F F tF∂∂+∂∂+∂∂+∂∂=∇⋅+∂∂VF t )(∇⋅+∂∂=V证其一: dt d V =V V V∇⋅+∂∂t 由 dt d V=tt ∆-→∆V V 'lim 0因 V=V (x ,y , z,t )V ’=V (x+Δx ,y+Δy ,z+Δz,t+Δt )所以 V ’=V++∆∂∂x x V +∆∂∂y y V z z∆∂∂V t t ∆∂∂+V 代入上式得dt d V==∆∆∂∂+∂∂∆+∂∂∆+∂∂∆→∆tt z z y x xt tV V y V V lim 0V V V z V y V x V t V ∇⋅+∂∂=∂∂+∂∂+∂∂+∂∂=tw v u 可见, 在欧拉法中质点速度的随体导数(即加速度)由两部分组成。

流体力学讲义2-2

流体力学讲义2-2
ij
ji
S ij S ji,有六个独立分量;

ji
只有三个独立分量。
0
V i V i 0 S ij x j ij x j
0
平移
变形
转动
Turbulence Research Laboratory
在直角坐标系中, S ij 、 ij
S ij V j 1 Vi 2 x j x i Vi 1 V j 2 x i x j S ji
右端第二项用ij表示,它是反对称张量,表示微团的转动
ij 1 Vi V j 2 x j x i 1 V j Vi 2 x i x j
S xx
的各个分量表示式如下:
w z
1 w u , 2 x z S x y S yx 1 u v 2 y x
u x
, S yy
v y
, S zz
S yz S zy
1 v w , 2 z y

,微团上参考点
x 0 处的速度为 V 0 V ( x 0 , t )
我们来考察微团上任意一点
x x 0 x
的速度
V ( x 0 x 0 , t )
V V V 2 V( x 1 x 2 x 3 o ( x ) Taylor展开,我们有:x 0 x 0 , t ) V ( x 0 , t ) 利用 x 1 x 2 x 3
O A O A A A O O V xex V0 x t V 0 t x 0 V t ex x x 0 V u v w ex ey ez 代入上式(为了简明起见取消下标0),得: x x x x

流体力学第二部分优秀课件

流体力学第二部分优秀课件
周传辉zhou3126com第三章一元流体动力学基础第一节描述流体运动的两种方法第二节恒定流动和非恒定流动第三节流线和迹线第四节一元流动模型第五节连续性方程第六节恒定元流能量方程第七节过流断面的能量方程第三章一元流体动力学基础第八节恒定总流的能量方程第九节能量方程的应用第十节总水头线和测压管水头线第十一节恒定气流能量方程第十二节总压线和全压线第十三节恒定流动量方程流场
1 2
Q1 Q2 v1 A1 v 2 A2
——不可压缩流体恒定总流的连续性方程
综合:
表明:在不可压缩流体一元流动中,平均流速与断面面积成反比.
第六节 恒定元流能量方程
假设条件:
不可压恒定流体、只有重力作用。在某一恒定流场中,任取一个元流,在元流上 任取两个过流断面1-1,2-2。
设:断面面积为dA1,dA2, 速度为u1,u2, 压强为p1,p2, 位置高度z1,z2
第二节 恒定流动与非恒定流动
非恒定流动:运动不平衡的流动,各点的流速随时间变化, 由流速决定的压强,粘性力和惯性力也随时间变化。
ux
u
x
(
a
,
b
,
c
,
t
)
uy uy (a,b,c, t)
uz
uz (a,b,c, t)
恒定流动:运动平衡的流动,各点的流速不随时间变化, 由流速决定的压强,粘性力和惯性力也不随时间变化。
的直线)。长而直的管道内的流动就是均匀流。 非均匀流:过流断面沿流程变化。比如:弯管、变径管,由于过流断面
的变化,引起流速的大小或方向发生变化。
第七节 过流断面的压强分布
2、渐变流与急变流 渐变流:流速沿流动方向变化极为缓慢地非均匀流。渐变流的流线趋近
于平行的直线,因此渐变流的过流断面可以近似的认为是平面 (过流断面有时是曲面)。 急变流:流速沿流程变化显著的流动

流体力学讲义

流体力学讲义

140第六章、 流體動量分析(Momentum analysis offlow systems )牛頓第二定率 – 動量守衡牛頓第二定律: ∑===F dtV m d dt V d m a m)({}⎭⎬⎫⎩⎨⎧++=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧)(viscous pressure forcces surface force body system the of m om entum the of change of rate tim e⎰∑=syssys F V d V Dt Dρ 假設系統與控制容積於時間 t 時互相重疊,如下圖所示:∑∑=CVcoincident the of contents sys F F則由雷諾轉換定理,∑∑⎰⎰⎰⎰-+∂∂=∙+∂∂=in in in in out out out out CVCS CVsys V A V V A V V d V t dA n V V V d V t V d V Dt D ρρρρρρ)(或141⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧+⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧....V C the of out m om entum of flow of rate net V C coincident the of m om entum the of change of rate tim e system coincident the of m om entum the of change rate tim e 故以控制容積而言∑∑∑⎰⎰⎰=-+∂∂=∙+∂∂CVtheof contents in in in in out out out out CV CS CVF V A V V A V V d V t dA n V V V d V t ρρρρρ)( (注意:上式中,每一項單位均為 kg.m/s ,並為一向量方程式,故有三分量。

) 此式可以下式表示之:∑=+-CVtheof contents F S I O∑=⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧⋅+⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⋅-⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⋅ii F s m kg CV the in m om entum of rate torage S s m kg m om entum of rate nflow I s m kg m om entum of rate utflow O )/()/()/(此為控制容積法表示之動量守衡定律。

流体力学完整讲义

流体力学完整讲义

流体力学一、流体静力学基础 包括内容三部分:01流体主要物理特性与牛顿内摩擦定律 02流体静压强 03流体总压力01流体主要物理特性与牛顿内摩擦定律 水银的密度13.6g/cm 3重度γ(也成为容重,N/m3),单位体积流体所具有的能量。

=g γρ流体的压缩系数:1=pa d dV V dp dpρρβ-=-(单位:) ,β值越大,流体的压缩性也越大。

压缩系数的倒数成为流体的弹性模量,用表示,21()dpdV V β=-k=单位:pa=N/m流体的体膨胀系数a :1=(:)d dVV a T dT dTρρ--=单位质量力:大小与流体的质量成正比(对于均质流体,质量与体积成正比,故又称为体积力)表面力:作用在流体表面的力,大小与面积成正比,它在隔离体表面呈连续分布,可分为垂直于作用面的压力和平行于作用面的切力。

流体的黏性:流体内部质点间或流层间因相对运动而产生内摩擦力以反抗相对运动的性质叫做黏性。

此内摩擦力成为黏制力。

du d T AA dy dtθμμ== 式中:T 流体的内摩擦力μ为流体的动力黏度,单位Pa s •。

A 为流体与管壁的接触面积dudy为速度梯度,表示速度沿垂直于速度y 轴方向的变化率 d dtθ为角变形速度 气体动力黏度随温度的升高而增加。

液体动力黏度随温度的升高而降低,例如:油。

运动黏度v (单位:2/m s )(相对黏性系数):v μρ=理想流体:假想的无黏性的流体,即理想流体流过任何管道均不会产生能量损失。

[推导过程]:tan()dudt d d dy θθ≈=,即:d dudt dyθ=。

02流体静压强流体净压强的特性:①流体静压强方向与作用面垂直;②各向等值性:静止或相对静止的流体中,任一点的静压强的大小与作用面方向无关,只于该点的位置有关。

帕斯卡定律:0P P gh ρ=+式中:P 为液体内某点的压强0P 为液面气体压强 h 为某点在液面下的深度等压面:流体中压强相等的点所组成的面成为等压面。

流体力学第二章ppt课件

流体力学第二章ppt课件
o
P ghC A 225kN
yC
4 sin 60
11
6.6m
IC
b 12
h3
4 3
1.33m4
4m
C D
60° y
yD
yC
IC yC A
6.6
1.33 6.6 4
6.6
0.05
6.65m
yC
图解法(求解矩形平面)
1 水静压强分布图 用一定比例的线段表示压强的大小。 与作用面垂直的箭头表示压强的方向。
(H 13.6103 kg/m 3, 1103 kg/m 3 )
解题步骤
解:
已知断面1上作用着大气压, 因此可以从点1开始,通过等 , 压面,并应用流体静力学基 本方程式,逐点推算,最后 便可求得A点压强。
, 因2-2、3-3、4-4为等压面,根据静压强公式可得
p2 H g(1 2 )
p3 p2 g(3 2 )
根据力的作用方式不同
质量力:指某种力场作用在流体的每一个质点上,大小 与受作用的流体质量成正比的力。
lim X
FBX
V M m
单位质量力轴向分力
lim Y
FBY
V M m
lim Z
FBZ
V M m
单位:N/kg
表面力:是指作用于流体表面上,大小与作用表面积成 正比的力。
P
法向分力
p lim A A A
➢与两流层间的速度差du及流层的接触面积A成正比,和流层间距dy成反比。 ➢与流体种类有关。 ➢与流体的压力大小无关。
T A du dy
T A du 或 du
dy
dy
牛顿内摩擦定律
§1.3 流体的力学模型

流体力学讲义 第二章 流体静力学

流体力学讲义 第二章 流体静力学

第二章流体静力学作用在流体上的力有面积力与质量力。

静止流体中,面积力只有压应力——压强。

流体静力学主要研究流体在静止状态下的力学规律:它以压强为中心,主要阐述流体静压强的特性,静压强的分布规律,欧拉平衡微分方程,等压面概念,作用在平面上或曲面上静水总压力的计算方法,以及应用流体静力学原理来解决潜体与浮体的稳定性问题等。

第一节作用于流体上的力一、分类1.按物理性质的不同分类:重力、摩擦力、惯性力、弹性力、表面张力等。

2.按作用方式分:质量力和面积力。

二、质量力1.质量力(mass force):是指作用于隔离体内每一流体质点上的力,它的大小与质量成正比。

对于均质流体(各点密度相同的流体),质量力与流体体积成正比,其质量力又称为体积力。

单位牛顿(N)。

2.单位质量力:单位质量流体所受到的质量力。

(2-1) 单位质量力的单位:m/s2 ,与加速度单位一致。

最常见的质量力有:重力、惯性力。

问题1:比较重力场(质量力只有重力)中,水和水银所受的单位质量力f水和f水银的大小?A. f水<f水银;B. f水=f水银;C. f水>f水银;D、不一定。

问题2:试问自由落体和加速度a向x方向运动状态下的液体所受的单位质量力大小(fX. fY. fZ)分别为多少?自由落体:X=Y=0,Z=0。

加速运动:X=-a,Y=0,Z=-g。

三、面积力1.面积力(surface force):又称表面力,是毗邻流体或其它物体作用在隔离体表面上的直接施加的接触力。

它的大小与作用面面积成正比。

表面力按作用方向可分为:压力:垂直于作用面。

切力:平行于作用面。

2.应力:单位面积上的表面力,单位:或图2-1压强(2-2)切应力(2-3) 考考你1.静止的流体受到哪几种力的作用?重力与压应力,无法承受剪切力。

2.理想流体受到哪几种力的作用?重力与压应力,因为无粘性,故无剪切力。

第二节流体静压强特性一、静止流体中任一点应力的特性1.静止流体表面应力只能是压应力或压强,且静水压强方向与作用面的内法线方向重合。

流体力学_lecture2_前言_流体性质(1次课)详解

流体力学_lecture2_前言_流体性质(1次课)详解

E
1
p
lim (
V 0
p V
/V
) V
lim
V 0
p V
V
dp dV
20
§1-2 前言
1 流体的性质
• 压缩系数及弹性模量的密度表达方式
– 质量守恒
V C
– 微分
dV Vd 0
P
d dp
1
dV Vd
E dp d
dV d
V
E 的单位:bar或Pa,与压强的单位相同
物理意义:相对变化率;E大p小不易压缩;
l 原油的动力粘度 0.021Pa s
轻油的动力粘度 0.0021Pa s
18
§1-2 前言
1 流体的性质
• 粘性
– 流体粘性的应用 是一切动力装置中不可缺少的
润滑——机床导轨 空气轴承——牙钻(20万转) 水润滑——冰块在冰上滑行 油轴承 汽轮机滑动轴承
铁路车辆滑动轴承 汽车轮胎的沟纹
三种温度的换算
Air
28.96 287
摄氏温度 t ,C
CO2 CO
44.01 188.9 28.01 296.5 开氏温度 T=273+t ,K
H2
2.016 4124 华氏温度 F=t9/5+32 ,F
O2
32.00 259.8
22
§1-2 前言
1 流体的性质
气体的弹性模量或体积压缩系数
等温压缩过程:T=c
上盘下表面切应力为 r
B点微元摩擦面积为
dA 2rdr
流体对微元表面的摩擦力
dF dA 2
r2dr
流体对微元表面的摩擦力矩 dT dF r 2 r3dr

清华大学任玉新---高等计算流体力学

清华大学任玉新---高等计算流体力学

主要内容任玉新1.Basics2.Methods for compressible flows1) The mathematical properties of Euler equations2) Shock wave and entropy conditions3) Riemann problem and the Godunov scheme4) Approximate Riemann solvers: HLL solver and Roe solver5) TVD scheme6) ENO/WENO scheme7) The compact scheme3.Methods for incompressible flows1) The staggered and the colocated grids2) The MAC method3) The SIMPLE method4) The projection method5) Other methods6) Solution of N-S equations on the nonstaggered gridReferences[1] E. F. Toro, Riemann solvers and numerical methods for fluid dynamics, Springer, 1997 (First edition)[2] J.D. Anderson, Computational fluid dynamics: basics with applications, Springer (清华大学出版社影印版)[3] Barth and Deconinck (eds.) High order method for computational physics, Lecture Notes in Computational Science and Engineering, 9. Springer, 1999[4] Ferziger and Peric, Computational method for fluid dynamics, Springer, 1996[5] T. J. Chung, Computational fluid dynamics, Cambridge University Press, 2002[6] J. W. Thomas, Numerical partial differential equations: conservation laws and elliptic equations. Texts in applied mathematics 33, Springer, 1999[7] 吴子牛,计算流体力学基本原理,科学出版社, 2002.[8] Sherrie L. Krist, Robert T. Biedron, Christopher L. Rumsey,CFL3D User's Manual, The NASA Langley Research Center,Hampton, VA[9] S. K. Lele, J. Comput. Phys. 103, 16 (1992)[10] S. Pirozzoli, J. Comput. Phys. 178 (2002)[11]Yu-Xin Ren, Miao'er Liu, Hanxin Zhang, J. Comput. Phys. 192 (2003)3FTP: 166.111.37.201Usr:cfdPasswd:cfd2005Email:****************.cn高等计算流体力学讲义(1)第一章计算流体力学基本原理第1节流体力学基本方程一、非定常可压缩Navier-Stokes方程5不计品质力的情况下,在直角坐标系中,守恒型N -S 方程可以写为下列向量形式:()()()0v v v t x y z∂∂-∂-∂-+++=∂∂∂∂U F F G G H H , (1) 其中u v w E ρρρρρ⎛⎫ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪⎝⎭U 2()u u p uv uw E p u ρρρρρ⎛⎫ ⎪+ ⎪ ⎪= ⎪ ⎪ ⎪+⎝⎭F 2()v vu v p vw E p v ρρρρρ⎛⎫ ⎪ ⎪ ⎪=+ ⎪ ⎪ ⎪+⎝⎭G 2()w uw vw w p E p w ρρρρρ⎛⎫ ⎪ ⎪⎪= ⎪+ ⎪ ⎪+⎝⎭H ,0xx xyv xzxx xy xz T u v w kx ττττττ⎛⎫⎪ ⎪⎪=⎪ ⎪⎪∂+++ ⎪∂⎝⎭F 0xy yy v yzxy yy yz T u v w k y ττττττ⎛⎫⎪ ⎪ ⎪= ⎪ ⎪ ⎪∂+++ ⎪∂⎝⎭G ,0xz zyv zzxz zy zz T u v w k z ττττττ⎛⎫⎪ ⎪ ⎪= ⎪ ⎪ ⎪∂+++ ⎪∂⎝⎭H 。

名师讲义【中国计量大学】工程流体力学第二章 流体静力学

名师讲义【中国计量大学】工程流体力学第二章 流体静力学

用dx、dy、dz除以上式,并化简得
X 1 p 0 (1) x
同理
Y 1 p 0 (2) —欧拉平衡微分方程(2.4)
y
Z 1 p 0 (3)
z
意义:平衡流体所受的质量力分量等于表面力分量。该
方程用于可压、不可压流体,理想和黏性流体。是流体静 力学最基本的方程。
9
现代设计制造研究所
18
现代设计制造研究所
静止液体中的压强计算和等压面
等压面
1、在重力作用下,不可压缩静止流体中的等
高面为等压面; 2、自由表面。
p p0 gz0 z p0 gh
静压强分布
19
现代设计制造研究所
静止液体中的压强计算和等压面
习题1:水池中盛水如图。已知液面压强 p0 98.07kN / m2,
解:圆柱体底面上各点所受到的计示 压强为:
F mg 100 5.1 9.807
pe d 2 / 4 0.7854 (0.12)2 13263(Pa)
pa F
H h
pe g(h H )
1
H pe h 0.8524(m)
g
w 1
d
24
现代设计制造研究所
流体静压强的测量
1. 静压强的单位
物理意义:在重力作用下的连续均质不可压静止流体
中,各点单位重量流体的总势能保持不变(能量守恒)。
16
现代设计制造研究所
静止液体中的压强计算和等压面
p gz C
p gz p0
C由边界条件确定。如果假定在液
面上,Z=0,p=p0,则C=p0。
p p0 gz
如果选取h的坐标方向与z轴相反,则: p p0 gh
积分 p gz c

流体力学实验讲义

流体力学实验讲义

流体力学实验讲义临沂大学建筑学院实验二雷诺实验实验三伯努利方程实验实验四、五流体力学综合实验实验二 雷诺实验一、实验目的1、观察流体在管内流动的两种不同流型。

2、测定临界雷诺数Re c 。

二、基本原理流体流动有两种不同型态,即层流(或称滞流,Laminar flow )和湍流(或称紊流,Turbulent flow ),这一现象最早是由雷诺(Reynolds )于1883年首先发现的。

流体作层流流动时,其流体质点作平行于管轴的直线运动,且在径向无脉动;流体作湍流流动时,其流体质点除沿管轴方向作向前运动外,还在径向作脉动,从而在宏观上显示出紊乱地向各个方向作不规则的运动。

流体流动型态可用雷诺准数(Re )来判断,这是一个由各影响变量组合而成的无因次数群,故其值不会因采用不同的单位制而不同。

但应当注意,数群中各物理量必须采用同一单位制。

若流体在圆管内流动,则雷诺准数可用下式表示:μρdu =Re (1-1)式中:Re —雷诺准数,无因次;d —管子内径,m ;u —流体在管内的平均流速,m /s ;ρ—流体密度,kg /m 3;μ—流体粘度;Pa ·s 。

层流转变为湍流时的雷诺数称为临界雷诺数,用Re c 表示。

工程上一般认为,流体在直圆管内流动时,当Re ≤2000时为层流;当Re ≥4000时为湍流;当Re 在2000至4000范围内,流动处于一种过渡状态,可能是层流,也可能是湍流,或者是二者交替出现,这要视外界干扰而定,一般称这一Re 数范围为过渡区。

式(1-1)表明,对于一定温度的流体,在特定的圆管内流动,雷诺准数仅与流体流速有关。

本实验即是通过改变流体在管内的速度,观察在不同雷诺准数下流体的流动型态。

三、实验装置及流程实验装置如图1-1所示。

主要由玻璃试验导管、流量计、流量调节阀、低位贮水槽、循环水泵、稳压溢流水槽等部分组成,演示主管路为220⨯φmm 硬质玻璃。

637891-红墨水储槽;2-溢流稳压槽;3-实验管;4-转子流量计;5-循环泵;6-上水管;7-溢流回水管;8-调节阀;9-储水槽图1-1 流体流型演示实验装置示意图实验前,先将水充满低位贮水槽,关闭流量计后的调节阀,然后启动循环水泵。

计算流体力学(中科院力学所)_第2讲-双曲型方程组

计算流体力学(中科院力学所)_第2讲-双曲型方程组

令:
R=u+
∫ ρ dρ
c
同理,沿特征线 : 同理,沿特征线2: 对于等熵完全气体
dx / dt = u c
2c R=u+ γ 1 2c S = u + γ 1
du c dρ + =0 沿特征线1: 沿特征线 : dα ρ dα u 1 c S = + dρ 2 2 ρ 保持不变 dR / dα =
A sin x 0 ≤ x ≤ 2π u ( x,0) = 0 others ρ ( x,0) = 1; p( x,0) = 1
考虑一维无粘流动( 方程), 考虑一维无粘流动(Euler方程),初始时 方程),初始时 刻(t=0)流动状态如下: )流动状态如下:
xa ≤ x ≤ xb u ′( x), ρ ′( x) u, ρ = 0, ρ 0 (= const ) others
(3) ) C (2) ) (1) )
x B
A
给定x3,t3 利用 (假设t3充分小) 给定
x3 x1 = (u1 + c1 )(t3 t1 ) x3 x2 = (u 2 c2 )(t3 t 2 )
区域( ),( ),(4) 区域(2),( ) 未扰动 区域( ) 区域(1)内的流动使用基本 方法计算
双曲型
Copyright by Li Xinliang
2
1) 一阶常系数偏微方程组
U U +A =0 x t U = (u1 , u 2 ,......u m )T
如果矩阵A 可以被对角化: 如果矩阵 可以被对角化: A = S 1 ΛS
U U + S 1 ΛS =0 t x S U U + ΛS =0 t x

高等流体力学第二部分讲义

高等流体力学第二部分讲义

p y
dxdydz
z方向,微元流体所受合压力
C
D.Βιβλιοθήκη NBp zdxdydz
微元流体所受合压力
A ZY
∂p ∂p ∂p
X
- ( ∂xi + ∂yj+ ∂zk)dxdydz
G
H
.M
.
OF
E
第二章 流体静力学
2、微元体所受的质量力:
F=F i +F j+F k=(Xi +Yj+Zk)ρdxdydz
绝对真空
则:绝对压强=相对压强+大气压强 p´=p+pa
第二章 流体静力学
绝对压强总是≥0,但相对压强不一定。若某流体
点处在B点,从图可知,B点相对压强为负。
pv=pa- p´
p
2、压强的度量单位
(1) 以压强的基本定义出
A.
. A点相对压强 大气压强pa
B
真空度
发即单位面积上的压力,单位 A点绝对压强 B点绝对压强 绝对真空
hD hC h
o α
a y
左侧受水压力,水面大 气压强为pa,在平板表面所在 y b 的平面上建立坐标,原点o取在平板
. .dA C
.
yC yD
x
D
表面与液面的交线上,ox轴与交线重合,oy轴沿平
板向下。
第二章 流体静力学
则微元面dA所受压强p=γh
压力dP=pdA=γhdA=γysinαdA
整个平面由无数dA组成, 则整个平板所受水静压力 由dP求和得到。
第二章 流体静力学
第五节 压强的计算基准和度量单位
1、 计算基准
(1) 绝对压强:
以无一点气体存在的绝对真空为零点起算的压

流体力学II教材讲解

流体力学II教材讲解

流体力学II(Viscous Fluid and Gas Dynamics)讲义第一章、粘性不可压缩流体运动基本方程组(学时数:6)1-1.绪论流体力学是力学的一个重要分支,主要研究流体介质(液体、气体、等离子体)的特性、状态,在各种力的作用下发生的对流、扩散、旋涡、波动现象和质量、动量、能量传输,以及同化学、生物等其他运动形式之间的相互作用。

它既是一门经典学科,又是一门现代学科,对自然科学和工程技术具有先导作用。

历史上,力学包括流体力学,曾经经历基于直观实践经验的古代力学、基于严密数学理论的经典力学、基于物理洞察能力的近代力学三个阶段。

在人类早期的生产活动过程中,力学即与数学、天文学一起发展。

17世纪,Newton基于前人的天文观测和力学实验,发明了微积分,并总结出机械运动三大定律和万有引力定律,发表了著名的《自然哲学的数学原理》一书。

由于原理是普适自然与工程领域的规律,从而使力学成为自然科学的先导。

从17世纪开始,人们逐步建立了流体力学的基本理论体系,从Pascal定律、Newton粘性定律、Pitot 管测速,到Euler方程和Bernoulli方程,标志着流体动力学正式成为力学的一个分支学科。

18世纪,人们着重发展无粘流体的位势理论。

到了19世纪,为了解决工程实际问题,开始注重粘性的影响,Navier-Stokes方程的建立为流体力学的进一步发展奠定了完整的理论基础,但该方程解的存在性与光滑性的证明至今仍是一大难题。

20世纪初,Prandtl凭借出色的物理洞察能力,提出边界层理论,从而开创了流体力学的近代发展阶段,使力学成为人类实现“飞天”梦想的重要理论先导。

60年代以来,由于超级计算机、先进测试技术的发展和应用,力学进一步凸显宏微观结合和学科交叉的特征,进入现代力学发展新阶段。

刚刚过去的2011年,人类遭遇了一系列极端事件:日本海底地震导致海啸和福岛核电站泄露事故;澳大利亚飓风;我国干旱洪水灾害等异常气候问题。

流体力学上课讲义第二章

流体力学上课讲义第二章

Figure 9-110 Cengel & Cimbala Laminar flat plate boundary layer
The no-slip condition 無滑動邊界條件
• No-slip condition: a fluid in direct contact with a solid stick to the surface due to viscous effects, i.e., the fluid has a zero velocity relative to the surface.
g
Specific gravity
= the ratio of the density of a substance to the density of some standard substance at a specified temperature (usually water at 4°C). 比重
Bulk modulus of elasticity
dp change in pressure Ev dV / V fractional change in volume
A minus sign is used in the definition to yield a positive Ev
• Kinematic viscosity

N s/m 2 2 m /s 3 kg/m
(2.8)
Effect of temperature on viscosity
• Viscosity is caused by the cohesive forces between the molecules in liquids and by the molecular collisions in gases. • In a liquid the molecules possess more energy at higher temperatures, so they can move more freely. • In a gas, the intermolecular forces are negligible, and the gas molecules at high temperatures move randomly at higher velocities. This results in more molecular collisions per unit volume per unit time and therefore in greater resistance to flow.

流体力学讲义2-1

流体力学讲义2-1

(4)例:Euler描述法中质点轨迹的计算 已知 x t A V x ,t 计算 x x A ,t
给定
V1 x 1 t , V2 x 2 t , V3 0
x1 t x 2 t x 3 t
求:t=0时位于(-1,-1,0)处的质点运动轨迹。 解:


x3 t V V
V3 V3 V3 V1 V2 V1 V2 e 3 e1 e 1 e2 e 2 e2 e 3 e2 e 1 e3 e 2 e3 e 3 e3 x1 x2 x2 x2 x3 x3 x3
w 0
求时 t=0,通过点(-1,-1,0)的流线。 解:
dx dy dz u v w dx dy ,dz 0 x y xy const., z=const.
代入起始点(-1,-1,0),得
xy 1, z=0
此流线是双曲线
(4) 流线的性质
①非定常流场中不同时刻通过同一点的流线一般不重合。 定常速度场和时间无关,
V1 x1 , x 2 , x3 ,t 0
dx1

V 2 x1 , x 2 , x3 ,t 0
dx 2

V3 x1 , x 2 , x3 ,t 0
dx3
注意: (1)在流线方程中时间是常参量
(2) 流线由不同质点组成
(3)计算流线举例
给定平面速度场
u x t,
v y t,
任何时刻通过同一点的流线 都是相同的。 ② 流线是向量线,是有走向的几何线,流线的走向由速度场给出。 ③一般情况下,同一时刻通过一点只有一根流线。般情况下,同一时刻流场中的流 线不能相交。 有两种例外: (i)流场中速度等于零的点,因为速度等于零的点其速度方向是任意的,不同方向 流线都可以通过该点; (ii)流场中速度为无限大的点,因为这种点上流体的速度方向也是不定的。 ④非定常流场中,一般情况下,流线和迹线不重合;定常流场中,通过同一点的流 线和迹线重合。 ⑤垂直流线方向的流体质点速度等于零。

流体力学讲义1-2

流体力学讲义1-2


;
运动粘性系数:



,
p
量纲 : [] [F][L]2[T ]
[ ] [L]2[T ]1
国际单位 NS/m2
m2/S
实验表明,在常温常E 压下:
液体:
气体:
3、压缩性 (可压缩流体、不可压缩流体)
3、一点上的应力及其性质
应力
Tn
是 向量,一般不垂直于它的作用面,
可分解为 Tnn Tns Tnt
约定
n,
t,
s
右手直角坐标系
nr Tn
n
Tns `sr
r t
Tns
定义:应力向量在作用面法线方向的分量称为正应力
x2
Tn
( 正值正应力是拉力;负的正应力指向作用面内是压力) 定义:应力向量在作用面切向的分量称为剪应力,
r
rr r
令V =, 求V 和gV .
Turbulence Research Laboratory
R1 R2 R
R
R d
2 cos
p1
pa

2g
R

pa
4g cos
d
静止液体中压强公式: p1 pa gH
H为液柱高,消去Pa后得
H 4g cos gd
空气和水银
毛细现象
Turbulence Research Laboratory
课后学习与练习
➢ 自学《附录》386-397页 ➢ 习题 ➢ (1)《补充习题》0.1-0.7
(3)
ar

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

只有重力作用下:
X 0、Y 0、Z g
dp (gdz) dz
const p z c
z1
p1
z2
p2
2020/10/16
液体的表面上
z 0, p p0 c p0
p p0 z p p0 h
24
2.2 欧拉平衡微分方程
❖重力作用下流体的压强分布规律
静止液体的基本方程(静止流体能量守恒方程),其适用条 件为(1)均质或不可压缩流体;(2)体积力为重力;(3)同 种流体的连通范围内。
2020/10/16
30
流体静力学
1
流体静压强的特性
2
欧拉平衡微分方程
3
液体压强的测量
4
静止流体对平面的作用力
5
静止流体对曲面的作用力
2020/10/16
31
2.3 液体压强的测量
❖绝对压强、相对压强、真空度
p p0 h p - pa h
p h
绝对压强(absolute pressure):以完全真空状态 为基准算起的压强
特性一 静止液体压强(压力)垂直指向作用面。
❖反证法:若存在拉应力或切应力,由流体易流动性质可知, 平衡状态破坏,产生流动。
2020/10/16
5
2.1 流体静压强的特性
❖流体静压强的两个特性:
特性二 静止液体中任意点的压强值
大小与作用面的方位无关。
❖证明:对右图四面体做受力分析 ❖如图所示,设在密度为ρ的静止 流体中有边长为 dx,dy ,dz 的四面体微元 PABC
❖等压面(isobaric surface):
在同一种液体中,如果各处的压强均相等,即 p const,
由各压强相等的点组成的面称为等压面。
因为 p C dp 0
由 dp dU 可知 dp dU 0
等压面性质
(1)等压面与等势面重合
dp dU 0
U C
2020/10/16
20
2.2 欧拉平衡微分方程
2020/10/16
22
2.2 欧拉平衡微分方程
❖等压面(isobaric surface):
等压面的选取 只有重力作用下,等压面选取的条件: (1)静止; (2)连通; (3)同种 (4)同一水平面。
连通容器
2020/10/16
连通容器
连通器被隔断
23
2.2 欧拉平衡微分方程
❖重力作用下流体的压强分布规律
作用在六面体微元体的面ABCD上的总压力为 pdydz ,
作用在从面ABCD开始计算起在 x方向上距离为 dx的面EF
GH上的总压力p 为 p / x dx dydz
由此可得,x方向上力 的平衡条件为
pdydz
X
dxdydz
p
p x
dx
dydz
0
由上式可得
p / x X
平衡流体微元体
2020/10/16
流体静力学—引言
❖本章学习要点: ⑴静止流体中应力的特性 ⑵流体平衡微分方程,等压面。 ⑶静止液体和相对静止液体压强的分布。绝对 压强、相对压强、真空度,测压管水头。 ⑷液体作用在平面上的总压力,压力中心,压 强分布图法。 ⑸液体作用在曲面上的总压力,压力体。 ⑹浮力,浮体的稳定性(介绍)。
2020/10/16
p dx p dy p dz (X dx Y dy Z dz)
x y z
2020/10/16
17
2.2 欧拉平衡微分方程
❖欧拉平衡方程(Euler′s equilibrium equation):
因为 p p(x, y, z)
所以 dp p dx p dy p dz x y z
即 dp (Xdx Ydy Zdz)
px d
Ax
1 2
px dydz
Py
py d
Ay
1 2
py
dxdz
Pz
pz
d
Az
1 2
pz dxdy
Pn pn d An
2020/10/16
8
2.1 流体静压强的特性
❖流体静压强的两个特性:
❖质量力 F, 分量为Fx、Fy、Fz
单位质量的质量力: X、Y、Z
Fx
1 6
dxdydzX
1 6
dxdydzX
相对压力(gage pressure) 以同高程大气压为为基准 起算的压强
2020/10/16
32
2.3 液体压强的测量
❖绝对压强、相对压强、真空度
pabs pa p
p pabs pa pabs pa
p为负值,我们说该处有真空存在
负压=真空 真空度 pV 绝对压强不足一 个大气压的不足部分
pV | pabs pa || p |
2020/10/16
33
2.3 液体压强的测量
❖绝对压强、相对压强、真空度
大气压强的测量(气压计) Evangelista Torricelli(1608-1647)
patm ρgh
一个标准大气压, 0℃海平面大气压(1atm=101325pa)
一个工程大气压, 1at=98070pa=1 kgf/cm2
p p0 h
上式说明: (1)静止液体中,任意点的压强由两部分组成,一部 分是表面压强(当液面与大气接触时,表面压强为当地大 气压强);另一部分是液重压强。液重压强与液面以下的 深度成线性关系。 (2)表面压强与液重压强无关。液面压强的增量同时 等值地传递到液体中每一点,这就是著名的巴斯卡原理。
2020/10/16
36
2.3 液体压强的测量
❖绝对压强、相对压强、真空度 P14,例题2.4 通过上面的例子可以看出,压强的大小可 用单位面积上的力或水柱高度、水银柱高 度表示,应用这些原理可以制成测量压强 的仪器,测压管、水银测压计、水银差压 计和金属测压计等。
2020/10/16
37
2.3 液体压强的测量
13
流体静力学
1
流体静压强的特性
2
欧拉平衡微分方程
3
液体压强的测量
4
静止流体对平面的作用力
5
静止流体对曲面的作用力来自2020/10/1614
2.2 欧拉平衡微分方程
❖欧拉平衡方程(Euler′s equilibrium equation): 在密度ρ为常数的静止流体中,有边长为 dx,dy,dz的平
16
2.2 欧拉平衡微分方程
❖欧拉平衡方程(Euler′s equilibrium equation): 同样地,根据y, z方向上的力平衡条件可得
p / y Y p / z Z
上述三式就是静止液体的平衡微分方程,也称为欧拉平 衡微分方程。
对上述三式的 x, y, z分量分别乘以dx, dy, dz ,相加合并 之后可得
❖等压面(isobaric surface): 等压面性质
(2)质量力与等压面正交
f Xi Yj Zk
ds dxi dyj dzk
等压面和质量力向量的关系
f ds Xdx Ydy Zdz 0
2020/10/16
21
2.2 欧拉平衡微分方程
❖等压面(isobaric surface): 等压面性质 (3)两种不相混合液体的交界面为等压面
F
y
1 6
dxdydzY
Fz
1 6
dxydzZ
该平衡四面体,应满足 F 0
2020/10/16
9
2.1 流体静压强的特性
❖流体静压强的两个特性: x轴方向
y轴方向
z轴方向
2020/10/16
10
2.1 流体静压强的特性
2020/10/16
11
2.1 流体静压强的特性
❖流体静压强的两个特性: 则有
2020/10/16
25
2.2 欧拉平衡微分方程
❖重力作用下流体的压强分布规律
2020/10/16
26
2.2 欧拉平衡微分方程
❖重力作用下流体的压强分布规律
2020/10/16
27
2.2 欧拉平衡微分方程
❖重力作用下流体的压强分布规律 几何意义 z——是研究点相对于某一 水平面的高度,称为位置高度或 z 位置水头。该水平面称为基准面 或参考面;
由于压强单位Pa的量值对于很多实际应用很小,我们通常也 引入1kPa=103Pa,1MPa=106Pa,1bar=105Pa,atm,mmHg
1 kgf/cm2=9.807 N/cm2=9.807× 104 N/m2 =9.807× 104 Pa=0.9807 bar=0.9679 at
m
2020/10/16
1
流体静力学
1
流体静压强的特性
2
欧拉平衡微分方程
3
液体压强的测量
4
静止流体对平面的作用力
5
静止流体对曲面的作用力
2020/10/16
2
2.1 流体静压强的特性
❖流体静压强:静止流体的压力强度,单位面积上的压力大小
❖固体力学中对应的概念“正应力”
❖由定义“单位面积上的压力大小”可知,压强的单位 N/m2 称为“帕斯卡(Pa)”1Pa=1 N/m2
行六面体微元(见图)。作用在该微元体上的表面力必须与 作用在微元体上的质量力达到静平衡状态。
现设流体微团上单位质量力为
Fx Xdxdydz Xdm
Fy Ydxdydz
Fz Zdxdydz
2020/10/16
平衡流体微元体
15
2.2 欧拉平衡微分方程
❖欧拉平衡方程(Euler′s equilibrium equation):
γ 为斜面ABC的法线与x , y , z 轴
相关文档
最新文档