第一讲变化方程组传递过程的理论基础PPT讲稿.ppt

《第一讲变化方程组传递过程的理论基础PPT讲稿.ppt》由会员分享，可在线阅读，更多相关《第一讲变化方程组传递过程的理论基础PPT讲稿.ppt（102页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第一讲变化方程组传递过程的理论基础第1页，共102页，编辑于2022年，星期一教 学 安 排 教学环节由三部分组成教学环节由三部分组成：1、授课：、授课：15 15 次课堂讲授次课堂讲授次课堂讲授次课堂讲授 2、作业与学习报告：、作业与学习报告：15 15 次课外作业，每周第一次课前交。次课外作业，每周第一次课前交。次课外作业，每周第一次课前交。次课外作业，每周第一次课前交。1 1 份学习报告，考试前交。份学习报告，考试前交。份学习报告，考试前交。份学习报告，考试前交。3、课程考试、课程考试分数分布分数分布 1 1 1 1、作业作业作业作业 15%15%15%15%2 2 2 2、学习报告学习

2、报告学习报告学习报告 15%15%15%15%3 3 3 3、考试考试考试考试 70%70%70%70%第2页，共102页，编辑于2022年，星期一课 程 网 站网址：ftp:/账号：学号学号s初始密码：12345第3页，共102页，编辑于2022年，星期一变化方程组传递过程的理论基础第4页，共102页，编辑于2022年，星期一变化方程组-传递过程的理论基础传递过程的传递过程的第5页，共102页，编辑于2022年，星期一第 一 个 问 题专业术语 “传递过程”代表的是什么自然现象?第6页，共102页，编辑于2022年，星期一答案：一个定义“传递过程”代表了 “由于物质所具有的三种基本性质：动量

3、、能量和质量在空间中的非均匀分布而导致其在时空中的迁移过程。”第7页，共102页，编辑于2022年，星期一观 点 1根据这个定义,“传递过程是发生在非平衡态下的动态过程暨速率过程。”第8页，共102页，编辑于2022年，星期一观 点 2根据这个定义,“传递过程进行的方向合理且自然地确定为：当其自发进行时，从非平衡态趋向于平衡态。”第9页，共102页，编辑于2022年，星期一关于能量的说明 在传递过程领域里,被传递的能量特指物质的内能,而引起传递过程的推动力所涉及的能量则可能包含多种形式的能量，如机械能、电能、磁能等等。第10页，共102页，编辑于2022年，星期一一 个 推 论 “当动量、能量

4、和质量密度在空间中的非均匀分布是唯一的推动力时，它们将自发地从高密度区域传递到低密度区域。”第11页，共102页，编辑于2022年，星期一动量和它的密度动量，M=mv，是一个矢量。速度,v=M/m,等于单位质量的动量,是动量密度的一种表征。因此速度在空间中的分布代表了动量密度在空间中的分布。第12页，共102页，编辑于2022年，星期一能量和它的密度物质的内能,U=mCVT,是一个标量。温度正比于单位质量的内能,是内能密度的一种表征。因此温度在空间中的分布代表了内能密度在空间中的分布。第13页，共102页，编辑于2022年，星期一质量和它的密度传递过程中被传递的质量特指混合物中化学组分的质量。

5、化学组分的质量密度可以用它的浓度、质量分率或摩尔分率来表征。因此浓度、质量分率或摩尔分率在空间中的分布代表了组分质量密度在空间中的分布。第14页，共102页，编辑于2022年，星期一常规传递过程“常规传递过程”表示 “由于速度、温度和组成在空间中的非均匀分布导致的动量、内能和化学组分在时空中的迁移过程。”第15页，共102页，编辑于2022年，星期一三 个 尺 度传递过程可以在三个尺度上进行描述:宏观尺度 Macroscopic Level 微观尺度 Microscopic Level 分子尺度 Molecular Level第16页，共102页，编辑于2022年，星期一宏 观 尺 度主要在工

