粮食颗粒流动特性的试验研究与数值模拟

1、内蒙古农业大学硕士学位论文粮食颗粒流动特性的试验研究与数值模拟姓名：赵会永申请学位级别：硕士专业：机械设计及理论指导教师：塔娜20090601摘要在粮食颗粒气力输送中，通过研究粮食颗粒的流场分布，得到悬浮速度，对于降低粮食颗粒气力输送中的动力消耗，提高输送的连续性和稳定性有重要意义。从流体力学湍流理论以及颗粒动力学理论出发，考虑气固两相的相互影响、颗粒与颗粒之间以及颗粒与输送管壁之间的相互作用，采用计算流体力学软件建立二维气固两相湍动双流体模型。对垂直管内的气固两相流流场进行数值模拟，获得颗粒浓度及气流流场的分布规律。通过与试验结果比较，表明模拟计算结果的气流速度与试验结果基本吻合。在本模型的

2、基础上，对影响流场的参数（输送气流速度、颗粒粒径、管径）进行数值模拟研究。结果表明，输送气流速度越大，颗粒运动越激烈，颗粒运动形式为中心处向上、近壁面处向下的内循环流动结构：达到稳定输送时，颗粒粒径越小，悬浮速度越大；管径越大，悬浮速度越大。粮食颗粒流动特性的研究，为粮食气力输送机械的设计、优化提供参考依据。关键词：气固两相流；双流体模型；流动特性：数值模拟一），）；九，、）】）（），。），；，；？，。：，；角，；：？：（）（）：）插图和附表清单图试验装置简图图测试系统流程图互童已吼图轴向测点分布图图径向测点分布图图不同加载量的压降轴向分布不同颗粒的压降轴向分布图不同风速的压降轴向分布图图不同

3、风速的压降径向分布图不同高度的压降径向分布加挖坞巧埔墙图结构简图图几何模型图图图求解器多相流模型图湍流模型图材料属性图基本相设置图操作条件图入口边界条件图均加殂毖船孔筋弱出口边界条件残差监视图图控制方程图图颗粒浓度分布图图轴向速度图不同风速的速度矢量图图不同管径的速度矢量图图不同粒径的速度矢量图表试验材料的物理特性表径向测点分布表模型中的系数表试验数据与计算结果比较表（速度）内蒙古农业大学研究生学位论文独创声明本人申明所呈交的学位论文是我本人在导师指导下进行的研究工作及取得的研究成果。据我所知，除了文中特别加以标注和致谢的地方外，论文中不包括其他人已经发表或撰写过的研究成果，也不包括为获得我校

4、或其他教育机构的学位或证书而使用过的材料，与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示谢意。申请学位论文与资料若有不实之处，本人承担一切相关责任。论文作者签名：蒸羔公箍日期：迦：么：尘导师指导研究生学位论文的承诺本人承诺：我的研究生蒸冬金基所呈交的学位论文是在我指导下独立开展研究工作及取得的研究结果，属于我现任岗位职务工作的结果，并严格按照国家及学校有关学术道德规范的规定要求而获得的研究结果。如果违反，我必须接受按国家及学校有关规定的处罚处理并承担相应导师连带责任。指导教师签名：日期：丝罗：受！内蒙古农业大学研究生学位论文版权使用授权书本人完全了解内蒙古农业大学有关

5、保护知识产权的规定，即：研究生在攻读学位期间论文工作的知识产权单位属内蒙古农业大学。本人保证毕业离校后，发表论文或使用论文工作成果时署名单位为内蒙古农业大学，且导师为通讯作者，通讯作者单位亦署名为内蒙古农业大学。学校有权保留并向国家有关部门或机构送交论文的复印件和电子文档，允许论文被查阅和借阅。学校可以公布学位论文的全部或部分内容（保密内容除外），采用影印、缩印或其他手段保存论文。论文作者签名：日垫金基指导教师签名：期：洫弘厶乒内蒙古农业大学硕士学位论文引言据中国国家统计局统计，年全年农业生产继续稳定发展，粮食实现连续五年增产。全年粮食总产量达到万吨，比上年增产万吨，增长，连续年增产。粮食产量

