旋流导叶式生物质热解反应器内气固两相涡旋流动特性.doc

《旋流导叶式生物质热解反应器内气固两相涡旋流动特性.doc》由会员分享，可在线阅读，更多相关《旋流导叶式生物质热解反应器内气固两相涡旋流动特性.doc（20页珍藏版）》请在第一文库网上搜索。

1、旋流导叶式生物质热解反应器内气固两相涡旋流动特性摘要：针对生物质现有热解反应器存在的传热效果差、产物与催化剂不能及时分离导致产率降低等问题，采用计算流体力学方法与脉冲示踪实验测量法对旋流导叶式反应器内的固相体积分数梯度分布、停留时间分布等进行了研究。结果表明气固两相在新型反应器不同区域不同强度涡旋流场内哥氏力、离心力的作用下达到了较好的混合接触和分离效果，为实现反应器内的反应分离一体化过程提供了前提。研究发现，当催化剂粒径为10m、剂气比为10.5时，反应器达到最佳混合与分离效果。催化剂颗粒停留时间分布曲线呈平滑的单峰分布，拖尾较小，平均停留时间为1.0551.235s，催化剂颗粒流动接近平推

2、流。引言生物质催化热解是改善生物油品质、高效利用生物质资源的有效方法之一，生物质催化热解技术操作简便，只需在常压条件下即可实现，具有广阔的发展前景13。目前研究者普遍认为催化热解过程可以分成两步：一是生物质原材料在高温下生成热解气，二是在催化剂作用下热解气（主要是有机蒸气）进一步裂解成较小的分子，即发生催化裂解4。基于上述两步法发展了生物质原位催化热解和非原位催化热解两种技术，非原位催化热解与原位催化热解相比，由于可以根据两步反应特点有针对性地设计不同的反应器类型，且对于工艺条件改进和结构优化的研究可以根据不同反应特性分别开展，使其更具针对性、更加具体有效，因此非原位催化热解成为目前具有潜力和

3、研究价值的新技术56。目前文献报道的生物质催化裂解的常用反应器类型主要包括固定床反应器和流化床反应器79。在热解蒸气催化过程中，固定床反应器因结构简单、运行费用低，且对催化剂的强度要求不高，应用较多。但其暴露的缺点也越来越受到国内外研究学者的重视，包括反应器内部温度分布不均匀，传热差以及原料处理能力弱和催化剂无法更换等问题10，因此目前常用流化床反应器代替固定床反应器。与固定床反应器相比，流化床反应器易实现连续操作，且生产能力强，便于放大，但因存在环-核非均匀流动、流化介质与催化剂的分离以及逆重力场输送运行成本较高等问题1112，在生物质催化热解中尚未广泛应用和推广。文献1314提出将涡旋流动

4、场内反应分离一体化技术应用于生物质非原位催化热解工艺，其原理是：在涡旋流动场内生物质热解蒸气可通过旋转获得持续、稳定和可控制的离心力场，其巨大的剪切力和快速更新的相界面强化了热解蒸气这一气相与催化剂固相之间的相对速度和相互接触，从而实现高效的传质传热过程和化学反应过程，能够有效提高催化反应转化率，且离心力场被用于相间分离也是常用技术1516。基于上述指导思想，本文在顺流式旋风分离器的基础上设计一种新型旋流导叶式下行床反应器17，以提高催化裂解反应效率，并采用实验研究与数值模拟1820相结合的方法，对旋流导叶式下行床反应器内气固两相涡旋流场进行分析，研究不同催化剂粒径、剂气比对反应器内气固梯度分

5、布特性、分离特性以及催化剂停留时间分布的影响规律。1实验装置该反应器结构如图1所示，入口到导向叶片上端为接触主反应区，在两个气相切向入口的引流作用下，该反应区形成弱旋流场，剪切力和快速更新的相界面强化了热解气与催化剂之间的相互接触，能够有效提高催化裂解转化率；导向叶片下部到排气管上端为反应分离耦合区，此时因存在大量气相热解产物，受到导向叶片的导流造旋作用，在强旋流场内主要进行的是各产物和催化剂因密度不同的分离过程，以及少量的催化反应和过裂解等副反应。实验装置主要由喂料系统、催化剂循环装置、反应器本体及测量系统4部分组成，如图2所示。实验装置采用有机玻璃材质进行冷态实验，可以对反应器内三维速度场

