颗粒流体两相流动课件

第四章颗粒流体两相流动 第四章颗粒流体两相流动 1流体与颗粒的相对运动流体与颗粒的相对运动 曳力与曳力系数（Drag and drag coefficient）流体与固体颗粒之间有相对运动时，将发生动量传递。颗粒表面对流体有阻力，流体则对颗粒表面有曳力。阻力与曳力是一对作用力与反作用力。由于颗粒表面几何形状和流体绕颗粒流动的流场这两个方面的复杂性，流体与颗粒表面之间的动量传递规律远比在固体壁面上要复杂得多。爬流(Creeping flow)：来流速度很小，流动很缓慢，颗粒迎流面与背流面的流线对称。流体与颗粒的相对运动 曳力与曳力系数（Drag and dr2曳力与曳力系数（Drag and drag coefficient）在球坐标系中用连续性方程和N-S方程可得到颗粒周围流体中剪应力 r 和静压强 p 的分布为 式中p0为来流压力。流体对单位面积球体表面的曳力（表面摩擦应力）为 曳力与曳力系数（Drag and drag coeffici3曳力与曳力系数（Drag and drag coefficient）r 在 z 轴的分量为 表面曳力(Wall drag)所以整个球体表面摩擦曳力在流动方向上的分量 F 为 z dd曳力与曳力系数（Drag and drag coeffici40 曳力与曳力系数（Drag and drag coefficient）流体静压强对整个球体表面的作用力在流动方向上的分量为 浮力 Fb与流体运动无关流体对颗粒的形体曳力 Fp正比于流速 u形体曳力(Form drag)0 曳力与曳力系数（Drag and drag coeffi5曳力与曳力系数（Drag and drag coefficient）流体流动对颗粒表面的总曳力为摩擦曳力与形体曳力之和 斯托克斯（Stockes）定律 严格说只有在 Rep 0.1 的爬流条件下才符合上式的求解条件 颗粒雷诺数颗粒雷诺数曳力与曳力系数（Drag and drag coeffici6曳力与曳力系数（Drag and drag coefficient）颗粒表面的总曳力 Fd(1)Rep2，层流区(斯托克斯定律区)(2)2Rep500，过渡区(阿仑定律区)(3)500Rep2105，湍流边界层区 边界层内的流动也转变为湍流，流体动能增大使边界层分离点向后移动，尾流收缩、形体曳力骤然下降，实验结果显示此时曳力系数下降且呈现不规则的现象，CD 0.1。曳力与曳力系数（Drag and drag coeffici7曳力与曳力系数（Drag and drag coefficient）曳力系数 CD 与颗粒雷诺数 Rep 的关系 流体绕球形颗粒流动时的边界层分离 曳力与曳力系数（Drag and drag coeffici8自由沉降与沉降速度（Free settling and settling velocity）单颗粒（或充分分散、互不干扰的颗粒群）在流体中自由沉降时在所受合力方向上产生加速度 合力为零时，颗粒与流体之间将保持一个稳定的相对速度。Fd Fg Fb ut 由颗粒与流体综合特性决定，包括待定的曳力系数CD重力场中的沉降速度自由沉降与沉降速度（Free settling and se9自由沉降与沉降速度（Free settling and settling velocity）颗粒-流体体系一定，ut一定，与之对应的Rep 也一定。根据对应的 Rep，可得到不同Rep范围内 ut 的计算式：(1)Rep2，层流区(斯托克斯公式)(2)2Rep500，过渡区(阿仑公式)(3)500Rep ut，up 0，颗粒向上运动；u ut，up 0，颗粒向下运动。自由沉降与沉降速度（Free settling and se11非球形颗粒的几何特征与曳力系数非球形颗粒的几何特征与曳力系数 一般采用与球形颗粒相对比的当量直径来表征非球形颗粒的主要几何特征。等体积当量直径 deV 等表面积当量直径 deA 等比表面积当量直径 dea 颗粒形状系数 非球形颗粒4个几何参数之间的关系 工程上多采用可以测量的等体积当量直径 deV 和具有直观意义的形状系数A。非球形颗粒的几何特征与曳力系数 一般采用与球形颗粒相对比的当12流体通过固定床的流动流体通过固定床的流动 固定床（Fixed bed）：固定不动的固体颗粒层例：固定床催化反应器、吸附分离器、离子交换器等。流体在固定床中的流动状态直接影响到传热、传质与化学反应。颗粒床层的几何特性颗粒床层的几何特性 粒度分布粒度分布 测量颗粒粒度有筛分法、光学法、电学法、流体力学法等。工业上常见固定床中的混合颗粒，粒度一般大于70mm，通常采用筛分的方法来分析颗粒群的粒度分布。标准筛：国际标准组织 ISO 规定制式是由一系列筛孔孔径递增（0.045 mm 4.0mm）的，筛孔为正方形的金属丝网筛组成，相邻两筛号筛孔尺寸之比约为2。由于历史的原因，各国还保留一些不同的筛孔制，例如常见的泰勒制，即是以筛网上每英寸长度的筛孔数为筛号，国内将其称之为目数。流体通过固定床的流动 固定床（Fixed bed）：固定不动13密度函数（频率函数）和分布函数密度函数（频率函数）和分布函数若筛孔直径为 di-1 和 di 相邻两筛的筛留质量为mi，质量分率为xi，则有粒度等于和小于 dpi 的颗粒占全部颗粒的质量分率 混合颗粒粒度分布函数 两函数可相互转换 密度函数（频率函数）和分布函数若筛孔直径为 di-1 和 d14混合颗粒的平均直径混合颗粒的平均直径由于颗粒的比表面对流体通过固定床的流动影响最大，通常以比表面积相等的原则定义混合颗粒的平均直径 dpm。若密度为p的单位质量混合球形颗粒中，粒径为dpi的颗粒的质量分率为xi，则混合颗粒的比表面为 比表面相等对于非球形颗粒，按同样的原则可得 也可用质量平均求混合颗粒的平均直径 混合颗粒的平均直径由于颗粒的比表面对流体通过固定床的流动影响15床层的空隙率、自由截面和比表面床层的空隙率、自由截面和比表面 床层空隙率 颗粒床层中空隙体积与床层总体积之比床层自由截面 颗粒床层横截面上可供流体流通的空隙面积床层比表面 单位体积床层具有的颗粒的表面积床层的空隙率、自由截面和比表面 床层空隙率 颗粒床层中16流体通过固定床的压降流体通过固定床的压降 流体在颗粒床层纵横交错的空隙通道中流动，流速的方向与大小时刻变化，一方面使流体在床层截面上的流速分布趋于均匀，另一方面使流体产生相当大的压降。困难：通道的细微几何结构十分复杂，即使是爬流时压降的理论计算也是十分困难的，解决方法：用简化模型通过实验数据关联。简化的机理模型简化的机理模型 流体通过固定床的压降 流体在颗粒床层纵横交错的空隙通道中流动17流体通过固定床的压降流体通过固定床的压降 把颗粒床层的不规则通道虚拟为一组长为 Le 的平行细管，其总的内表面积等于床层中颗粒的全部表面积、总的流动空间等于床层的全部空隙体积。该管组（即床层）的当量直径可表达为 将流体通过颗粒床层的流动简化为在长为 Le、当量直径 deb 的管内流动，床层的压降 p 表达为 u1 流体在虚拟细管内的流速，等价于流体在床层颗粒空隙间的实际(平均)流速。u1 与空床流速(又称表观流速)u、空隙率 的关系流体通过固定床的压降 把颗粒床层的不规则通道虚拟为一组长为 18流体通过固定床的压降流体通过固定床的压降 工程上为了直观对比的方便而将流体通过颗粒床层的阻力损失表达为单位床层高度上的压降 固定床流动摩擦系数康采尼（Kozeny）式：Reb 2床层雷诺数K 康采尼常数，=5.0 康采尼（Kozeny）方程流体通过固定床的压降 工程上为了直观对比的方便而将流体通过颗19流体通过固定床的压降流体通过固定床的压降 欧根（Ergun）关联式：Reb=（0.17420）可用 A 与 deV 的乘积(A deV)代替 dea。欧根(Ergun)方程流体通过固定床的压降 欧根（Ergun）关联式：Reb=（020流体通过固定床的压降流体通过固定床的压降 当 Reb 2.8(Rep 280(Rep 1000)时，欧根方程右侧第一项可忽略。即流动为湍流时，压降与流速的平方成正比而与粘度无关。流体通过固定床的压降 当 Reb 2.8(Rep 121固体颗粒流态化固体颗粒流态化(Fluidization)Fluidization)流态化(流化床)：颗粒在流体中悬浮或随其一起流动。强化颗粒与流体间的传热、传质与化学反应特性。流态化过程及流化床操作范围流态化过程及流化床操作范围 流态化过程床层压降及床高变化曲线 初始流态化：临界流化速度 umf 临界空隙率mf 颗粒被气流带出：带出速度 u（=ut）流化床操作范围：临界流化速度 umf 与带出速度之间固体颗粒流态化(Fluidization)流态化(流化床22临界流化速度临界流化速度umf是流化床的特性，是固定床变为流化床的一个转折点。可由实验测定的pb u 曲线得到较准确的值。初始流化时，床层内颗粒群（注意不是单颗粒）所受的曳力、浮力与重力相平衡，即流体通过床层的阻力 pb 等于单位床层面积上颗粒所受的重力与浮力之差 因该状态下床层压降也符合欧根方程，将其与上式联立并用（AdeV）代替 dea，可得 临界流化速度umf是流化床的特性，是固定床变为流化床的一个转23临界流化速度临界流化速度当 deV 较大，umf 对应的Rep1000 时，左侧第一项可忽略，注意：计算 umf 的准确程度及可靠范围取决于关联式本身。应充分估计 umf 计算值的误差。最好以实验测定为准。颗粒几何性质及床层 mf 可用经验式估算 当 deV 较小，umf 对应的 Rep10%pb，且0.35mmH2O)可使气体初始分布均匀，以抑制气泡的生成和沟流的发生。多孔板风帽管式内部构件：阻止气泡合并或破碎大气泡。宽分布粒度：宽分布粒度的细颗粒可提高床层的均化程度。床层振动：气流脉动：改善聚式流化质量的措施气体分布板：高阻分布板(10%p27聚式与散式流态化的判断聚式与散式流态化的判断散式流态化（散式流态化（Particulate fluidizationParticulate fluidization）：）：特征：颗粒分散均匀，随着流速增加床层均匀膨胀，床内空隙率均匀增加，床层上界面平稳，压降稳定、波动很小。散式流态化是较理想的流化状态。一般流-固两相密度差较小的体系呈现散式流态化特征，如液-固流化床。聚式流态化（聚式流态化（Aggregative fluidizationAggregative fluidization）：）：特征：颗粒分布不均匀，床层呈现两相结构。即颗粒浓度与空隙率分布较均匀且接近初始流化状态的连续相(乳化相)和以气泡形式夹带着少量颗粒穿过床层向上运动的不连续相(气泡相)。又称鼓泡流态化。一般出现在流-固两相密度差较大的体系，如气-固流化床。聚式与散式流态化的判断散式流态化（Particulate f28聚式与散式流态化的判断聚式与散式流态化的判断气-固流态化与液-固流态化并不是区分聚式与散式流态化的唯一依据，在一定的条件下气-固床可以呈现散式流态化(密度小的颗粒在高压气体中流化)或者液-固床呈现聚式流态化(重金属颗粒在水中流化)行为。散式流态化聚式流态化临界流化条件下的弗鲁德数，D为床径 根据流-固两相的性质及流化床稳定性理论，B.Bomero 和I.N.Johanson 提出了如下的准数群判据：聚式与散式流态化的判断气-固流态化与液-固流态化并不是区分聚29流化床床层高度及分离高度流化床床层高度及分离高度流化数流化数散式流化具有空隙率随流化数均匀变化的规律聚式流化乳化相的空隙率几乎不变，床层膨胀主要由气泡相的膨胀所引起。聚式流化床膨胀比是一个较难确定的参数。床层的流化状态和流化质量与流化数有很大关系实际操作流速与临界流化速度之比 u/umf膨胀比膨胀比 R R流化床的膨胀高度 L 与临界流化高度之比流化床床层高度及分离高度流化数散式流化具有空隙率随流化数均匀30流化床床层高度及分离高度流化床床层高度及分离高度分离高度 H 或 TDH（Transport Disengaging Height）：流化床膨胀高度以上颗粒可以依靠重力沉降回落的高度。超过这一高度后颗粒将被带出。TDH 的确定对流化床气体出口位置的设计具有重要意义。流化床床层高度及分离高度分离高度 H 或 TDH（Trans31广义流态化体系广义流态化体系对高流化数(数百)下的操作，可在床顶设置旋风分离器将随气流带出的颗粒(ut u)回收并返回床内。广义流态化体系：包括密相层、稀相段和颗粒输送段。例：流态化催化裂化装置：原料油高温气化后与催化剂颗粒在提升管内形成高速并流向上的稀相输送，57秒即可完成原料油的催化裂解反应。催化剂经旋风离器分离后由下行管进入再生器，被从底部送入的空气流化再生，停留时间约为712分钟。广义流态化体系对高流化数(数百)下的操作，可在床顶设置旋风分32气力输送气力输送(Pneumatic transport)Pneumatic transport)气力输送：在密闭的管道中借用气体(最常用的是空气)动力使固体颗粒悬浮并进行输送。输送对象：从微米量级的粉体到数毫米大小的颗粒。优点：效率高；全密闭式的输送既可保证产品质量、又可避免粉体对环境的污染；容易实现管网化和自动化；可在输送过程中同步进行气固两相的物理和化学加工(颗粒干燥、表面包裹、气固反应等)。缺点：能耗高，设计和操作不当易使颗粒过度碰撞而磨蚀、破碎，同时造成管道和设备的磨损。气力输送(Pneumatic transport)气力输33气力输送气力输送(Pneumatic transport)Pneumatic transport)气源颗粒进料与加速段稳定输送段气固分离装置气力输送(Pneumatic transport)气源颗34垂直气力输送管内流型 气力输送颗粒气力输送颗粒-流体两相流流动特性与流型图流体两相流流动特性与流型图 颗粒-流体两相的流体动力学特征常表现为流型转变垂直气力输送垂直气力输送 影响参数：气体流速敏感参数：输送管内的压降系统动力消耗评价指标用来表征流型稀相输送与密相输送 均相垂直气力输送流型图 压降最低曲线 密相区稀相区“哽噎”速度 轻微团聚 聚团 节涌 垂直气力输送管内流型 气力输送颗粒-流体两相流流动特性与流型35气力输送颗粒气力输送颗粒-流体两相流流动特性与流型图流体两相流流动特性与流型图 水平气力输送水平气力输送 输送中重力的作用方向与流动方向垂直，使颗粒保持悬浮的不再是曳力、而是水平流动的气流对颗粒产生的升力，因此管内流型（主要是密相）也有所不同。垂直气力输送管内流型 均匀稀相 颗粒堆积“沉寂”速度“沙丘”流 水平“拴塞”水平气力输送流型图“沉寂”速度 气力输送颗粒-流体两相流流动特性与流型图 水平气力输送 36气力输送颗粒气力输送颗粒-流体两相流流动特性与流型图流体两相流流动特性与流型图 气力输送装置的压降包括输送段压降、除尘装置压降和系统内各管件、阀件压降。pf 气体与管壁的摩擦损失pa 颗粒加速所需的惯性压降pr 使颗粒悬浮并上升的重力压降pi 颗粒自身及与管壁的碰撞与摩擦压降直管输送段压降 p气力输送颗粒-流体两相流流动特性与流型图 气力输送装置的压降37两相模型两相模型流体相 把流体和颗粒看作具有相互作用的两相，在微元长度 L 内，分别以流体相和颗粒相为控制体进行动量衡算，得到 颗粒相u,c 气相与颗粒相在管内的平均流速mg,mp 气相和颗粒相在控制体内的质量 若微元管段内的空隙率为，则 颗粒相流体相 两相模型流体相 把流体和颗粒看作具有相互作用的两相，在微元长38两相模型两相模型气相对颗粒相的曳力 Fd：对粒径为 dp 的颗粒 两相滑移速度 流体相摩擦阻力Ff,g：假定管内自由截面分率与 相等，则 而将颗粒相的摩擦阻力 Ff,p 表达为 压降梯度对两相的作用力 Fp,g 和 Fp,p 分别表达为 两相模型气相对颗粒相的曳力 Fd：对粒径为 dp 的颗粒 39两相模型两相模型以上各式中所有动力学参数及颗粒相摩擦系数 p 直接与管内空隙率有关。气力输送中固体加料速率和两相的流速都直接影响空隙率的大小。气力输送加料比颗粒质量流率为G，流体质量流率为w，则粗略估算时常以加料比判断流型，例如有人将 =15 作为密相输送与稀相输送的分界线。实际上，即使加料比相同，两相的物性或流速不同，气力输送管道中固体颗粒的真实体积密度并不一样。或以体积流率之比 来表达则为 两相模型以上各式中所有动力学参数及颗粒相摩擦系数 p 直接40两相模型两相模型在均匀分布条件下，空隙率与颗粒流速 c、气体实际流速 u1的关系为 气固两相间的相互作用力 Fd 是两相模型的核心，目前要预测其大小尚有许多困难，因此限制了两相模型的实际应用。两相模型在均匀分布条件下，空隙率与颗粒流速 c、气体实际流速41拟均相法拟均相法将两相视为一体来考察则不用考虑两相间的相互作用 Fd，使问题得到简化(在稀相输送条件下与实际情况相符)。稳定输送段，颗粒、流体两相的加速度都为零时：水平输送前两项可略去颗粒流速 c 可用 IGT(Institute of Gas Technology)修正式计算：Konno&Saito 计算公式 式中：流体相的摩擦系数 p 颗粒相的摩擦系数拟均相法将两相视为一体来考察则不用考虑两相间的相互作用 Fd42气力输送的类型及装置气力输送的类型及装置 稀相输送（15）负压体系：一般为稀相输送 气力输送的类型及装置 稀相输送（15）和密相输送（43气力输送的类型及装置气力输送的类型及装置 正压体系：低压：100 kPa；中压：300 kPa；高压：1000 kPa 气力输送的类型及装置 正压体系：44气力输送的类型及装置气力输送的类型及装置 组合体系：气力输送的类型及装置 组合体系：45气力输送的类型及装置气力输送的类型及装置 气力输送系统风机的选用气力输送系统风机的选用风机特性曲线与流型图上的压降特性曲线适配。如图，加料量为 G1 时，风机 II 和风机 I 均能满足稀相输送的操作条件。加料量增加到 G2 时，风机 I 已不能满足稀相输送的条件。加料量增加到 G3 时，两台风机都不能满足稀相输送的条件。风机特性曲线越陡峭（如正位移式风机），在稀相输送区操作范围内颗粒加料量的调节余地越大。气力输送的类型及装置 气力输送系统风机的选用风机特性曲线与流46气力输送的类型及装置气力输送的类型及装置 固体颗粒加料器的选用固体颗粒加料器的选用 需考虑颗粒特性（流动性、粘附性、易碎性、大小、形状、温度等），操作压力，是否连续加料以及加料量的控制精度等。旋转阀 文丘里管 气力输送的类型及装置 固体颗粒加料器的选用 47固体颗粒加料器的选用固体颗粒加料器的选用 螺旋加料器板阀 密相脉冲加料固体颗粒加料器的选用 螺旋加料器板阀 密相脉冲加料48