流体力学PPT幻灯片课件.ppt

流体力学 二 主讲教师 宗智孙雷 船舶与海洋工程 船舶工程学院 1 B3 1微分形式的连续性方程B3 2作用在流体微元上的力B3 3微分形式的动量方程B3 4纳维 斯托克斯 N S 方程B3 5边界条件和初始条件B3 6压强场 B3微分形式的基本方程 2 B3微分形式的基本方程 本章讨论流体力学三要素中第三要素 力 微分形式的流体力学基本方程描述空间点邻域内的物理量关系 求解这些方程可得到物理量在空间分布的细节主要内容 微分形式的连续性方程和动量方程 作用在流体微元上的体积力和表面力 重力场 应力场 压强场 边界条件和初始条件等 重点 1 不可压缩流体连续性方程 2 纳维 斯托克斯方程 3 压强的表达方式和单位 4 静止和运动流体中压强分布特征 3 B3 1微分形式的连续性方程 B3 1 1流体运动的连续性17世纪初 英 哈维 W Harvey 运用伽利略倡导的定量研究原则 测量出人的心脏每小时泵出约540磅 245Kg 的血 相当于人体重的两倍多 这么多血来自何方流向何方呢 哈维通过实验和逻辑思维否定了统治人类1400多年的陈旧观念 大胆提出从动脉到静脉的血液循环理论 虽然当时还不知道毛细血管的存在 直至45年后从发明的显微镜里首次观察到毛细血管 证实了哈维的理论 4 根据质量守恒定律 不可压缩流体流进控制体的质量应等于流出控制体的质量 称其为流体运动的连续性原理 B3 1 1流体运动的连续性 流体运动的连续性是物质质量守恒定律在流体运动中的特殊体现 血液循环理论是流体连续性原理的胜利 在科学史上有里程碑的意义 5 B3 1 2微分形式的连续性方程 如图所示 设流体流过以M x y z 为基点 以dx dy dz为边长的控制体元 在 t时间内沿x方向净流出控制体 流出质量减去流入质量 的质量为 按质量守恒定律 在 t时间内沿三个方向净流出控制体的总质量应等于控制体内减少的质量 6 B3 1 2微分形式的连续性方程 取极限后可得 利用质点导数概念 可改写为 1 2 式均为微分形式的三维流动连续性方程 7 不可压缩流动对于不可压缩流体 由于密度恒为常数 则不可压缩流体的连续性方程为 在直角坐标系中为 B3 1 2微分形式的连续性方程 在柱坐标系中为 在不同条件下连续方程有不同形式 速度散度为零意味着在空间一点邻域内流体的体积相对膨胀率恒为零 这是保证流体密度恒等于常数的运动学条件 8 可压缩流体定常流动对定常流动 可压缩流体定常流动的连续性方程为 B3 1 2微分形式的连续性方程 在直角坐标系中为 9 思考题 连续性方程适用于 连续性方程适用于 A 不可压缩流体 B 不定常流体 C 定常流体 D 任何流体 B3 1 2微分形式的连续性方程 10 例题B3 1 2 不可压缩流动连续性方程 1 已知 一不可压缩平面流动的x方向速度分量为求 y方向的速度分量v 解 由不可压缩流动连续性方程 y方向的速度分量为 c为常数 式中f x 为任何仅包含x变量的函数 讨论 本例说明对不可压缩流动 任一点的速度分量不能是任意的 而是受到不可压缩流动连续性方程的约束 若设f x 0 该流场代表位于原点的点涡流 若f x v 代表位于原点的点涡流叠加y方向速度为v的均匀流等等 他们均满足不可压缩流动条件 11 例题B3 1 2A 不可压缩流动连续性方程 2 已知 在收缩喷管流场中 设A1截面附近的a1点的轴向速度为u 10 38m s 速度梯度为 a1 点在a1点的上方5mm处 求 a1点y方向的速度分量v 解 由不可压缩流动连续性方程 在a1点v 0 在a1 点v va v 方向如图示 y 5mm 0 005m 讨论 本例说明a1点x方向正的速度梯度引起y方向负的速度梯度 两侧质点向轴心流动 12 B3 2作用在流体元上的力 B3 2 1体积力和表面力 体积力 穿越空间作用在所有流体元上的非接触力 如 重力 惯性力 电磁力等 作用在流体元上的体积力 Fb 大小一般与流体元体积 成正比 故名体积力 重力和惯性力正比于流体元的质量 又称质量力 体积力可表示为空间位置和时间的分布函数 作用在M x y z 点邻域内单位质量流体元上的体积力f为 13 B3 2 1体积力和表面力 作用在有限体积域 的流体团上的体积力合力为 作用在单位体积流体元上的体积力为 f 对被考察的流体团 体积力一般当作外力 当体积力仅为重力时 流体可称为重力流体 14 表面力 表面力为流场中假想面一侧的流体 或固体 对另一侧流体的接触力 如压强 粘性切应力等 作用在流体元上的表面力 Fs 除了与空间位置 时间有关外 还与面积元的方位有关 作用在过M x y z 点 外法线单位矢量为n的面积元A上的单位面积表面力为 B3 2 1体积力和表面力 称为表面应力 脚标n代表面积元的方位 作用在有限表面域A上的表面力合力为 15 B3 2 1体积力和表面力 16 思考题 表面应力Pn的脚标n代表面积元的方位 指该面积元的单位外法矢量 Pn的方向为 A 垂直于面积元 方向与一致 B 垂直于面积元 方向与相反 C 不一定与面积元垂直 B3 2 1体积力和表面力 17 B3 2 2重力场 重力场 在Z轴垂直向上的直角坐标系中 作用在单位质量流体之上的重力构成重力场 g为重力加速度 重力是有势力 设 简称为重力势 是单位质量流体元具有的重力势能 重力势梯度的负值即为单位质量流体元的重力 18 思考题 在重力场中设Z轴垂直向上 1 单位体积的流体元受到的体积力大小为 相应的重力势能为 2 单位质量的流体元受到的体积力大小为 相应的重力势能为 3 单位重量的流体元受到的体积力大小为 相应的重力势能为 回答组合 a g gz b 1 z c g gz A 1 a 2 c 3 b B 1 c 2 a 3 b C 1 b 2 c 3 a B3 2 2重力场 19 在静止流体中没有切向应力 只有法向应力 静止流体中的表面应力始终与作用面垂直 在静止流体中一点的法向应力在各个方向均相等 B3 2 3流体应力场 静止流体中的应力状态 称p为静压强 就是热力学中的平衡压强 负号表示流体只受压 运动的无粘性流体中也没有切向应力 应力状态与静止流体相似 20 思考题 如下图 一圆柱可绕轴心转动 左半圆浸没于水中 右半圆暴露于空气中 判断下面的说法是否正确 由于左半圆受到水的浮力产生力矩使圆柱做顺时针转动 A 这种说法是正确的 B 这种说法是错误的 B3 2 3流体应力场 所有的表面应力均垂直于圆柱面 因此都通过轴心 无力矩 21 运动流体的应力状态 B3 2 3流体应力场 在运动粘性流体中 一点的表面应力与作用面不垂直 即有法向分量又有切向分量 而且这些分量的大小与作用面的方位有关 称其为应力状态 一点的应力状态可用通过该点三个互相垂直的面积之上三组表面应力分量完全确定 如外法矢沿x轴正向的面积元dAx上一组应力分量为 pxx x法向 xy y切向 xz x切向 上式中表面应力分量的第一个脚标代表面积元的方位 即外法矢的指向 第二个脚标代表表面应力作用方向 称为应力表示约定 22 B3 2 3流体应力场 同另外两个正交面积元上的两组应力分量共九个分量构成应力矩阵 张量 可以证明九个分量中只有六个是独立的 通常约定 当法向应力与外法矢n方向一致时为正 被作用的流体元受拉伸 方向相反时为负 被作用的流体元受压缩 23 应力矩阵的常用表达式 运动的可压缩粘性流体各方面的法向压应力可以不相等 引入平均压强 并认为它也等于热力学中的平衡压强 简称为压强p B3 2 3流体应力场 把压强从法向应力中分离出来 式中 x y z是运动粘性流体偏离平均压强的附加法向应力 与流体元线应变率有关 24 B3 2 3流体应力场 应力矩阵可写成 上式右边第一项称为静压强项 第二项称为 偏应力 项 由流体运动产生 静止时为零 25 B3 3微分形式的动量方程 微分形式的动量方程 流体运动微分方程 用牛顿第二定律描述流体运动 可得在直角坐标系中微分形式的动量方程如下 上式表明 单位体积流体元上的体积力及三个方向的表面应力梯度造成了单位体积流体元的加速度 26 如下图所示 在正方体微元三组平面上x方向的表面应力梯度构成表面应力合力 B3 3微分形式的动量方程 流体运动微分方程适用于任何流体 对不同类型的流体将具有不同的形式 27 B3 4纳维 斯托克斯 N S 方程 不可压缩牛顿流体本构关系对于不可压缩牛顿粘性流体 将牛顿粘性定律从一维推广到三维 法向应力和切向应力分别与线应变率和角变形率成线性关系 Stokes假设 28 N S方程将不可压缩牛顿流体的本构关系代入直角坐标系中微分形式的动量方程可得 B3 4纳维 斯托克斯 N S 方程 29 上式称为均质不可压缩牛顿流体的纳维 斯托克斯方程 习惯上简称为N S方程 B3 4纳维 斯托克斯 N S 方程 N S方程是本课程中占主导地位的控制方程 在不同条件下 对不同流体模型可化为不同形式 N S方程加上连续性方程构成封闭的方程组 可在适当的边界条件和初始条件下求解 矢量形式 30 思考题 A 体积力 压强 粘性应力 B 体积力 压强梯度 粘性应力 C 体积力 压强梯度 粘性应力散度 B3 4纳维 斯托克斯 N S 方程 N S方程是牛顿第二定律应用于牛顿粘性流体流动中的表达式 由N S方程可看到 引起单位体积流体元加速度的作用力是 压强和粘性应力是表面力 当它们作用在流体元某一方向上处于平衡状态时不引起该方向的加速度 只有存在梯度 粘性应力在各个方向上的作用合力是粘性应力的散度 时才引起加速度 31 内流问题 出入口的速度和压强分布已知 一般由实验测得 外流问题 无穷远处的速度和压强分布已知 两种流体交界面 界面上的速度 压强和粘性切应力应连续 边界条件固体边界粘性流体 必须满足固壁面不滑移条件 或速度连续条件 无粘流体 无需满足不滑移条件 但法向速度仍应连续 B3 5边界条件与初始条件 32 两种流体交界面应满足的边界条件为 初始条件对定常流动 无初始条件 对于非定常流动应知道初始时刻的速度和压强分布 B3 5边界条件与初始条件 33 已知 牛顿流体 在重力作用下沿斜坡 倾角为 做定常层流流动 液面上方为大气压 流层深h 设图中坐标系中速度 体积力 压强分别为 解 平面流动的N S方程为 例题B3 5 1 沿斜坡的定常层流 N S方程与边界条件 求 验证是否满足N S方程及边界条件 34 例题B3 5 1 沿斜坡的定常层流 N S方程与边界条件 本例中 1 式左边 0右边 2 式左边 0右边 满足N S方程 在斜坡上 y 0 u 0在液面上 y h 压强 满足不滑移条件 满足切应力为零 y h 0 为大气压强 满足边界条件 35 B3 6压强场 压强在流体运动 流体与固体相互作用中扮演重要角色 如机翼升力 高尔夫球及汽车的尾流阻力都与压强有关 龙卷风产生强大的负压强作用 液压泵和压缩机推动流体做功是正压强作用的结果 36 B3 6 1静止重力流体中的压强分布 均质流体压强一般表达式 静止流体中无惯性力和粘性力 体积力为重力 由N S方程可得 前两式表明p与xy无关 对均质流体 常数 由第三式积分可得 上式表明静止流体中的压强沿垂直坐标为线性分布 常数c由边界条件决定 37 公式常用来表示具有自由液面的液体内的压强分布 均质液体压强公式 静止液体中的压强分布示意图 B3 6 1静止重力流体中的压强分布 设自由液面的坐标为z0 压强为p0 可得 在工程上通常用自由液面下的深度 称为淹深 h z0 z 表示一点的垂直位置 右图 则上式可改写为 上式称为匀质静止液体中的压强公式 它表明 在垂直方向 压强与淹深成线性关系 在水平方向 h 常数 压强为常数 水平面是等压强面 简称等压面 38 思考题 如下图所示的一U形管 管内有两种液体处于平衡状态 试指出图中所画断面中的等压面 A 1 1面 B 2 2面 C 3 3面 B3 6 1静止重力流体中的压强分布 判断等压面的条件是 连通的同种流体 39 例题B3 6 1 静压强分布图 已知 静止液体的自由表面上方为大气压强 求 定性的画出液体中斜面和曲面上的压强分布 解 40 B3 6 2压强计示方式与单位 压强计示方式 压强公式可作为压强计算的基础 其中为基准压强 两个基准 绝对真空 和当地大气压 三种计示方式 绝对压强 相对于绝对真空计量之值 标注为 ab 表压强 相对于当地大气压计量之值 当低于时为负 标注为 g 真空度 当表压强为负时 取其绝对值 标注为 v 约定 除特别说明外 压强均以表压强计算 41 思考题 pab pg pv分别表示绝对压强 表压强和真空度 pa表示大气压强 试判断下列表达式哪个是对的 A B C B3 6 2压强计示方式与单位 42 压强单位 B3 6 2压强计示方式与单位 国际单位制 SI 帕斯卡 Pa 1Pa 1N m2 1kPa 103N m2 1MPa 106N m2 物理单位制 cgs 毫米汞柱 mmHg 单位制 43 例题B3 6 2 单管与U形管测压计 已知 一封闭容器中充满密度为 的液体 求 用单管和U形管测压计测量内壁面上任一点A的压强 解 在A点处壁面上开一小孔 接液柱式测压计 1 若pA 0 接单管测压计 如图 液体在压力作用下上升至h高度 液面上为大气 按下式 pA pa gh 绝 gh 表 h pA g m a 称h为A点的测压管高度 还可以表示为能量形式 gh pA gh表示重力势能 pA 称为压强势能 b 44 例题B3 6 2 单管与U形管测压计 2 若pA 0 接U型管测压计 如图 U形管内有一段重液体 如汞 密度为 设其液面差为 h A点离左支管液面距离为h1 pA gh1 mg h 0 U形管测压计也适用于测量气体压强 1 1 由等压面1 1列压强平衡方程 pA gh1 mg h 0 用被测液体的测压管高度表示 U形管液面差折合成测压管高度 45 例题B3 6 2A U形管压差计 已知 二个封闭容器A B中分别充满密度为 的流体 气体或液体 求 用U形管测量A B两点的压强差 p pA pA 解 将U形管两支接到A B两点 U形管内有一段重液体 密度为 m 液体差为 h 取0 0线为基准面 A B的位置为zA zB pA g zA h pB gzB mg h p pA pB g zB zA m g h 用被测流体的测压管高度表示 由等压面1 1列压强平衡方程 U形管液面压强差 位置差 46 例题B3 6 2B 压强计示与单位 已知 设水泵吸水管的绝对压强为p 8N cm2 大气压强为pa 1 013 105Pa 试用 国际单位制表示其绝对压强 表压强 真空压强和真空度 解 pv pg 2 13 104Pa 21 3kPa pab 8 104N m2 Pa 80kPa 或表示为 p 80kPa ab p pg pab pa 8 104 1 013 105 Pa 2 13 104Pa 21 3kPa 或 p 2 13 104Pa 1 013 105Pa 21 v 绝对压强 表压强 真空压强 真空度 p 21 3kPa v 47 B3 6 3运动流体中的压强分布 运动流体中 影响压强分布的因素除体积力外 还有惯性力和粘性力等 例一 圆柱绕流 惯性力和粘性力的影响 设流体对圆柱作定常平面绕流 圆柱表面的压强分布在无粘性流体和粘性流体中有不同的概念 设压强系数为 式中p为圆柱面上压强 p0 v0为无穷远处压强和速度 48 图b为粘性流体绕流时 Re 105 由于边界层分离在圆柱后部形成尾流区 见动画 前后压强分布不对称 作用在圆柱上的压强合力不为零 形成压差阻力 图a为无粘性流体绕流的压强系数分布图 为前后对称分布 B D点是最大正压强点 驻点 C E点是最大负压强点 作用在圆柱上的压强合力为零 达朗贝尔佯谬 B3 6 3运动流体中的压强分布 49 机翼上下表面压强分布示意图 下表面以正压强为主 上表面以负压强为主 压强合力形成升力 NACA标准翼型 2412 在攻角分别为7 4度和2 8度时的压强系数分布图 可见主要以上表面负压强为主 例二 机翼绕流 B3 6 3运动流体中的压强分布 50 在风洞里沿轿车中剖面测量的压强系数分布图 可见除迎风面为正压强外 其他部位大多是负压强 例三 汽车绕流 B3 6 3运动流体中的压强分布 51 普通型轿车在车速很高时将产生升力 使轮胎与地面咬合力减小 造成驱动效率降低 稳定性差 为了克服这些缺点 可采取如下改进措施 A 增加轮胎的表面粗糙度 B 改变车身形线 使高速时升力减少 C 在轿车车身上安装产生负升力的辅助装置 B3 6 3运动流体中的压强分布 思考题 答案 b c 目前流行的楔形车身在高速运动时不仅不产生升力 反而产生向下的压力 另外在轿车后部安装倒置的翼形片 产生的升力向下 可抵消车身的升力 52