6、及浓度分布进行测量，测量系统利用粒子图像测速法（Particle image velocimetry，PIV）技术2122并配置七孔球探针测速系统和PV-6A型颗粒浓度测量仪23。生物质催化裂解这种非均相反应体系中，两相接触效果越好，反应进行得越完全，但是时间过长，副反应则增多。而催化剂的停留时间分布（esidence time distribution，TD）是评价反应器内气固相间接触时间的重要参考依据，对于正确把握物料在反应器内的混合情况和反应时间、增加产物收率、提高反应选择性、反应器过程强化与放大具有重要意义2426。本文采用高速摄影、脉冲示踪法结合图像灰度化处理技术对停留时间进行实验测

7、量27。停留时间的具体测量方法为：使用染色剂将部分催化剂进行染色作为示踪剂，在装置运行平稳后由入口处一次性注入少量示踪剂，同时在反应器出口利用高速摄像机获取图像监测示踪剂浓度，然后将其随时间变化处理为每隔0.025s的逐帧图像，因采用黑色染色剂，故可采用Matlab将其处理为灰度图像，从而得到各图像的像素特征值，即对应于出口处各时间点的示踪剂浓度变化情况，进而最终得到TD曲线。停留时间分布函数指的是在某一停留时间下出口处示踪剂含量占示踪剂总量的比值，定义为2几何与数学模型2.1几何模型与网格划分模拟采用的下行床反应器的几何结构如图1所示，计算坐标系采用直角坐标系，坐标原点设在接触主反应区顶端中

8、心，Z轴沿轴向向下。根据结构化网格优于非结构化网格的原则，除入口和导流锥附近区域结构复杂采用四面体网格外，反应器其余区域均采用六面体网格，在验证了网格无关性基础上，最终确定网格节点数为33万，图3是反应器网格示意图。2.2湍流模型采用NGk-28湍流模型，NGk-模型与标准k-模型相比，修正了湍动黏度，考虑了平均流动中的旋转及旋转流动情况，故NGk-模型可以更好地处理高应变率及流线弯曲程度较大的流动。该模型中k和的输运方程分别为2.3数值求解和边界条件以计算流体力学软件Fluent为平台，使用有限体积法建立离散方程，除压力梯度项的离散采用PESTO29格式外，控制方程其他各项的离散均采用QUI

9、CK差分格式，并采用SIMPLE30算法进行数值求解。在下行床反应器内气体与催化剂形成气固两相流，并将催化剂视为拟流体，两相模型选择混合模型，反应器入口设为均匀速度入口，出口为充分发展的管流条件，气固无滑移边界条件31，其他参数设置如表1所示，其中入口速度和体积分数由实验测得的物料质量流量换算得到。3计算结果及分析3.1模拟结果准确性验证将空气与催化剂两相冷态流场模拟中的固相体积分数计算结果与利用PV-6A型颗粒浓度测量仪测量结果进行对比。如图4所示，图中r为同一截面以中心位置为圆心的径向距离，为半径，r/表示径向不同位置。由图可见，在z=300mm和z=800mm截面，实验值与模拟值基本重合

10、，吻合较好。总体来看，模拟结果与实验结果有较好的一致性，证明采用该数值模拟方法对旋流导叶式生物质热解反应器内流场有良好的预测精度。因此在以下的数值模拟中仅将气相物质改为生物质热解蒸气，气相入口温度为500，其他边界条件等设置均不改变。3.2涡旋流动场特性图5分别是反应器接触主反应区和反应分离耦合区横截面的速度矢量分布图。相同点是速度方向均为顺时针，表明反应器内呈现涡旋流动形态，不同点是z=300mm截面的速度明显小于z=800mm截面，这与设计时的初衷也是相符的，接触主反应区的弱旋流场能够增强两相接触效果，而反应分离耦合区的强旋流场则是为了提高分离效率。图6是不同位置截面的催化剂固相体积分数分

