1、液面高度H0/cm12为了简化模型,将纸杯近似为圆柱形,底面半径取r1和r2的平均值r=7cm。忽略温度变化对于静置的碳酸饮料的影响,认为该体系处于293K(约20C)的恒温环境中,并且碳酸饮料为均匀溶液,CO2浓度在任意时刻不随空间点的变化而变化。同时,为了便于查找数据,不妨设碳酸饮料的液相主体物化性质与水相同,此时二氧化碳的流失主要以如下两种传质扩散的形式进行:(i)液面扩散传质:碳酸饮料在液面出的传质过程相当于二氧化碳在水和空气中气液两相间的扩散传质过程,亦即CO2水溶液在空气环境中的解吸过程,该解析过程属于CO2的由液相主体向气相主体的单向扩散。,以A代表CO2,B代表空气,C代表水,
2、设pA为气相主体中CO2分压,cA为液相主体中CO2浓度,pAi、cAi分别为相界面处CO2的气相一侧分压与液相一侧浓度,(NA)G、(NA)L分别为溶质通过气、液膜的传质通量,根据双模理论模型可知:已知20C时,水的密度=998.2kg/m3,E=1.44105kPa,Ms=18g/mol,则溶解度系数为:同时有:通过摩尔平衡可得,因此:其中,pA0恒定为空气中CO2分压(CO2体积分数0.03%),即pA0=0.03%101325Pa=30.3975Pa,但是液相主体内CO2浓度随着CO2的扩散会逐渐减少,因此这是一个非稳态的传质过程。同时,注意到;且在CO2在气液两相中含量均较少的情况下
3、,漂流因子,则可以认为气液相传质系数在此过程中均恒定。设(NA)1为液面扩散传质部分的传质通量(mol/(m2s),联立以上各式可得:若设液面扩散传质引起的液相CO2物质的量变化为,浓度变化为,同时液面面积为S1恒定,液相体积为V0,则:(ii)杯壁扩散传质:由于初始碳酸浓度过高,碳酸饮料会在纸杯壁附近形成诸多气泡,这种方式形成的气泡仍可按照单向扩散模型处理,膨胀阶段时气泡所受到的约束力大于上升力,但随着气泡的增长气泡所受上升力逐渐增大,直至上升力等于约束力时气泡脱离壁面,下一个新的气泡又开始产生。同时,随着传质推动力的减小,二氧化碳进入气泡的流率逐渐减小,认为在二氧化碳扩散过程中壁面各处气泡
4、的排布密度相同(近似均匀气膜,厚度随时间周期性变化),只是单个气泡生成所用的时间逐渐变长。设(NA)2为某时刻CO2相气泡内壁面扩散的传质通量(mol/(m2s),S为某时刻全部气泡的表面积,设此部分扩散传质引起的液相浓度变化为,同时液相体积为V0,则:上式中,为了简化模型可以认为S为一常量,可以认为在CO2扩散过程中在杯壁表面形成了一个表面积稳定但厚度不断波动的气膜对气泡生成及其消散过程进行模拟,因此气泡总表面积S为杯底面积S1与杯侧面积S2之和,分别对应于杯底气泡总表面积与。但是(NA)2为时间和位置的函数,随着时间的延续和深度的改变传质通量也由于传质推动力的变化而有所不同,下面分别就杯底
5、和杯侧壁的气泡生成情况对其进行分析。对于杯底气泡,认为其在生成过程中其内气体的压力等于其所在位置所受到的液体压强,设液相水密度为,纸杯高度为H,底面半径为r,气泡内CO2压强为p,则根据单向扩散传质模型有:其中,设此部分扩散传质引起的液相浓度变化为,则:对于杯侧壁气泡,认为其在生成过程中其内气体的压力等于其所在位置所受到的液体压强,设液相水密度为,纸杯侧壁某液膜在液面以下深度为h(0hg,一般情况下,k=r0/re10时,可以认为气泡形状为球型(忽略球冠部分)水的表面张力/(N/m)0.07197设k=r0/re,利用MATLAB求解方程作出rek的关系,源程序如下:%求解气泡最大半径源程序:
6、k=0:0.001:0.1;dt=0.07197;r=1000*sqrt(3*k*dt/2/g/p);plot(k,r) r0/re re /mm图三 rek的关系曲线变化图通过实际观察可知图三中点(0.091, 1.00161)与实际情况(每个脱离底面的气泡半径大约为1mm)较为相符,因此气膜的最大厚度为Dmax=2.00322mm。利用碳酸饮料中二氧化碳流失过程模拟的结果,对杯底气膜的厚度周期性变化模拟,设(NA)1为液面扩散传质部分的传质通量(mol/(m2假设杯底气膜中气体为理想气体,设杯底面积为S0,则杯底气膜厚度为:当DDmax时,则气膜厚度清零重新增长,自此周期性往复变化。3.3
7、程序编制上述模型的求解要利用碳酸饮料中二氧化碳流失过程模拟的结果,MATLAB原程序如下:%插值求解初始时刻二氧化碳浓度随时间变化关系:25;t0=t(1:26);cA00=cA(1:cA=lagr(t0,cA00,ts);%插值求解5h附近二氧化碳浓度随时间变化关系:ts=17975:18000;t0=t(17976:18001);cA00=cA(17976:%求解初始时刻/5h附近杯底气膜厚度随时间变化关系源程序:NA=kG*(pA0+a*cA)/(1+a*H)-pA0);Dmax=0.00200322;j=1;for i=1:25001t=ts(j:i);Na=NA(j:DD=simp(
8、t,Na,0.001)*R*T/p/H0/g;if DD=Dmax D(i)=DD;else D(i)=0; j=i;end%作图:plot(ts,D)时间t/s杯底气膜厚度D/m3.4结果讨论程序输出结果:图四 初始时刻杯底气膜厚度D随时间变化关系图五 5h附近杯底气膜厚度D随时间变化关系观察图四、图五,可知初始时刻每层杯底气膜形成所需时间约为1s,而5h附近每层杯底气膜厚度形成所需时间约为10s,可以看出随着时间的推移,形成单层气膜所需时间逐渐变长,这与实际情况是相符的。4.总结由于期中之前主要学习了传动量部分,而第一次大作业主要针对传热部分,因此本次大作业主要针对传质部分。在解释现象时由于第一次大作业已经运用数值解法求解非稳态过程,因此本次主要应用了双膜理论对非稳态模型进行了大量简化后进行模拟,但是双膜理论中膜厚的具体数值难以确定,只能凭借文献值和自己的经验给出气、液相传质系数等未知参数较为准确的量级,通过后续工作地分析去验证校核这些参数最终得到较为符合生活实际的模型。在这个过程中,很好地锻炼了自己的工程经验。在模拟生活现象的同时也充分锻炼了自己的编程能力。通过此次作业发现随处可见的生活现象几乎都与传动、传热、传质相关,因此均可以通过三传原理进行解释或模拟,深刻体会到了传递过程原理在生活中应用的普遍性及其重要意义。
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1