6、程应用中涉及考虑一个系统与环境发生动量、能量和组分质量的相互传递时这些物理量的总量的变化。第17页，共102页，编辑于2022年，星期一微 观 尺 度考虑动量、能量和质量在时空中的分布，导出相关物理量在一个系统内部的分布图象。采用物理量场的表述方法采用物理量场的表述方法 本课程的重点研究领域本课程的重点研究领域第18页，共102页，编辑于2022年，星期一分 子 尺 度采用统计物理的研究方法主要在传递性质的计算中涉及 从物质结构和分子间相互作用的角度寻求对动量、能量和质量传递现象机理的基本理解。第19页，共102页，编辑于2022年，星期一传递过程的三种机制（1 1）对流传递 物理量由流动的流

7、体携带着进入或离开一个区域，从而导致该物理量在空间中的传递。分子传递 分子布朗运动所引起的物理量在空间中的传递.这两种机制被分类为近程传递，物理量在空间中逐点传递。第20页，共102页，编辑于2022年，星期一传递过程的三种机制（2 2）远程传递 被传递的物理量在空间某一处被转化为一种势能场，例如重力场、电磁场、电场和磁场，这种势能场被传播到空间另一处，然后被物质分子所接收并再转化为被传递的物理量。从表观上看，该物理量在空间中实现的是点对点传递。第21页，共102页，编辑于2022年，星期一动动 量量 与与 力力(1)牛顿第一定律(惯性定律)每一个物体都保持静止或匀速直线运动状态，除非被施加于

8、它的力强迫改变这种状态。经典力学的大厦是建立在牛顿运动三定律的基础上的。第22页，共102页，编辑于2022年，星期一动动 量量 与与 力力(2)牛顿第二定律(加速度定律)一个物体的加速度等于所有作用于它的外力的合力除以它的质量。牛顿第三定律(作用力与反作用力定律)对于每一个作用力，总是存在一个大小相等而方向相反的反作用力。第23页，共102页，编辑于2022年，星期一动动 量量 与与 力力(3)在相对论效应可以忽略的情况下，物体的质量不随其运动状态而改变。于是把它代入牛顿第二定律，得到此式表明，作用于一个物体的合外力等于传递动量给该物体的速率。即，该式给出了力的定义。第24页，共102页，编

9、辑于2022年，星期一动动 量量 与与 力力(4)根据力的这个定义，我们可以对其它两条牛顿运动定律给出另一种解释：动量守恒。牛顿第一定律 一个物体的动量保持恒定，除非有外加的动量传递给它。牛顿第三定律 发生相互作用时，一个物体获得的动量总是等于另一个物体失去的动量。第25页，共102页，编辑于2022年，星期一动量通量与应力 应力被定义为作用于单位面积上的力，具有两个特征方向：力的方向和受力表面的法线方向。xy 表示 y 方向的力作用在垂直于 x 的单位面积上。根据牛顿第二定律，xy 可以被解释为单位时间里通过垂直于 x 的单位面积传递的 y 方向的动量，即动量通量。第26页，共102页，编辑

10、于2022年，星期一变化方程组-传递过程的理论基础分子传递现象的分子传递现象的第27页，共102页，编辑于2022年，星期一分子传递现象的本构方程分子传递现象的本构方程 在常规传递过程中，动量、内能和组分质量的传递是由于速度、温度和组分浓度在空间中的非均匀分布，而描述分子传递通量与速度梯度、温度梯度和组分浓度梯度之间的定量关系的数学表达式即称为分子传递现象的本构方程。分子传递本构方程的最初形式都源于对一维传递过程的实验观察和分析，发现在多数情况下均可表示为线性齐次函数。第28页，共102页，编辑于2022年，星期一牛顿粘性定律牛顿粘性定律 描述一维动量分子传递规律的本构方程称为牛顿粘性定律：该

11、定律表明：动量通量是动量密度梯度的线性齐次函数。式中线性项的系数称为流体的粘度。如果流体的粘度值不依赖于速度梯度，则归类为牛顿流体；否则归类为非牛顿流体。(1.？-？)第29页，共102页，编辑于2022年，星期一傅立叶导热定律 内能的分子传递通常被称作热传导。描述一维热传导规律的本构方程称为傅立叶定律：(9.1-2)该定律表明：内能通量是内能密度梯度的线性齐次函数，式中线性项的系数称为材料的导热系数。第30页，共102页，编辑于2022年，星期一费克扩散定律 化学组分的分子传递通常被称为扩散。描述二元体系(A+B)中的一维扩散规律的本构方程称为费克定律：(17.1-2)该定律表明：组分 A

