降温除尘洗涤塔流动与传热性能的数值模拟及结构优化
传热实验实验报告-传热实验报告
传热实验 一、实验目的 1、了解换热器的结结构及用途。 2、学习换热器的操作方法。 3、了解传热系数的测定方法。 4、测定所给换热器的传热系数K。 5、学习应用传热学的概念和原理去分析和强化传热过程，并实验之。 二、实验原理 根据传热方程 Q=KA△ tm，只要测得传热速率 Q，冷热流体进出口温度和传 热面积 A，即可算出传热系数 K。在该实验中，利用加热空气和自来水通过列管 式换热器来测定 K, 只要测出空气的进出口温度、自来水进出口温度以及水和空 气的流量即可。Q 与自来水在工作过程中，如不考虑热量损失，则加热空气释放出的热量 1Q 得到的热量 Q 应相等，但实际上因热损失的存在，此两热量不等，实验中以 22为准。 三、实验流程和设备 实验装置由列管换热器、风机、空气电加热器、管路、转子流量计、温度计 等组成。空气走管程，水走壳程。列管式换热器的传热面积由管径、管数和管长 进行计算。 实验流程图： 空气进口水进口温度计 温度计列管式 转子流 换热器 转子流量计量计 风机温度计温度计 空气电 调节阀 加热器 传热系数K 测定实验流程图
四、实验步骤及操作要领 1、熟悉设备流程，掌握各阀门、转子流量计和温度计的作用。 2、实验开始时，先开水路，再开气路，最后再开加热器。 3、控制所需的气体和水的流量。 4、待系统稳定后，记录水的流量、进出口温度，记录空气的流量和进出 口温度，记录设备的有关参数。重复一次。 5、保持空气的流量不变，改变自来水的流量，重复第四步。 6、保持第 4 步水的流量，改变空气的流量，重复第四步。 7、实验结束后，关闭加热器、风机和自来水阀门。 五、实验数据记录和整理 1、设备参数和有关常数 换热流型错流；换热面积 0.4 ㎡ 2、实验数据记录 序号风机出口空气流量空气进口温空气出口温度℃水流量水进口温度℃水出口温度℃2 度℃L/h 压强 mHO 读数 m3/h 1 1.61611029.28018.921.9 2 1.61611029.48018.921.9 1 1.61611029.96018.922.4 2 1.61611029.96018.922.3 1 1.61611031.92019.024.8 2 1.61611032.02019.024.9 1 1.61111029.62019.123.0 2 1.61111029.62019.023.0 1 1.6611027.82019.021.3 2 1.6611027.82019.021.3 3、数据处理 空气流量水流量水的算术水的比热 传热速对数平均换热面传热系数K 的平均 序号平均温容 J/ m3/s kg/s率 J/s温度△ t m积 m2K W/m2K值 W/m2K 度℃（ kg·℃） 10.00440.022220.404183278.86736.24790.419.2333 19.1717 20.00440.022220.404183278.86736.48160.419.1101
袋式除尘器的选型
袋式除尘器的选型1．处理气体量的计算 计算袋式除尘器的处理气体量时，首先要求出工况条件下的气体量，并且还要考虑除尘器本身的漏风量。Q=Qs－(273+Tc)×101.324/273Pa×﹙1＋K﹚，Q——通过除尘器含尘气体量，m3/h;Qs——生产过程产生的气体量，m/h;Tc——除尘器内气体的温度，摄氏度；Pa——环境大气压，kpa;K——除尘器前漏风系数。 2.过滤风速的选取 多数反吹风袋式除尘器的过滤风速在0.6-1.3m/s之间。脉冲袋式除尘器的过滤风速在1.2-2.0m/s左右，玻璃纤维袋式除尘器的过滤风速约为0.5-0.8m/s。 3.总过滤面积 S=S1+S2=Q／60V+S2 式中：S——总过滤面积，m2; S1——滤袋工作过滤面积，m2 S2——滤袋清灰部分的过滤面积，m2; Q——通过除尘器的总气体量，m3/h; V——过滤速度，m/min 4.滤袋直可取150—250mm，长度以2-3米。由于清灰强度不大，滤袋寿命较长，一般可达1-3年。过滤风速一般为0.5-0.8m/min,阻力约500-1000pa,除尘器的入口含尘质量浓度
通常不超过3-5m/min.本体阻力大体在50-500pa之间。使用脉冲除尘器时如果滤袋上端带有文氏管导流器，则需要加上50-150pa的阻力损失。这部分阻力是不可忽视的。 5.由于机械振打的振打加速度分布均匀，因此设计中，通常需要停风清灰。但是机械振打袋式除尘器通常是小型设备，不停风清灰的场合也很多。 6.脉冲阀。脉冲阀是脉冲喷吹清灰装置的执行机构和关键部件，主要分直角式和淹没式两类，每类有6个规格接口从20-76(3/4英寸至3英寸)。每个阀一次喷吹耗气量30-600m3/min(0.2-0.6MPa).值得注意的是国产脉冲阀的工作压力直角式阀是0.4-0.6Mpa,淹没式阀是0.2-0.6Mpa.进口产品不管哪一种阀，工作压力范围均是0.06-0.86MPa，两类阀没有承受压力和应用压力高低之区分。 直角式脉冲阀构造的工作原理。直角式脉冲阀的构造。阀内的膜片把脉冲分成前后两个气室，当接通压缩空气时，压缩空气通过节流孔进入后气室，此时后气室压力将膜片紧贴阀的输出口，脉冲阀处于关闭状态。 脉冲喷吹控制仪的电信号使电磁脉冲阀衔铁移动，阀后气室放气孔打开，后气室迅速失压膜片后移，压缩空气通过阀输出口喷吹，脉冲阀处于开启状态。压缩空气瞬间从阀内喷出，形成喷吹气流。当脉冲控制仪电信号消失，脉冲阀衔铁复位，后气是放气孔关闭，后气室压力升高使膜片紧贴阀的出口，
传热系数计算方法
第四章循环流化床锅炉炉内传热计算 循环流化床锅炉炉膛中的传热是一个复杂的过程，传热系数的计算精度直接影响了受热面设计时的布置数量，从而影响锅炉的实际出力、蒸汽参数和燃烧温度。正确计算燃烧室受热面传热系数是循环流化床锅炉设计的关键之一，也是区别于煤粉炉的重要方面。 随着循环流化床燃烧技术的日益成熟，有关循环流化床锅炉的炉膛传热计算思想和方法的研究也在迅速发展。许多著名的循环流化床制造公司和研究部门在此方面也做了大量的工作，有的已经形成商业化产品使用的设计导则。 但由于技术保密的原因，目前国内外还没有公开的可以用于工程使用的循环流化床锅炉炉膛传热计算方法，因此对它的研究具有重要的学术价值和实践意义。 清华大学对CFB锅炉炉膛传热作了深入的研究，长江动力公司、华中理工大学、浙江大学等单位也对CFB锅炉炉膛中的传热过程进行了有益的探索。根据已公开发表的文献报导，考虑工程上的方便和可行，本章根椐清华大学提出的方法，进一步分析整理，作为我们研究的基础。为了了解CFB锅炉传热计算发展过程，也参看了巴苏的传热理论和计算方法，浙江大学和华中理工大学的传热计算与巴苏的相近似。 4.1 清华的传热理论及计算方法 4.1.1 循环流化床传热分析 CFB锅炉与煤粉锅炉的显著不同是CFB锅炉中的物料(包括煤灰、脱硫添加剂等)浓度C p 大大高于煤粉炉，而且炉内各处的浓度也不一样，它对炉内传热起着重要作用。为此首先需要计算出炉膛出口处的物料浓度C p，此处浓度可由外循环倍率求出。而炉膛不同高度的物料浓度则由内循环流率决定，它沿炉膛高度是逐渐变化的，底部高、上部低。近壁区贴壁下降流的温度比中心区温度低的趋势，使边壁下降流减少了辐射换热系数；水平截面方向上的横向搅混形成良好的近壁区物料与中心区物料的质交换，同时近壁区与中心区的对流和辐射的热交换使截面方向的温度趋于一致，综合作用的结果近壁区物料向壁面的辐射加强，总辐射换热系数明显提高。在计算水冷壁、双面水冷壁、屏式过热器和屏式再热器时需采用不同的计算式。物料浓度C p对辐射传热和对流传热都有显著影响。燃烧室的平均温度是床对受热面换热系数的另一个重要影响因素。床温的升高增加了烟气辐射换热并提高烟气的导热系数。虽然粒径的减小会提高颗粒对受热面的对流换热系数，在循环流化床锅炉条件下，燃烧室内部的物料颗粒粒径变化较小，在较小范围内的粒径变化时换热系数的变化不大，在进行满负荷传热计算时可以忽略，但在低负荷传热计算时，应该考虑小的颗粒有提高传热系数的能力。 炉内受热面的结构尺寸，如鳍片的净宽度、厚度等，对平均换热系数的影响也是非常明显的。鳍片宽度对物料颗粒的团聚产生影响；另一方面，宽度与扩展受热面的利用系数有关。根
热物理过程的数值模拟-计算传热学3汇总
四、非线笥问题迭代式解法的收敛性 每一层次上满足迭代法求解的收敛条件+相邻次间代数方程的系数变化不太大（亦即未知量 的变化不太大J多数情形下非线性问题迭代式解法是可以收敛的）。 使相邻两层次间未知量变化不太大的措施： 1欠松弛迭代常用逐次欠弛线迭法（SLUR）：一组临时系数下逐线迭代求解+对所得的解 施以欠松弛，再用欠松弛后的解去计算新的系数，常数，以进入下一层次的迭代。 实施：常把欠松弛处理纳入迭代过程，而不是在一个层次迭代完成后再行欠松弛。 .（n 1）川）.'a n bt n b t p =t p （t p ） a p （先）t p n1） = 7a n b t n b b （1一?）屯t p n） co o a'p t p n 9 、a n bt n b b' a'p -a^ ■, b' = b （^ ）（a p ）t p n），用交替方向线迭代法求解这一方程，就实现了SLUR 的迭代求解。为一般化起见，上式中t n b上没有标以迭代层次的符号（J, GS时不相同）。 2、采用拟非稳态法 前面已指出，稳态问题的迭代解法与非稳态问题的步进法十分相似。对于非线性稳态问题， 从代数方程的一组临时系数进入到另一组临时系数亦好象非稳态问题前进了一个时间层，非稳态问题的物理特性：系数热惯性越大（a； = PM v/也I ），温度变化越慢，仿此，对稳态非线性 问题，可在离散方程中加入拟非稳态项，以减小未知量托两个层次间的变化，即 由 （=a n b -S p：V）t p n。= Ua n bt n b b=（3a n b - S p：V a；）t p n。=二a n bt n b b a p tf Za n bt n b - b - a；t p n） （n 1） t p o Ea n b -S p心V +a p 一直进行到t p,t n b收敛，虚拟时间步的大小通过计算实践确定。 3、采用Jacobi点迭代法 中止迭代的判据（该层次迭代）除前述变化率判据外，还可以规定迭代的轮数，例如规定进 行4-6次ADI线迭代就结束该层次上的计算。此时，用收敛速度低的丁迭代也就起到了欠松弛的作用。 五、迭代法的收敛速度 1收敛速度 对给定的代数方程组（包括是临时系数的情形），采用不同的迭代方法求解时，使一定的初始误差缩小成：?倍所需要的迭代轮数K是不相的
导热系数实验报告
一、【实验目的】 用稳态法测定金属、空气、橡皮的导热系数。 二、【实验仪器】 导热系数测定仪、铜-康导热电偶、游标卡尺、数字毫伏表、台秤(公用)、杜瓦瓶、秒表、待测样品（橡胶盘、铝芯）、冰块 三、【实验原理】 1、良导体（金属、空气）导热系数的测定 根据傅里叶导热方程式，在物体内部，取两个垂直于热传导方向、彼此间相距为h 、温度分别为θ1、θ2的平行平面（设θ1>θ2），若平面面积均为S ，在t ?时间内通过面积S 的热量Q ?免租下述表达式： h S t Q ) (21θθλ-=?? （3-26-1） 式中， t Q ??为热流量；λ即为该物质的导热系数，λ在数值上等于相距单位长度的两平面的温度相差1个单位时，单位时间内通过单位面积的热量，其单位是)(K m W ?。 在支架上先放上圆铜盘P ，在P 的上面放上待测样品B ，再把带发热器的圆铜盘A 放在B 上，发热器通电后，热量从A 盘传到B 盘，再传到P 盘，由于A,P 都是良导体，其温度即可以代表B 盘上、下表面的温度θ1、θ2，θ1、θ2分别插入A 、P 盘边缘小孔的热电偶E 来测量。热电偶的冷端则浸在杜瓦瓶中的冰水混合物中，通过“传感器切换”开关G ，切换A 、P 盘中的热电偶与数字电压表的连接回路。由式（3-26-1）可以知道，单位时间内通过待测样品B 任一圆截面的热流量为 冰水混合物 电源 输入 调零 数字电压表 FD-TX-FPZ-II 导热系数电压表 T 2 T 1 220V 110V 导热系数测定仪 测1 测1 测2 测2 表 风扇 A B C 图4-9-1 稳态法测定导热系数实验装置
2 21)(B B R h t Q πθθλ-=?? (3-26-2) 式中，R B 为样品的半径，h B 为样品的厚度。当热传导达到稳定状态时，θ1和θ2的值不变， 遇事通过B 盘上表面的热流量与由铜盘P 向周围环境散热的速率相等，因此，可通过铜盘P 在稳定温度T 2的散热速率来求出热流量 t Q ??。实验中，在读得稳定时θ1和θ2后，即可将B 盘移去，而使A 盘的底面与铜盘P 直接接触。当铜盘P 的温度上升到高于稳定时的θ2值若干摄氏度后，在将A 移开，让P 自然冷却。观察其温度θ随时间t 变化情况，然后由此求出铜盘在θ2的冷却速率 2 θθθ=??t ,而2 θθθ=??t mc ,就是铜盘P 在温度为θ2时的散热速率。 2、不良导体（橡皮）的测定 导热系数是表征物质热传导性质的物理量。材料结构的变化与所含杂质的不同对材料导热系数数值都有明显的影响，因此材料的导热系数常常需要由实验去具体测定。 测量导热系数在这里我们用的是稳态法，在稳态法中，先利用热源对样品加热，样品内部的温差使热量从高温向低温处传导，样品内部各点的温度将随加热快慢和传热快慢的影响而变动；适当控制实验条件和实验参数可使加热和传热的过程达到平衡状态，则待测样品内部可能形成稳定的温度分布，根据这一温度分布就可以计算出导热系数。而在动态法中，最终在样品内部所形成的温度分布是随时间变化的，如呈周期性的变化，变化的周期和幅度亦受实验条件和加热快慢的影响，与导热系数的大小有关。 本实验应用稳态法测量不良导体(橡皮样品)的导热系数，学习用物体散热速率求传导速率的实验方法。 1898年C ．H ．Le e s ．首先使用平板法测量不良导体的导热系数，这是一种稳态法，实验中，样品制成平板状，其上端面与一个稳定的均匀发热体充分接触，下端面与一均匀散热体相接触。由于平板样品的侧面积比平板平面小很多，可以认为热量只沿着上下方向垂直传递，横向由侧面散去的热量可以忽略不计，即可以认为，样品内只有在垂直样品平面的方向上有温度梯度，在同一平面内，各处的温度相同。 设稳态时，样品的上下平面温度分别为 12θθ，根据傅立叶传导方程，在t ?时间内通过 样品的热量Q ?满足下式：S h t Q B 21θθλ-=?? （1） 式中λ为样品的导热系数，B h 为样品的厚度，S 为样品的平面面积，实验中样品为圆盘状。设圆盘样品的直径为B d ，则半径为B R ，则由（1)式得： 2 21B B R h t Q πθθλ-=?? （2） 实验装置如图1所示、固定于底座的三个支架上，支撑着一个铜散热盘P ，散热盘P 可以借助底座内的风扇，达到稳定有效的散热。散热盘上安放面积相同的圆盘样品B ，样品B 上放置一个圆盘状加热盘C ，其面积也与样品B 的面积相同，加热盘C 是由单片机控制的自适应电加热，可以设定加热盘的温度。
传热学数值计算大作业2014011673
数值计算大作业 一、用数值方法求解尺度为100mm×100mm 的二维矩形物体的稳态导热问题。物体的导热系数λ为1.0w/m·K。边界条件分别为： 1、上壁恒热流q=1000w/m2; 2、下壁温度t1=100℃； 3、右侧壁温度t2=0℃; 4、左侧壁与流体对流换热，流体温度tf=0℃，表面传热系数 h 分别为1w/m2·K、10 w/m2·K、100w/m2·K 和1000 w/m2·K; 要求： 1、写出问题的数学描述； 2、写出内部节点和边界节点的差分方程； 3、给出求解方法； 4、编写计算程序(自选程序语言)； 5、画出4个工况下的温度分布图及左、右、下三个边界的热流密度分布图； 6、就一个工况下（自选）对不同网格数下的计算结果进行讨论； 7、就一个工况下（自选）分别采用高斯迭代、高斯——赛德尔迭代及松弛法（亚松弛和超松弛）求解的收敛性（cpu 时间，迭代次数）进行讨论； 8、对4个不同表面传热系数的计算结果进行分析和讨论。 9、自选一种商业软件（fluent 、ansys 等）对问题进行分析，并与自己编程计算结果进行比较验证（一个工况）。（自选项） 1、写出问题的数学描述 设H=0.1m 微分方程 22220t t x y ??+=?? x=0，0 y=H ，0 1、已知：一块厚度为0.1mm 的无限大平板，具有均匀内热源，q ＝50×103W/m 3，，导热系数K ＝10W/m.℃，一侧边界给定温度为75℃，另一侧对流换热，T f ＝25℃，，h=50W/m 2.℃,求解稳态分布。（边界条件用差分代替微分和能量平衡法），画图。(内,外节点) 2、试以下述一维非稳态导热问题为模型，编写求解一维非稳态扩散型问题的通用程序： 00 00000()()()() L L f x x x x L fL L x x x x T T k s c x x T k h T T W x T k h T T W x T T x τρτ =====???+=????=-+??-=-+?= 其中，x 是空间坐标变量，τ是时间坐标变量，T 是温度（分布），k 是材料的导热系数，s 是内热源强度，ρ是材料的密度，c 是材料的比热，h 0和h L 分别是x 0和x L 处流体与固体壁面间的换热系数，而T f0和T fL 分别是固体壁两侧流体的温度，W 0和W L 是x 0和x L 处(非对流换热)热流密度，T 0(x )是固体壁内初始温度分布。注意k 、ρ、c 、s 、h 0 、h L 、W 0和W L 均可以是温度T 和/或空间坐标x 的函数。 具体要求： 1) 将数学模型无量纲化； 2) 考虑各种可能的边界条件和初始条件组合 3) 提供完整的程序设计说明，包括数学推导过程和程序使用说明 3、对于有源项的一维稳态方程， s dx d T dx d u dx d +=)()(φφρ 已知 x=0,φ=0,x=1, φ=1.源项S=0.5-X 利用迎风格式、混合格式、乘方格式求解φ的分布. 实验三 对流传热实验 一、实验目的 1.掌握套管对流传热系数i α的测定方法，加深对其概念和影响因素的理解，应用线性回归法，确定关联式4.0Pr Re m A Nu =中常数A 、m 的值； 2.掌握对流传热系数i α随雷诺准数的变化规律； 3．掌握列管传热系数Ko 的测定方法。 二、实验原理 ㈠ 套管换热器传热系数及其准数关联式的测定 ⒈ 对流传热系数i α的测定 在该传热实验中，冷水走内管，热水走外管。 对流传热系数i α可以根据牛顿冷却定律，用实验来测定 i i i S t Q ??= α （1） 式中：i α—管内流体对流传热系数，W/(m 2?℃)； Q i —管内传热速率，W ； S i —管内换热面积，m 2； t ?—内壁面与流体间的温差，℃。 t ?由下式确定： 2 2 1t t T t w +- =? （2） 式中：t 1，t 2 —冷流体的入口、出口温度，℃； T w —壁面平均温度，℃； 因为换热器内管为紫铜管，其导热系数很大，且管壁很薄，故认为内壁温度、外壁温度和壁面平均温度近似相等，用t w 来表示。 管内换热面积： i i i L d S π= （3） 式中：d i —内管管内径，m ； L i —传热管测量段的实际长度，m 。 由热量衡算式： )(12t t Cp W Q m m i -= （4） 其中质量流量由下式求得： 3600 m m m V W ρ= （5） 式中：m V —冷流体在套管内的平均体积流量，m 3 / h ； m Cp —冷流体的定压比热，kJ / (kg ·℃)； m ρ—冷流体的密度，kg /m 3。 m Cp 和m ρ可根据定性温度t m 查得，2 2 1t t t m +=为冷流体进出口平均温度。t 1，t 2， T w ， m V 可采取一定的测量手段得到。 ⒉ 对流传热系数准数关联式的实验确定 流体在管内作强制湍流，被加热状态，准数关联式的形式为 n m A Nu Pr Re =. （6） 其中： i i i d Nu λα= ， m m i m d u μρ=Re ， m m m Cp λμ=Pr 物性数据m λ、m Cp 、m ρ、m μ可根据定性温度t m 查得。经过计算可知，对于管内被加热的空气，普兰特准数Pr 变化不大，可以认为是常数，则关联式的形式简化为： 4.0Pr Re m A Nu = （7） 这样通过实验确定不同流量下的Re 与Nu ，然后用线性回归方法确定A 和m 的值。 ㈡ 列管换热器传热系数的测定 管壳式换热器又称列管式换热器。是以封闭在壳体中管束的壁面作为传热面的间壁式换热器。这种换热器结构较简单，操作可靠，可用各种结构材料（主要是金属材料）制造，能在高温、高压下使用，是目前应用最广的类型。由壳体、传热管束、管板、折流板（挡板）和管箱等部件组成。壳体多为圆筒形， 实验三袋式除尘器性能测定 1、实验目的及意义 袋式除尘器又名过滤式除尘器，是使含尘气流通过过滤材料将粉尘分离捕集的装置。采用纤维织物作滤料的袋式除尘器，在工业废气除尘方面应用广泛。本实验主要研究这类除尘器的性能。 袋式除尘器的性能与其结构型式、滤料种类、清灰方式、粉尘特性及其运行参数等因素有关。袋式除尘器性能的测定和计算是袋式除尘器选择、设计和运行管理的基础，是本科学生必须具备的基本能力。本实验要求学生在认真了解实验原理、装置、方法、内容和实验要求的基础上，综合应用已掌握的基本知识和基本技能，自行完成实验方案步骤设计和实验测定记录表设计，独立完成本实验。 通过本实验使学生进一步提高对袋式除尘器结构形式和除尘机理的认识；掌握袋式除尘器主要性能的实验研究方法；了解过滤速度对袋式除尘器压力损失及除尘效率的影响；提高对除尘技术基本知识和实验技能的综合应用能力；并通过实验方案的设计和实验结果分析，加强创新能力的培养。 2、实验原理和方法 本实验是在除尘器结构型式、滤料种类、清灰方式和粉尘特性一定的条件下，测定袋式除尘器主要的技术性能指标，并在此基础上，测定处理气体量Q、过滤速度U f对袋式除尘器压力损失△P和除尘效率η的影响。 2.1处理气量和过滤速度的测定与计算 (a) 动压法测定：测定袋式除尘器处理气体量Q，同时测出除尘器进出口连接管道肿的气体流量，取其平均值作为除尘器的处理气体量，取其平均值作为除尘器的处理气量。 漏风率的计算见《大气污染控制工程》5－44式，pp141。 一般要求除尘器的漏风率小于±5％。 (b) 静压法测定：采用静压法测定袋式除尘器进口气体流量（Q1N），根据在除尘器入口管道系统的测口测得系统入口管道处的平均静压P s，具体计算公式见流 一维稳态导热的数值计算 1.1物理问题 一个等截面直肋，处于温度t ∞=80 的流体中。肋表面与流体之间的对流换热系数为 h =45W/(m 2?℃),肋基处温度t w =300℃，肋端绝热。肋片由铝合金制成，其导热系数为λ=110W/(m ?℃)，肋片厚度为δ=0.01m ，高度为H=0.1m 。试计算肋内的温度分布及肋的总换热量。 1.2数学描述及其解析解 引入无量纲过余温度θ = t?t ∞t w ?t ∞ ，则无量纲温度描述的肋片导热微分方程及其边界条件： 22 20d m dx θθ-= x=0,θ=θw =1 x=H, 0x θ?=? 其中m = 上述数学模型的解析解为：[()] ()() w ch m x H t t t t ch mH ∞∞--=-? ()()w hp t t th mH m ∞?= - 1.3数值离散 1.3.1区域离散 计算区域总节点数取N 。 1.3.2微分方程的离散 对任一借点i 有：22 2 0i d m dx θ θ??-= ??? 用θ在节点i 的二阶差分代替θ在节点i 的二阶导数，得：211 2 20i i i i m x θθθθ+--+-= 整理成迭代形式：()1122 1 2i i i m x θθθ+-=++ (i=2,3……,N-1) 1.3.3边界条件离散 补充方程为：11w θθ== 右边界为第二类边界条件，边界节点N 的向后差分得：1 0N N x θθ--= ，将此式整理为 迭代形式，得：N 1N θθ-= 1.3.4最终离散格式 11w θθ== ()1122 1 2i i i m x θθθ+-= ++ (i=2,3……,N-1) N 1N θθ-= 1.3.5代数方程组的求解及其程序 假定一个温度场的初始发布，给出各节点的温度初值：01θ，02θ，….，0 N θ。将这些初值代 入离散格式方程组进行迭代计算，直至收敛。假设第K 步迭代完成，则K+1次迭代计算式为： K 11w θθ+= () 11 11 2212i i K K K i m x θθθ+-++= ++ (i=2,3……,N-1) 1 11N K K N θθ-++= #include 化工原理实验报告 实验三 传热膜系数测定实验 实验日期：2015年12月30日 班级： 学生姓名： 学号： 同组人： 报告摘要 本实验选用牛顿冷却定律作为对流传热实验的测试原理,通过建立不同体系的传热系统,即水蒸汽—空气传热系统、分别对普通管换热器和强化管换热器进行了强制对流传热实验研究。确定了在相应条件下冷流体对流传热膜系数的关联式。此实验方法可以测出蒸汽冷凝膜系数和管内对流传热系数。采用由风机、孔板流量计、蒸汽发生器等组成的自动化程度较高的装置，让空气走内管，蒸汽走环隙，用计算机在线采集与控制系统测量了孔板压降、进出口温度和两个壁温，计算了传热膜系数α，并通过作图确定了传热膜系数准数关系式中的系数A 和指数m （n 取0.4），得到了半经验关联式。实验还通过在内管中加入混合器的办法强化了传热，并重新测定了α、A 和m 。 二、 目的及任务 1.掌握传热膜系数α及传热系数K 的测定方法； 2.通过实验掌握确定传热膜系数准数关系式中的系数A 和指数m 的方法； 3.了解工程上强化传热的措施。 三、基本原理 对流传热的核心问题是求算传热膜系数α，当流体无相变时对流传热准数关 系式的一般形式为：p n m Gr A Nu Pr Re 对于强制湍流而言。Gr 数可忽略，即 n m A Nu Pr Re = 本实验中，可用图解法和最小二乘法计算上述准数关系式中的指数m 、n 和系数A 。 用图解法对多变量方程进行关联时，要对不同变量Re 和Pr 分别回归。本实验可简化上式，即取n=0.4(流体被加热)。这样，上式即变为单变量方程，在两边取对数，得到直线方程为 Re lg lg Pr lg 4.0m A Nu += 在双对数坐标中作图，求出直线斜率，即为方程的指数m 。在直线上任取一点函数值带入方程中，则可得系数A ，即 m Nu A Re Pr 4.0= 用图解法，根据实验点确定直线位置有一定人为性。而用最小二乘法回归，可得到最佳关联结果。应用计算机辅助手段，对多变量方程进行一次回归，就能的道道A 、m 、n 。 对于方程的关联，首先要有Nu 、Re 、Pr 的数据组。其特征数定义式分别为 μρ du = Re ， λμ Cp = Pr ， λαd Nu = 实验中改变空气的流量，以改变Re 值。根据定性温度(空气进、出口温度的算数平均值)计算对应的Pr 值。同时，由牛顿冷却定律，求出不同流速下的传热膜系数值，进而求得Nu 值。 牛顿冷却定律为 Q=αA △t m 式中α——传热膜系数，W/(m 2.℃)； 第四章 复习题 1、 试简要说明对导热问题进行有限差分数值计算的基本思想与步骤。 2、 试说明用热平衡法建立节点温度离散方程的基本思想。 3、 推导导热微分方程的步骤和过程与用热平衡法建立节点温度离散方程的过程十分相似， 为什么前者得到的是精确描述，而后者解出的确实近似解。 4、 第三类边界条件边界节点的离散那方程，也可用将第三类边界条件表达式中的一阶导数 用差分公式表示来建立。试比较这样建立起来的离散方程与用热平衡建立起来的离散方程的异同与优劣。 5．对绝热边界条件的数值处理本章采用了哪些方法？试分析比较之． 6．什么是非稳态导热问题的显示格式？什么是显示格式计算中的稳定性问题？ 7．用高斯－塞德尔迭代法求解代数方程时是否一定可以得到收敛德解？不能得出收敛的解时是否因为初场的假设不合适而造成？ 8．有人对一阶导数()()()2 21,253x t t t x t i n i n i n i n ?-+-≈ ??++ 你能否判断这一表达式是否正确，为什么？ 一般性数值计算 4-1、采用计算机进行数值计算不仅是求解偏微分方程的有力工具，而且对一些复杂的经验公式及用无穷级数表示的分析解，也常用计算机来获得数值结果。试用数值方法对Bi=0.1,1,10的三种情况计算下列特征方程的根:)6,2,1( =n n μ 3,2,1,tan == n Bi n n μμ 并用计算机查明，当2 .02≥=δτ a Fo 时用式（3-19）表示的级数的第一项代替整个级数（计 算中用前六项之和来替代）可能引起的误差。 解：Bi n n =μμtan ，不同Bi 下前六个根如下表所示： Bi μ 1 μ2 μ3 μ 4 μ 5 μ 6 0.1 0.3111 3.1731 6.2991 9.4354 12.5743 15.7143 1.0 0.8603 3.4256 6.4373 9.5293 12.6453 15.7713 10 1.4289 4.3058 7.2281 10.2003 13.2142 16.2594 Fo=0.2及0.24时计算结果的对比列于下表： Fo=0.2 δ=x Bi=0.1 Bi=1 Bi=10 第一项的值 0.94879 0.62945 0.11866 前六和的值 0.95142 0.64339 0.12248 比值 0.99724 0.97833 0.96881 Fo=0.2 0=x Bi=0.1 Bi=1 Bi=10 第一项的值 0.99662 0.96514 0.83889 前六项和的值 0.994 0.95064 0.82925 比值 1.002 1.01525 1.01163 Fo=0.24 δ=x 布袋除尘器 简介及设计说明 限公司 2016年3月 1 / 13 概述 一、设计要求 设计原则 选择行之有效的适应性强、操作灵活、效果稳定、管理简便、节约能耗的工艺处理流程，减少投资和运行费用； 平面布置要求分区明确，便于管理； 管理控制采用集中监测管理、集中控制的方式，对整个烟气处理过程进行监测和控制； 整个工程做到卫生安全、无扰民危害及有效控制和妥善处理为原则，避免造成二次污染。 二、设计依据 设计依据相关标准的规定 JB/T8532-1997 《反吹喷吹袋式除尘器技术标准》 ZBJ88002.2-88 《除尘器性能测试方法》 ZBJ88002.3-88 《除尘器效率测试》 JB/T5917-91 《袋式除尘器应用滤袋框架技术条件》 QB700-88 《国产碳素钢Q235材质标准》 布袋除尘器概述及方案 脉冲袋式除尘器离线清灰方式的有关技术，并借鉴国外先进技术，根据业主要求选择推出的长袋低压脉冲袋式除尘器,该类型除尘器是一种处理风量大、过滤风速低、清灰效果好、除尘效率高、运行可靠、维护方便、占地面积小的除尘设备。广泛应用于电力、冶金、建材、化工等行业的锅炉、烟气除尘及物料回收、粉尘治理。 1、除尘器阻力控制技术：通过烟气流通途径的设计保证除尘器的阻力损失水平。 2、滤料运用技术：根据设备的不同运用场合选用性能价格比较好的滤料。 3、喷吹技术：低压、高效、长寿命膜片电磁脉冲阀的运用，加上喷吹独到设计 和加工手段，使布袋除尘器的清灰方式得到了彻底的改变。 4、检测、监控技术的运用：针对除尘器使用特点，设置了除尘器温度、运行压 力检测、料位检测、运行设备故障检测等先进了在线检测、监控设备。 中国石油大学(华东) 储建学院热能与动力工程系 《计算传热学程序设计》 设计报告 1引言 有关墙体传热量计算的方法是随着人们对房间负荷计算精度要求的不断提高而不断发展的.考虑辐射强度和周围空气温度综合作用，当外界温度发生周期性的变化时，屋顶内部的温度和热流密度也会发生周期性的变化。 计算题目 有一个用砖墙砌成的长方形截面的冷空气通道，其截面尺寸如图1所示。假设在垂直于纸面方向上冷空气及砖墙的温度变化相对较小，可近似地予以忽略。试计算稳态时砖墙截面的温度分布及垂直于纸面方向1米长度的冷量损失。设砖墙的导热系数为(m·℃)。内、外壁面均为第三类边界条件，外壁面：t f1=30℃,h1=10W(m2·℃);内壁面：t f2=10℃, h2=4W(m2·℃)。 图1 砖墙截面 已知参数 砖墙的基本尺寸，砖墙的导热系数，外壁面的表面传热系数，对应的流体温度，内壁面的表面传热系数，对应的流体温度。 2 物理与数学模型 物理模型 由题知垂直于纸面方向上冷空气及砖墙的温度变化相对较小，可近似予以忽略，墙面为常物性，可以假设： 1）砖墙在垂直于纸面方向上没有导热。 2）由于系统是几何形状与边界条件是对称的，它的中心对称面就是一个绝热边界，这时只需求解1/4个对称区域就可以得到整个区域的解。 数学模型 考虑到对称性，取右下的1/4为研究对象，建立如图2的坐标系。 a 图2 砖墙的稳态导热计算区域 由上述的物理模型与上面的坐标系，该问题的数学模型可直接由导热微分方程简化而来，即 22220T T x y ??+=?? (1) 相应的边界条件是： 1.1 0y T y =?=? 1.5 0x T x =?=? (2) 110 ()f x x T h T T x λ ==?-=-? (3) 111.1 1.1 ()f y y T h T T y λ ==?-=-? (4) 22(0.5,00.6)(0.5,00.6) ()f x y x y T h T T x λ =<<=< 套管换热器传热实验实验报告数据处理 我们组做的是实验I : 1， Q=m s1c 1 △t 1 求K 得先求Q Q=m s 1C 1△t 1 ,其中，C 1=所以得先求m s 1 , C 1， △t 1, ◇ 1m s1 =V s1 ρ 要得求V s1，V s1=u 1A ，V s1 =C 0A 0ρρρ/o (2）-gR C 0为空流系数，C 0=0.855,A 0为空口面积，A 0的计算方法如下：A 0 =π4 d 02 , d 0=20.32 mm,故 A 0= π4 ×（20.32 1000 ）2=3.243293×10-4 m 2 R 为压计差读数 A=π4 d 2 ,d 为内管内径=20mm ， 用内插法求解空气密度 ρ 值 这样求得m s 1， ◇ 2 C 1 的求法为先查表的相近温度下空气的C 值，然后用内插法求得对应平均温 度对应的的C 1值 ◇ 3 求△t 1= t △ t 1 ,= t = t 1 + t 2 2 t 1 为进口温度 t 2 为出口温度 进口温度t 1的求解方法 由热电偶中的电位Vt ，按照公式求[]2 000000402.00394645.0t t V E t t ++=得 Et ，再由852.4901004.810608.1105574.15 43-??+?=---t E t 求得t 1值 出口温度t 2的求解方法 由热电偶中的电位Vt ，按照公式[]2 000000402.00394645.0t t V E t t ++=求得 Et ，再由852.49010 04.810608.1105574.15 43-??+?=---t E t 求得t 2值 由以上步骤求出 Q 2 ，由Q=KA △t m 求出K 值 K= Q A △t m Q 由第一步已经求出，A 为内管内径对应的面积，A=2π rL ,r=17.8mm=0.0178 m, A=2×3.14×0.0178×1.224=0.13682362 m 2 3 ,求Re ，Nu 流体无相变强制湍流经圆形直管与管壁稳定对流传热时，对流传热准数关联式的函数关系为： (,,)l Nu f Re Pr d = 对于空气，在实验范围内，Pr 准数基本上为一常数；当管长与管径的比值大于50 时，其值对 Nu 的影响很小；则 Nu 仅为 Re 的函数，故上述函数关系一般可以处理成： m Nu aRe = 式中，a 和 m 为待定常数。 Re=du ρ μ d=2×0.0178 m =0.0356 m , u=Vs/（π×0.01782 ）μ 和ρ用内插法，先查表 的相近温度的μ，ρ，再用线性关系计算求得。 测量空气一侧管壁的中区壁温T W ，由热电偶按前面公式求得；由下式可以计算空气与管壁 传热膜系数测定实验（第四组） 一、实验目的 1、了解套管换热器的结构和壁温的测量方法 2、了解影响给热系数的因素和强化传热的途径 3、体会计算机采集与控制软件对提高实验效率的作用 4、学会给热系数的实验测定和数据处理方法 二、实验内容 1、测定空气在圆管内作强制湍流时的给热系数α1 2、测定加入静态混合器后空气的强制湍流给热系数α1’ 3、回归α1和α1’联式4 .0Pr Re ??=a A Nu 中的参数A 、a *4、测定两个条件下铜管内空气的能量损失 二、实验原理 间壁式传热过程是由热流体对固体壁面的对流传热，固体壁面的热传导和固体壁面对冷流体的对流传热三个传热过程所组成。由于过程复杂，影响因素多，机理不清楚，所以采用量纲分析法来确定给热系数。 1）寻找影响因素 物性：ρ，μ ，λ，c p 设备特征尺寸：l 操作：u ，βgΔT 则：α=f （ρ，μ，λ，c p ，l ，u ，βgΔT ） 2）量纲分析 ρ[ML -3]，μ[ML -1 T -1]，λ[ML T -3 Q -1]，c p [L 2 T -2 Q -1]，l [L] ，u [LT -1]， βg ΔT [L T -2]， α[MT -3 Q -1]] 3）选基本变量（独立，含M ，L ，T ，Q-热力学温度） ρ，l ，μ, λ 4）无量纲化非基本变量 α：Nu ＝αl/λ u: Re ＝ρlu/μ c p : Pr ＝c p μ/λ βgΔT : Gr ＝βgΔT l 3ρ2/μ2 5）原函数无量纲化 ??? ? ???=223,,μρβλμμρλαtl g c lu F l p 6）实验 Nu ＝ARe a Pr b Gr c 强制对流圆管内表面加热：Nu ＝ARe a 圆管传热基本方程： m t A K t T t T t T t T A K Q ???=-----?=111 22112211 1ln ) ()( 热量衡算方程： )()(12322111t t c q T T c q Q p m p m -=-= 圆管传热牛顿冷却定律： 2 2112211 22211221121 1ln ) ()(ln )()(w w w w w w w w T T T T T T T T A t t t t t t t t A Q -----?=-----?=αα 圆筒壁传导热流量：)] /()ln[)()()/ln(11221122121 2w w w w w w w w t T t T t T t T A A A A Q -----?-?=δλ 空气流量由孔板流量测量：54 .02.26P q v ??= [m 3h -1，kPa] 空气的定性温度：t=(t 1+t 2)/2 [℃]计算传热学
对流传热实验实验报告
实验三 袋式除尘器性能测定
传热学上机C程序源答案之一维稳态导热的数值计算
传热膜系数实验报告
传热学第四版课后题答案第四章
布袋除尘器介绍
计算传热学程序设计
套管换热器传热实验实验报告数据处理
化工原理实验传热实验报告