6、的大幅增长有效解决了我国的粮食问题，但同时突出了粮食流通领域的一些问题。实现粮食“四散”（即粮食的散装、散运、散卸、散存）流通技术，可以解决我国粮食流通中存在的一些突出问题，推进我国粮食物流的现代化。农业现代化的重要组成部分之一是粮食流通领域的现代化，而粮食“四散”流通技术是实现粮食物流现代化的关键方法和技术。因此适合我国国情的散装瞳及散卸专用配套设施的实现可以尽快且有效的解决粮食流通中的散运问题。相关研究表明，输送的气流速度是影响粮食气力输送能耗的最主要因素。输送气流速度越大，提升物料的压力损失就越高，从而能耗越高；同时，较高的气流速度也会造成输料管的磨损和物料的破碎。但是气流速度过低时，粮

7、食的输送出现输料管掉料、堵塞现象，而且输送会变得不稳定。气流速度与能耗呈平方关系，气流速度的增加引起能耗的迅速增加。在生产实践中，选择合理的输送气流速度成为气力输送过程中的重点和难点。气力输送设备的广泛应用，已在我国的国民经济建设中占有十分重要的地位。目前气力输送技术已朝着气流斜槽输送、柱塞流静压输送、高压吸送等方向发展。随着气力输送设备的广泛使用，也促成了气力输送理论的产生和发展，而气力输送理论研究，则又是以气固两相颗粒动力学和流体动力学的基本原理为基础。但由于气力输送设备本身具有噪声大、磨损大、动力消耗大等缺点，也限制了它的进一步的应用。为了避免和克服这些不足，在流体力学中两相流被作为一门

8、已有高度发展的分支学科来研究。随着计算机的快速发展，计算流体力学在多相流中的应用越来越受到人们的关注，研究人员对这类反应器中多相流动过程己经进行了大量的试验和模型化研究。计算流体力学模型基于流场中质量、动量和能量守恒规律，建立反映气固两相流动的基本流体力学方程组，加上与具体过程相关的反应过程，最终构成对整个反应器的理论描述。经过试验检验的计算流体力学模型，有利于对反应器更详细的了解，进而指导反应器的设计、优化操作和放大。早在十九世纪，气力输送提供了原始的理论基础。前人采用欧拉坐标和拉格朗日坐标，提出了气固两相流流动状态的公式和方程。进入二十世纪中叶以来，随着颗粒动力学和流体动力学的出现和快速发

9、展，又产生了多种基于湍流动能方程、能量守恒方程、动量守恒方程、连续性方程等守恒定律。这些方程是用以更精确描述气固两相流流动状态的数学模型。在计算流体力学和计算机技术出现后，气力输送理论的发展进入了一个新的阶段。粮食颗粒流动特性的试验研究与数值模拟研究现状气力输送技术是一项综合技术学科，涉及流体力学、材料科学、自动化技术、制造技术等领域，属高新技术项目。其广泛应用于铸造、电力、石油、化工、冶金、建材、。粮食等部门，是一种适合于散粒状物料输送的先进的运输设备。气力输送是全封闭型管道输送系统、布置灵活、无二次污染、高效节能、便于物料输送和回收、为无泄漏输送、气力输送系统以强大的优势，将取代传统的各种

10、机械输送。气固两相流的普遍性和复杂性使它成为一门前沿科学。国内外的气固两相流的相关计算和研究工作较多，其中有化学反应堆研究、流化床技术船、气固两相流风机设计和气力输送佰等，气固两相流在工农业中得到广泛应用，但人们对气固两相流运动规律的认识仍处于不成熟阶段。随着计算机技术的发展，数值模拟在颗粒流动的研究方面得到广泛地应用，目前所采用的方法主要分为两种：一种是连续介质模型，另一种是非连续介质模型。国外研究现状国外众多科学家对于颗粒流动特性的研究相对较早，他们在试验测试的基础上进行建模分析，有的还考虑了流动中颗粒的旋转、碰撞等问题，取得了较大成就。等阳、等睁用颗粒动力学双流体模型对循环流化床提升管内

