微分形式的基本方程流体力学幻灯片.ppt

《微分形式的基本方程流体力学幻灯片.ppt》由会员分享，可在线阅读，更多相关《微分形式的基本方程流体力学幻灯片.ppt（39页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、微分形式的基本方程流体力学第1页，共39页，编辑于2022年，星期六B3.1.1 B3.1.1 流体运动的连续性流体运动的连续性(2-2)(2-2)1717世纪哈维：血液循环理论世纪哈维：血液循环理论 解剖发现解剖发现：从心脏到动脉末端血液单向：从心脏到动脉末端血液单向 流动，从静脉末端到心脏也流动，从静脉末端到心脏也 是单向流动是单向流动 定量测量定量测量：每小时流出心脏血液：每小时流出心脏血液245kg 大胆预言大胆预言：从动脉到静脉再回心脏：从动脉到静脉再回心脏 45年后发现年后发现：毛细血管的存在毛细血管的存在血液循环理论血液循环理论流体连续性原理流体连续性原理的胜利的胜利血液循环图血

2、液循环图第2页，共39页，编辑于2022年，星期六B3.1.2 B3.1.2 微分形式的连续性方程微分形式的连续性方程 x，y，z方向净流出质量为方向净流出质量为因密度变化引起的质量减少为因密度变化引起的质量减少为由由质量守恒定律质量守恒定律单位时间单位体积内单位时间单位体积内边长为边长为 ,的长方体控制体元，的长方体控制体元，内内x方向净流出的质量方向净流出的质量B3.1.2 微分形式的连续性方程微分形式的连续性方程(2-1)第3页，共39页，编辑于2022年，星期六B3.1.2 B3.1.2 微分形式的连续性方程微分形式的连续性方程(2-2)(2-2)用场量公式并运用质点导数概念，微分形式

3、用场量公式并运用质点导数概念，微分形式连续性方程连续性方程为为或改写为：或改写为：左左边边代代表表一一点点邻邻域域内内流流体体体体积积的的相相对对膨膨胀胀速速率率，右右边边代代表表密密度度相相对减少率。连续性方程适用于任何同种流体。对减少率。连续性方程适用于任何同种流体。不可压缩流体不可压缩流体连续性方程连续性方程第4页，共39页，编辑于2022年，星期六 例例B3.1.2B3.1.2 不可压缩流动连续性方程不可压缩流动连续性方程 已知：已知：不可压缩流体平面流动不可压缩流体平面流动（C为常数）为常数）求：求：v 解：解：由不可压缩流动连续性方程的二维形式由不可压缩流动连续性方程的二维形式可得

4、可得（B3.1.113.1.11）当当f(x)=0，表示位于原点的点涡流动；，表示位于原点的点涡流动；当当f(x)=U，表示点涡流叠加，表示点涡流叠加y方向速度为方向速度为U的均流；的均流；讨论：讨论：本例说明对不可压缩流动，任一点的各速度分量不能是任意的本例说明对不可压缩流动，任一点的各速度分量不能是任意的,而是受到（而是受到（B3.1.113.1.11）式制约的。）式制约的。第5页，共39页，编辑于2022年，星期六B3.2 B3.2 作用在流体元上的力作用在流体元上的力B3.2.1 B3.2.1 体积力和表面力体积力和表面力1.1.体积力体积力长程力长程力穿越空间作用穿越空间作用到流体元

5、上到流体元上万有引力万有引力电磁力电磁力惯性力惯性力与流体元体与流体元体积成正比积成正比体积力体积力单位质量流体上的体积力单位质量流体上的体积力 单位体积流体上的体积力单位体积流体上的体积力 B3.2.1 体积力和表面力体积力和表面力(2-1)第6页，共39页，编辑于2022年，星期六B3.2.1 B3.2.1 体积力和表面力体积力和表面力(2-2)(2-2)2.2.表面力表面力短程力短程力通过接触面通过接触面作用作用压强压强粘性切应力粘性切应力与表面面积与表面面积和方位有关和方位有关表面力表面力表面力定义：作用在单位平面面积元上的短程力。表面力定义：作用在单位平面面积元上的短程力。n面积元外