12、的质量通量是 A 的质量密度梯度的线性齐次函数，式中线性项的系数称为A对B的二元扩散系数。第31页，共102页，编辑于2022年，星期一三维情况下的分子传递(1)工程应用中的传递过程大多数在三维空间中发生。为了处理三维空间中的传递问题，我们必须给出三维形式的分子传递通量的表达式。既然实验已经证明由相应物理量密度的非均匀分布引起的一维分子传递通量是物理量密度梯度的线性齐次函数，因而把这种函数关系推广到三维情况就是很自然的选择。第32页，共102页，编辑于2022年，星期一三维情况下的分子传递(2)首先考虑热传导。我们假设热通量仍然是温度梯度 的线性齐次函数，称为广义傅里叶定律。不过，由于热通量和

13、温度梯度两者都是矢量，而两个矢量变量之间的线性齐次函数的一般形式是(9.1-7)第33页，共102页，编辑于2022年，星期一三维情况下的分子传递(3)对于各向同性材料，导热系数应该是各向同性二阶张量，其一般形式为 因此一般而言，导热系数是一个对称二阶张量，具有6个独立分量。第34页，共102页，编辑于2022年，星期一三维情况下的分子传递(4)于是各向同性材料的三维热传导方程为(9.1-6)(9.1-6)式在直角坐标系中的分量方程为(B.2-3)(B.2-2)(B.2-1)第35页，共102页，编辑于2022年，星期一各向异性材料的实例木材中的热传导过程 木材就是一种各向异性材料，顺着纹理方

14、向的导热性能与垂直于纹理方向的导热性能有所不同。xyz第36页，共102页，编辑于2022年，星期一三维情况下的分子传递(4)类似地，在广义费克扩散定律中，扩散系数也是一个对称二阶张量，其一般形式为(17.7-3)对于各向同性材料才简化为(17.1-10)第37页，共102页，编辑于2022年，星期一三维情况下的分子传递(5)在直角坐标系下，(B.3-3)(B.3-2)(B.3-1)第38页，共102页，编辑于2022年，星期一三维情况下的分子传递(6)动量传递的情况更为复杂，由于动量通量和速度梯度都是二阶张量，所以线性齐次函数关系中的系数是一个四阶张量：此式表明粘度具有54个独立分量!由于

15、是对称的，所以 应该对前两个下标对称：(1.2-3)即第39页，共102页，编辑于2022年，星期一三维情况下的分子传递(7)它的非零分量可以用三个标量表达。除了液晶、长链高分子液体以及高度取向的纤维悬浮液以外，大多数流体都是均相的各向同性材料。我们规定牛顿流体也必须是各向同性流体，其粘度必然是各向同性张量。而四阶各向同性张量的一般形式为第40页，共102页，编辑于2022年，星期一三维情况下的分子传递(8)因为粘度张量对于前两个下标是对称的，所以必然有第41页，共102页，编辑于2022年，星期一三维情况下的分子传递(9)将此结果代入(1.2-3)式中，我们得到：右侧第二项可以拆分为各向同性

16、部分和各向异性部分：第42页，共102页，编辑于2022年，星期一三维情况下的分子传递(10)于是，(1.2-3)式简化为它表征流体的膨胀或压缩对动量通量的贡献。其 中 各 向 同 性 项 的 系 数 被 称 为 膨 胀 粘 度(dilatational viscosity)，记为第43页，共102页，编辑于2022年，星期一三维情况下的分子传递(11)把上述方程改写为张量符号的形式，我们就得到了广义牛顿粘性定律：膨胀粘度仅仅对于高温或高频率振动下的多原子气体才有实际意义，在其它其况下的效应通常可以忽略不计。对于不可压缩流体，速度的散度等于零，因此各向同性项自然消失。(1.2-7)第44页，共