11、的两相流动进行了二维数值模拟，得到了固体体积分数分布、颗粒速度、浓度的径向分布。等町采用同样的模型对二维流化床进行了模拟，得到了固体分数沿床层变化、气泡长大等详细信息，其预测结果和试验结果相符。等幻用模型和离散单元法对大颗粒流化床的气泡行为进行了模拟，模拟所得到的气泡形成、长大、运动及在床层表面的破碎等情况与试验结果定性地一致。长期以来，用轨道模型研究流体一固体两相流动主要集中在假定颗粒间无相互碰撞的稀相流动方面。在稠密颗粒流体一固体两相流中，颗粒间、颗粒与管壁间的相互作用不能忽视。口和钔讨论了颗粒间碰撞的机理及其重要性。和根据局部分布而推导出的碰撞概率引入了假设环境粒子的串行轨道跟踪法来模拟

12、研究水平管内非稀相流中颗粒间的碰撞。但应用得较多的还是矛叼提出的黏接一滑移碰撞模型，由冲量来计算碰撞后颗粒的速度。¨刀等对气相采用七一双方程模型，并对颗粒相采用描述颗粒间相互作用的硬球碰撞模型，模拟了流化床内气泡的长大、聚并、破碎等现象。用同样的模型对水平管中的气固两相流动特性进行了研究，得到了两相在管中的平均速度分布规律。等¨”对气相采用改进的氐数七一模型，考虑颗粒间的碰撞，用双向耦合真实地模拟了循环流化床上升管中气固两相流动的过程。内蒙古农业大学硕士学位论文国内研究现状我国的散料颗粒学及气力输送技术的研究起步较晚，年代在中科院化冶所郭慕孙院士的倡议下，我国成立了“中国颗

13、粒学会”，我国从此开始研究散料颗粒学及气力输送技术。西安交通大学、同济大学、大连理工大学、浙江大学、中科院化冶所、清华大学等单位在散料颗粒学及气力输送技术方面作了大量的工作啪一明。谢灼利等研究了颗粒粒径和密度对悬浮颗粒的浓度分布的影响心。主要论述对气体相湍动能采用修正的方程模型，建立了水平气力输送的数学模型和相应的计算方法。陆慧林等泓用理论预测得到循环流化床内的“环一核”流动结构，速度的模拟计算结果和颗粒相浓度与试验结果基本吻合。其方法是采用气固两相流动模型，运用气相湍流模型和颗粒湍流模型，模拟计算管内气固两相流动。李晓光等对影响颗粒运动无规则性的主要因素进行了分析幢舶。主要研究在密相大颗粒流

14、化床气固两相流的数值模拟中，借助数值模拟手段，将气相当作连续介质，颗粒相作为离散体系通过提取不同操作条件下大颗粒流化床内颗粒受力的瞬间信息。自世纪年代中期以来，周力行等心研究了湍流两相燃烧的双流体模型。双流体模型的优点是易于给出三维空间颗粒的速度、浓度和温度分布，可以用统一的数值方法求解气固两相方程组，处理颗粒进入和离开回流区没有困难。为了解决用双流体模型模拟颗粒反应经历的问题，周力行等卫提出了湍流两相燃烧的双流体一轨道模型（或颗粒相的连续介质一轨道模型，（，）和全双流体模型（）。两者都用同样的气相欧拉方程组、颗粒相欧拉连续方程和动量方程，包括同样的两相湍流模型。在用双流体一轨道模型模拟多相湍

15、流反应流动时，周力行等妲提出了气相欧拉和颗粒的欧拉一拉氏（）模型算法，称为算法。在这个算法中，除了气相内部的迭代之外，还有气相欧拉计算一颗粒欧拉计算一颗粒拉氏计算三者之间的相互藕合和相互迭代。存在问题分析粮食颗粒的流动特性分析属于气固两相流问题。多年来，针对气固流化系统的基本现象，围绕着了解颗粒流体系统的内在规律，许多研究者开展了大量的研究工作。但由于颗粒流体两相流的复杂结构，使得没有统一的理论或模型能对其流动特性进行完全的描述。实际的颗粒流问题看似简单，但许多机理问题至今还不清楚。颗粒材料范围很广，物理性质差别较大，运动方式不同，其机理也各不相同，虽然对气力输送的研究取得了很大的进展，但是研

