基于fluent玉米烘干塔内温度仿真研究.docVIP

基于fluent玉米烘干塔内温度仿真研究 　　摘 要：本文通过分析玉米烘干塔的结构及工作原理，利用fluent软件对烘干塔内传热过程进行仿真分析，利用fluent完成建模、网格化及边界值的设定，从而生成了烘干塔内部温度分布图，它表明高温热流从入料口流失较大，低温空气流可有效降低玉米表面温度，模拟结果可为进一步优化烘干塔的设计提供依据。 　　关键词：玉米烘干塔；fluent；有限元 　　中图分类号：TS210 文献标识码：A DOI：10.11974/nyyjs.20150932013 　　引言 　　吉林省是全国的粮食大省，同时也是世界三大黄金玉米带之一，其玉米产量占全省粮食产量的70%以上， 由于我省玉米田间收获时间一般从9月中旬到10月中旬，除部分小水分玉米可直接脱粒外，大部分是在温度达到冰点下才开始脱粒[1]。这就造成玉米的含水量偏高，在粮食的存储运输过程中易造成霉变，从而影响玉米的品质，所以玉米的烘干变得极为重要。同时，因面临的成本压力，多数企业会人为提高粮食烘干温度，以加快生产速度.但玉米籽粒的温度超过333K时会造成玉米的蛋白质变性和淀粉糊化等问题[2]。因此有必要对烘干塔内的温度分布进行仿真研究。 　　Fluent是目前国际上比较流行的商用CFD软件包，它具有丰富的物理模型、先进的数值方法和强大的前后处理功能，可进行多耦合场分析，可直观的表达玉米烘干塔内热质交换过程，为进一步改进烘干设备提供相应的支持。 　　1 仿真模型设计方案 　　1.1 烘干塔模型搭建 　　根据长春吉大科学仪器有限公司提供的常规8吨玉米烘干塔工艺参数如表1所示，烘干塔设计高度为20m，宽度为7m，而单个玉米籽粒与玉米芯的形状多呈不规则圆锥型，为了准确模拟玉米籽粒及玉米芯的几何形状，利用数字工业相机对多个玉米剖面进行实体拍摄，并且利用0.02mm游标卡尺进行玉米剖面的测量，间接计算获得玉米籽粒及玉米芯的平均三维尺寸[3]，其长度约为20cm，直径为7cm，玉米与烘干塔之间数量级差距约为10-3，其差距过大，因此将几何模型成比例缩小。按缩放比例1：10进行建模。建设模型如图1所示。 　　其中玉米烘干塔高温入口温度为65℃，中温入口温度为45℃，低温入口温度为20℃，冷却温度为10℃。玉米从玉米塔顶部投入，玉米塔底部流出。 　　1.2 烘干玉米模型搭建 　　因玉米模型与玉米塔一样按缩放比例为1：10进行缩放，图2为简化玉米进料模型，其为图1烘干塔模型顶部的长方形平面放大图。在一个长方形平面上排列着多个玉米模型，玉米排列顺序为乱序。在分析当中，长方形平面将会从烘干塔的顶端竖直平移到烘干塔底部，进而模拟进料过程。 　　2 fluent分析 　　2.1 几何构建及网格划分 　　因烘干塔的3D几何模型为平面对称模型，其主要观察结果在对称平面切面上，并且采用2D几何模型可以大大缩减网格数目从而加快计算时间。因此2D几何模型能满足本有限元分析的要求。本次分析中利用Gambit软件进行几何构建及网格划分。按CAD中每个点的坐标在Gambit中建立点，通过点建立线，通过线建立面，划分的几何单元为Tri，划分网格的Spacing数值为9，划分好的网格模型如图3所示。 　　2.2 边界条件设定 　　根据玉米烘干塔的技术模型，高温入口边界条件设为inlet1，中温入口边界条件为inlet2，低温入口边界条件为inlet3，冷空气入口边界条件为inlet4，其中：Inlet1、inlet2、inlet3、inlet4均设为质量流量入口（MASS_FLOW_INLET）；玉米入口设为outlet1，玉米出口为outlet2，其中：Outlet1、outlet2均设置为压力出口（PRESSURE_OUTLET）；玉米模型边界设置为airwall；烘干塔模型其余边界一次设置为wall。 　　2.3 CFD基本模型算法设定 　　选择fluent基于压力求解器求解矢量形式的连续方程、动量方程、能量方程、压力方程。其中，平面设置选择2D平面非稳态时间分析，其余默认。打开能量方程（Energy Equation），流体模型选择标准k-ε模型，选择该模型是因为，该模型在模拟近壁面边界层、自由剪切和低雷诺数流动时性能更好、在强流线弯曲、旋涡和旋转有好的表现。 　　2.4 材料选择 　　因玉米烘干塔内部流体选择是air，设置空气为不可压缩流体（ideal-gas），热值为1006J/Kg，热传导系数为0.0242w/m-k，粘度为1.7984e-5kg/m-s。分子量为28.966。玉米模型固体选择Wood，密度为700kg/m3，热值为2310 J/Kg，热传导系数为0.5 J/Kg。 　　2.5 边界条件设置 　　打开重力设置，将边界wall