《流体力学课件丁祖荣孔珑.ppt》由会员分享，可在线阅读，更多相关《流体力学课件丁祖荣孔珑.ppt（70页珍藏版）》请在三一办公上搜索。
1、工程流体力学,李宝宽,流体流动现象,自然界固有的流动现象,自然的流动现象,黄河壶口,济南趵突泉,波滑,人类的利用,绪论 人类：海洋里50亿年，空气：700万年。流体运动与流体力学三个问题（1）高尔夫球，飞得远应表面光滑还是粗慥？40-200米（2）汽车，阻力来自前部还是后部？Cd=1-0.2,未来型，0.137。（3）机翼：升力来自下部还是上部？数百吨飞机，数万吨巨轮：机理不一样流体力学：研究流体运动规律，流体之间、流固之间相互作用、流动过程中动量、质量和能量的传输等。,流体力学与科学纪元前：大禹治水、古罗马城市供水系统、阿基米德浮力定律18世纪：理想流体力学：L.Euler,D.Bernou
2、lli,J.DAlembert,J.Lagrange,P.Laplace 实验流体力学：G.Hagen,J.Poiseuille,A.Chezy19世纪：模型实验法则：W.Froude 量纲分析法：L.Reyleigh 两种流态：O.Reynolds 粘性流体的运动方程：C.Navier,C.stokes现代流体力学：以普朗特（L.Prandtl）边界层理论为标志。中国：周培源，钱学森,汉灵帝,流体力学方法推动其他学科领域的发展流体力学边界层理论导致应用数学中渐进展开匹配法的形成 1904年，边界层理论流体力学孤立波理论成为新学科光通信的基石 1844年，J.Russell:孤立波，20世纪6
3、0年代，光通信从流体力学劳伦兹方程发现混沌 1995年，E.Lorenz,求解大气对流，相同输入数据，第四位小数四舍五入 输出结果无法预测。混沌理论是现代非线性科学基础，解释生命和社会现象。,流体力学与工程技术：三分之二以上自然与工程涉及流体流动。冶金,流体力学研究水平最能体现国家的综合实力,Electric Arc Furnace,Continuous Casting,Argon Oxygen Decarburization,RH Degassing,Magnet,钢精炼与连铸工艺流程,专业应用例,转炉氧枪钢包精炼轧钢加热炉,流体力学研究方方法：理论分析实验分析数值计算上述方法的结合,单位制
4、：用于定量地表示物理量B的参考量B（i），类别与相应的物理量相同，大小人为规定。物理量的大小由相应的单位表示：B=k(i)B(i)k(i)为数值，单位与物理量之间仅存在数量的差别。流体力学问题中四个基本量：质量，长度，时间和温度。辅助单位：平面角：弧度导出单位：力，压强，密度，黏度，能量，功率。,流体及其物理性质,物质三态：固体、液体和气体。前者不易变形，后两者易变形。高温和高压时还有等离子态和超固态。气体在温度和压力影响下体积变化大，液体则变化小，且有自由表面。物态差别取决于内部分子结构、分子热运动和分子之间的相互作用力。分子之间距离较近 10-9 m：吸引力，分子之间距离非常近：排斥力平衡
5、距离：10-10 m,阿夫加德罗常数：1 mol 物质包含6.022136710-23 分子。分子之间平均距离为分子线尺度的10倍。3.310-7 cm.气体 液体和固体三态可以互相转化。主要是因温度和压力的变化导致。,流体力学将气体和液体的宏观力学行为用连续介质模型描述。连续介质假定：物质结构：流体由分子组成，分子之间有空隙。流体力学 研究宏观运动而非分子本身。连续介质模型认为物质连续无间隙分布于物质所占有的整个空间，流体宏观物理量是空间及时间的连续函数。连续介质模型方法不针对单个粒子，而考虑大量粒子的平均运动和统计特性，速度、密度和温度等。涉及到平均尺度问题。某一单元体内密度值随单元尺度的
6、变化，如图其中：l 分子间平均距离，(气体10-7 cm,液体 10-8 cm)流体宏观物理量统计平均的尺度 l a L,流体质点：微观上充分大，宏观上充分小。质点具有的物理量充分均匀的，它是大量粒子的统计平均值。1）无线尺度，平移、无变形。2）不做热运动。3）质点中心统计平均值为质点的物理量。连续介质是由这些连续分布的质点组成。流体的密度、比重和比容是流体的最基本物理量，流体密度是流体中某空间点上单位体积的平均质量，流体的比重是该流体的重量与同体积水在4摄氏度时的重量之比。流体的可压缩性和热膨胀性在外力作用下，流体体积或密度可以改变的性质，称之为流体的可压缩性；在温度改变时，流体体积或密度可
7、以改变的性质，称之为流体的热膨胀性。这部分在工程热力学和传热学中将详细介绍。,流体的输运性质流体由非平衡态转向平衡态时物理量的传递性质，统称为流体的输运性质。流体的输运，包括动量输运、能量输运和质量输运。动量输运-粘滞现象1687年，牛顿平行平板实验：或写成微分形式 称为牛顿切应力公式。的单位Pa.s或1N.s/m2,亦即1kg/(m.s).运动黏度：的单位m2/s.牛顿流体：符合牛顿切应力公式的流体，空气、水等。非牛顿流体：不符合牛顿切应力公式的流体，如奶油，果酱，血液等。,U,y,x,粘度或动力粘度：表示产生单位切变率所需要的切应力大小。与温度关系：问题：1550钢液的动力粘度和运动粘度，
8、并与25水的相应粘度比较。热能输运-热传导现象1822年傅里叶实验得到傅里叶定律导热系数单位w/(m.K)质量输运-扩散现象1855年，菲克实验 扩散系数单位m2/s.,表面张力：一种液体与气体、另一种不相溶液体或固体接触时，在交界面表层内表现出来的张力。,习题：,流体力学三大要素,流体（密度、粘度、压缩性。）运动（拉格朗日，欧拉）力（质量力、表面力和惯性力）,描述流体运动的方法例（1）交通流量(2)小河、船和桥。1标记质点法流体物理量=流体物理量（标记，时间）2流体物理量=流体物理量（空间点，时间）拉格朗日坐标标记(a,b,c)-拉格朗日坐标随体坐标方法：常常采用初始时刻t=t0的位置坐标(
9、x0,y0,z0)来表示。但拉格朗日坐标不随时间变化。欧拉坐标固定于空间的坐标系的一组坐标（q1,q2,q3）表示不同时刻流体质点的空间位置。流体质点连续布满流体空间。故流体质点与空间点也与一组坐标（q1,q2,q3）对应。这是欧拉坐标。可以用直角坐标、柱坐标和球坐标等表示。不同时刻可有不同坐标。拉格朗日坐标与欧拉坐标不同但相关。,若拉格朗日坐标为(a,b,c)，则某一物理量的拉格朗日数学表达式为：,若位置以r表示，则,质点加速度,欧拉描述/空间描述流体物理量=流体物理量（空间点，时间）欧拉描述是场描述。某一时刻在空间点上必有一流体质点占据，故欧拉描述的物理量也就是占据该空间点的流体质点的物理
10、量。,拉格朗日描述与欧拉描述之间的关系拉格朗日描述流体质点 物理量=物理量（随体坐标，时间）欧拉描述空间点 物理量=物理量（空间坐标，时间）,随体导数流体质点物理量随时间的变化率称为随体导数，或物质导数，质点导数。对拉格朗日描述，随体导数就是偏导数。例如：,对欧拉描述，因x,y,z也是时间的函数。因此，物理量的随体导数是,流体质点的加速度,速度场,标量矢量张量,流量平均速度,维数定常流与非定常流,流动的几何描述迹线:流体质点的迹线就是流体质点的运动的轨迹,也就是该流体质点在不同时刻运动位置的连线.如:即为参数形式的迹线方程,消去t,得这是一条平面双曲线。对于欧拉描述，迹线微分方程为：流线:用来
11、描述流场中各点流动方向的曲线.是矢量线,即在该点的切线方向与该点在该时刻的速度矢量方向一致.在某一时刻流线上取一微小弧元素,由流线定义,应有:直角坐标系中,它是流线:某一时刻的迹线:某一质点的,习题：收缩喷管流动：迁移加速度,Helmholtz速度分解定理：,流场中一点邻域的相对运动分析速度分解:流场中 的邻域,设M点速度为v,由泰勒级数展开得：,写成分解式：,写成矩阵式,第一个矩阵S是对称的，第二个矩阵是反对称的。反对称的矩阵A的九个分量中只有三个独立分量，即,因而有：,上式等价于,它是相对速度的一部分。因而：,Helmholtz速度分解定理：某一点邻点速度等于相对该点的平移速度和相对该点旋
12、转速度及相对该点的变形速度之和。同刚体比较，流体多一个变形速度：包括相对伸长率，体积膨胀率和角变形率。,流体的变形线应变速率,2.体膨胀率或速度散度,3 角变形速率,流体的旋转,层流:阻力与速度成正比,湍流:阻力与速度的二次方成正比.雷诺数:,雷诺实验,内流:管道,发动机,炉窑等.外流:汽车外形,飞机机翼等,无旋与有旋,涡量是否为零,涡线为任意一点的切线与涡量方向的瞬时矢量线,流场中涡量是否处处为零.,作用在流体单元上的力,质量力：作用于体积内每一质量微元上的力称为质量力，或称体积力，如重力、电磁力等。，是空间点和时间的函数。作用于整个流体体积上的质量力为。表面力：作用于流体表面上的力称为表面
13、力，大气压力和摩擦力等。，表示以n为法向的面积元上的表面力，而作用于整个表面S上的表面力为:是空间点和时间的函数，同时还与面积的取向有关。如果Pn与n的方向不一致，可以分解成n方向和切向方向。,流体中任一点的应力、应力张量：推导Pn与n的关系,做四面体，设,ABC的法向单位矢量为n:则：,或简写成：,设ABC面积为S,于是MBC,MCA,MAB的面积分别用Sx,sy,Sz表示为,，MABC四面体体积为,四面体受有质量力，表面力和惯性力，按照达朗贝尔受力平衡原理，三种力平衡。因质量力和惯性力与四面体体积成正比，故为三阶小量，而表面力与面积成正比，故为二阶小量。四面体收缩时，忽略三阶小量，则有：,
14、故得：,该式在直角坐标系中的投影为,写成张量形式,写成矩阵形式,其中：,应力张量与应变率张量之间的关系本构方程,1848年stokes假定（1）应力张量是应变率张量的线性函数。（2）流体是各向同性的。（3）流体静止时，应变率为零，流体中应力就是静压强。,应用牛顿切应力公式，可得：,假定：,微分形式的基本方程,连续性方程,设t时刻A点流体密度为(x,y,z,t),速度为v(x,y,z,t),其分量是u、v、w.考虑六面体每个元上质量的流入或流出，由于每个面只与一个坐标轴垂直，故每个面上只有一个速度分量使相应的质量流入或流出该六面体。先计算与x轴相垂直的两个面上的质量流量。,流体运动的连续性原理,
15、在ABCD面上，,时间将有,的流体质量进入六面体，而在EFGH面上，,时间内则将有,的流体质量流出该六面体。这样，通过这两个面，,时间内就有,时间内，通过DAEH及BFGC这两个面净流出的流体质量为,这样，,时间内通过六面体的全部六个面净流出的质量就为,与此同时，此六面体内的质量将发生变化，因t时刻，在六面体内的流体质量为,故经过,时间，即在t+t,时刻，六面体内的质量将是,在,时间内，六面体的质量增加了,或减少了,于是，依据质量守恒定律，,时间内，六面体内所减少的质量一定与同一时间内从六面体中流出的质量相等：,这是直角坐标系中流体运动的微分形式的连续方程。上式第一项代表单位时间内，单位体积的
16、质量增量，第二、三、四项代表单位时间内、单位体积内质量的净流出。利用散度公式：,于是,再由,及,也可写为,可压缩与不可压缩流体,对于不可压缩流体,pxx,pxy,pxz,微分形式的运动方程,动量变化包括两部分 1 微元体内的动量变化率2 经微元体控制面迁移的动量,t 时刻控制体内动量,t+t 时刻控制体内动量,于是单位时间内净变化,经微元体控制面迁移的动量：六个面其中垂直于 x 轴两个面流进和流出的动量为,流进：,流出：,净流出：,其他面：,全部控制面净流出：,系统动量变化率,利用连续性方程，上式中小括号内之值为零。于是微元体内流动系统变化率为,计算质量力和表面力质量力,表面力,相对面上的应力
17、为,六个面上的力在x，y，z轴上的投影分别是,上式是矢量,的三个投影。,依动量平衡定律：,运动方程：,惯性力=质量力+表面力,分量形式为：,Naiver-stokes方程,最终形式：,n,n,n,n,A5,A4,A2,A3,CS=A2+A3+A4+A5,雷诺守恒方程,连续性方程,不可压缩流体流动n-s方程,矢量式为：,初始条件边界条件,一类边界:给定变量的值二类边界：给定变量的导数值。三类边界：部分变量值，部分变量导数值。,流体的平衡 V=0平衡是相对于某一坐标系静止不动。特点是没有切应力法向为n的面积元，作用在此面积元上的应力,单位：Pa=N/m2.,平衡时，体积力和表面力相等,亦即,A,B
18、,两个封闭容器A、B分别充满密度为density的流体，用U形管测量A，B两点的压差，推倒压差与液位差的关系。,zA,h,zB,D,流体压强全微分式流体压强在邻点的空间增量,等压面 dp=0,习题两道,静压强，静水中物体所受的压强。,绝对压力与表压的关系真空，真空度,0,p,pa,绝对压,绝对压,表压,表压（-）,真空度,单位：1 pa=1N/m2=1kg/m.s2,1kpa=103N/m2,1Mpa=106 N/m2,1 atm=1.013X105 pa=760 mm Hg=10.33 mH20,积分形式的连续性方程,雷诺输运公式,“总”的物理量系统的广延量,是质量，动量或能量等。,连续性方
19、程,不可压缩流体,对于流管内,n1,n2,v1,v2,多个进出口情况,例题,伯努利积分,忽略Navier-Stokes方程的粘性作用，则有：,可写成：,假定体积力有势，并定义体积力势,流体为正压，即密度为常数或仅为压强的函数，即,则,于是为：,伯努利积分,流动为定常，用流线微元点乘各项后有,伯努利方程及其应用伯努利，1738，瑞士。动能与压强势能相互转换。沿流线的伯努利方程,将牛顿第二定律应用于控制体内的流体元，沿流线切线方向整理后因为,则导出,此式为一维欧拉方程，使用下述关系将方程沿流线积分。两边乘以ds,得：,沿流线积分,此式为欧拉方程的积分式，适合于可压、无粘不定常运动。对于不可压定常流动，则可简化为,此式为伯努利方程，三项分别表示单位质量流体具有的动能、位置势能和压强势能。即总机械能守恒。,伯努利方程使用的限制条件（1）无粘性流体，（2）不可压流体（3）定常流（4）沿流线。,加入能量损失就可适应粘性流体。,将伯努利方程三项机械能在有效截面A上按质量流量积分，总机械能沿流束仍保持守恒，即,以截面平均流速V代替不均匀的速度分布，引入动能修正因子。有,