16、究气力输送的基础理论等很多方面时，常常会遇到下面这些问题：（）气力输送过程是一个非常复杂的物理过程，管内流场会受到许多不确定性粮食颗粒流动特性的试验研究与数值模拟因素的干扰和影响，管内流场分布需进一步研究。同时，因为其中许多状态量之间存在极为繁琐的非线性关系，所以目前整个气力输送过程的完整运动状态还无法用一个数学上解析的模型来描述和实现。（）国内外众多学者就气固两相流建立起了许多数学模型，并将这些模型应用于输送对象十分广泛的气力输送中，但是就具有一定的局限性。尤其是对其流动状态醮大量研究表明这些模型大都是对应于直径在微米级，或更小范围内的颗粒所建立的。因此有必要对直径在微米级或更大范围的颗粒迸

17、一步研究。选题意义及研究署的气力输送技术是符合现代物料输送和物流工程要求的输送方式之一，与其他机械输送方式螺旋式、刮板式、如斗式输送相比具有许多独特的优点，如结构简单，管路布置灵活方便，造价低，能保证良好的工作环境，连续性好，易于实现自动化等，因此在工农业许多部门得到广泛使目，然而粮食气力输送还存在着一些制约其进一步推广及应用的因素，其中最主要的因素是动力消耗大。长期以来人们已在气力输送领域进行了大量的研究，但仍有许多问题没有得到很好解决。从应用基础研究方面来看，气力输送是一个复杂的多相流动过程，输送管道内粮食颗粒，动状态既有悬浮又有滚动，同时还发生粮食颗粒与壁面、粮食颗粒之间的碰撞，粮食颗粒

18、的旋转还产生举力，完全考虑这些问题是相当复杂的。气力输送越来越多地应用于工农业生产中，同时作为各种散装物料的一种连续输送方式，在整个物料输运中占有很大的比重。对输送风速、不同物种、风压相互关系以及管道压力损失未能进行深入研究，以至无法合理确定气力输送的主要参数。为了有效的控制粮食颗粒气力输送系统，根据生产要求而得到适当的设计参数并制定合理的输送工艺是十分必要的。因此，在粮食颗粒气力输送中，研究粮食颗粒的流动特性，了解流场分布，得到其悬浮速度，对于降低粮食颗粒气力输送中的动力消耗，提高输送的连续性和稳定性有重要意义。研究粮食颗粒的流动特性，不仅可为气力输送装置的设计和粮食颗粒加工过程中干燥提供参

19、考依据；而且，有助于粮食颗粒加工、储存、运输设备向精密化、自动化、大型化的方向发展，所以有必要对粮食颗粒流动特性进行深入的研究。课题研究内容本文以试验为基础，对垂直管道中粮食颗粒的流动特性进行研究，利用。软件对垂直管道中粮食颗粒的流动特性进行模拟。试验和数值模拟相辅相成，成为研究气固两相流动的重要手段。内蒙古农业大学硕士学位论文试验研究的主要内容（）测量所研究物料的物理性质。包括粮食颗粒的散粒体密度、千粒重、真密度，以及测量单颗粒的三轴尺寸用以计算颗粒的体积当量直径。（）通过试验来研究不同加载量、不同风速、不同物料的情况下，压降在轴向分布，分析得出一般的变化规律；研究不同风速、不同截面的情况下

20、，压降的径向分布，分析得出般的变化规律。数值模拟研究的主要内容随着计算机的发展，数值模拟研究得到了广泛的应用，为研究者带来了极大的方便。本文将在的平台上，选择欧拉双流体模型，对垂直管内的粮食颗粒进行简单的气固两相流动模拟试验，通过模拟试验来研究管内气固两相的流动规律及流场特性，其主要内容：（）建立以双流体模型为基础，运用软件气固两相流进行数值模拟。试验数据与模拟结果进行对比，确定模型的准确性。（）分析不同参数对粮食颗粒流动特性的影响。两相流和软件的介绍两相流两相流是两相物质（至少一相为流体）所组成的流动系统。若流动系统中物质的相态多于两个，则称为多相流。两相或多相流是化工生产中为完成相与相传质