6、法线单位矢面积元外法线单位矢n面积元内法线单位矢面积元内法线单位矢（注意：（注意：和和 不一定与不一定与 垂直）垂直）第7页，共39页，编辑于2022年，星期六B3.2.2 B3.2.2 重力场重力场在直角坐标系的重力场中在直角坐标系的重力场中称为重力势，代表单位质量流体具有的重力势能称为重力势，代表单位质量流体具有的重力势能B3.2.2 重力场重力场第8页，共39页，编辑于2022年，星期六B3.2.3 B3.2.3 应力场应力场1.1.运动粘性流体中的应力状态运动粘性流体中的应力状态一点的表面一点的表面应力应力用过该点三个坐标用过该点三个坐标面上三组表面力分面上三组表面力分量唯一确定量唯一

7、确定应力状态应力状态与作用力的大小、方向、作用面方位有关与作用力的大小、方向、作用面方位有关上的应力分量为上的应力分量为上的应力分量为上的应力分量为上的应力分量为上的应力分量为B3.2.3 应力场应力场(4-1)应力矩阵应力矩阵第9页，共39页，编辑于2022年，星期六作用在任意方位作用在任意方位面元上的面元上的表面应力表面应力表面应力的分量式表面应力的分量式B3.2.3 应力场应力场(4-2)作用在作用在外法矢沿外法矢沿x轴向的面积元轴向的面积元dAx上三个应力分量如图示上三个应力分量如图示第10页，共39页，编辑于2022年，星期六B3.2.3 B3.2.3 应力场应力场(4-3)(4-3

8、)2.2.静止流体中的应力状态静止流体中的应力状态静止流体的应力状态静止流体的应力状态结论：静止流体中一点的应力状态只用一个标量静压强结论：静止流体中一点的应力状态只用一个标量静压强p p表示表示.只有法向应力只有法向应力无切应力无切应力第11页，共39页，编辑于2022年，星期六B3.2.3 B3.2.3 应力场应力场(4-4)(4-4)3.3.应力的常用表达式应力的常用表达式运动粘性流体中的运动粘性流体中的(平均平均)压强压强在法向应力中把压强分离出来在法向应力中把压强分离出来为附加法向应力分量（与流体元线应变率有关）为附加法向应力分量（与流体元线应变率有关）压强矩阵压强矩阵 偏应力矩阵偏

9、应力矩阵 应力矩阵表示为应力矩阵表示为 第12页，共39页，编辑于2022年，星期六 例例B3.2.3B3.2.3 平面线性剪切流中的应力状态平面线性剪切流中的应力状态 已知：已知：平面线性剪切流平面线性剪切流求：求：应力状态应力状态 解解:附加法向应力附加法向应力切应力切应力讨论：讨论：附加法向应力与该方向的线应变率有关，平面线性剪切流中任附加法向应力与该方向的线应变率有关，平面线性剪切流中任一点处在一点处在x、y方向的线应变率均为零，因此相应的附加法向应力方向的线应变率均为零，因此相应的附加法向应力也均为零，也均为零，x,y方向的法向应力均等于平衡压强；粘性切应力则方向的法向应力均等于平衡

10、压强；粘性切应力则在全流场保持常数。在全流场保持常数。法向应力法向应力（k为常数）为常数）第13页，共39页，编辑于2022年，星期六 例例B3.2.3AB3.2.3A 刚体旋转流动刚体旋转流动:纯旋转纯旋转(2-1)(2-1)已知：已知：二维不可压缩平面流场为二维不可压缩平面流场为求：求：试分析该流场中的试分析该流场中的应力状态应力状态（k为常数）为常数）解：解：附加法向应力附加法向应力第14页，共39页，编辑于2022年，星期六流体中任一点的法向流体中任一点的法向应力为应力为 切向应力为切向应力为讨论：讨论：（1 1）线应变率处处为零，附加法向应力为零，全流场）线应变率处处为零，附加法向应