11、布云图。从图上可以看出，z=300mm截面各部位的固相体积分数梯度较小，说明催化剂从入口进入反应器后在气相湍动以及固相颗粒之间的碰撞接触作用下迅速扩散到各处，对气相起到较好的催化作用。到达z=800mm截面时，强旋产生的离心力使各组分因密度不同而逐渐分离，图上可见催化剂的固相体积分数梯度显著增大，中心大部分区域基本为零，具体数值将在下文中详细论述。如图7所示，气相由切向入口进入反应器后形成一定的螺旋下行流，经过导向叶片后旋流强度增大且最终由排气口离开反应器，也有少量由排剂口流出。固相催化剂颗粒由轴向入口竖直向下进入后在气流作用下扩散并下移，经过导向叶片后在较强离心力场作用下沿边壁螺旋下行，最终

12、由排剂口离开反应器，实现了催化剂与气体产物的分离。3.3颗粒粒径对固相体积分数分布的影响催化剂颗粒粒径对于其能否有效扩散到反应器内各处以及反应后的分离效果影响非常大。因此将颗粒粒径为5、10、20m的算例分别定义为情况1、情况2、情况3，考察颗粒粒径对固相体积分数分布的影响规律。图8a是5m粒径颗粒在反应器内不同轴向位置截面的固相体积分数分布曲线。由图可见，催化剂由轴向入口进入反应器后z=50mm截面梯度较大，r/=0.4位置处出现高峰，r/=0.5位置处出现区域小高峰，这是由于该处位于竖直催化剂入口正下方，催化剂聚集程度大，而后跟随气流在下移过程中逐渐扩散，梯度有所减小，在弱旋流作用下两侧边

13、壁附近浓度高于中间区域。到达z=500mm截面，即接触主反应区末端时浓度分布明显改善，说明颗粒达到了有效扩散效果。当颗粒粒径为10m时，由图8b可见沿轴向梯度固相体积分数逐渐减小，在z=100mm截面已有明显改善，到达z=500mm截面时催化剂体积分数集中分布在4%8%之间，梯度明显减小，达到了较好的均匀扩散效果。20m颗粒因粒径较大，受到的曳力作用也较大，因此从图8c可见在z=50mm截面催化剂明显集中在轴向竖直入口正下方以及边壁附近，与中间部位浓度差别较大，而后催化剂颗粒的分布情况基本相同，都是边壁附近颗粒集中，中间部位基本为零，这是因为颗粒粒径最大，所受离心力也最大，即使在弱旋流动下也已

14、呈现出较好的分离效果，且由于颗粒粒径较大、聚集程度较高，颗粒间的相互碰撞和磨损也最为剧烈。这种现象对于提高反应器的整体催化热解效率是不利的。选取各轴向截面催化剂体积分数曲线上具有代表性的10个不同位置点计算其不均匀度Mes，公式为如图9所示，整体来看，3种催化剂粒径条件下催化剂的分布都是沿轴向下移不均匀度逐渐减小，但是，从数值上来看，情况2在各截面的不均匀度分布都低于情况1和情况3，优势明显。图10是3种粒径催化剂z=800mm截面的固相体积分数分布对比曲线。情况1中催化剂分离效果最差，中间大部分区域的固相体积分数在2%4%之间，明显高于其他两种粒径。情况2中在z=800mm截面中心大部分区域

15、固相体积分数小于2%，仅在r/0.8区域数值明显增大。情况3与情况2的不同之处在于其中心大部分区域的固相体积分数在04%之间。综合图8、9的分析结果，当颗粒粒径dp=10m时催化剂在接触主反应区达到了较好的均匀扩散效果，保证了与生物质热解蒸气的充分接触，最终达到增强催化作用效果的目的，不均匀度在3种粒径条件下最低；在反应分离耦合区强旋造成的离心力作用下颗粒显示出较好的分离效果，保证了催化剂与反应后的产物快速分离，减少了二次反应的发生，从而能有效提高反应转化率。因此dp=10m为最优催化剂粒径。3.4剂气比对固相体积分数分布的影响通过改变入口气固相速度的方式来实现剂气比（催化剂循环量与进料质量之比，Catalyst to gasratio，CTG）的变化，在情况2的基础上增加情况4和情况5，剂气比以及入口气固相速度如表2所示。如图11a所示，由于剂气比最大，在z=50mm截面可见催化剂明显集中在轴向竖直入口正下方，在z=100mm和z=300mm截面可见反应器内由上到下各截面颗粒体积分数梯度不断减小，在接触主反应区的弱旋流作用下，分布曲线基本为“U”形。到z=500mm截面时，中心部位的固相体积分数均大于2%，边壁处约为9.5%。减小剂气比时，催化剂