21、和反应过程所涉及的最普遍的粘性流体流动。自然界和其他工程领域中两相流也广泛存在，如雨、雪、云、雾的飘流，生物体中的血液循环，水利工程中的泥沙运动和高速掺气水流，环境工程中烟尘对空气的污染等。两相流的流动形态有多种。除了同单相流动那样区分为层流和湍流外，还可以依据两相相对含量（常称为相比）、相界面的分布特性、运动速度、流场几何条件（管内、多孔板上、沿壁面等）划分流动形态。对于管内气液或固系统，随两相速度的变化，可产生气泡流、塞状流、层状流、波状流、冲击流、环状流、雾状流等形态；对于多孔板上气液系可以产生自由分散的气泡、蜂窝状泡沫、活动泡沫、喷雾等形态。两相流研究的一个基本课题是判断流动形态及其相

22、互转变，流动形态不同，则热量传递和质量传递的机理和影响因素也不同。两相流研究的另一个基本课题，是关于分散相在连续相中的运动规律及其对传递和反应过程的影响。当分散相液滴或气泡时，有很多特点。例如液滴和气泡在运动中会变形，在液滴或气泡内出现环流，界面上有波动，表面张力梯度会造成复杂的表面运动等。这些都会影响传质通量，进而影响设备的性能。两相流研究的课题，还有两相流系统的摩擦阻力，系统的振荡和稳定性等。粮食颗粒流动特性的试验研究与数值模拟两相流的实验研究，是掌握两相流规律的基本方法。目前广泛应用光学法（包括光吸收、散射、干涉、折射等），其他辐射吸收和散射法，示踪法，以及电容和电导法等测定两相流中的重

23、要参数，如压力、空隙率、平均膜厚、液滴直径、运动速度等。在某种意义上说，对两相流规律更深入的了解，有赖于实验技术的进步。试验研究、理论分析和数值模拟是当前研究流体运动规律的三种基本方法。试验研究以试验为研究手段，耗费巨大；理论分析是利用简单流动模型假设，给出某些问题的解析解，但对于复杂的非线性流动现象目前还有些无能为力。数值模拟以计算流体力学为基础，采用数值计算方法，通过计算机求解流体流动的数学方程，对具体问题进行数值求解，己经成为研究流体流动的一个重要研究方向和方法。这方面，已经开发出了计算准确、界面友好、使用简单叉能解决问题的大型商业软件，目前比较著名的有，等。本课题将借助于常用的计算流体

24、力学软件一来进行多种工况的数值模拟试验，以此来研究各种因素对粮食颗粒流动特性的影响规律。软件简介年，美国的流体技术服务公司公司的软件部（￥公司的前身），推出了第一个商用软件包。自软件面世以来，以其丰富的物理模型、先进的数值方法及技术人员高质量的技术支持和服务，软件很快成为市场的领先者。年公司正式成立。年，公司并入了在数字和能源电子等领域具有主导地位的公司。的软件设计基于软件群的思想，采用不同的离散格式和数值方法，从用户需求角度出发，针对各种复杂流动的物理现象，在特定的领域内使稳定性、精度和计算速度等方面达到最佳组合，从而高效率地解决各个领域的复杂流动计算问题。基于上述思想开发了适用于各个领域的

25、流动模拟软件，软件之间采用了统一的共同的图形界面及网格生成技术，这些软件能够模拟传热传质、流体流动、化学反应和其它复杂的物理现象。由于应用的工业背景不同使各软件之间有所区别因此大大方便了用户。软件是计算流体力学（）软件的简称，是专门用来进行流场预测、流场计算、流场分析的软件。的使用方法是：首先建立装置的计算模型：然后将流体流动的物理特性应用到虚拟的计算模型；最后软件将输出流体动力特性。是一种高级的分析技术软件，它不仅可以预测流体的行为，同时还可以得到相变（如凝固和沸腾）、化学反应（如燃烧）、机械运动（涡轮机）、传质（如分离和溶解）、传热，以及相关结构的变形（如风中桅杆的弯曲）和压力等的性质。内