11、力为零，全流场 的法向应力均等于平衡压强。的法向应力均等于平衡压强。（2 2）角变形率也处处为零，全流场的粘性切应力为）角变形率也处处为零，全流场的粘性切应力为零，流体和刚体一样作定轴旋转运动。零，流体和刚体一样作定轴旋转运动。例例B3.2.3AB3.2.3A 刚体旋转流动刚体旋转流动:纯旋转纯旋转(2-2)(2-2)第15页，共39页，编辑于2022年，星期六B3.3 B3.3 微分形式的动量方程微分形式的动量方程按牛顿第二定律，长方体流体元的按牛顿第二定律，长方体流体元的运动方程运动方程为为 各面元上各面元上 x 方向表面应力的分量方向表面应力的分量如图示。如图示。B3.3 微分形式的动量

12、方程微分形式的动量方程(2-1)表面力合力表面力合力 dFsx 由应力梯度造成由应力梯度造成第16页，共39页，编辑于2022年，星期六x方向的体积力分量为方向的体积力分量为 将将dFsx和和dFbx代入运动方程，并利用代入运动方程，并利用 和质点导数和质点导数概念，可化为概念，可化为 同理可得同理可得 上式称为上式称为粘性流体运动一般微分方程粘性流体运动一般微分方程，适用于任何流体。，适用于任何流体。B3.3 微分形式的动量方程微分形式的动量方程(2-2)第17页，共39页，编辑于2022年，星期六B3.4 B3.4 纳维斯托克斯方程纳维斯托克斯方程 斯托克斯假设：斯托克斯假设：1.1.将牛

13、顿粘性定律从一维推广到三维；将牛顿粘性定律从一维推广到三维；2.2.流体各向同性；流体各向同性；3.3.静止时法向应力等于静压强。静止时法向应力等于静压强。均代入粘性流体运动一般微分方程均代入粘性流体运动一般微分方程对牛顿流体（对牛顿流体（常数）常数）B3.4 纳维斯托克斯方程纳维斯托克斯方程(4-1)不可压缩条件（不可压缩条件（常数）常数）第18页，共39页，编辑于2022年，星期六B3.4 B3.4 纳维斯托克斯方程纳维斯托克斯方程(4-2)(4-2)可得均质不可压缩牛顿流体的可得均质不可压缩牛顿流体的纳维纳维-斯托克斯方程斯托克斯方程(NS S方程方程)NS S方程的适用条件是：方程的适

14、用条件是：常数常数常数常数，第19页，共39页，编辑于2022年，星期六B3.4 B3.4 纳维斯托克斯方程纳维斯托克斯方程(4-3)(4-3)NS S方程的矢量式为方程的矢量式为N NS S方程方程的意义和求解：的意义和求解：物理意义是：物理意义是：惯性力与体积力、压力、粘性力平衡惯性力与体积力、压力、粘性力平衡 u、v、w、p，方程组是封闭的；，方程组是封闭的；加上连续性方程加上连续性方程 ，四个方程求解四个未知数，四个方程求解四个未知数 在边界条件较简单时可求解析解在边界条件较简单时可求解析解;在边界条件较复杂时在边界条件较复杂时可求数值解可求数值解;对不同的流动专题可作不同程度的简化（

15、见专题篇）。对不同的流动专题可作不同程度的简化（见专题篇）。第20页，共39页，编辑于2022年，星期六B3.4 B3.4 纳维斯托克斯方程纳维斯托克斯方程(4-4)(4-4)N NS S方程方程平衡方程平衡方程相对平衡方程相对平衡方程欧拉方程欧拉方程惯性力惯性力体积力体积力粘性力粘性力 压力压力00第21页，共39页，编辑于2022年，星期六B3.5 B3.5 边界条件与初始条件边界条件与初始条件 1.1.常见边界条件常见边界条件(1)(1)固体壁面固体壁面粘性流体：不滑移条件粘性流体：不滑移条件(图图a)a)无粘性流体：法向速度连续无粘性流体：法向速度连续(图图b)b)v=v固 vn=v

