1、1北京工业大学市政学科部马长明 高等流体(水)力学讲稿一、课程名称一、课程名称:高等流体力学高等流体力学水利水电工程水利水电工程高等水力学高等水力学给排水工程(土木工程)给排水工程(土木工程)学什么学什么?二、教材:二、教材:1、高等流体力学?高等流体力学?天津大学天津大学新世纪研究生适用教材新世纪研究生适用教材相对于本科相对于本科“水力学水力学”或或“流体力学流体力学”,在相关问题上进行更深入的理论分析,在相关问题上进行更深入的理论分析和论述,以满足现代水力工程对流体力学的要求,有助于提高理论修养,深入理解现代和论述,以满足现代水力工程对流体力学的要求,有助于提高理论修养,深入理解现代流体力
2、学的内容。是水力工程以及学科各硕士专业的学位课。流体力学的内容。是水力工程以及学科各硕士专业的学位课。2、高等水力学?环境流体力学导论、高等水力学?环境流体力学导论余常昭编著余常昭编著1980年开始作为清华大学环境系和水利系研究生课的教材,年开始作为清华大学环境系和水利系研究生课的教材,1991年清华大学年清华大学出版社出版。利用流体力学理论以研究、解决环境问题出版社出版。利用流体力学理论以研究、解决环境问题污染源的扩散与弥散。污染源的扩散与弥散。3、环境水力学、环境水力学董志勇董志勇 编著编著课程简介课程简介2北京工业大学市政学科部马长明 高等流体(水)力学讲稿课程简介课程简介三、内容三、内
3、容高等流体力学(共高等流体力学(共6章)章)1、流体力学的基本概念;、流体力学的基本概念;2、流体运动的基本方程;、流体运动的基本方程;3、势流运动;、势流运动; 4、粘性流体运动;、粘性流体运动;5、紊流运动;、紊流运动;6、涡旋运动、涡旋运动附录附录 张量与场论基础知识张量与场论基础知识环境流体力学导论(共环境流体力学导论(共8章)章)1、流体运动的基本概念与基本方程;、流体运动的基本概念与基本方程;2、紊流基础(一);、紊流基础(一);3、紊流基础(二);、紊流基础(二);4、扩散理论;、扩散理论;5、剪切流中的离散;、剪切流中的离散;6、射流、羽流和浮射流;、射流、羽流和浮射流;7、分
4、层流;、分层流;8、地下水中的弥散、地下水中的弥散3北京工业大学市政学科部马长明 高等流体(水)力学讲稿课程简介课程简介三、内容三、内容环境流体力学环境流体力学董志勇(共董志勇(共8章)章)1、绪论;、绪论;2、迁移扩散理论;、迁移扩散理论;3、剪切流离散;、剪切流离散;4、射流、羽流和浮射流;、射流、羽流和浮射流;5、水质模型;、水质模型;6、地下水污染模型;、地下水污染模型;7、分层流;、分层流;8、生态水力学引论。、生态水力学引论。4北京工业大学市政学科部马长明 高等流体(水)力学讲稿课程简介课程简介讲授内容讲授内容第一讲第一讲 数学基础知识、流体力学的基本概念数学基础知识、流体力学的基
5、本概念第二讲第二讲 流体运动基本方程与小雷诺数下流体运动基本方程与小雷诺数下N-S方程的基本解方程的基本解第三讲第三讲 大雷诺数问题、势流理论与边界层的解析解大雷诺数问题、势流理论与边界层的解析解第四讲第四讲 紊流运动基础紊流运动基础第五讲第五讲 扩散理论基础扩散理论基础第六讲第六讲 剪切流中的离散剪切流中的离散第七讲第七讲 地下水中的弥撒地下水中的弥撒第八讲第八讲 射流、羽流和浮射流射流、羽流和浮射流5北京工业大学市政学科部马长明 高等流体(水)力学讲稿课程简介课程简介四、英文参考教材:四、英文参考教材:1、环境流体力学(、环境流体力学(Environmental Fluid Mechani
6、cs- Theories and Applications)2、高等环境流体力学(、高等环境流体力学(MIT开放教案开放教案MIT OcwWeb)Advanced Fluid Dynamics of the Environment, Fall 2002 五、教、学与评价探讨五、教、学与评价探讨课程特点:课程特点:1)要求数学知识多;方程、公式多,推演论证繁琐;解题)要求数学知识多;方程、公式多,推演论证繁琐;解题难度大。难度大。 2)学时少()学时少(32),所留自学时间也少,而教学内容多。),所留自学时间也少,而教学内容多。6北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数
7、学基础知识一、正交曲线坐标系一、正交曲线坐标系1、直角坐标系、柱坐标系与球坐标系、直角坐标系、柱坐标系与球坐标系1)坐标线与坐标面)坐标线与坐标面2)坐标系间的转化)坐标系间的转化2、一般正交坐标系的表示、一般正交坐标系的表示1)正交坐标系与直角坐标系的表示)正交坐标系与直角坐标系的表示2)正交坐标系下弧元素表示)正交坐标系下弧元素表示拉梅(拉梅(Lame)系数;单位矢量;弧元素;面积元素;体积元素)系数;单位矢量;弧元素;面积元素;体积元素3 , 2 , 1 ),(3 , 2 , 1 ),(321jqqqxxizyxqqjjii7北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知
8、识数学基础知识(1)拉梅系数()拉梅系数(Lame Coefficient)(2)单位矢量)单位矢量(3)弧元素)弧元素(4)面积元素)面积元素(5)体积元素)体积元素)3 , 2 , 1( 321iqxqxqxHiiii)3 , 2 , 1( 1iqrHeiii)3 , 2 , 1( idqHdsiii212131313232321 ; ;dqdqHHddqdqHHddqdqHHd 321321dqdqdqHHHdV 8北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识3)对正交曲线坐系单位向量的导数)对正交曲线坐系单位向量的导数(详细证明可参考(详细证明可参考“流
9、体力学流体力学”,周光坰等,高等教育出版社,周光坰等,高等教育出版社,2001))( 1jieqHHqejijiji 1 1kkikjjijiieqHHeqHHqe9北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识二、梯度、散度、旋度二、梯度、散度、旋度1、定义、定义梯度(梯度(gradient):标量场中一点沿其等值面法线方向的变化率;):标量场中一点沿其等值面法线方向的变化率;散度(散度(divergence):向量场中围绕点):向量场中围绕点p的封闭曲面的通量与该封闭的封闭曲面的通量与该封闭曲面所围体积之比的极限值;曲面所围体积之比的极限值;旋度(旋度(cur
10、l or rotation):向量场中围绕一点的封闭曲线的环量与该):向量场中围绕一点的封闭曲线的环量与该封闭曲线所张的面积比值的极限;封闭曲线所张的面积比值的极限;nngrad VdsnaadivSV0lim Sl daarotS0lim 10北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识2、直角坐标下的计算、直角坐标下的计算梯度:梯度:散度:散度:旋度:旋度:kzjyixgrad zayaxaadivzyx zyxaaazyxkjiarot 11北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识3、正交曲线坐标下的计算、正交曲线坐标下
11、的计算梯度:梯度:散度:散度:旋度:旋度:333222111111 eqHeqHeqHgrad 3321221311323211 qaHHqaHHqaHHHHHadiv3322113213322113211 aHaHaHqqqeHeHeHHHHarot12北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识4、哈密顿(、哈密顿(Hamilton)(或那勃拉)(或那勃拉Nabla)算子)算子(1)直角坐标系)直角坐标系(2)正交曲线坐标系)正交曲线坐标系(3)用算子表示梯度、散度、旋度)用算子表示梯度、散度、旋度梯度:梯度:散度:散度:旋度:旋度:zkyjxi332211
12、1111qHeqHeqHe gradadiva arota 13北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识三、张量基础知识三、张量基础知识正交坐标系正交坐标系(一)张量定义(一)张量定义1、求和约定、求和约定哑指标:两相同的指标,规定哑指标:两相同的指标,规定遍历求和遍历求和;自由指标:可任意取下标值,自由指标:可任意取下标值,aij ,不求和不求和。2、坐标变换关系、坐标变换关系坐标变换系数坐标变换系数ij旋转变换的新旧新两个坐标系(旋转变换的新旧新两个坐标系(ox1x2x3)和()和(ox1x2x3),其基向量),其基向量分别为:则两坐标轴的方向余玄为:分
13、别为:则两坐标轴的方向余玄为:jjiizyxeaeaeaeaeakajaiaa332211333231232221131211 jijiee333231232221131211aaaaaaaaaaij14北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识(1)标量的坐标变换)标量的坐标变换对点对点P在旧、新坐标系中的位置(在旧、新坐标系中的位置(x1,x2,x3)()(x1,x2,x3),),则有:则有:(2)矢量的坐标变换关系)矢量的坐标变换关系各分量间满足:各分量间满足:及:及:),(),()(321321xxxxxxP jjiiexxxaexxxaPa),(),
14、()(321321ijijiijaeeaa jjiijjiaeeaa 15北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识(3)应力量)应力量1)应力定义:)应力定义:2)应力张量分量变换关系:)应力张量分量变换关系:即有即有同理可得:同理可得:neenSFjijinS 0limmlmljijieeee ijmjlimjijillmeeee )()(lmmjlimjlmilijeeee )()(16北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识3、张量的定义、张量的定义保证物理量在不同坐标系表示下量不变,坐标转换应具有保证物理量在不同坐
15、标系表示下量不变,坐标转换应具有相同的数学结构,将具有此类特性的量称为张量(相同的数学结构,将具有此类特性的量称为张量(Tensor)。)。设设A的分量表示为的分量表示为其其3n个分量在新坐标下的分量满足:个分量在新坐标下的分量满足:则称则称A为为n阶张量。阶张量。按此定义,按此定义,标量为标量为0阶张量阶张量;矢量为矢量为1阶张量阶张量;应力量为应力量为2阶张量阶张量。niiiA.21nnnniiiijijijjjjAA.21221121 17北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识(二)张量的基本运算规则(二)张量的基本运算规则1、两个基本运算符、两个基
16、本运算符(1)克罗内克()克罗内克(Kroneker)符号)符号ijij是二阶单位张量。是二阶单位张量。两矢量的点积可表示为:两矢量的点积可表示为:ji 01ij ji jjiiijjijjiibababaebeaba jiijee 18北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识(2 2)里奇()里奇(RicciRicci)置换符号)置换符号ijkijk 三指标不同,顺序排列;三指标不同,顺序排列; 三指标不同,逆序排列;三指标不同,逆序排列; 具有相同指标时。具有相同指标时。ijkijk 是三阶张量。是三阶张量。两矢量的矢量积两矢量的矢量积: :同理,单位矢
17、量的矢量积可表示为:同理,单位矢量的矢量积可表示为:kjiijkjjiiebaebeaba 011ijk kijkjieee 19北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识(3)ijkijk与与ijij的关系恒等式的关系恒等式当当k=nk=n时,有:时,有:knkmkljnjmjlinimillmnijk jlimjmiljmjlimillmkijk1000020北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识2、张量的基本运算、张量的基本运算(1)张量的加减运算)张量的加减运算对同阶张量对同阶张量(2)张量乘法)张量乘法1)张量的
18、并积)张量的并积C为为r+s阶张量。阶张量。如:应力张量是力作用面方向矢量与力矢量的并积。如:应力张量是力作用面方向矢量与力矢量的并积。nnnniiiiiiiiiiiibabBaA.21212121BAC ;则:及srsrjjjiiijjjiiibabBaA.21212121ABC ;则:及21北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识2)张量的缩并)张量的缩并一一n阶张量具有两相同的指标阶张量具有两相同的指标 时,经求和运算,张量时,经求和运算,张量A缩并为缩并为n - 2阶张量。阶张量。3)张量内(点)积)张量内(点)积r阶张量阶张量s阶张量阶张量内积所得
19、张量内积所得张量C为为r+s-2阶张量阶张量nijjiiaA.21srsrssrrjjiiijjiiijjjjjjiiiiiieeeeecbeeeeeeaBAC. .2121212121212121.ssjjjjjjeeebB.2121.rriiiiiieeeaA.2121.22北京工业大学市政学科部马长明 高等流体(水)力学讲稿数学基础知识数学基础知识3、二阶张量的分解、二阶张量的分解对称张量与反对称张量对称张量与反对称张量对任意二阶张量对任意二阶张量C对称部分:对称部分:有有6个独立分量。个独立分量。反对称部分:反对称部分:有三个独立分量。有三个独立分量。4、梯度、散度及旋度的张量表示、梯
20、度、散度及旋度的张量表示)(21)(21jiijjiijijcccccC)(21jiijijccsS)(21jiijijccaAiiexgradiixaadivkijijkexaarot23北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念一、连续介质模型(假设)一、连续介质模型(假设)尺度:尺度: 分子尺度分子尺度 l;研究问题的宏观计量尺度研究问题的宏观计量尺度L;流体质点尺度流体质点尺度a。1、假设成立条件:、假设成立条件:即流体质点在宏观上无限小,在微观上无限大。即流体质点在宏观上无限小,在微观上无限大。2、连续介质模型含义、连续介质
21、模型含义(1)按流体质点尺度,流体在时间和空间上在域内连续分布;)按流体质点尺度,流体在时间和空间上在域内连续分布;(2)描述流体特性和运动的物理量是定义在流体质点上的,几乎处)描述流体特性和运动的物理量是定义在流体质点上的,几乎处处连续可微;处连续可微;(3)不发生化学反应,是热力平衡过程,服从热力学关系。)不发生化学反应,是热力平衡过程,服从热力学关系。Lal24北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一讲第一讲 流体力学的基本概念流体力学的基本概念二、描述流体运动的两种方法二、描述流体运动的两种方法1、拉格朗日法(、拉格朗日法(Lagrangian Method)(1)质点
22、运动方程:)质点运动方程:a,b,c:拉格朗日变量,为:拉格朗日变量,为t=0时,流体质点的坐标值。时,流体质点的坐标值。(2)特点:质点运动学的研究方法,难以形成对流体域整体)特点:质点运动学的研究方法,难以形成对流体域整体运动特性的描述。运动特性的描述。(3)流体质点的运动速度:)流体质点的运动速度:(4)流体质点的运动加速度:)流体质点的运动加速度:)3 , 2 , 1( ),(itcbaxxii)3 , 2 , 1( itxvii)3 , 2 , 1( 22itxtvaiii25北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念2、欧
23、拉法(、欧拉法(Eulerian Method)(1)速度场:在流体场中的每个空间点上设想有测速仪,则可记录:)速度场:在流体场中的每个空间点上设想有测速仪,则可记录:(2)特点:关注流体运动场的研究,而不是单个流体质点的运动。)特点:关注流体运动场的研究,而不是单个流体质点的运动。(3)欧拉法对流体质点参数变化的描述)欧拉法对流体质点参数变化的描述质点(随体)导数(质点(随体)导数(mterial derivative, substantial derivative)对定义在质点上的物理量对定义在质点上的物理量(x,yx,y, ,z,tz,t) ),其随质点运动的时间变化,其随质点运动的时间
24、变化率表示为:率表示为:引入随体导数算子:引入随体导数算子: ) 3 , 2 , 1( ),(itzyxuuiivtxvtDtDii iixvtDtD26北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念3、欧拉法与拉格朗日法之间的相互转化、欧拉法与拉格朗日法之间的相互转化(1)基本关系)基本关系其中:其中: 欧拉坐标;欧拉坐标; 拉格朗日流体质点运动坐标。拉格朗日流体质点运动坐标。对所描述的物理量对所描述的物理量F,有:,有:(2)变换过程)变换过程1)运动方程()运动方程(L)与流动场()与流动场(E)间的变换)间的变换1,2,3)(i
25、),(tcbaxxLiEiEixLix ),(),(tcbaFtzyxFLE1,2,3)(i ),(),(ttcbaxvtzyxvLiLiEi1,2,3)(i ),(tzyxaaii; 1,2,3)(i ),(tcbaxxii27北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念变换成立条件,变换成立条件,Jacobi行列式行列式2)物理量间的变换)物理量间的变换0),(),(McbazyxJ ),(),(tcbaFtzyxFLE1,2,3)(i ),(tzyxaaii; 1,2,3)(i ),(tcbaxxii28北京工业大学市政学科部马长
26、明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念4、流动的基和描述与显示、流动的基和描述与显示(1)迹线()迹线(Pathline)(2)流线()流线(Streamline)(3)脉线()脉线(Streak line,染色线),染色线)(4)时间线()时间线(time line, 流体线)流体线)流管的概念流管的概念流体面保持定理流体面保持定理29北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念三、流体微团的运动分析三、流体微团的运动分析1、亥姆霍兹(、亥姆霍兹(Helmholtz)运动分解)运动分解(1)速度梯
27、度张量)速度梯度张量 其中对称张量:其中对称张量: 反对称张量:反对称张量:) 3 , 2 , 1( )(200jOxxvvviMijMjMj ijijijASxvjiijijxvxvS21jiijijxvxvA2130北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念(2)速度梯度张量的对称张量)速度梯度张量的对称张量线变形率与角变形率线变形率与角变形率1)线变形率)线变形率2)体积变形率)体积变形率3)角变形率)角变形率yyyvs 22zzzws 33xxxus 11yuxvssxy212112 zvywssyz213223 xwzuss
28、zx211331 vzwyvxudtd 1lim031北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念(3)速度梯度张量的反对称张量)速度梯度张量的反对称张量转动角速度转动角速度1)流体微团的转动角速度)流体微团的转动角速度2)旋度(涡量)旋度(涡量Vorticity)与旋转角速度的关系)与旋转角速度的关系有:有:)(211zvywx )(213xvyuz )(212xwzuy kkkijijkeexvwvuzyxkjiv kk21 32北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念流
29、体微团的运动:流体微团的运动:平动平动 + 转动转动 + 膨胀变形膨胀变形 + 前切变形前切变形2、按运动学角度分类流体运动、按运动学角度分类流体运动(1)有旋运动与无旋运动)有旋运动与无旋运动有旋:有旋:00;无旋:无旋: = = 0 0。(2 2)可压缩流动与不可压缩流动)可压缩流动与不可压缩流动可压缩:可压缩:不可压缩:不可压缩:0 v0 v33北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念四、作用在流体上的力、应力张量及牛顿本构方程四、作用在流体上的力、应力张量及牛顿本构方程1、作用于流体上的力、作用于流体上的力(1)质量力(亦称
30、体力)质量力(亦称体力body forces,远程力,远程力long range forces)1)定义;)定义;2)实例;)实例;3)数学表示)数学表示单位质量力:单位质量力:量纲;量纲;质量力合力:质量力合力: bbFf1lim0VbbdfF 34北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念(2)表面力()表面力(surface forces, 近程力近程力Short range forces)1)定义;)定义;2)实例;)实例;3)表述形式)表述形式单位面积上的应力:单位面积上的应力:表面力的合力:表面力的合力:SPSn0lim
31、SndsP 35北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念(3)表面应力张量及其特性)表面应力张量及其特性1)应力张量)应力张量2 2)应力张量的特性)应力张量的特性i i)应力张量是二阶对称张量;)应力张量是二阶对称张量;iiii)静止流体、无粘性流体的应力张量为:)静止流体、无粘性流体的应力张量为:zzzyzxyzyyyxxzxyxxij nn ijijp 36北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念2、牛顿流体的本构方程、牛顿流体的本构方程应力张量与变形率张量的关系应
32、力张量与变形率张量的关系(1)基本假设()基本假设(Stokes 在在1848年提出)年提出)1)ijij是是s sijij的线性函数的线性函数2 2)流体是各向同性的,满足牛顿切应力定律)流体是各向同性的,满足牛顿切应力定律3 3)流体静止时,应变率张量为零,流体中的应力等于流体静)流体静止时,应变率张量为零,流体中的应力等于流体静压强压强p pijijijbsa jixuyujixy ij 0ijijijps 37北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念(2)理论推导结果)理论推导结果其中:其中:为第二粘性系数为第二粘性系数体积膨胀粘性系数,与体积变化的能量体积膨胀粘性系数,与体积变化的能量耗散机制有关。耗散机制有关。当流体压缩性不大时,可取:当流体压缩性不大时,可取:ijij 2 kkijijxusp32 32 ijij 32 2kkijijxusp38北京工业大学市政学科部马长明 高等流体(水)力学讲稿第一章第一章 流体力学的基本概念流体力学的基本概念本讲结束本讲结束