26、蒙古农业大学硕士学位论文试验方法与数据处理试验装置图试验装置简图本研究采用的试验装置如图所示。气流经布风板进入垂直管道，粮食颗粒由加料装置倒入，气流作用到颗粒床层产生气固两相流，对粮食颗粒流动特性进行信号采集，空气排入大气，粮食颗粒由回收器收回。以的风机作为风力源，电动机型号（）。布风板的设计、安装有助于更真实、准确的认识管道中颗粒的流动特性。布风板使作用在颗粒上的气流均匀稳定，它放置在风机出口处，与管道相连，可有效改善气流在管道径向的均匀性，滤除输入管道中气流自身所携带的波动。为了使风机的振动不会影响到所研究颗粒的流动，风机和管道使用帆布连接，帆布连接是一种柔式连接方式，它可将风机的振动与管

27、道相隔离，通过帆布连接还可以防止接口处气流的泄漏。为了使承载粮食颗粒的物料网有较高的强度和刚度，在颗粒放置较多或风速较大时不会产生明显变形，选用了钢丝网（目）作为物料网，物料网夹置在法兰之间，安装时用法兰将其夹紧。为避免颗粒从管口顶部飞出，。在管道的顶端安装了挡料装置。粮食颗粒流动特性的试验研究与数值模拟管径的选择工程上输送流体的管道，大多为圆管，设圆管的内径为，则管道的截面积为（）：（除特别指明之外，均指内径），流速为：“七口蒯二式中：一一管道内径，：流体体积流量，：流体在管道内的流速，；拈、匿：糕叫等。（）根据流量和流速，可用式（）算得管道内径，其中流量通常是为生产任务所决定，所以关键在于

28、选择合适的流速。由式（）可知，当流量一定时，流速越大，管径越小，设备费用可减小，但此时流体流速相应增大，其在管道中流动阻力也越大，使操作费用（动力消耗）增加。反之，流速减小，阻力降低，操作费用减少，但管径增大，设备费用增加。设计管道时，尤其是输送距离较长时，需要综合考虑这两个相互矛盾的因素，确定适宜的流速，使操作费用与设备费用之和为最低。在管径的选择时，先根据情况选定流速，再经式（）算出后，从而选用标准管的规格。本试验管径取为来逆行试验。试验介质气体：空气，取室温。，空气的密度为运动粘度为一颗粒：大豆、黑豆、绿豆，颗粒物理特性参数如表所示。表试验材料的物理特性，内蒙古农业大学硕士学位论文当量直

29、径的测量对绿豆、大豆、黑豆进行随机取样，对每一百粒用电子天平称其重量、称十组得其千粒重。并随机在每一百粒中取十粒用游标卡尺测量其三个垂直方向的尺寸长（）、宽（）、厚（），通常；每一种粮食颗粒测量组（即粒）。体积当量直径是指如物体体积与某个球体的体积相等，则这个球体的直径就称为该物体的体积当量直径。将粮食颗粒看作椭球形计算其体积，则：竺（）体积当量直径：拓瓦密度的测量（）粮食颗粒在贮存、加工、输送、流动等过程中，它的流动性、摩擦特性、贮存特性都与密度有密切关系。因此，密度是物料的一个重要物理量，在很多研究过程中都与其有关。由于研究的对象和研究过程的不同，表示密度的方法也不同：（）研究单个颗粒时，

30、其密度可由其质量和占有的体积之比求得，称为颗粒密度。（）研究散粒体群体时，用其所占的容积以及其质量，计算出的密度称为散粒体密度。（）真密度，：颗粒或散粒体的质量除以其固体部分的体积。（）表观密度。：颗粒或散粒体的质量除以除去表面开口，缝隙后固体部分的体积。（）有效密度。：颗粒或散粒体的质量除以包括表面开，缝隙和内部空洞后固体部分的体积。显然，真密度表观密度有效密度。散粒体密度的测量：采用的环刀量取颗粒，然后用电子天平称其重量，多次测量取平均值（本试验中测量次），即可得物料的散粒体密度。真密度的测量：采用浸液法测量物料的真密度，首先在量程为的量筒中加入约的水，读取水的精确体积。：然后，称取的颗粒

31、物料，将其倒入量筒中，再次读取量筒中液面处的刻度：，则：一，即为物料的体积。重复上述操作，多次测量取平均值，通过计算便可得出颗粒物料的真密度。粮食颗粒流动特性的试验研究与数值模拟试验仪器（）直径为的型皮托管（公司，）：（）微电脑数字压力计（压力范围：卜，分辩率：预热时间：）测静压和风速；、（）微差压变送器（公司系列，量程：一）：（）数据采集卡（，分辨率位，多通道采样频率）。（）热球式风速仪（量程：）由以上仪器构成测试系统，系统流程如图所示：通过型皮托管（公司，）获取压力波动信号，由数据采集平台进行数据显示、分析、记录。圈测试系统流程图舄试验的测量点分布垂直管轴向测点分布图轴向测点分布图内蒙古农

