商品详情
书名:高等流体力学(一)
定价:36.0
ISBN:9787030563200
作者:无
版次:1
出版时间:2018-01
内容提要:
本书主要围绕黏性对流体流动影响的物理本质及数学描述展开讨论。全书共六章,主要介绍黏性流体运动的基本方程及层流运动的流动特性。
第1章介绍黏性流体流动的基本概念。第2章讲述黏性流体力学的基本方程和黏性流体运动的基本性质。第3章举例介绍圆管内、平行平板间和同轴旋转圆筒间这三种能够获得解析解的黏性流体的层流流动。第4章给出在小Re数下,黏性流体绕流小圆球蠕流流动的斯托克斯解和奥森解,以及轴承润滑理论。第5、6章为边界层理论的基础,主要介绍普朗特对边界层流动的论述、层流绕流平板的布拉修斯相似性解、绕楔形体流动的弗克纳-斯肯解,以及绕流边界层的卡门动量积分关系式解法。湍流运动及数值模拟将在后续出版的《高等流体力学(二)》中介绍。
目录:
目录
前言
第1章 绪论 1
1.1 概述 1
1.1.1 理想流体与黏性流体 1
1.1.2 可压缩流体与不可压缩流体 3
1.2 应力张量 5
1.2.1 理想流体中的应力 5
1.2.2 黏性流体中的应力 5
1.2.3 张量知识简介 6
1.2.4 应力张量 11
1.3 变形率张量 13
1.3.1 速度分解定理 13
1.3.2 流体微团的运动分析 14
1.3.3 变形率张量 17
1.4 本构方程(广义牛顿内摩擦定律) 18
1.4.1 牛顿内摩擦定律 18
1.4.2 斯托克斯的三点假设 19
1.4.3 广义牛顿内摩擦定律 19
第2章 黏性流体力学的基本方程 23
2.1 连续性方程 23
2.1.1 连续性方程的导出 23
2.1.2 随体导数 24
2.2 运动方程 26
2.2.1 用应力张量[τ]表示的运动方程 26
2.2.2 用变形率张量[ε]表示的运动方程 27
2.3 能量方程 28
2.3.1 用总能表示的能量方程 28
2.3.2 用内能表示的能量方程 30
2.3.3 耗散函数 30
2.3.4 用温度T表示的能量方程 31
2.4 状态方程 32
2.4.1 完全气体状态方程 32
2.4.2 其他热力状态参数间的关系 33
2.5 黏性流体运动方程组的封闭性和定解条件 34
2.5.1 方程组的封闭性 34
2.5.2 定解条件 34
2.6 黏性流体运动的基本性质 35
2.6.1 黏性流体运动的有旋性 35
2.6.2 黏性流体中涡旋的扩散性 36
2.6.3 黏性流体运动能量的耗散性 39
第3章 特殊条件下的黏性流体运动方程解 40
3.1 圆管内层流 40
3.1.1 圆管内层流流动的速度分布和流量表达式 40
3.1.2 圆管内层流流动的沿程阻力公式 42
3.1.3 入口段与充分发展的管内流动 42
3.2 平板间的层流 43
3.2.1 平行平板间层流流动的微分方程和速度分布 43
3.2.2 泊肃叶流动与库埃特剪切流 45
3.3 同轴旋转圆筒间黏性流体的定常流动 49
第4章 黏性流体绕固体物面的缓慢流动 52
4.1 黏性流体绕小圆球的蠕流流动 52
4.1.1 斯托克斯阻力系数 52
4.1.2 奥森解及其修正 56
4.2 颗粒在静止流体中的自由沉降 59
4.3 流体润滑 61
第5章 边界层层流流动及其相似性解 64
5.1 边界层流动的基本概念与基本特征 64
5.2 边界层的各种厚度 65
5.2.1 边界层的名义厚度δ 65
5.2.2 排挤厚度(位移厚度)δ* 66
5.2.3 动量损失厚度θ 67
5.2.4 δ,δ*,θ的图解 67
5.3 边界层微分方程 68
5.4 绕曲面流动和边界层的分离 71
5.4.1 绕曲面流动边界层的分离 71
5.4.2 边界层分离的原因和后果 74
5.4.3 卡门涡街 75
5.5 层流边界层的相似性方程 76
5.5.1 边界层相似的概念 76
5.5.2 相似性方程 79
5.5.3 存在相似性解的物面条件 81
5.6 绕平板层流流动边界层方程的布拉修斯解 88
5.6.1 布拉修斯解 88
5.6.2 布拉修斯解的应用 92
5.7 绕楔形体流动的弗克纳-斯肯解 96
5.7.1 弗克纳-斯肯方程的解 96
5.7.2 弗克纳-斯肯解的应用 100
第6章 层流边界层积分关系式解法 104
6.1 卡门边界层动量积分关系式 104
6.2 单参数速度剖面和相容边界条件 106
6.2.1 单参数速度剖面 106
6.2.2 相容边界条件 106
6.3 绕曲面流动的边界层动量积分关系式解法 108
6.3.1 卡门-波尔豪森单参数方法 108
6.3.2 霍斯汀的改进 113
6.3.3 斯韦茨解法 116
6.4 绕平板流动的边界层动量积分关系式解法 119
参考文献 123
附录 常用正交坐标系中基本量和基本方程的表达式 124
在线试读:
第1章 绪论
1.1 概述
1.1.1 理想流体与黏性流体
我们在“工程流体力学”等课程中已学过,流体有三大特性:易流动性、黏性和可压缩性。其中,易流动性,即流体是一种受任何微小剪切力作用都能发生连续变形的物质,是流体区别于固体的*本质的属性。黏性是流体的第二大属性,任何流体都具有黏性。但是,黏性的存在给流体运动的数学描述和处理带来很大困难,因此在流体的黏性作用反映不出的场合,用不考虑黏性的理想流体代替黏性流体,从而可简化数学模型及其求解。黏性作用由黏性切应力τ反映。根据牛顿内摩擦定律
(1-1)
式中,τ为黏性切应力;μ为流体的动力黏度;为流体运动的速度梯度。可见,黏性切应力由两个因素决定,一为流体的动力黏度,二为流体运动的速度梯度。因此,当流体处于静止状态或以相同的速度流动(即速度梯度为零)时,流体的黏性作用反映不出,此时就可用理想流体代替黏性流体。另外,对于一些动力黏度较小的流体,如水和空气,当其运动的速度梯度较小时,由于黏性的作用较弱,可将其视为理想流体处理,再对黏性的影响进行修正,使问题由繁变简。在流体力学的发展史上,不考虑黏性的理想流体力学理论*先得到发展,其中*著名的就是不可压缩理想流体的欧拉(Euler)运动方程
(1-2)①
式中,ρ为流体密度;V为流体的速度矢量;f为流体所受的质量力;p为流体所受压强。
在重力场中,对定常流动的欧拉运动方程沿流线积分,就可得到著名的伯努利(Bernoulli)方程
(1-3)
式中,V为流体运动速度;g为重力加速度;z为高度。伯努利方程的物理意义为:不可压缩理想流体在重力场中做定常流动时,沿流线单位质量流体的动能、位势能和压强势能之和等于常数。伯努利方程的应用解决了许多流体流动的问题。但是应该看到,理想流体力学理论是建立在不考虑流体黏性基础上的。由于不考虑黏性,就出现了以下两个与事实不符合的简化假定:①流体不承受流层与流层之间的切向力(即τ=0),只承受法向力(压力);②在固体壁面上存在滑动。而这两个简化假定使理想流体力学理论在求解流动阻力时与实验结果不符。
考虑如图1-1所示的理想流体平行流绕流半径为R的圆柱体的二维流动。按理想流体理论,在极坐标中,其径向和切向速度分量分别为
(1-4)
式中,为来流速度。在圆柱体表面上,即径向速度为零,但切向速度除在前、后驻点处为零外,其余各点不为零,且按正弦规律分布,即壁面上存在切向滑移速度,这与黏性流体绕流圆柱体时在壁面上无滑移的实际流动情况不一致。对于二维不可压缩理想流体,可按伯努利方程求圆柱体表面上的压强分布
(1-5)
式中,为来流压强。上式所示压强分布也可用量纲为一的量(简称为量纲一量)——压强系数Cp表示
(1-6)
图1-1 理想流体平行流绕流半径为R的圆柱体的二维流动
图1-2给出了按式(1-6)计算得到的理想流体的压强系数与实际测量值的比较。由图1-2可见,理想流体的压强系数前后对称,表明理想流体平行流绕圆柱体流动时并未受到任何阻力,这与达朗贝尔对理想流体进行一系列严谨研究后得出的结论一致,即“一个任意固体在无穷大静止的理想流体中做匀速直线运动时,在其运动方向不承受作用力”。这就是历史上著名的达朗贝尔佯谬。产生达朗贝尔佯谬的原因在于理想流体忽略了流体的黏性作用,使得按理想流体力学理论计算流动阻力时与实验结果不符。由图1-2可知,在圆柱的前缘(θ=0,θ=2π)附近,理想流体的计算结果与实际符合较好,但在此外的其他各处,两者差别较大。对于理想流体,圆柱前后的流动是完全对称的,所以按理想流体力学理论求得的阻力为零。但是,实测的压强分布前后不对称,圆柱后部的实测压强系数低于前部对应点处的值,且随Re数不同而不同,这样就使圆柱受向后作用的力,即压差阻力。理想流体力学理论计算不出这些阻力,从而产生了与实验不一致的结果。
图1-2 理想流体的压强系数与实际测量值的比较
因此,如果要研究阻力,就要研究壁面影响下的流体流动,必须考虑黏性。此时,理想流体的两个简化假定必须抛弃,必须承认流体既承受法向力(压力),也承受流层之间的切向力;必须承认流体在固体壁面上无滑移,即壁面处切向速度处处为零。这样,描写流体运动的不再是不可压缩理想流体的欧拉方程,而是不可压缩黏性流体的纳维-斯托克斯(Navier-Stokes)方程。
(1-7)
式中,μ为流体的动力黏度;Δ为拉普拉斯算符。与式(1-2)欧拉方程相比,增加了黏性项μΔV。
1.1.2 可压缩流体与不可压缩流体
可压缩性是流体的第三大特性。实际中的任何流体,无论是液体还是气体,都是可以压缩的,不可压缩流体并不存在。同样为了研究方便,人们提出了不可压缩流体的概念。当流体受压体积不减小,受热体积不增加,其密度保持为常数,这种流体称为不可压缩流体。液体的压缩性很小,随着压强和温度的变化,液体的密度仅有微小的变化。在绝大多数情况下,可以忽略压缩性的影响,认为液体的密度是常数。于是通常把液体看成不可压缩流体。
气体的压缩性很大。由热力学知,当温度不变时,完全气体的体积与压强成反比,压强增加一倍,体积减小为原来的1/2;当压强不变时,温度升高1℃,体积增加0℃时体积的1/273。所以,通常把气体看成可压缩流体,即它的密度不能作为常数,而是随着压强和温度的变化而变化。
在工程实际中,是否要考虑流体的压缩性,要视具体情况而定。若流体密度的相对变化较大,需考虑流体的压缩性;若流体密度的相对变化较小,则可不考虑。在流体流动时,这个密度相对变化值可用流速v来反映。
物质的压缩特性可用弹性模数反映,它定义为
(1-8)
式中,Δp为压强的变化量,Δρ为密度的变化量。
考察同一水平面内的A1,A2两点,A1点的压强为p1,流速为0,A2点的压强为p2,流速为v,则根据伯努利方程,有
即
可见与相当,代入式(1-8),可得
因为声速
(1-9)
马赫数
(1-10)
所以
不可压缩流体,也就是若将作为是否考虑压缩性的阈值,则可将作为是否考虑压缩性的分界线。对常温下的空气,C=334m/s,则由M2=0.1,求得v≈100m/s。所以,对空气,当流速小于100m/s时,可不考虑压缩性的影响;对液体,弹性模量大大增加,所对应的声速大大增加。所以,可不考虑压缩性的流速范围也大大增加。
1.2 应力张量
1.2.1 理想流体中的应力
在工程流体力学中已经阐明,作用在静止流体单位面积上的表面力(应力)永远沿着作用面的内法线方向,而且其大小与作用面所处的方向无关,只与截面的位置有关,即一点的静压力各方向相等。
对于运动的理想流体,由于忽略黏性,所以没有切向力,只有法向力(压力),因而单位面积上的表面力(压力)垂直于作用面,而且各个方向相等。所以,对理想流体,无需特别关注应力(压力)的方向,即使考虑,也仅需考虑作用面方向即可。
1.2.2 黏性流体中的应力
对于黏性流体,由于存在黏性,除法向应力外,还有切向应力存在,因此单位面积上的表面力(应力)就不一定垂直于作用面,而且各个方向的大小也不一定相等。
1.应力的表示
在运动流体中选取一微元曲面dA,若微元曲面dA是闭曲面的一部分,则取外法线方向为dA的正方向;若dA所在的曲面不封闭,则可以规定某一法线方向为正。在dA正方向的面上作用的表面应力用τn表示,它表示的是微元曲面dA正方向所指流体作用在dA上的应力。同样,用表示dA负方向面上作用的表面力,即表示微元曲面dA负方向所指流体作用在dA上的应力,如图1-3所示。可见
(1-11)
图1-3 作用在微元面dA上的应力
2.应力分量
若应力作用面垂直于某坐标轴,则应力可以分解成三个分量。如图1-4(b)所示,应力的作用面的法线方向沿z轴方向,则可将分解为三个分量:一个为垂直于作用面的法向应力,另外两个是与作用面平行的切向应力的分量和。可见,应力分量的第*个下标表示作用面的法线方向,第二个下标表示与应力分量平行的坐标轴,即应力投影方向。同理,若作用面法线方向沿x轴,如图1-4(c)所示,则作用在该面上的应力用表示,且可分解成;若作用面法线方向沿y轴,如图1-4(a)所示,则作用在该面上的应力用表示,且可分解成。我们把称为黏性流体的法向应力,简称为法向应力或正应力,把称为黏性流体的切向应力,简称切应力。
图1-4 作用在各坐标轴面上的应力
设直角坐标系(x,y,z)三个坐标轴上的单位矢量为e1、e2、e3,则上述应力分解关系可表示成
(1-12)
式(1-12)中有9个分量,可以排成矩阵,如式(1-13)所示,称为二阶应力张量[τ]。
(1-13)
这样,我们就引入了一个新的量,称为张量。它和我们以前学过的矢量、标量一样,也是为了表示一定的物理量而引入的。例如,对时间t、长度l、面积A等量,仅需表示大小,无需表示方向,用标量表示。对速度v、加速度a、力f等量,既要表示大小,又要表示方向,用矢量表示。现在对应力,既要表示大小、方向,还要表示作用面的方向,矢量也不够用了,所以引入二阶张量[τ]。所以,黏性流体中一点的应力既有大小、方向,又与作用面方向有关,用二阶应力张量表示。
为统一起见,将既表示大小、方向,又表示作用面方向的量称为二阶张量,将矢量(表示大小、方向的量)称为一阶张量,将标量(仅表示大小的量)称为零阶张量。
1.2.3 张量知识简介
1.张量的基本特征
张量是一种用以表示物理量的数学工具,而物理量本身是不随坐标系的选取而变化
定价:36.0
ISBN:9787030563200
作者:无
版次:1
出版时间:2018-01
内容提要:
本书主要围绕黏性对流体流动影响的物理本质及数学描述展开讨论。全书共六章,主要介绍黏性流体运动的基本方程及层流运动的流动特性。
第1章介绍黏性流体流动的基本概念。第2章讲述黏性流体力学的基本方程和黏性流体运动的基本性质。第3章举例介绍圆管内、平行平板间和同轴旋转圆筒间这三种能够获得解析解的黏性流体的层流流动。第4章给出在小Re数下,黏性流体绕流小圆球蠕流流动的斯托克斯解和奥森解,以及轴承润滑理论。第5、6章为边界层理论的基础,主要介绍普朗特对边界层流动的论述、层流绕流平板的布拉修斯相似性解、绕楔形体流动的弗克纳-斯肯解,以及绕流边界层的卡门动量积分关系式解法。湍流运动及数值模拟将在后续出版的《高等流体力学(二)》中介绍。
目录:
目录
前言
第1章 绪论 1
1.1 概述 1
1.1.1 理想流体与黏性流体 1
1.1.2 可压缩流体与不可压缩流体 3
1.2 应力张量 5
1.2.1 理想流体中的应力 5
1.2.2 黏性流体中的应力 5
1.2.3 张量知识简介 6
1.2.4 应力张量 11
1.3 变形率张量 13
1.3.1 速度分解定理 13
1.3.2 流体微团的运动分析 14
1.3.3 变形率张量 17
1.4 本构方程(广义牛顿内摩擦定律) 18
1.4.1 牛顿内摩擦定律 18
1.4.2 斯托克斯的三点假设 19
1.4.3 广义牛顿内摩擦定律 19
第2章 黏性流体力学的基本方程 23
2.1 连续性方程 23
2.1.1 连续性方程的导出 23
2.1.2 随体导数 24
2.2 运动方程 26
2.2.1 用应力张量[τ]表示的运动方程 26
2.2.2 用变形率张量[ε]表示的运动方程 27
2.3 能量方程 28
2.3.1 用总能表示的能量方程 28
2.3.2 用内能表示的能量方程 30
2.3.3 耗散函数 30
2.3.4 用温度T表示的能量方程 31
2.4 状态方程 32
2.4.1 完全气体状态方程 32
2.4.2 其他热力状态参数间的关系 33
2.5 黏性流体运动方程组的封闭性和定解条件 34
2.5.1 方程组的封闭性 34
2.5.2 定解条件 34
2.6 黏性流体运动的基本性质 35
2.6.1 黏性流体运动的有旋性 35
2.6.2 黏性流体中涡旋的扩散性 36
2.6.3 黏性流体运动能量的耗散性 39
第3章 特殊条件下的黏性流体运动方程解 40
3.1 圆管内层流 40
3.1.1 圆管内层流流动的速度分布和流量表达式 40
3.1.2 圆管内层流流动的沿程阻力公式 42
3.1.3 入口段与充分发展的管内流动 42
3.2 平板间的层流 43
3.2.1 平行平板间层流流动的微分方程和速度分布 43
3.2.2 泊肃叶流动与库埃特剪切流 45
3.3 同轴旋转圆筒间黏性流体的定常流动 49
第4章 黏性流体绕固体物面的缓慢流动 52
4.1 黏性流体绕小圆球的蠕流流动 52
4.1.1 斯托克斯阻力系数 52
4.1.2 奥森解及其修正 56
4.2 颗粒在静止流体中的自由沉降 59
4.3 流体润滑 61
第5章 边界层层流流动及其相似性解 64
5.1 边界层流动的基本概念与基本特征 64
5.2 边界层的各种厚度 65
5.2.1 边界层的名义厚度δ 65
5.2.2 排挤厚度(位移厚度)δ* 66
5.2.3 动量损失厚度θ 67
5.2.4 δ,δ*,θ的图解 67
5.3 边界层微分方程 68
5.4 绕曲面流动和边界层的分离 71
5.4.1 绕曲面流动边界层的分离 71
5.4.2 边界层分离的原因和后果 74
5.4.3 卡门涡街 75
5.5 层流边界层的相似性方程 76
5.5.1 边界层相似的概念 76
5.5.2 相似性方程 79
5.5.3 存在相似性解的物面条件 81
5.6 绕平板层流流动边界层方程的布拉修斯解 88
5.6.1 布拉修斯解 88
5.6.2 布拉修斯解的应用 92
5.7 绕楔形体流动的弗克纳-斯肯解 96
5.7.1 弗克纳-斯肯方程的解 96
5.7.2 弗克纳-斯肯解的应用 100
第6章 层流边界层积分关系式解法 104
6.1 卡门边界层动量积分关系式 104
6.2 单参数速度剖面和相容边界条件 106
6.2.1 单参数速度剖面 106
6.2.2 相容边界条件 106
6.3 绕曲面流动的边界层动量积分关系式解法 108
6.3.1 卡门-波尔豪森单参数方法 108
6.3.2 霍斯汀的改进 113
6.3.3 斯韦茨解法 116
6.4 绕平板流动的边界层动量积分关系式解法 119
参考文献 123
附录 常用正交坐标系中基本量和基本方程的表达式 124
在线试读:
第1章 绪论
1.1 概述
1.1.1 理想流体与黏性流体
我们在“工程流体力学”等课程中已学过,流体有三大特性:易流动性、黏性和可压缩性。其中,易流动性,即流体是一种受任何微小剪切力作用都能发生连续变形的物质,是流体区别于固体的*本质的属性。黏性是流体的第二大属性,任何流体都具有黏性。但是,黏性的存在给流体运动的数学描述和处理带来很大困难,因此在流体的黏性作用反映不出的场合,用不考虑黏性的理想流体代替黏性流体,从而可简化数学模型及其求解。黏性作用由黏性切应力τ反映。根据牛顿内摩擦定律
(1-1)
式中,τ为黏性切应力;μ为流体的动力黏度;为流体运动的速度梯度。可见,黏性切应力由两个因素决定,一为流体的动力黏度,二为流体运动的速度梯度。因此,当流体处于静止状态或以相同的速度流动(即速度梯度为零)时,流体的黏性作用反映不出,此时就可用理想流体代替黏性流体。另外,对于一些动力黏度较小的流体,如水和空气,当其运动的速度梯度较小时,由于黏性的作用较弱,可将其视为理想流体处理,再对黏性的影响进行修正,使问题由繁变简。在流体力学的发展史上,不考虑黏性的理想流体力学理论*先得到发展,其中*著名的就是不可压缩理想流体的欧拉(Euler)运动方程
(1-2)①
式中,ρ为流体密度;V为流体的速度矢量;f为流体所受的质量力;p为流体所受压强。
在重力场中,对定常流动的欧拉运动方程沿流线积分,就可得到著名的伯努利(Bernoulli)方程
(1-3)
式中,V为流体运动速度;g为重力加速度;z为高度。伯努利方程的物理意义为:不可压缩理想流体在重力场中做定常流动时,沿流线单位质量流体的动能、位势能和压强势能之和等于常数。伯努利方程的应用解决了许多流体流动的问题。但是应该看到,理想流体力学理论是建立在不考虑流体黏性基础上的。由于不考虑黏性,就出现了以下两个与事实不符合的简化假定:①流体不承受流层与流层之间的切向力(即τ=0),只承受法向力(压力);②在固体壁面上存在滑动。而这两个简化假定使理想流体力学理论在求解流动阻力时与实验结果不符。
考虑如图1-1所示的理想流体平行流绕流半径为R的圆柱体的二维流动。按理想流体理论,在极坐标中,其径向和切向速度分量分别为
(1-4)
式中,为来流速度。在圆柱体表面上,即径向速度为零,但切向速度除在前、后驻点处为零外,其余各点不为零,且按正弦规律分布,即壁面上存在切向滑移速度,这与黏性流体绕流圆柱体时在壁面上无滑移的实际流动情况不一致。对于二维不可压缩理想流体,可按伯努利方程求圆柱体表面上的压强分布
(1-5)
式中,为来流压强。上式所示压强分布也可用量纲为一的量(简称为量纲一量)——压强系数Cp表示
(1-6)
图1-1 理想流体平行流绕流半径为R的圆柱体的二维流动
图1-2给出了按式(1-6)计算得到的理想流体的压强系数与实际测量值的比较。由图1-2可见,理想流体的压强系数前后对称,表明理想流体平行流绕圆柱体流动时并未受到任何阻力,这与达朗贝尔对理想流体进行一系列严谨研究后得出的结论一致,即“一个任意固体在无穷大静止的理想流体中做匀速直线运动时,在其运动方向不承受作用力”。这就是历史上著名的达朗贝尔佯谬。产生达朗贝尔佯谬的原因在于理想流体忽略了流体的黏性作用,使得按理想流体力学理论计算流动阻力时与实验结果不符。由图1-2可知,在圆柱的前缘(θ=0,θ=2π)附近,理想流体的计算结果与实际符合较好,但在此外的其他各处,两者差别较大。对于理想流体,圆柱前后的流动是完全对称的,所以按理想流体力学理论求得的阻力为零。但是,实测的压强分布前后不对称,圆柱后部的实测压强系数低于前部对应点处的值,且随Re数不同而不同,这样就使圆柱受向后作用的力,即压差阻力。理想流体力学理论计算不出这些阻力,从而产生了与实验不一致的结果。
图1-2 理想流体的压强系数与实际测量值的比较
因此,如果要研究阻力,就要研究壁面影响下的流体流动,必须考虑黏性。此时,理想流体的两个简化假定必须抛弃,必须承认流体既承受法向力(压力),也承受流层之间的切向力;必须承认流体在固体壁面上无滑移,即壁面处切向速度处处为零。这样,描写流体运动的不再是不可压缩理想流体的欧拉方程,而是不可压缩黏性流体的纳维-斯托克斯(Navier-Stokes)方程。
(1-7)
式中,μ为流体的动力黏度;Δ为拉普拉斯算符。与式(1-2)欧拉方程相比,增加了黏性项μΔV。
1.1.2 可压缩流体与不可压缩流体
可压缩性是流体的第三大特性。实际中的任何流体,无论是液体还是气体,都是可以压缩的,不可压缩流体并不存在。同样为了研究方便,人们提出了不可压缩流体的概念。当流体受压体积不减小,受热体积不增加,其密度保持为常数,这种流体称为不可压缩流体。液体的压缩性很小,随着压强和温度的变化,液体的密度仅有微小的变化。在绝大多数情况下,可以忽略压缩性的影响,认为液体的密度是常数。于是通常把液体看成不可压缩流体。
气体的压缩性很大。由热力学知,当温度不变时,完全气体的体积与压强成反比,压强增加一倍,体积减小为原来的1/2;当压强不变时,温度升高1℃,体积增加0℃时体积的1/273。所以,通常把气体看成可压缩流体,即它的密度不能作为常数,而是随着压强和温度的变化而变化。
在工程实际中,是否要考虑流体的压缩性,要视具体情况而定。若流体密度的相对变化较大,需考虑流体的压缩性;若流体密度的相对变化较小,则可不考虑。在流体流动时,这个密度相对变化值可用流速v来反映。
物质的压缩特性可用弹性模数反映,它定义为
(1-8)
式中,Δp为压强的变化量,Δρ为密度的变化量。
考察同一水平面内的A1,A2两点,A1点的压强为p1,流速为0,A2点的压强为p2,流速为v,则根据伯努利方程,有
即
可见与相当,代入式(1-8),可得
因为声速
(1-9)
马赫数
(1-10)
所以
不可压缩流体,也就是若将作为是否考虑压缩性的阈值,则可将作为是否考虑压缩性的分界线。对常温下的空气,C=334m/s,则由M2=0.1,求得v≈100m/s。所以,对空气,当流速小于100m/s时,可不考虑压缩性的影响;对液体,弹性模量大大增加,所对应的声速大大增加。所以,可不考虑压缩性的流速范围也大大增加。
1.2 应力张量
1.2.1 理想流体中的应力
在工程流体力学中已经阐明,作用在静止流体单位面积上的表面力(应力)永远沿着作用面的内法线方向,而且其大小与作用面所处的方向无关,只与截面的位置有关,即一点的静压力各方向相等。
对于运动的理想流体,由于忽略黏性,所以没有切向力,只有法向力(压力),因而单位面积上的表面力(压力)垂直于作用面,而且各个方向相等。所以,对理想流体,无需特别关注应力(压力)的方向,即使考虑,也仅需考虑作用面方向即可。
1.2.2 黏性流体中的应力
对于黏性流体,由于存在黏性,除法向应力外,还有切向应力存在,因此单位面积上的表面力(应力)就不一定垂直于作用面,而且各个方向的大小也不一定相等。
1.应力的表示
在运动流体中选取一微元曲面dA,若微元曲面dA是闭曲面的一部分,则取外法线方向为dA的正方向;若dA所在的曲面不封闭,则可以规定某一法线方向为正。在dA正方向的面上作用的表面应力用τn表示,它表示的是微元曲面dA正方向所指流体作用在dA上的应力。同样,用表示dA负方向面上作用的表面力,即表示微元曲面dA负方向所指流体作用在dA上的应力,如图1-3所示。可见
(1-11)
图1-3 作用在微元面dA上的应力
2.应力分量
若应力作用面垂直于某坐标轴,则应力可以分解成三个分量。如图1-4(b)所示,应力的作用面的法线方向沿z轴方向,则可将分解为三个分量:一个为垂直于作用面的法向应力,另外两个是与作用面平行的切向应力的分量和。可见,应力分量的第*个下标表示作用面的法线方向,第二个下标表示与应力分量平行的坐标轴,即应力投影方向。同理,若作用面法线方向沿x轴,如图1-4(c)所示,则作用在该面上的应力用表示,且可分解成;若作用面法线方向沿y轴,如图1-4(a)所示,则作用在该面上的应力用表示,且可分解成。我们把称为黏性流体的法向应力,简称为法向应力或正应力,把称为黏性流体的切向应力,简称切应力。
图1-4 作用在各坐标轴面上的应力
设直角坐标系(x,y,z)三个坐标轴上的单位矢量为e1、e2、e3,则上述应力分解关系可表示成
(1-12)
式(1-12)中有9个分量,可以排成矩阵,如式(1-13)所示,称为二阶应力张量[τ]。
(1-13)
这样,我们就引入了一个新的量,称为张量。它和我们以前学过的矢量、标量一样,也是为了表示一定的物理量而引入的。例如,对时间t、长度l、面积A等量,仅需表示大小,无需表示方向,用标量表示。对速度v、加速度a、力f等量,既要表示大小,又要表示方向,用矢量表示。现在对应力,既要表示大小、方向,还要表示作用面的方向,矢量也不够用了,所以引入二阶张量[τ]。所以,黏性流体中一点的应力既有大小、方向,又与作用面方向有关,用二阶应力张量表示。
为统一起见,将既表示大小、方向,又表示作用面方向的量称为二阶张量,将矢量(表示大小、方向的量)称为一阶张量,将标量(仅表示大小的量)称为零阶张量。
1.2.3 张量知识简介
1.张量的基本特征
张量是一种用以表示物理量的数学工具,而物理量本身是不随坐标系的选取而变化