化工原理指导书-05.doc

化工原理上机指导书中国矿业大学 化工系丁曙光2002.5.14学生上机守则化工学院计算机实验室是提高学生政治思想素质和文化素质的重要阵地，为建立一个健康、安全的上机环境，根据国家的有关法律、法规和学校的有关规定，制定本规定。1、凡进入机房的人员除须遵守校规校纪外，还必须遵守机房的各项管理规定。上机学生必须出示学生证，认真登记上下机时间、机号及上机内容。不得随意更换机位。2、严禁任何人利用计算机网络观看、传播、拷贝、制作淫秽、反动、迷信等不健康内容。3、注意个人卫生，保持机房安静，不要随意走动，以免影响他人上机。4、机房内严禁吸烟，严禁吃零食和乱丢杂物，严禁乱写乱画。5、严禁携带一次性水杯和潮湿物品进入机房，以免发生电故障。6、严禁私接电源、拉线路，严禁乱动电闸和消防器材，严禁使用明火。7、爱护机房内的设备，不得私自拆卸、搬移设备，更不得将机房物品带出。8、未经允许，不准私自携带各种存储盘进入机房，严防计算机病毒。9、严禁在计算机上玩游戏。对违反规定者，机房管理人员有权予以制止。10、上机者必须遵守安全操作程序，不得破坏计算机系统。 11、尊重机房管理员的安排与建议，如有意见，可向有关领导提出。本规定自公布之日起执行，由学院有关部门负责解释。未尽事宜，参照国家及学校相关规定。 化工学院2002年1月1日上机实验说明本课程的上机内容，可以使用任何一种或几种程序语言(包括Excel等)编程实现。其目的：编程训练；加深课程的理解。上机实验报告由源程序(调试前初稿及调试好的补充内容)、实际运行结果、结果分析讨论及教师签字共四部分组成。统一用报告纸，并装订成册后上交评分。 目 录 化工原理上机(1)实验一 粘度的插值求取 (1学时)实验二 绘制速度分布及剪应力分布的曲线 (1学时)实验三 试差法求管径 (2学时)实验四 求平均速度及方次 (2学时)实验五 选择泵型号 (4学时)实验六 K判据法求颗粒沉降速度ut (2学时)实验七 沉降法测定液体粘度 (2学时)实验八 求沉降槽的最大沉降面积 (2学时)实验九 求导热通量及两层砖间界面温度 (2学时) 化工原理上机(2)实验十 画出不同a值下的气液平衡线 (2学时)实验十一 用安托尼公式求苯-甲苯数据 (2学时)实验十二 比较不同q及不同R下的理论板数 (2学时)实验十三 求理想物系精馏塔的理论板数 (2学时)实验十四 求非理想物系精馏塔的理论板数 (2学时)实验十五 操作型精馏塔的计算 (2学时)实验十六 求吸收中的传质单元数 (2学时)实验十七 气体流量对塔径的影响 (2学时)实验十八 多组分吸收的两步图解法 (选做，2学时)实验十九 求泡点温度 (2学时)实验二十 多组分精馏气化率 (2学时)实验一 粘度的插值求取 (1学时)一、 实验目的、原理学习在给定的数据组中，进行插值。可以采用“两点一线的线性插值法”和“三次样条插值法”等。线性插值公式：二、 实验内容某水溶液在20C下，浓度(重量%)与粘度mPa.s的关系如下：浓度，重量%02030404550粘度，mPa.s1.0051.7762.4803.6524.6215.921三、 实验任务 试求出粘度在1.15.9之间，任意一个粘度值下的浓度值。保留三位小数。四、 参考程序框图x=x+0.5Y.输入给定数据组xi，yi和组数n输入待求粘度初值x i=1 to nENDY.Y. xxn ?or: xxi ? xxi+1 ? 输出x,y i=i+1五、 实验结果记录六、 思考题经过实践，用结果表述线性插值方法和三次样条插值法相比较，哪种方法精确度更高？还可以用其它的那些方法实现本题。实验二 绘制速度分布及剪应力分布的曲线 (1学时)一、 实验目的、原理 学习用程序语言或其它计算机软件绘制曲线。二、 实验内容粘度为50mPa.s的油在内径100mm的直圆管内流动，其速度分布表示为：u = 20 y 200 y 2 剪应力分布表示为：三、 实验任务 以你认为最完整的方式，在同一个坐标上绘制这两条曲线。四、 实验结果记录五、 思考题 速度分布与剪应力分布在管中心和管壁处分别有何表现？实验三 试差法求管径 (2学时)一、 实验目的、原理 练习用化工原理中最常用的试差法求取管径。二、 实验内容 10C的水以500L/min的流量流过一根300m长的水平管,管壁的绝对粗糙度为0.05mm。有6m的压头可供克服流动的摩擦阻力(含进、出口阻力)。三、 实验任务 用试差法求管径的最小尺寸。给出试差中的若干组“摩擦系数， 雷诺数， 管径”值。四、 实验结果记录五、 思考题如何考虑直管阻力hf与压头Hf的关系。可以借适宜的流速范围帮助求解d。六、参考程序框图 图中e为能够接受的相对误差值，建议为1%。输入流量、管长等已知数据设摩擦系数的初值l设由Hf求d算Re及e/d输入查得的摩擦系数或算其值|l计-l设|/l设e？l=l计END输出结果END实验四 求平均速度及方次 (2学时)一、 实验目的、原理 学习将实验测得的不同半径点处的“点速度”数据，处理成平均速度值。 R管半径, mm；Vs流量, m3/s 二、 实验内容 某种流体以稳态在直径为75mm的管中流动，其中各点的点速度ur与该点距离管中心的半径r之间的关系如下表：r，mmur，m/sr，mmur，m/s01.04222.500.9193.751.03326.250.8647.501.01930.000.80911.250.99633.750.69915.000.97833.6250.50718.750.95537.50三、 实验任务 求管内的平均流速及n值。四、 实验结果记录五、 思考题 如何正确地解决积分项的问题？实验五 选择泵型号 (4学时)一、 实验目的、原理 本题是一个工程性很强的训练题，它包括泵的设计计算和操作计算，还有资料查询和经济核算的训练。二、 实验内容 有40C的水按140m3/h的流量从水池送到某设备上方。出水口比池中水面高11.5m，设备内为常压。所有管长及管件的当量长度总共为880m(进、出口阻力略)。三、 实验任务 由适宜的流速范围，以经济核算为准，确定管径(按无缝钢管系列圆整)并选择适宜的泵型号。泵和电机的价格去图书馆或网上查取；管子的价格按2250元/吨管材计算。在结果记录中列出清单后，指明最佳方案。 在所选泵允许的压头下(最佳效率范围内)，求该管路可提供的最大流量。四、 实验结果记录udd圆H泵型号泵/电机价格管路价格总价格五、 思考题1. 管路特性曲线和泵的特性曲线对你解这道题有何帮助?2. 还有哪些因素影响泵使用的价格？实验六 K判据法求颗粒沉降速度ut (2学时)一、 实验目的、原理 学习使用“多条件”处理问题。二、 实验内容 (1)粒径为30mm,颗粒密度rs=2000kg/m3, 求其在r=1.2kg/m3, m=0.0185mPa.s的空气中的的沉降速度。 (2)求d=0.08mm, rs=7500kg/m3的方铅矿在r=1000kg/m3水中的沉降速度。三、 实验任务 K判据以 z = b / Re n为通式，代入沉降速度计算式中，成为： 令： 不同的K值对应着不同的Re范围，相应的b，n值如下：ut区域ReKbn斯托克斯 1 69.10.440四、 实验结果记录五、 思考题 K判据可以免去试差求解的繁琐，但是有什么情况未考虑进去？实验七 沉降法测定液体粘度 (2学时)一、 实验目的、原理 1.了解如何用沉降法测定液体粘度。 2.学习输出表格式的清单。二、 实验内容 一个d=0.5mm的铜球, rs=8.9103kg/m3。在某种油中沉降, 测得沉降速度ut=1.5cm/s, 油的密度r=850kg/m3, 求油的粘度m Pa.s三、 实验任务 连同其它几组数据一起列出清单:drsrutm1.07.81030.92.01.02.61030.820.03.02.61030.810.01.00.91030.830.0四、 实验结果记录五、 思考题 在不同的Re范围，有着不同的沉降速度计算式，式中又含有待求的粘度项。求取的方法是什么？实验八 求沉降槽的最大沉降面积 (2学时)一、 实验目的、原理 练习数值的排序和比较；选择安全系数。二、 实验内容 将固相浓度为236kg/m3的石灰料浆在连续沉降槽增稠至700kg/m3。由间歇沉降试验求得颗粒的表面沉降速度与悬浮液中固相浓度的关系如表C, kg/m3 236 358 425 525 600 714U0, m/h 0.157 0.05 0.0278 0.0127 0.00646 0.00158三、 实验任务 按每天处理的料浆体积为3785m3，确定沉降槽的面积。四、 实验结果记录C, kg/m3236 358 425 525 600 714A, m2 安全系数： 沉降槽的面积：五、 思考题料浆的浓度与沉降槽的面积有无关系？实验九 求导热通量及两层砖间界面温度 (2学时)一、 实验目的、原理 掌握界面温度的求法。二、 实验内容 燃烧炉的内层为460mm厚的耐火砖，外层为230mm厚的绝缘砖。若炉的内表面温度为1400C，外表面温度为100C。设两层砖的接触良好，已知耐火砖的导热系数 ，绝缘砖的导热系数 。式中的t为该层壁的平均温度。 三、 实验任务 试求导热过程的热通量q，以及两层砖之间的界面温度。输入内、外侧温度t0，t2设界面温度t取两壁的平均温度以及两壁的导热系数分别求两壁的导热通量q1及q2|q1-q2|/q1q2?输出t值, ENDt=t-Dtt=t+DtY.Y.四、 参考程序框图五、 实验结果记录 六、 思考题 对于不同材料或不同厚度的平壁层，阻力与动力之间有何比例关系。实验十 画出不同a值下的气液平衡线 (2学时)一、 实验目的、原理 自作CAI演示件。二、 实验内容 以a值分别等于2，4，6，8，10时的气液平衡线。三、 实验任务 要给出坐标值，绘出每一条平衡线时要同时显示该线的a值，每一条平衡线和相应的a值同色，不同值不同色。四、 实验结果记录五、 思考题 为什么说，某物系的a值越大，分离越容易？ 实验十一 用安托尼公式求苯-甲苯数据 (2学时)一、 实验目的、原理 掌握t与x，y间的运算关系，为求灵敏板数打下基础。二、 实验内容 安托尼公式 ，p饱和蒸汽压，kPa； t物系温度，C；A、B、C安托尼常数，分别为：组分ABC苯6.0231206.35220.24甲苯6.0781343.94219.58a三、实验任务 求常压下，(1)对应x=0和x=1的温度；(2)同(1)范围的 ；(3)任意一个x值下的温度。四、 实验结果记录五、 思考题由拉乌尔定律考虑x=0和x=1时，组分蒸汽与操作压强的关系。实验十二 比较不同q及不同R下的理论板数 (2学时)一、 实验目的、原理 (1)比较不同进料热状况及不同回流比下理论板数的变化情况。 (2)掌握Rmin、xd的求法.二、 实验内容 常压连续精馏, a=2.46, xD=0.9, xF=0.4, xW=0.02(摩尔分率)。三、 实验任务 (1)填写结果记录所给各项内容。 (2)求理论板数，精确到保留小数3位。四、 实验结果记录q1.21.00.60-0.2R=3xdNTR=2RminRxdNT五、 思考题在不同的R的计算基准下，xd 及NT 分别有何变化规律?实验十三 求理想物系精馏塔的理论板数 (2学时)一、 实验目的、原理 练习从进料板处向塔的两端用逐板计算法作理论板。为求灵敏板奠定基础。二、 实验内容 常压精馏苯-甲苯物系，进料热状况分别为：0.5、1.0、1.2，进料含苯0.4(摩尔分率, 下同)塔顶馏出物含苯0.9，釜液含苯0.1，R=1.5Rmin。三、 实验任务求该塔的理论板数；加料板位置及各板的温度、组成分布四、 实验结果记录q0.51.01.2RminNT（ 釜）NF序号备注xitxitxit五、 思考题 进料热状况对理论板数有何影响，为什么？实验十四 求非理想物系精馏塔的理论板数 (2学时)一、 实验目的、原理 练习用下列四种计算方法之一，处理非理想物系气、液平衡。二、 实验内容 常压精馏乙醇-水物系，xD=0.82, xF=0.2, xW=0.004(摩尔分率)。进料q=1，R=2.0。三、 实验任务 求该塔的理论板数；加料板位置及各板的温度、组成分布。 气、液平衡关系处理的计算方法：(1)三次样条插值法； (2)最小二乘法回归曲线； 乙醇-水物系t-x-y数据txy95.585.582.781.579.878.7478.4178.150.0190.12380.23370.32730.50790.67630.74720.89400.17000.47040.54450.58260.65640.73850.78150.8940 (3)直接利用xy数据； (4)两点线性插值法四、 实验结果记录NT序号备注xit五、 思考题处理气、液平衡关系的计算方法还有哪些?你所选用的计算方法有何优点?实验十五 操作型精馏塔的计算 (2学时)一、 实验目的、原理 掌握“操作型”计算的方法。二、 实验内容 某精馏塔有10块理论板(不含釜)，用于分离苯-甲苯物系。xF=0.4, xD0.95, xW 0.06(摩尔分率)，进料q=1, (设a=2.46)。三、 实验任务 求：R至少为多大？NF=？理论板数保留两位小数时，xW=？四、 实验结果记录五、 思考题 当理论板数相对误差的允许值由0.1增至0.01时, R值的变化有多大? 建议用“二分法”求。即：取R1=Rmin，R2=4 Rmin；若以R=R2求出一个NT满足 题目给定的理论板数NT，则初始R2确定；然后以R=(R1+R2)/2求出NT，与任务给定值相比较，若NT，则再在R与R1之间确定新的中间值R；否则在R与R2之间确定新的R，直至所求NT满足要求为止。实验十六 求吸收中的传质单元数 (2学时)一、实验目的、原理练习用积分图解法或近似图解法计算。二、实验内容 填料吸收塔，用水吸收SO2气体，原料气中SO2含量为5%(体积%)，要求吸收率达95%。塔径0.8m，操作温度303K，压力100kPa；入塔气流量1000m3/h(操作状态)，用水量为最小值的1.2倍。三、 实验任务 求传质单元数。 SO2与水的平衡关系：a,kg SO2/100kgH2O0.10.20.30.50.71.01.5P,mmHg4.711.819.736.052.079.0125四、 实验结果记录五、 思考题 在什么情况下, NT=NOG?实验十七 气体流量对塔径的影响 (2学时)一、 实验目的、原理 了解气体的体积流量对气、液分离塔的塔径尺寸的影响。二、 实验内容 连续精馏塔分离乙醇-水物系。已知xD=0.85, xF=0.4, xW=0.002(摩尔分率)，F=100kmol