17、102页，编辑于2022年，星期一三维情况下的分子传递(12)在直角坐标系下，式(1.2-7)的分量方程为：(B.1-1)(B.1-2)(B.1-3)第45页，共102页，编辑于2022年，星期一三维情况下的分子传递(13)(B.1-4)(B.1-5)(B.1-6)在柱坐标系和球坐标系下的展开式要复杂得多，可参看附录B选用。第46页，共102页，编辑于2022年，星期一变化方程组传递过程的理论基础第47页，共102页，编辑于2022年，星期一物理量的通用衡算方程 (1)(1)考虑直角坐标系下的一个微元六面体，物理量 E 的通用衡算方程可表达为来自对流传递的 E 收入来自分子传递的 E 收入来自

18、远程传递的 E 收入+E 在该体积中的积累速率=第48页，共102页，编辑于2022年，星期一物理量的通用衡算方程 (2)(2)上式左侧的积累项可以表示为体积分式中 E 是物理量 E 的空间密度，即，单位体积中所包含的 E 的量。第49页，共102页，编辑于2022年，星期一物理量的通用衡算方程 (3)(3)根据积分中值定理，上述体积分等于被积函数在积分体中某一点处的值乘以积分体的体积：第50页，共102页，编辑于2022年，星期一物理量的通用衡算方程 (4)(4)E进入立方体中的对流传递项等于第51页，共102页，编辑于2022年，星期一物理量的通用衡算方程 (5)(5)根据积分中值定理，上

19、述面积分等于被积函数在积分面上某一点处的值乘以积分面的面积：第52页，共102页，编辑于2022年，星期一以上各式中的各个、的值不一定相同，但均满足物理量的通用衡算方程 (6)(6)第53页，共102页，编辑于2022年，星期一物理量的通用衡算方程(7)(7)E进入立方体中的分子传递项等于第54页，共102页，编辑于2022年，星期一物理量的通用衡算方程 (8)(8)根据积分中值定理，上述面积分等于被积函数在积分面上某一点处的值乘以积分面的面积：第55页，共102页，编辑于2022年，星期一以上各式中的各个、的值不一定相同，但均满足物理量的通用衡算方程 (9)(9)第56页，共102页，编辑于

20、2022年，星期一物理量的通用衡算方程(10)(10)物理量E的远程传递项（即局部产生项）可以表示为 V 上的体积分式中 RE 是 E 的局部产生速率，即，单位时间里在单位体积中产生的 E 的量。第57页，共102页，编辑于2022年，星期一物理量的通用衡算方程 (11)(11)根据积分中值定理，上述体积分等于被积函数在积分体中某一点处的值乘以积分体的体积：第58页，共102页，编辑于2022年，星期一物理量的通用衡算方程 (12)(12)把以上结果代入原方程并用 x y z除以整个方程，然后取 x、y、z 0的极限，得到第59页，共102页，编辑于2022年，星期一物理量的通用衡算方程 (1

21、3)(13)第60页，共102页，编辑于2022年，星期一物理量的通用衡算方程 (14)(14)根据偏导数的定义，上式取极限的结果是以下微分方程第61页，共102页，编辑于2022年，星期一物理量的通用衡算方程 (15)(15)将左侧方括号里的项展开并移项整理，得到物理量通用衡算方程的微分形式(直角坐标系)：第62页，共102页，编辑于2022年，星期一随 体 导 数(1)令 E(x,y,z,t)是定义在时空中的一个连续函数，对于一个观察者而言，E(x,y,z,t)的变化速率是它对时间的全导数：式中vob=(vob,x,vob,y,vob,z)是观察者的运动速度。第63页，共102页，编辑于2

22、022年，星期一随 体 导 数(2)特别地，如果 vob 等于该位置处流体流动的速度 v，我们则得到一个特别的全导数我们随流体一起运动时观察到的E随时间的变化率，称为E的随体导数，记为 随体导数也称为物质导数，因为其物理意义亦为一个流体质点的物理性质在流动过程中随时间的变化率。第64页，共102页，编辑于2022年，星期一物理量的通用变化方程 (1)(1)将随体导数引入物理量通用衡算方程，得到该方程的随体导数表达形式(直角坐标系)：物理量通用衡算方程亦称为物理量通用变化方程(Equation of Change)。第65页，共102页，编辑于2022年，星期一连续性方程(1)取通用变化方程中的

