化工传递 3动量传递方程的若干解ppt课件.ppt

Ch3：动量传递方程的若干解,本章讨论重点流体作简单层流流动时，动量传递方程的典型求解。主要包括：,1 两平壁间的稳态层流；,2 圆管与套管环隙间的稳态层流；,4 极慢黏性流动（爬流）；,5 势函数与理想流体的流动。,3 非稳态流动的求解；,课后学习与作业：,第三章的概念和例题；第三章作业：3-1，3-5，3-8，3-11，3-14,动量传递方程的分析,动量传递方程组：,当流体不可压缩时，=常数,(2-46),(2-41b),变量数：ux，uy，uz，p；方程数：4,动量传递方程组的特点：,（1）非线性偏微分方程；,方程组的求解目的 获得速度与压力分布,动量传递系数 CD（或 f ）等。,（2）质点上的力平衡，仅能用于规则的层流求解。,方程组求解的分类：,（1）对于非常简单的层流，方程经简化后，其形式非常简单，可直接积分求解解析解；,（2）对于某些简单层流，可根据流动问题的物理特征进行化简。简化后，积分求解物理近似解；,（3）对于复杂层流，可采用数值法求解；将方程离散化，然后求差分解；,（4）对于湍流，可先进行适当转换，再根据问题的特点，结合实验，求半理论解。,1 两平壁间的稳态层流 P50,平壁间的轴向平行层流,应用场合：板式热交换器，各种平板式膜分离装置等；,特点：平壁无限宽，忽略平壁宽度方向流动的变化，可认为是一维流动；,一维流动：,不可压缩流体：,平壁无限宽：,连续性方程,物理模型：流体在两平壁间作平行稳态层流流动，例如板式热交换器、各种平板式膜分离装置等。,设=常数；稳态；远离流道进、出口；流体仅沿 x方向流动：,一、方程的简化,平壁间的轴向平行层流,（1）连续性方程的简化,（2）运动方程的简化,x 方向：,z 方向：,y 方向：,(b),(c),（a),（b）对 y 积分得,对x 微分得,因,仅是 y 的函数,(3-18),二、方程的求解,边界条件（B.C.）：,（1）,（2）,速度分布为,抛物线形,当,时速度最大,平壁间不可压缩流体作稳态层流的速度分布,忽略流道进、出口处的影响，流体速度分布呈抛物线形状,0,三、平均流速与流动压降,在流动方向上，取单位宽度的流通截面,则通过该截面的体积流率为,平均流速：,压降：,范宁摩擦因子（推导过程？）：,圆管壁面处的剪应力,流体的平均流速,流动阻力,(3-30),与例题 3-3不同之处？,6,例： 10摄氏度的水以4m3/h的流率流过以宽1m，高0.1m的矩形水平管道。假定流动已经充分发展，流动为一维，试求截面上的速度分布及通过每米长管道的压力降。已知10摄氏度水的粘度为1.307mN*s/m2,解：主体流速,为了判断此情况下流体的流型，需计算Re，流道为矩形，故Re中的几何尺寸应采用当量直接de替代，de的值为,故流动为层流，可确定速度分布方程，即,每米长管道的压力降为,应用场合：膜状冷凝，湿壁塔吸收等；,特点：稳态层流，一维流动；一侧紧贴壁面，另一侧为自由表面；,不可压缩流体在液膜内速度分布方程：,主体流速：,液膜厚度：,平壁面上的降落液膜流动 P54,(3-36),(3-37),(3-38),例： 某流体的运动粘度为2*10-4m2/s，密度为800kg/m3，欲使该流体沿宽为1m的垂直平壁下降的液膜厚度达到2.5mm，则液膜下降的质量流率应为多少？,解：由,因此，单位宽度的质量流量为,上述计算结果仅当液膜内流动为层流时才是正确的，因此，需要验算流动的Re数。当量直径,故,由此可知，流动确为层流，上述计算结果是正确的。,一、圆管中的轴向稳态层流,二、套管环隙中的轴向稳态层流,三、旋转黏度计的测量原理,2 圆管与套管环隙间的稳态层流 P55,一、圆管中的轴向稳态层流,流体在圆管中的流动问题许多工程科学中遇到。,设：不可压缩流体在水平圆管中作稳态层流流动，所考察的部位远离管道进、出口，流动为沿轴向的一维流动。,柱坐标连续性方程的简化,N-S方程简化,r 分量:,z 分量:, 分量:,(3-44),速度分布,管中心最大流速,平均流速,一、圆管中的轴向稳态层流,不可压缩流体在水平圆管中作稳态层流流动：,(3-48),(3-49),压力降,范宁摩擦因子,动阻力：,圆管壁面处的剪应力：,摩擦系数= 64/Re,(3-51),例3-5 毛细管粘度计测量流体粘度的原理是使被测流体在一细长的圆管（毛细管）中作稳态层流流动，测定流体流过整个圆管的压力降，从而求出流体的粘度。已知甘油在299.6K下流过长度为0.3048m，内径为0.00254m的水平圆管。在体积流率为1.878*10-6m3/s时，测得压降为2.76*105pa。在299.6K时甘油的密度为1261kg/m3。试求甘油的粘度？,解：由式（3-51）得,式中,L = 0.3048 m将以上各值代入上式中，得,校核流动的雷诺数,因此流动为层流，计算是正确的。,二、套管环隙中的轴向稳态层流 P58,流体在两根同心套管环隙空间沿轴向的流动在物料的加热或冷却时经常遇到，如套管换热器。,设：不可压缩流体在两管环隙间沿轴向流过。设所考察的部位远离进、出口，求解套管环隙内的速度分布、主体流速以及压力降的表达式。,套管环隙中层流的变化方程与圆管相同，即,B.C. 为,(3-44),速度分布,由B.C.+,由B.C.+,联立二式,P58, (3-55) ？推导,(3-55),(3-56),主体流速,压力降,范宁摩擦因子,(3-58),三、旋转黏度计的测量原理 P61,两垂直的同轴圆筒，内筒的直径为a, 外筒的直径为b, 在两筒的环隙间充满不可压缩流体。当内筒以角速度 、外筒以角速度 旋转时，将带动流体沿圆周方向绕轴线作层流流动。若圆筒足够长，端效应可以忽略。,连续性方程简化,运动方程简化,B.C.,通解：,方向上的剪应力与形变速率的关系：,代入速度分布方程,通常，旋转粘度计的内筒固定不动（ ），故作用于外圆筒内壁上的剪应力为：,(3-70),已知旋转粘度计圆筒长为L，则作用于外筒内壁上的摩擦力为,外筒绕轴旋转的力矩为,(3-74),测定某未知粘度的液体时，规定外圆筒转速 ，测定相应的转动力矩 ，可由上式计算待测液体的粘度。,(3-74),一、方程的简化,二、方程的求解,3 非稳态流动的求解（略）,非稳态流动的求解的例靠近平板的无限大流体当平板突然运动时的流动分布。,一、方程的化简,考察位于平板（ xz 平面）上方的半无限大流体，流体初始静止。在 t= 0 时刻，平板以恒速 u0 沿 x 正方向运动，求流速分布。,设 x 方向压力梯度为零，流动为层流。,由于,由连续性方程和运动方程化简可得,I.C. 0，ux0；(所有 y),B.C. y = 0，uxu0；(所有 0),y = ，ux0；(所有 0),令,B.C.（1）,，,B.C.（2）,两次积分得,速度分布, 误差函数。,二、方程的求解,剪应力,任一瞬时，壁面上的剪应力为,4 极慢黏性流动（爬流）P62,奈维-斯托克斯方程,Case 1：粘性力 惯性力，则可忽略惯性力,爬流（蠕动流）：流速非常低的流动,e.g. 细粒子在流体中的自由沉降、气溶胶粒子的运动以及某些润滑问题,雷诺数,Case 2：惯性力 粘性力，则可忽略粘性力,势流：理想流体的无旋流动,e.g. 流体绕过沉浸物体流动,Re1, CD,雷诺数非常大,(2-46),一、爬流的概念与爬流运动方程 P62,惯性力,粘性力,重力,压力,本节求运动方程的物理近似解。,流体流动时，起支配作用的是惯性力和黏性力。,Re =惯性力/粘性力,本节讨论 Re 1的流动，称为爬流 ( Creeping flow) 。,（1）当流动的Re很大，惯性力黏性力，惯性力起主导作用，黏性力是次要因素。,（2）当流动的Re很小，惯性力黏性力，粘性力起主导作用，惯性力是次要因素。,略去运动方程中的惯性力和重力项：,方程的特点：线性偏微分方程组，4个方程，4个未知量，可直接求出解析解。,(3-76),(2-21),球粒子在流体中的沉降,半径为r0的球粒子在静止无界不可压缩流体中以u0运动。设流动 Re很小。,二、斯托克斯定律,1.流体的速度分布；,2.压力分布；,3.粒子沉降的阻力。,（1）方程简化,连续性方程：,稳态,轴对称,运动方程,B.C.,（2）速度及压力分布,(3-79a),(3-79b),(3-79c),球坐标本构方程,（3）流动阻力,故,球面上,则,球面上,代入本构方程,代入速度和压力分布方程,本构方程用于球面上,流动阻力（斯托克斯方程）,阻力系数,适用条件,Re 1,(3-83),例：根据粒子沉降的终端速度测定流体粘度。,解：,球粒子在粘性流体中降落时，它将一直加速，直至获得一恒定的终端速度为止。当达到这一状态时，作用于粒子上的合力必为零。因此,式中，r0为球粒子半径，s 为球粒子密度，为流体密度。,则可解得,当且仅当Re 1时适用。,(3-83),曳力：,流体对物体施加的总曳力,远离物体表面的流体速度,物体表面的受力面积,曳力系数,绕流流动与曳力系数 P48,一、理想流体的运动方程,二、流体的旋度与速度势函数,三、势流的速度与压力分布,5 理想流体的势流流动 P65,当流动 Re 很大，粘滞力惯性力，惯性力起主导作用，除壁面附近的流体层外，大部分区域可按理想流体考虑。,研究理想流体流动的学科称为理论流体动力学，在航空、航天、水利工程等领域应用广泛。如研究流体绕过沉浸物体的流动的问题时，理想流体的理论可以用来解决压力分布等问题。,一、理想流体的运动方程,N-S方程,连续性方程,一、理想流体的运动方程,方程特点：非线性偏微分方程组，4个方程，4个变量。为求解析解，可将方程转化为线性方程。为此，引入势函数的概念。,二、流体的旋度与速度势函数 P66,流体的旋度,，为有旋流动。,当 ，为无旋流动；,以二维流动（x，y方向）为例讨论：,无旋流动时，,0,0,速度势函数,令,积分,令C=0,定义,二、流体的旋度与速度势函数,引入势函数的目的是将速度变量 用一个变量 代替，从而使方程的求解得以简化。,速度势函数存在的唯一条件是：流动无旋。因此，在三维流动中，也存在相应的速度势函数,速度势函数(续),关系：,三、势流的速度与压力分布,势流 理想流体的无旋流动。,势流的求解：,B.C.,1.速度分布,2.压力分布,x 方向的欧拉方程：,令,x 方向:,y 方向:,z 方向:,无旋流动,取坐标 x，y 水平，z垂直向下,重力场为有势场，令单位质量流体所具有的势能为,积分,Bernoulli equation,理想流体作无旋流动时，动能、位能与压力能之和为常数，不产生流动阻力。,6 平面流与流函数的概念 P69,一、平面流,许多流动体系，其一个方向的尺度要比另外两个方向的尺度大得多，例如矩形管道、流体在一个很宽的平壁面的流动等。对于这类问题，由于流体的物理量在一个方向上无变化或变化很小，可将其按二维流动处理。,稳态不可压缩流体的平面流的变化方程：,二、流函数,不可压缩流体的平面流动的连续性方程为,令,积分,令C = 0,定义, 流函数,流函数满足连续性方程，证明？,思考题 试证明：不可压缩流体作二维平面无旋流动时，流函数满足拉普拉斯方程：,二、流函数,习 题,1. 不可压缩流体在相距为 2y0 的两平行平壁间作一维稳态层流流动。若上板以 u0 的速度移动 ，下板静止不动，试导出其速度分布、剪应力、体积流率和平均速度的表达式。,习 题,2. 粘性流体沿垂直圆柱体的外表面以稳态的层流液膜向下流动，如本题附图所示。试求该流动的速度分布。该液体的密度和粘度分别为 和 。,习 题,3. 如本题附图所示，两平行的水平平板间有两层互不相溶的不可压缩流体，这两层流体的密度、动力粘度和厚度分别为 和为 ，设两板静止，流体在常压力梯度作用下发生层流运动，试求流体的速度分布。,4. 半径为 r0 的无限长圆柱体以恒定角速度 在无限流体中绕自身轴作旋转运动。设流体不可压缩，试从一般柱坐标系的运动方程出发，导出本流动问题的运动方程，并求速度分布与压力分布的表达式。,习 题,5. 试求与速度势 相对应的流函数 ，并求流场中点（2，5）的压力梯度（忽略质量力）。,