流体力学 丁祖荣 孔珑

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