23、物理量为流体的总质量，则 E=，FE=0，RE=0，我们得到变化方程的一个特定形式连续性方程：该方程的矢量形式为：(3.1-4)(3.1-3)第66页，共102页，编辑于2022年，星期一连续性方程(2)重新排列方程中的各项并采用用随体导数表达形式，得到其矢量形式为：(Tab.3.5-1.A)第67页，共102页，编辑于2022年，星期一连续性方程(3)对于不可压缩流体，流体微元的密度在流动过程中保持不变，即于是，不可压缩流体的连续性方程简化为(3.1-5)第68页，共102页，编辑于2022年，星期一组分连续性方程 (1)令物理量E为化学组分 的质量，则E的密度就是的浓度E的分子传递通量就是

24、的扩散通量E的局部产生速率就是在化学反应中的生成速率第69页，共102页，编辑于2022年，星期一组分连续性方程 (2)将这些量代入通用变化方程，我们得到关于化学组分 的组分连续性方程：(19.1-7)组分连续性方程也可以表达为随体导数形式：第70页，共102页，编辑于2022年，星期一组分连续性方程 (3)将上两式代入式(19.1-7)，移项整理后得到 因为的质量浓度等于流体密度乘以的质量分率 =w ,所以式(19.1-7)的左侧和右侧第一项可展开为第71页，共102页，编辑于2022年，星期一组分连续性方程 (4)根据式(3.1-4)，方括号的值恒为零。圆括号正好是质量分率w 的随体导数，

25、于是我们得到用质量分率表达的组分连续性方程的随体导数形式：(19.1-7a)在直角坐标系下为第72页，共102页，编辑于2022年，星期一组分连续性方程 (5)加和所有化学组分的组分连续续性方程，就得到总连续性方程，其结果与第三章所导出的连续性方程完全一致。第73页，共102页，编辑于2022年，星期一组分连续性方程 (6)对于仅含A和B两种组分的二元体系，分子传递通量可以用广义费克扩散定律(式17.7-3)表达，将其代入组分连续性方程(式19.1-7)，得如果 DAB为常数，则有(19.1-16)第74页，共102页，编辑于2022年，星期一组分连续性方程 (7)在直角坐标系下展开上式，有或

26、表达成随体导数形式：第75页，共102页，编辑于2022年，星期一运 动 方 程(1)取通用变化方程中的物理量为动量。因为动量是质量，包含三个分量，所以动量密度亦为矢量，动量密度在直角坐标下的表达式为第76页，共102页，编辑于2022年，星期一运 动 方 程(2)动量的每个分量的通量亦为一矢量，在直角坐标系下可表示为当传递方向与动量方向相同时，动量通量除粘性传递通量外还包含流体的静压强。第77页，共102页，编辑于2022年，星期一运 动 方 程(3)动量的局部产生速率亦为一矢量，在直角坐标系下可表示为式中 代表导致远程动量传递的场强矢量。如果只有重力场，则等于重力加速度矢量；如果还有电场、

27、磁场等其它远程传递场，则为各作用场强矢量的和。第78页，共102页，编辑于2022年，星期一运 动 方 程(4)把与动量的x分量相关的量带入通用变化方程，得到第79页，共102页，编辑于2022年，星期一运 动 方 程(5)把圆括号展开并移项整理，得到根据连续性方程，方括号的值等于零，于是第80页，共102页，编辑于2022年，星期一运 动 方 程(6)类似地，对动量的y分量和z分量，有第81页，共102页，编辑于2022年，星期一运 动 方 程(7)三个分量方程也可以表达成矢量和张量形式：上述方程被称为运动方程。第82页，共102页，编辑于2022年，星期一运 动 方 程(8)运动方程也可表

28、达为随体导数形式：在直角坐标系下的展开式为：(Tab.3.5-1.B)第83页，共102页，编辑于2022年，星期一运 动 方 程(9)将广义牛顿粘性定律式(式1.2-7)代入运动方程(式Tab.3.5-1.B)，对于常物性流体，得：(3.5-6)引入修正压力(亦称为动压头)P=p+gh，得：(3.5-7)上两式即著名的纳维-斯托克斯(Navier-Stokes)方程。第84页，共102页，编辑于2022年，星期一运 动 方 程(10)N-S方程在直角坐标系下的展开式为：第85页，共102页，编辑于2022年，星期一能 量 方 程 (1)对于能量传递，我们关注的是内能传递。但由于在流体流动过程