32、业大学硕士学位论文在管壁高度上布测个点，原则上布测点越多越精确，但布测点太多，可能会使玻璃破裂或造成裂纹，而且管道漏风更多导致误差增大。因此，考虑有机玻璃管的强度及试验精度根据实际情况，以物料网为参考面，布置个测量点（，）。布测点的高度分别为，单位为米（如图所示）。截面测量选择原理及测点分布由于管道内壁附近附面层的形成和风机的结构，将导致管中风速的分布的不均匀。应用在工程实践中测量平均风速的方法来布置测量点。对于圆形管道，通常采用等环面积法（中华人民共和国机械行业标准）。这就是将圆形断面分成个面积相等的同心圆环，再将每个圆环分为两个面积相等的部分，测点就放在这两部分的分界线上。测点所在圆中的半

33、径分别为，乃，则根据等面积划分的原则可求得：罢吲，例怯一：生，。一，陲（，丛：立，：型一（）其中：一垂直管道半径，卜一圆环圈数：一分割圆环数；一第测量点所在的同心环中心圆半径。根据以上方法，确定本试验垂直管某截面的测点，将截面分成个等面积的环形，取为从而计算得到测点的位置如图所示和表所示，表数据的负值代表方向，负值说明在零点左侧。粮食颗粒流动特性的试验研究与数值模拟图径向测点分布图表径向测点分布表（续）径向测点分布试验结果与分析轴向压力分布（）试验颗粒：绿豆，风速：，加载量分别为，、轴向压力分布如图所示。内蒙古农业大学硕士学位论文一一生量醛避床层商度【）图不同加载量的压降轴向分布通过图可以看出

34、当把风速固定在，在同一高度，压降是随加载量的增大而增大的。分析原因，风速固定，加载量增加，颗粒被流化的几率就增加，曳力就增加，所以压降也增加。无论何种加载量，压降都随高度的增加而减小。由于主直管下半部分颗粒浓度比较大，上半部分浓度逐渐减小，压降也减小。压降在下半部分变化率剧烈，而在下半部分则较平缓。（）试验颗粒：绿豆，加载量：，风速分别为、。轴向测试结果如图所示。床层高度（图不同风速的压降轴向分布由图可知，当加载量不变，风速从增加到勺过程中，同一高度截面上的压降是增加的，而增加到过程压降反而减小。原因是随着风速增加流化床处于由鼓泡流态化向湍流流态化转变过程中，流化的颗粒量增加，同一截面粮食颗粒

35、流动特性的试验研究与数值模拟上的颗粒重量增加，所托起颗粒的曳力也加增加，压降就会增加。风速增大到时，流化床转化为快速流态化，颗粒的速度也会增加，颗粒变稀浓度变小，所托起的颗粒曳力也减小，所以压降也减小。无论何种风速，压降都随高度的增加而减小。颗粒在垂直管下半部分压降变化率剧烈，上半部分比较平缓。（）将大豆，黑豆，绿豆分别为在同一风速（）下轴向压降分布测试的结果如图所示。床层高度（）图不同颗粒的压隆轴向分布：通过图可知，压降都随着床层高度的增大而减小。在同一截面上，颗粒密度越大，压降也增加。在同样加载量下，在较小气流速度下，密度越大，其颗粒数密度越低，因此颗粒与壁面、颗粒与气体及颗粒之间的摩擦减

36、弱，从而消耗的能量越低，导致压降降低。在较大气流速度（）下，压降增加的原因可能是颗粒重力引起压降的显著增大，颗粒密度越大，其惯性越大。粒径增加，压降减小。在悬浮输送情况下，对同样加载量的情况下，颗粒粒径减小，颗粒数密度将增加，因此颗粒与气体相互作用强烈，消耗能量增多，压降自然增大：颗粒粒径减小，颗粒数密度增大，颗粒与壁面之间以及颗粒之间的相互作用随之增大，从而引起阻力增加。径向压力分布径向测试点由等面积同心圆环法（中华人民共和国机械行业标准）选取，垂直管半径，测量。，测点与管径之比为径向相对距离。内蒙古农业大学硕士学位论文颗粒：绿豆，加载量：，截面：高度处，风速分别为，。径向测试结果如图所示。

