高等流体力学第3讲
- 格式:doc
- 大小:256.00 KB
- 文档页数:13
第三章流体动力学基础本章是流体动力学的基础。
主要阐述了流体运动的两种描述方法,运动流体的基本类别与基本概念,用欧拉法解决运动流体的连续性微分方程、欧拉运动微分方程及N-S方程。
此外,还阐述了无旋流与有旋流的判别,引出了流函数与势函数的概念,并且说明利用流网与势流叠加原理可解决流体的诸多复杂问题。
第一节流体流动的基本概念1.流线(1)流线的定义流线(stream line)是表示某一瞬时流体各点流动趋势的曲线,曲线上任一点的切线方向与该点的流速方向重合。
图3-1为流线谱中显示的流线形状。
(2)流线的作法:在流场中任取一点(如图3-2),绘出某时刻通过该点的流体质点的流速矢量u1,再画出距1点很近的2点在同一时刻通过该处的流体质点的流速矢量u2…,如此继续下去,得一折线1234 …,若各点无限接近,其极限就是某时刻的流线。
流线是欧拉法分析流动的重要概念。
图3-1 图3-2(3)流线的性质(图3-3)a.同一时刻的不同流线,不能相交。
图3-3因为根据流线定义,在交点的液体质点的流速向量应同时与这两条流线相切,即一个质点不可能同时有两个速度向量。
b.流线不能是折线,而是一条光滑的曲线。
因为流体是连续介质,各运动要素是空间的连续函数。
c.流线簇的疏密反映了速度的大小(流线密集的地方流速大,稀疏的地方流速小)。
因为对不可压缩流体,元流的流速与其过水断面面积成反比。
(4)流线的方程(图3-4)根据流线的定义,可以求得流线的微分方程:图3-4设d s为流线上A处的一微元弧长:u为流体质点在A点的流速:因为流速向量与流线相切,即没有垂直于流线的流速分量,u和d s重合。
所以即展开后得到:——流线方程(3-1)(或用它们余弦相等推得)2.迹线(1)迹线的定义迹线(path line)某一质点在某一时段内的运动轨迹线。
图3-5中烟火的轨迹为迹线。
(2)迹线的微分方程(3-2)式中,u x,u y,u z均为时空t,x,y,z的函数,且t是自变量。
Chapter 3 流体运动的基本方程组本章任务:建立控制流动的基本方程组,确定边界条件。
§3.1系统和控制体系统(sys )指给定流体质点组成的流体团,相当于质点或刚体力学中的研究对象——物体;系统在流动过程中可以不断改变自己的位置和形状,但维持其连续性,始终由固定的那些流体质点组成。
系统与外界可以有力的相互作用,可以有动量和能量交换,但是没有物质交换。
控制体(CV )指流动空间内的一个给定空间区域(子空间),其边界面称为控制面(CS )。
控制体一旦选定,其大小、形状和位置都是确定的,有流体不断出入。
物质体元即流体微团。
物质面元可以看成由连续分布的流体质点(看成是没有体积的几何点)构成的面元,物质面元在流动过程中可以变形,但始终由这些流体质点组成。
物质线元可以看成连续分布的流体质点(看成是没有体积的几何点)构成的线元,或者说是连续分布的流体质点的连线线元,物质线元在流动过程中可以变形,但始终由这些流体质点组成。
时间线就是物质线。
(三者如同面团、薄饼和面条) §3.2雷诺输运定理设(),f r t 代表流动的某物理量场(可以是密度场、温度场、动量密度分量场、能量密度场等),t 时刻某流体团(即系统)占据空间τ,取该空间为控制体。
t 时刻该流体团的总f 为()(),I t f r t d ττ=⎰。
(3-1)此I 也是t 时刻控制体内的总f 。
设t t δ+时刻(0t δ→)该系统运动到如图所示位置,占据空间τ',此时系统的总f 为()(),I t t f r t t d τδδτ'+=+⎰。
(3-2)该系统总f 的随体导数()()()0lim t I t t I t DI t Dt tδδδ→+-=。
(3-3)将空间II τ分为与空间I τ重合的部分2τ和其余部分1τ,空间I τ去除2τ后剩余部分记为3τ,于是13ττττ'=+-,(3-4)进而()()()()13I t t I t t I t t I t t τττδδδδ+=+++-+,(3-5)可得()()()()()130lim t I t t I t t I t t I t DI t Dt tττττδδδδδ→+++-+-=()()()()31000lim lim lim t t t I t t I t t I t t I t t t tττττδδδδδδδδδ→→→+++-=+-, (3-6)其中第一项()()()0limt I t t I t I t t t ττδδδ→+-∂=∂。
第三讲 流体静力学一、 静止流体中的应力特性静止流体中,流体质点之间没有相对运动,切应力必然为0,又由于流体分子之间的引力很小,流体质点之间几乎不能承受拉力。
因此,在静止流体中,只能存在指向作用面的法向应力。
即n p =-p n (3-1)式中的p n 就是工程流体力学中的流体静压力。
上式也可以写成张量形式P ==000000p p p -⎡⎤⎢⎥-⎢⎥⎢⎥-⎣⎦=p 00000011⎡⎤⎢⎥1⎢⎥⎢⎥⎣⎦= p I (3-2) 式中I 为单位张量。
静止流体中任意一点处的应力无论来自何方均相等,即任意一点处的静压力与作用方向无关。
二、 欧拉平衡方程惯性坐标系中,任何流体处于静止状态的必要条件是:作用在物体上的合外力为0,即0∑=F (4-3)在静止流场中任取一个流体团作为研究对象,作用在其上的质量力可表示为d ρττ⎰⎰⎰f (a ) 表面力可表示为d d AAp A p A -=-⎰⎰⎰⎰n n (b )根据第一个平衡条件(3-3)可得d d =0Aρτp A τ-⎰⎰⎰⎰⎰f n (c ) 根据高斯定理可知,若物理量p 在封闭空间τ中连续且存在连续的一阶导数,则有d =d Ap A p ττ∇⎰⎰⎰⎰⎰n (d )将(d)式代入(c)式则可得d 0ρp ττ-∇=⎰⎰⎰()f 由于流体团是任意选取的,所以要使上式成立,则被积函数在该体积内任意点上的数值必须为0,于是有=0ρp -∇f或1=p ρ∇f (3-4)这就是欧拉平衡微分方程式,其在直角坐标系中可写为111x yzp f ρx pf ρy p f ρz ⎧∂=⎪∂⎪⎪∂=⎨∂⎪⎪∂=⎪∂⎩(3-5) 同时,合力矩为0是自动满足的。
三、 静压流场的质量力条件(自学)对于所有的静止流体,(3-4)式均成立,现对其两端同时取旋度可得1111==+=p p p p ρρρρ⎛⎫⎛⎫⎛⎫∇⨯∇⨯∇∇⨯∇∇⨯∇∇⨯∇ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭()f上式中应用了标量函数梯度的旋度为0这一结论,现证明之p ∇⨯∇()=p p p xy z ⎛⎫∂∂∂∇⨯++ ⎪∂∂∂⎝⎭i j k=x y z p p p xy z∂∂∂∂∂∂∂∂∂∂∂∂ij k =p p p p p p y z z y x z z x x y y x ⎛⎫⎛⎫∂∂∂∂∂∂∂∂∂∂∂∂⎛⎫---+-⎪ ⎪ ⎪∂∂∂∂∂∂∂∂∂∂∂∂⎝⎭⎝⎭⎝⎭i j k =0(矢量) 将上式与(3-4)式进行点乘则有()1=p p ρρ⎡⎤⎛⎫∇∇⨯∇⨯∇⎢⎥ ⎪⎝⎭⎣⎦f f 上式右端为矢量的混合积,由混合积的定义可知由于三个矢量中有两个同名,所以其值为0,可得()=0∇⨯f f (3-6)由此可以得出结论:流体静止的必要条件是质量力必须满足()=0∇⨯ff 。
(此式自然成立啊,为什么还要证明呢?)对于不可压缩流体,由于密度为常数,则平衡方程(3-4)式可写成=p ρ⎛⎫∇ ⎪⎝⎭f (3-7) 对上式两端取旋度则有=p ρ⎛⎫∇⨯∇⨯∇ ⎪⎝⎭f =0 (3-8)这是不可压缩流体静止时对质量力所附加的限制条件,即质量力必须无旋。
(矢量场的四等价定理:无旋必有势;有势必无旋;环量为零;线积分与路径无关。
) 四、 等压面等压面是指由静压力相等的点组成的面,等压面的性质包括以下几点: 1. 等压面与质量力相垂直; 2. 等压面是质量力的等势面;【证明】设质量力f =X i +Y j +Z k 所代表的质量力场中,若存在标量函数U 满足U Xdx Ydy Zdz =-++⎰或()d d d d U X x Y y Z z =-++则称之为质量的势函数,其中的负号表示质量力做正功时质量力的是函数减小。
在静止流场中任取一个微元矢量力d r =d x i +d y j +d z k ,()d d =d U U U U dx dy dz U Xdx Ydy Zdz xy z ⎛⎫∂∂∂∇=++==-++- ⎪∂∂∂⎝⎭r f r由于d r 是任意的微元矢量,所以有U =-∇f将其代入(3-4)式可得1∇=-∇p U ρ上式表明等压面方程p ∇=0与等势面U ∇=0等价。
【证毕】3. 等压面是不相容流体间的分界面;4. 等压面是等密度面。
五、 静力学基本公式——静压力分布规律21p p ρgh =±它表明:○1重力作用下的均质流体内部的静压力,与深度h 呈线性关系,因此,水坝都设计成上窄下宽的形状;○2静止流体内部任意点的静压力由液面上的静压力p 0与液柱所形成的静压力ρgh 两部分组成,深度h 相同的点静压力相等;○3静止流体边界上压力的变化将均匀地传递到流体中的每一点,这就是著名的帕斯卡定律。
六、 压力标准1. 流体静压力的计量标准流体力学中,静压力的计量有两个标准,一个是以物理真空为零点的标准,称为绝对标准,按照绝对标准计量的压力称为绝对压力;另一个是以当地大气压力为零点的标准,称为相对标准,按照相对标准计量的压力称为相对压力。
2. 流体静压力的表示方法 绝对压力用abp 表示,对敞口容器中液面以下深度为h 的点来讲,其绝对压力可表示为=ab a p p +ρgh (3-9) 绝对压力不能小于零。
工程中绝对压力的数值可以大于当地大气压力,也可以小于当地大气压力。
因此,相对压力便有了正负之分。
当绝对压力大于当地大气压力时,相对压力大于零,称为表压,用pM 来表示,即==M ab a p p p ρgh -。
(3-10)之所以称之为表压是因为压力表所显示的压力就是这个压力。
当绝对压力小于当地大气压力时,相对压力小于零,称为真空压力或真空度,用pv 来表示,即=v a abp p p - (3-11)真空压力采用真空表测量。
绝对压力、表压和真空度之间的关系可用图3-1来表述。
从中可以看出,表压的含义是比当地大气压力大多少,真空度的含义是比当地大气压力小多少。
表压越大压力越大,真空度越大压力越小,反之亦然。
工程流体力学中所说的压力如不特殊说明指的就是表压,为了简便通常将其下标M 省略,仅以p 表示。
至此,可以归纳出以下的关系=ab a p p +p(3-12)=ab a vp p p - (3-13)==ab a vp p p p -- (3-14)真空度与表压的符号相反,真空度不能直接参与计算,计算工程中必须将真空度转换为负的表压代入方程中,这一点在以后的计算中要予以注意。
七、 静力学基本方程式+=p z c ρg(3-15a )1212=p p z z ρgρg++(3-15b )物理意义:z 称为比位能,p /ρg 称为比压能,而z +p /ρg 称为总比能。
因此,静力学基本方程的物理意义是:静止流体中总比能为常数。
几何意义:z 称为位置水头,p /ρg 称为压力水头,而z +p /ρg 称为测压管水头。
因此,静力学基本方程的几何意义是:静止流体中测压管水头为常数。
八、 作用在平面上的总压力sin C C C P =ρgy αA =ρgh A =p A (3-16)式中pC 代表形心C 处的压力。
它表明:作用在任意形状平面上的总压力大小等于该平面的面积与其形心C 处的压力的乘积。
图3-1 绝对压力与相对压力p总压力的作用点称为压力中心。
压力中心的yD 坐标的计算公式为C D C C J y =y +y A(3-17)因为J C /y C A 恒为正值,所以y D >y C ,说明压力中心D 永远低于平面形心C 。
但是,这一结论对水平放置的平面不适用,此时的压力中心与形心重合。
在应用上述计算公式时应该注意以下两点:(1)没有考虑大气压的影响。
这主要是因为工程实际中容器外也会受到大气压的作用,两者形成的总压力相互抵消,所以在计算总压力时不考虑大气压力的影响,而仅仅考虑液体形成的总压力。
(2)在压力中心的计算公式中y 坐标原点的取法。
式(3-17)只是适用于液面压力为大气压时的情形。
即y 的坐标的原点位于自由液面与平面延长线的交点处。
但是,当自由液面上的压力不是大气压时,式(3-17)中y坐标的原点只能在等效自由液面与平面延长线的交点处。
那么什么是等效自由液面呢?现在考察图3-2所示的两种情形,左图为液面压力大于大气压的情形,其液面绝对压力为p 0’=p a +ρgh ,右图为将原有液面升高了h =(p 0’-p a )/ρg 后,且液面绝对压力等于大气压时的情形,两者对平面AB 形成了完全相同的压力分布,同时两者作用在平面上的总压力是完全相同的。
因此,称右图中的自由液面为左图中液面的等效自由液面。
在计算过程中绝对不可以将左图中的O ’点作为y 坐标的原点,而应取右图中的O 点作为y 坐标的原点。
也就是说,在进行压力中心位置计算时,应该将液面压力不是大气压的液面转换成等效自由液面,然后重新找出y 的原点进行计算。
对液面绝对压力低于大气压的情形应该用类似的方法来处理,具体如何处理请读者自己思考。
九、 作用在曲面上的总压力1. 总压力的水平分量x C x P =ρgh A y C y P =ρgh Az C zP =ρgh A图3-2 等效自由液面2. 总压力的铅直分量z P =ρgV3. 压力体压力体可以表述为:压力体是由受力曲面、液体的自由表面(或其延长面)以及两者间的铅垂面所围成的封闭体积。
压力体是从积分式A z hdA ⎰得到的一个体积,它是一个纯数学的概念,与这一体积内是否充满液体无关,图3-3是两个典型的压力体。
比较图中(a )和(b )的压力体,不难发现两者有着明显的不同,○1压力体所形成的总压力方向不同,图3-3a 中压力体形成的总压力方向向下,图3-3b 中的压力体所形成的总压力方向向上;○2是两者压力体与液体所处的位置不同,图3-3a 中的压力体与液体位于曲面的同一侧,图3-3b中的压力体与液体则不在曲面的同一侧。
由此可以引入定义:如果压力体与形成压力的液体在曲面的同侧,则称这样的压力体为实压力体,用(+)来表示;如果压力体与形成压力的液体在曲面的异侧,则称这样的压力体为虚压力体,用(-)来表示。
图3-3a 中的压力体是实压力体,它对曲面形成向下的压力。
图3-3b 中的压力体是虚压力体,它对曲面形成向上的浮力。
需要注意的是:以上的两个压力体给人的感觉是实压力体就是内部充满液体的压力体,虚压力体就是内部没有液体的压力体。
其实压力体的虚实与其内部是否充满液体无关,∇AB CD EF GHIJ 图3-4 压力体的合成(+)(-)(+)(-)图3-3 压力体的虚实(a)(b)图3-4中FG 的压力体是实压力体,但压力体的上半部分却没有充满液体。
再一个问题就是压力的合成,图3-4中CDE 曲面的压力体的画法对初学者来讲有一定的难度,首先我们把曲面划分为CD 和DE 两部分,先画出CD 部分的压力体,即图中的画右斜线部分,这部分压力体为虚压力体;后画出DE 部分的压力体,即图中的左斜线部分,这部分压力体为实压力体;最后将两者合成,交叉部分的压力体虚实相抵后剩下的凸出部分便是CDE 曲面的压力体,其压力体为实压力体,压力体对曲面的作用力是向下压力。