29、中，流体的机械能与内能之间会发生相互转化，所以我们首先令物理量 E 为动能与内能的和(称为总能量)，则 总能量空间密度为 总能量通量为 总能量局部产生速率为 第86页，共102页，编辑于2022年，星期一能 量 方 程 (2)把以上表达式代入通用变化方程，我们就得到了总能量方程(11.1-7)第87页，共102页，编辑于2022年，星期一机械能方程(1)为了获得机械能的表达式，用速度 v 与运动方程(Tab.3.5-1.B)式进行点积，可得此式的左侧可以展开为(*)第88页，共102页，编辑于2022年，星期一机械能方程(2)式中的负号项可以合并为下式的方括号项，根据连续性方程，方括号等于零。

30、对于(*)式的右侧，我们有第89页，共102页，编辑于2022年，星期一机械能方程(3)于是，(*)式可以改写为此式被称为机械能方程。(3.3-1)第90页，共102页，编辑于2022年，星期一能 量 方 程 (3)从(11.1-7)式中减去机械能方程(3.3-1)式，就得到了内能方程：借助于连续性方程，上式可以被表示成随体导数的形式：(11.2-1)(11.2-2)第91页，共102页，编辑于2022年，星期一能 量 方 程 (4)此式清楚地表明，流体微元内能的增加有四个来源：1、热传导；2、压缩功；3、粘性耗散；4、其它能量转化为内能。(11.2-2)第92页，共102页，编辑于2022年

31、，星期一能 量 方 程 (5)假定平衡态下的热力学关系在传递过程中局部适用，我们有将此式代入内能的随体导数中，我们可以把(11.2-2)式的左侧改写为第93页，共102页，编辑于2022年，星期一能 量 方 程 (6)比容的随体导数可表为 于是，我们得到了用温度场表达的能量方程：因为所以(Tab.11.4-1.I)第94页，共102页，编辑于2022年，星期一能量方程(7)把广义傅里叶定律(9.1-6)(9.1-6)式代入能量方程(Tab.11.4-1.I)(Tab.11.4-1.I)中，得到 如果希望使用等压热容，可将(9.1-6)(9.1-6)式代入(Tab.11.4-1.J)(Tab.1

32、1.4-1.J)式中，得到式中 QQp p 代表热量的局部产生速率。第95页，共102页，编辑于2022年，星期一能量方程(8)以上能量方程适用于常密度和均匀组成的系统。对于非均匀组成的系统，热通量中还应该计入扩散分子所携带的内能：(Tab.19.2-4.F)式中 是组分 的偏摩尔焓。第96页，共102页，编辑于2022年，星期一耗散函数 (1)根据广义牛顿粘性定律，粘性耗散可以被表示为：(3.3-3)式中 被称为耗散函数，被称为膨胀耗散函数。第97页，共102页，编辑于2022年，星期一耗散函数 (2)在直角坐标系下，耗散函数可以被展开为：(B.7-1)第98页，共102页，编辑于2022年

33、，星期一能量方程(9)对于常密度和均匀组成的系统，能量方程在直角坐标系下的展开式为：(B.9-1*)第99页，共102页，编辑于2022年，星期一变化方程小结(1)解决多组分系统中传递过程的基本出发点是一组变化方程：1、连续性方程(密度场方程)；2、运动方程(速度场方程)；3、能量方程(温度场方程)；4、组分连续性方程(浓度场方程)。第100页，共102页，编辑于2022年，星期一变化方程小结(2)1.2.3.4.(19.1-7)(Tab.11.4-1.I)(Tab.3.5-1.B)(Tab.3.5-1.A)第101页，共102页，编辑于2022年，星期一变化方程小结(3)1、通过包含在随体导数中的对流项，速度场显式地影响温度场和浓度场。2、通过与流体密度和分子传递性质的关系，温度场和浓度场隐式地影响速度场。3、对于常物性流体，速度场独立于温度场和浓度场。第102页，共102页，编辑于2022年，星期一