37、颗粒：绿豆，加载量：，在同一风速（）、不同高度径向测试的结果如图所示蹙向唁对琵盈限营；甸糟对距高限图不同风速的压降径向分布图不同高度的压降径向分布在不同的风速下，颗粒在垂直管某一高度截面（）径向分布如图所示，管中心压降略高于管壁处，风机出口气体风速并不均匀，经整流装置改善了气流分布。由于在壁面处颗粒浓度较小，壁面与颗粒的相互摩擦作用下，近壁面位置处运动的颗粒被减速，故在垂直管壁面附近，对颗粒的加速作用非常弱，颗粒不能被充分地带走。在管中心处颗粒浓度较大，气体能量主要用于颗粒的加速和输送，其它方面的能量损耗很少，所以管中心处压降高于管壁处。如图所示，不同高度截面的压降径向分布规律也是如此，管中心

38、处压降高于管壁面处。粮食颗粒流动特性的数值模拟及分析两相流模拟假设由于气固两相体系流动的复杂性，需要对两相进行一些合理假设。本文基本假设为：（）管内温度为常温，且恒定不变：（）气、固两相均连续，气相为牛顿流体，两相均有确定的物理特性参数；（）粮食颗粒由相同尺寸组成；（）流场内气固两相各占据不同的体积，气体相、颗粒相有自己占据的体积；（）仅考虑颗粒在流场中所受的阻力和重力。粮食颗粒流动特性的试验研究与数值模拟由以上假定，根据牛顿第二定律和守恒定律，可以得到气力输送过程中气固两相流动的守恒方程。建立几何模型数值计算几何模型根据现实中的试验设备尺寸进行合理的简化建立。粮食颗粒流动特性试验的结构简图如

39、图所示。垂直管中的气力悬浮输送，因气体速度一般较高，曳力方向与重力方向平行，气体以及颗粒的运动在切向和径向方向上不受重力的影响，所以可以将其简化为二维来处理，认为切向速度为零。数值计算的流场几何模型如图所示。压办出口壁面速度凡口风图结构简图图”几何模型网格划分利用前处理软件进行网格划分，网格是求解控制方程的基本单元，网格的形状及划分直接影响最终的求解结果。结构化网格节点排列有序、邻点间的关系明确，具有正交性，适用于具有规则几何外形的计算区域。由于几何外形相对简单，而且流动和几何外形很符合（长管），可以使用大比率的四边形单元。这种网格可能会比三角形网格少很多单元。几何外形不是很复杂时，四边形网格

40、可能会实现流动和网格成一条线。只有在简内蒙古农业大学硕士学位论文单的流动（长管流动中），才可以使用四边形网格来减少数值耗散，而且在这种情况下使用四边形有很多优点，因为与三角形网格相比可以用更少的单元得到更好的解。因此采用结构化四边形（长方形或正方形）网格。一般来说，对于同一差分格式，计算精度越高需要的网格越密。但是，网格越密计算量就越大，因此需要在计算精度内台理的细化网格。非均匀分布的网格体系无论从资源占用情况还是从计算结果来看都要优于均匀分布的网格体系。但是还没有一个通用的法则可以说明相邻网格百之间最大（最小）比值究竟应该保持多大。一般在工程计算上，在所求解的变量变化比较剧烈的地区网格分布应该稠密一些。本课题对计算区域进行网格划分，在下部由于颗粒相的存在网格划分比较周密，而在上部和出口处网格划分比较稀疏。读入网格及网格检查（）读入网格（术）（）检查网格在中对模型首先进行网格检查，这一步的目的是为了保证不会出现负体积网格，保证了网格体积为正值时，才可以进行下一步。以下为网格的检查结果。：（）：（），（）（）（）（）：（）：（）：一（）：粮食颗粒流动特性的试验研究与数值模拟掣血，