16、n固 (2)(2)外流无穷远条件外流无穷远条件v=v，p=p B3.5 边界条件与初始条件边界条件与初始条件(2-1)第22页，共39页，编辑于2022年，星期六(3)(3)内流出入口条件内流出入口条件v=vin(out)，p=p in(out)(4)(4)自由面条件自由面条件2.2.初始条件初始条件定常流时无初始条件定常流时无初始条件不定常流时给出某时刻的参数值：不定常流时给出某时刻的参数值：v(t0),p(t0),(t0)等等B3.5 边界条件与初始条件边界条件与初始条件(2-2)第23页，共39页，编辑于2022年，星期六 例例B3.5.1AB3.5.1A 沿斜坡的重力粘性层流沿斜坡的重

17、力粘性层流(3-1)(3-1)已知：已知：不可压牛顿流体在重力作用下沿斜坡不可压牛顿流体在重力作用下沿斜坡()作定常层流流动，流层作定常层流流动，流层深深h，自，自由面上为大气压（由面上为大气压（p0 0）。）。(a(a)求：求：(1)速度分布速度分布 (2)压强分布压强分布 (3)切应力分布切应力分布 (4)流量流量 解解：在图示坐标系中连续性方程在图示坐标系中连续性方程和和NS S方程方程为为(b(b)(c(c)第24页，共39页，编辑于2022年，星期六 例例B3.5.1AB3.5.1A 沿斜坡的重力粘性层流沿斜坡的重力粘性层流(3-2)(3-2)因因v0，由（，由（a a）式）式由（由

18、（c c）式）式由边界条件由边界条件(1):y=h,p=0,C(x)=,压强分布为压强分布为且且，由，由(b)式式积分两次积分两次第25页，共39页，编辑于2022年，星期六流量流量 速度分布为速度分布为讨论：讨论：压强和切应力为线性分布，速度分布为压强和切应力为线性分布，速度分布为y的二次函数，流量为的二次函数，流量为h 的三次函数。的三次函数。切应力分布切应力分布 例例B3.5.1AB3.5.1A 沿斜坡的重力粘性层流沿斜坡的重力粘性层流(3-3)(3-3)由边界条件由边界条件(2):y=0,u=0 可得可得 C2=0由边界条件由边界条件(3):y=b,第26页，共39页，编辑于2022年

19、，星期六B3.6B3.6压强场压强场 由由N NS S方程方程粘性流动粘性流动绝对平衡绝对平衡相对平衡相对平衡无粘性流动无粘性流动B3.6 压强场压强场 第27页，共39页，编辑于2022年，星期六B3.6.1 B3.6.1 静止重力流体中的压强分布静止重力流体中的压强分布 均质静止流体均质静止流体 =常数，常数，uvw0在重力场中在重力场中上式说明：上式说明：z方向压强梯度由单位体积流体的重力决定。方向压强梯度由单位体积流体的重力决定。积分可得积分可得B3.6.1 静止重力流体中的压强分布静止重力流体中的压强分布(3-1)1.1.压强分布一般表达式压强分布一般表达式由由N-SN-S方程可得方

20、程可得第28页，共39页，编辑于2022年，星期六B3.6.1 B3.6.1 静止重力流体中的压强分布静止重力流体中的压强分布(3-2)(3-2)2.2.具有自由液面的重力液体具有自由液面的重力液体 压强公式压强公式为自由面上的压强，为自由面上的压强，h为为淹深淹深(1)(1)在垂直方向压强与淹深成线性关系在垂直方向压强与淹深成线性关系 (2)(2)在水平方向压强保持常数在水平方向压强保持常数 第29页，共39页，编辑于2022年，星期六B3.6.1 B3.6.1 静止重力流体中的压强分布静止重力流体中的压强分布(3-3)(3-3)3.3.等压面等压面在连通的同种流体中的等压强面称为在连通的同

21、种流体中的等压强面称为等压面等压面。在静止重力流体中的等压面为水平面在静止重力流体中的等压面为水平面h常数常数右图中右图中3 33 3 为等压面为等压面非等非等压面压面1 11 1 为不连通液体为不连通液体2 22 2 为不同液体为不同液体第30页，共39页，编辑于2022年，星期六 例例B3.6.1B3.6.1 静压强分布图静压强分布图第31页，共39页，编辑于2022年，星期六B3.6.2 压强计示方式与单位压强计示方式与单位1.压强计示方式压强计示方式习惯上取习惯上取压强基准压强基准真空度真空度 完全真空完全真空绝对压强绝对压强表压强表压强大气压强大气压强B3.6.2 压强计算方法与单位

22、压强计算方法与单位(2-1)由压强公式由压强公式p0提供压强基准提供压强基准第32页，共39页，编辑于2022年，星期六B3.6.2 B3.6.2 压强计算方法与单位压强计算方法与单位(2-2)(2-2)2.压强单位压强单位标准大气压标准大气压atm(标准国际大气模型标准国际大气模型)液柱高：液柱高：国际单位制（国际单位制（SI）：帕斯卡）：帕斯卡Pa 毫米汞柱毫米汞柱mmHg（血压计）（血压计）米水柱米水柱mH2O（水头高）（水头高）测压管高度测压管高度 h=pA/g第33页，共39页，编辑于2022年，星期六 例例B3.6.2B3.6.2 单管测压计（单管测压计（2 21 1）已知：已知：

23、图示密封容器中液体图示密封容器中液体(),),在在A点接上单管测压计点接上单管测压计求：求：与测压管高度与测压管高度h 的关系的关系解：解：(表压强表压强)h为被测点的淹深，称为测压管高度为被测点的淹深，称为测压管高度.讨论：讨论：液面在压强液面在压强 推动下上升至推动下上升至 h 高度，压强势能转化为重力势能。高度，压强势能转化为重力势能。压强势能压强势能重力势能重力势能第34页，共39页，编辑于2022年，星期六 例例B3.6.2B3.6.2 U形管测压计（形管测压计（2 22 2）解解：沿沿U 形管右支液面取等压面，列平衡方程形管右支液面取等压面，列平衡方程已知已知：图示封闭容器中为水图

24、示封闭容器中为水,U形管水银测压计形管水银测压计中中h=10cm求：求：(,(,表压强表压强 真空压强真空压强 绝对压强）绝对压强）第35页，共39页，编辑于2022年，星期六 例例B3.6.2AB3.6.2A U形管差压计形管差压计解：解：沿沿U 形管左支液面取等压面形管左支液面取等压面11已知已知：图示盛满水封闭容器高差图示盛满水封闭容器高差 ,U形管水银测压计中液面差形管水银测压计中液面差h=10cm求：求：(,(,表压强表压强 绝对压强）绝对压强）第36页，共39页，编辑于2022年，星期六B3.6.3 B3.6.3 运动流场中的压强分布运动流场中的压强分布 压强系数压强系数1.1.惯

25、性力对压强分布的影响惯性力对压强分布的影响 p 0，v 0为参考值，对外流场取为参考值，对外流场取p，v B3.6.3 运动流场中的压强分布运动流场中的压强分布(3-1)文丘里管流动文丘里管流动 第37页，共39页，编辑于2022年，星期六B3.6.3 B3.6.3 运动流场中的压强分布运动流场中的压强分布 无粘流场无粘流场压强分布压强分布 静止流场压强分布静止流场压强分布2.粘性力对压强分布粘性力对压强分布的影响的影响B3.6.3 运动流场中的压强分布运动流场中的压强分布(3-2)第38页，共39页，编辑于2022年，星期六B3.6.3 B3.6.3 运动流场中的压强分布运动流场中的压强分布 汽车与飞机绕流汽车与飞机绕流B3.6.3 运动流场中的压强分布运动流场中的压强分布(3-3)3.复杂物面的压强分布复杂物面的压强分布第39页，共39页，编辑于2022年，星期六