化工数学答案(全)教学文案.DOC

1、此文档仅供收集于网络，如有侵权请联系网站删除化工数学各章习题选解（仅供参考）第一章习题1. （） 在一个有效容积为 V 的半连续式搅拌反应器中，由原料生产物质，若浓度为 c0 流量为 Q 的溶液加入空反应器，反应遵循以下连串 -可逆步骤k1Bk3A Ck2且所有的反应均为一级，证明在反应器中的克分子数 NB 是以下微分方程的解d2N2dtBPdNdtBRNBC式中Pk k1k2 3Rk k1 3C Qc0k1证明 ：对 A 、B 分别作质量衡算，有dNAA ： c0Q k1N k2 N (1)A BdtdNBB： k1N k2 N k3 N (2)A B Bdt 由（2）得到：dNAk1N c

2、0Q k2N (3)A Bdt （3）代入（ 2），得：2dN d NB Bk1c0Q k1k3N (k1 k2 k3 ) 2 (4)Bdt dt令P k1 k2 k3, R k1k3, C c0Q得2d N dNB B2 (5)P RN CBdt dt证毕。2. 冬天的池塘水面上结了一层厚度为 l 的冰层，冰层上方与温度为 Tw 的空气接触，下方与温度为 0的池水接触。当 Tw0时，水的热量将通过冰层向空气中散发，散发的热量转化为冰层增加的厚度。已知水结冰的相变潜热为 Lf，冰的密度为，导热系数为 k，导温系数为，求：只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除1) 当气温 T

3、w 不随时间变化时，给出冰层厚度随时间变化的关系，若 Lf3.35105J/kg，913kg/m3，k2.22W/m K，Tw 10，问冰冻三尺，需几日之寒？2）当气温随时间变化时，设 TwTw (t) 已知，导出冰层厚度变化的完整数学模型。解：(1) 冰层的温度为 0，水通过冰层向空气散发热量，记为 Q，该热量用于水结成冰。假设冰层面积为 s，厚度为 l 根据导热方程，可得：sQ k(Tw 0) dt Lf lsdl代入数值， Lf3.35105J/kg， 913kg/m3，k2.22W/m K，Tw10，l1m，求解积分上式得：22.2t 10 5 dt 1d t 913 3 .35 10

4、0t79.7 天80 天若冰冻三尺，在 Tw10时，需要约 80 天。(2) 若 Tw Tw (t)，冰层厚度为 l根据热量守恒：sQ k(Tw 0) dt Lf lsdlL ldlfkT dtw两边积分：l0L ldlft0kTwdt0.5 Lt2f l k Twdt0厚度变化与 Tw 的关系为：l2kLft0Twdt3. （） 在一个半分批式搅拌反应器中进行着一级放热化学反应，反应速率常数由Arrhenius 关系式给出， 反应热由釜内的冷却盘管移出， 请自行设定有关的参数， 导出该反应器的数学模型。解：设物料以恒定的体积流量 F 加入，则反应器中反应物浓度 CA 与温度 T 由以下物料衡

5、算与热量衡算方程给出1 物料衡算方程只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除FC Vr0 Ad (VC )Adtr kCA AV V Ft0(1)k k0Eexp( )RT2 能量衡算方程d(C TV )pF C T K (T T ) A kC V ( H ) (2)p 0 c A rdt1 2 合并，得数学模型为d (VC )AdtkVC FCA 0F C T K (T T )A kC V( H )p 0 c A rd(C TV)pdtV V Ft0(3)k k0Eexp( )RTc (0) 0, T (0) TA 0式中 K( TTc)A 为冷却移热， kCAV(-Hr)

6、为反应热。4（）采用微元分析法推导出柱坐标系中的不定常热传导方程。解： 考虑柱坐标系中热传导方程的形式。柱坐标系下的三个空间变量：向径 r，经度角 ，高度 z。在这三个方向上，与自变量的微分变化所对应的线段微元长度分别是dr (dr ,rd , dz) （4）由偏导数的定义，温度梯度 T 在三个方向的分量即温度在每个方向上的微元增量除以相应的线元长度，即TT T T( , , )r r Z（5）于是 Fourier 热传导定律在柱坐标系中的分量形式为T T T（q ,q ,q） -k( , , ) （6）r zr r Z 接着考虑各方向输入和输出的微元通量，首先考虑 r 方向T输入项： （7）

7、k (rd )( dz)rT T输出项： （8）k (rd )( dz) k (r )drd dzr r r只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除 T于是 r 方向的净输入通量为： k (r )drd dz r r对 方向作同样的分析，（9）T输入项： （10）k (dr )( dz)rT T输出项： （11）k (dr )( dz) k ( )drd dzr r T净输入通量： （12）k ( ) drd dzrT输入项： （13）k (rd )( dr )zz 方向的分析，T T输出项： （14）k (rd )( dr ) k (r ) drd dzz z zTk (r

8、 ) drd dzz z净输入通量： （15）T微元体内的积累项： C dr( rd )dzpt（16）将三个方向输入微元的热流净增量加和并令其等于积累项，就得到2 2T 1 T 1 T T (r ) 2 2 2t r r r r z（17）5. 风吹过皮肤表面时，人会有干燥凉爽的感觉，这是因为风的吹拂强化皮肤表面的对流传热与传质，形成一个速度，温度，浓度（含水量）的边界层，设流动为层流（微风） ，考虑出汗的蒸发潜热，求：1）列出皮肤表面的三传问题的边界层方程，根据实际情况适当简化并给出问题的边界条件；2）将上述问题无量纲化，并解释所得到的各无量纲参数的物理意义；3）试分析速度分布，温度分布，

9、含水量分布分别与哪些无量纲参数有关，并用简单的函数关系示意；4）根据所得结果定性的解释一些经验常识：为什么风越大越感觉到冷？为什么出汗后擦了汗感觉更凉快？当空气中湿度变化时， 对表面散热会带来哪些影响？在冬天和夏天， 人体对空气湿度的增加会有什么样的感觉？解：1）同时考虑流动传热传质时的边界层传递方程是2u u p u) -(u g (T T ) g (C C )2 x xx y x y2T T T u2C (u k ( ) VHr) p 2 ix y y y2c c ci i i ) (u D i 2x y yg 表示重力在 x 方向的分量， 为热膨胀系数， 为密度变化系数x只供学习与交流此

10、文档仅供收集于网络，如有侵权请联系网站删除VH 水汽化潜热 ri 水蒸发速度由于px可忽略， 0g ，xu2( )y可忽略，化简后(u2u u u) 2x y y2T T TC (u k VHr) p 2 ix y y2c c ci i) i (u Di 2x y y边界条件y0， u0（皮肤表面气流速度） TT0（皮肤表面温度） cc0（皮肤表面的含水量）y1 uu（速度边界层外气流速度）y2 TT（温度边界层外气流温度）y3 cc（浓度边界层外气流中含水量浓度）1，2， 3分别为速度边界层，温度边界层，浓度边界层的厚度。2）无量纲化uxuT T0TT T0c c0cc c0无量纲物理性质的

11、比值1 vPr TcDiSc无量纲化后2u v v vv v 2x y v yv2Vk Hru T T T iT T 2x y C T y Cp T p2Dc c i cu c c 2x y c y只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除边界条件在yy0, 0, 1对于较大的 Pr 或 S c，热传导与扩散效应与黏性比较相对较弱，热边界层和扩散边界层位于速度边界层内部，反之，对于较小的 Pr 或 Sc，热传导与扩散速率大于黏性传递速率，热和扩散边界层就有可能扩展到速度边界层之外。3）速度分布，温度分布，含水量分布的简单函数关系式uxu( 1)uxu1( 1)T T1 ( 1)

12、 ( 1)T TT T0 0TT01( 1)Tc c1 ( 1) ( 1)c cc c0 0cc01( 1)c4）风越大，皮肤表面的气体更新速度越快，水的蒸发速度变快，传热越快，感觉到冷出汗后感觉更凉快，是因为减小了汗水层的厚度，蒸发速度加快当空气中湿度变大时，皮肤表面水的蒸发速度变慢，不利于传热夏天空气湿度增加，汗水蒸发困难，人感觉闷热冬天空气湿度增加，少量的汗水在皮肤表面使人感觉温暖。6（）在管式反应器模型（ 1.4.15）中，当 Pe 0 时，相当于完全返混的情况。试从方程（4.15）出发，通过适当的体积积分和取极限 Pe 0，导出均相釜式反应器模型。解：当 Pe 0 时，由原方程（ 4

13、.15）及边界条件可知， c=const，说明在完全返混的情况下，反应器内具有均匀的浓度。对于任意的 Peclet 数，对方程 4.15 进行体积积分得到2c 1 c c1 1 1 1dz Dacdz dz dz20 0 0 0Pe z z 1c 1 c Dac cPe z010（31）式中 c 为反应器内的平均浓度。将边界条件（ 4.15）代入（ 31），得到只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除dcdDac cz1 1（32）= ， (0)=0 c 0上式对任意 Peclet 数均成立，仅当 Pe 0 时，反应器内浓度均匀，cc ，上式成为无量z 1纲的理想混合釜式反应

14、器数学模型。7. （） 烯烃在 ZieglarNatta 催化剂颗粒上的气相聚合过程可用最简单的固体核模型来描述，如附图所示。气相中的烯烃单体在催化剂颗粒（图中阴影部分）表面聚合后生成一多孔的固体聚合物壳层并将催化剂包裹在内部， 外部的气相烯烃单体只有扩散穿过此固体聚合物壳层后才能到达催化剂表面参与反应。试求：（i）证明单体在壳层中的扩散及聚合物粒子的生长由以下方程描述M 1 M2D ( r )2t r r rdRdtM DWsMr|rR式中为单体浓度 (mol/m 3)，s 为聚合物壳层的密度 (kg/m3),为单体在壳层中的扩散系数(m2/s), M为单体的分子量， R 为聚合物颗粒的半径

15、。（ii）设催化剂核半径为 rc，单体在外部气相本体中的浓度为 ，以上述参量为 r 和 M 的特征尺度，并引入适当的时间尺度，将上述方程无量纲化。然后根据气相单体与固体聚合物密度之间的巨大差别 (s/g103)将问题进一步简化。（iii ）设单体在催化剂核表面的浓度恒为 0（瞬时反应） ，R 的初始值为 R0(R0 rc)，求解上述简化后的模型并给出聚合物粒子半径 R 随时间的变化关系。提示：对单体的浓度分布可采用拟稳态假定。解：为简化计算，令单体分子量 M w 的单位是 kg（1） 问题建模如图 1 所示，对微元 dr 作物料衡算( MA )dr ( AJ )dr 0t rMJ D , A

16、4 rr2(18)得M 1 M2D (r ) (19)2t r r r如图 2.对微元 dR 作物料衡算42R dR M2s4 R M D (20)wdt rr RdR M D M得 w (21)dt rs r R只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除R RR+dRr+ drr图 1 图 2（2） 无量纲化与简化分析：本问题存在着两个特征时间尺度，一个是单体组分内扩散通过聚合物壳层的时间尺度12rcD，该尺度可以从内扩散方程 （19）中得出， 为此， 近似取M1代替Mt，M/r 代替Mr，M/r2 代替2M2r，就可估算出12rcD；另一个是聚合物颗粒生长的特征时间尺度 2

17、，可以从方s程（20）中用类似的比值代替微分的办法估算出 2 1 ? 1 。在对问题进行无量纲化时，g不同时间尺度的选择代表着所关注的不同过程。1 如果选取12rcD为时间尺度，式（ 20）和（ 21）可分别无量纲化为（仍然用当前变量表示无量纲变量） ：M 1 M2(r ) (22) 2t r r rdR Mdt rrRrc(23)gM M , = 1 g B ws此时式（ 23）中出现一个小参数 。时间尺度 1 称为快时间尺度，选择这一尺度所得到的方程（22）中不含有小参数 ，表示我们关注的是单体 M 通过聚合物壳层的不定常扩散而不是粒子的生长。略去（ 23）中的小参数项后得到 R cons

18、t，说明在考虑单体内扩散时，由于时间较短，可以将粒子半径作为常数考虑。因此，选择 1 为时间尺度显然不妥，得到的不是我们希望关注的问题。2 如果选取22rc sDg为时间尺度， （20）、（21）式可无量纲化为：只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除M 1 M2(r ) (25) 2t r r rdR Mdt rrRrc(26)此时粒子生长方程（ 26）不含小参数，粒径将随时间变化，表示我们关注的是颗粒的生长。而单体内扩散方程（ 25）中的时间导数项含小参数 ，可以略去，说明在慢时间尺度 2 上考虑粒子生长时，单体的内扩散过程可以忽略时间变化项，内扩散可以作为拟稳态过程来考

19、虑。从1 和2 中得到的不同简化模型说明时间尺度的选择需要根据建模目的来考虑， 使简化后的模型能够代表所关注的过程的主要特征。（3）对单体的浓度分布作拟稳态（时间导数项为零）假设，即（ 25）中的 M 0t，得1 M2(r ) 0 2r r rr 1,M 0 (27)Rr , M 1rc解得M1 R(1 ) (28)r R rc代入（ 26），得dR2rcdt R(R r )cR(0)R0rc(29)解得3 2 3 2R R r R Rc 2 0 0r t3 2 30 c3 2 3r 2rc c8. 在缺乏数学模型的某些情况下， 仅仅根据量纲分析或尺度比较也可以获得一些很有价值的结果，考虑以下

20、例子：1） 对于固体颗粒在黏性流体中的 Stock 流动问题，颗粒受到的阻力 f 仅仅与颗粒尺度 d，动力学粘度和速度有关，即ff（d，）根据量纲齐次化的要求，物理方程等式两边的量纲应该相同，而有参数 d，组成的具有离地量纲的参量只可能是 d，因此上述函数关系只可能取一下形式，fAd式中 A 是一个只与颗粒形状有关的常数，上式即为 Stock 定律。现根据上述量纲分析方法分析湍流的消磁度运动。 湍流中存在一系列大小不同的涡旋，能量从大尺度涡旋顺序传递给消磁度涡旋，同时将机械能耗散为热能， 其中最小的涡旋尺度称为Kolmogorov 尺度，在这个尺度上，黏性和能量耗散占优，因此只有运动学粘度 v

21、 (m2/s)和能量只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除耗散速率 (W/Kg ）两个产量起作用，其他物理量都可以用这两个量表示。试根据量纲齐次化原理推导出 Kolmogorov 尺度及局部速度 与 v，的关系（可相差一个常数）2）流体在自由空间中的射流形成一个夹角为的圆锥型区域，如图所示，设 UU（x）为距喷口 x 处的平均流速， RR（x）为 x 处的射流半径，试根据总动量2 2U R 沿 x 方向守恒的要求确定速度 U 和射流区总流量2Q UR 沿 x 的变化关系 （可相差一个常习题 8：湍动射流数）3）对于放热反应，当反应器尺寸增大时，其体积按长度的三次方增长，而表

22、面积却按平方增长，因此体积增大有利于热量的增加，而体积减小有利于冷却散热。这是化学工程中说明“放大效应” 的一个典型例子。根据类似的道理解释为什么生活在寒冷地区的动物一般体型较大（例如北方人就比南方人高大） ，而且形状趋于圆滑，而热带地区的动物体型较小且趋于瘦长（例如南方人比北方人相对较瘦，且身体凸出部分的轮廓更为明显） 。解：1）分析：题中出现的符号意义如下Kolmogorov 尺度 m能量耗散速率 W/Kg运动学粘度 v m2/s局部速度 m/s从 fAd可知 f 量纲与 d相同，f又的量纲 W/Kg 与的量纲相同V（f 力N，速度 m/s，密度 Kg/m 3，V 体积 m3）可知 V 与

23、3 量纲相同， f 与v 量纲相同，与量纲相同，f带入 中Vv=322v2）分析：总动量2 2U R 沿 x 方向守恒，设2 2U R C（C 为常数）。又 R tg x ，将 R 带入总动量表达式中，只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除U tg x C2 ( ) 2UC(tg x)C tg x 2 2Q UR gtg x C2 ( )1(tg x)3）分析：将人体看成一个反应器， 食物在体内消化放出热量，除供人体正常活动所需的能量和储存在体内外，以热量的形式通过体表释放到体外。人体散热与人体表面积和外界温度和人体温度的差值成正比， 表面积与人体尺寸的平方成正比。可记为散

24、热 A* 温度差 *表面积。人体放热与人体体积成正比，体积与人体尺寸的立方成正比，可记为放热 B* 体积（A，B为关于人体散热放热的常数） 。为保持人体温度一定，放热与散热需要平衡，散热放热。寒冷地区的外界温度和人体温度的差值比热带地区大，所以需要的体积与表面积比也大，即表现为人体尺寸大，也就是人的体型较大，而形状圆滑是为了减小表面积所致。热带地区的人体型瘦长，身体凸出部分轮廓明显，可以增大表面积，方便散热。9. 在水平液面上垂直插入一个半径为 R 的毛细管，此时液体将在表面张力的拉动下沿着管中上升。弯曲液面形成的毛细压强可以用以下 YongLaplace 方程计算P2CosR式中为气液表面张

25、力， 为气液界面与固壁之间的接触角，管中流体一方面受到毛细压强的驱动而上升，一方面又受到重力和粘性阻力的作用，设流动速度遵从粘性管流的 Poiseuille 分布，求：1) 对于两端开口的毛细管，证明液位高度 H 随时间 t 的变化满足以下方程dHdt21 2R gRCos8 H式中为液体的动力学粘度， g 为重力。2）对于上端封闭的毛细管，设总管长为 l，管内气体满足理想气体状态方程，试推导相应的液位高度 H 的变化方程。3）从上述方程中求出最大液位高度 Ho 和时间变化关系 H (t),据此讨论 H 变化的趋势。解：(1) 弯曲液面形成的毛细压强可以用以下 YongLaplace 方程计算

26、，同时又受到重力的作用产生压强，总的 P 为：P2CosRgH式中为气液表面张力，为气液界面与固壁之间的接触角， R 为毛细管半径。只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除Poiseuille 分布，体积流率的表达式为：dV2RdH4RPdt dt 8 H式中为液体的动力学粘度， g 为重力。dHdt2 2 2R P R 2 Cos 1 RgH8H 8H R 8 HCos2gR(2) 若毛细管上端是封闭的，则 P 由三部分组成，还有一部分是液面上端产生的压强。P2CosRgHHPoL H式中 Po 为大气压强， L 为毛细管长。2 2 2 2dH R P R 2 Cos Po

27、 H 1 R gR PogH Cosdt 8H 8H R L H 8 H L2RHdH(3) 当达到最高液位 Ho 时， 0dt，则：dHdt2 2 2R P R 2 Cos 1 RgH8H 8H R 8 HCos2gR0H o2CosgRdHdt82RHCos2gR2RCos 8H2gRH8 HH 20 2RCos gRo toH0dt积分得：8 2 2t0 gR H o 2 RCos 1 ln 2 RCos gR H o 2 RCos ln 2 RCos2 2 4g R1从上式可以看出， H 先增高，到最大值2CosgR后开始下降。10. 气液两相的传质过程与色谱过程有许多类似之处，例如，

28、气相通过反应器（鼓泡塔、板式只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除G塔、填料塔等）的流动可以看成是溶质通过固定相的运动，气液传质阻力可类比于气固传质阻力， 气液两相的逆流操作模式也与移动床相似。此外，气液两相在界面上处于平衡状L液位态，由 Henry 定律表述，与 7.2 节考虑的微孔分子筛的内扩散过程类似。与色谱问题不同的是，许多气液反应器（鼓泡塔与搅拌釜）中的液相或液固两相一般都处于全混流状l 态，而色谱柱中固定相是静止的， 移动床中固体接近平推流。试根据与移动床的类比建立如图所示的鼓泡塔反应器的稳态数学模型，图中气体从塔底加入，经分布器之后形成分散气体分布器的气泡并在

29、液体中浮升，最后从容器的上部输出；液体则从G塔顶加入， 从底部流出。 气相中的组分 A 被液体吸收后在液相中发生一级化学反应。 鼓泡塔中气相的流动可考虑为平推L流，液相考虑为全混流。其它已知的参数为：鼓泡塔液位高习题 10：鼓泡塔反度 l，气含率 g，空塔气速 U，加入液体的质量流率 F，单应器位体积气液传质系数 kLa，一级反应动力学常数 kA，Henry系数 HA。所建立的数学模型要求包括以下内容：1）设 cg 和 cL 分别为反应组分在气相和液相中的浓度，给出其方程和边界条件；2）如果是强放热反应，反应热通过溶剂蒸发和气液相的连续流动移出，请自行设定有关物性参数，给出温度 T 满足的方程

30、。解：数学模型一般包括物料衡算、热量衡算和动量衡算，对于鼓泡塔的气液反应体系，气液两相的温度场计算是不必要，因为塔内混合良好，温差很小。动量衡算也非必要。物料衡算是要考虑到对流、相间传质、轴向分散及化学反应等影响因素。(1) 假设塔内等温，作气相及液相的物料衡算：F F (F /U ) (Fig /U g )F 1ig ig ig gil气相： k l(H ) E l l g Zg gZ U U l Z Z tg l流入项 传质项 返混项 积蓄项液项：FF F 1 (F /U ) (Fil /U ligil il il lk l(H ) (1 )E kc (1 )ll g Zl l gZ U

31、U l Z Z tg l)流入项 传质项 返混项 反应项 积蓄项式中：Fig c U Ag g RUg UgFU 1c g AR只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除假设 气相平推流 E 0 液相全混流 EZl 1Zg边界条件：气相： Z0：0F Fig lgZ1： dFig / dZ 0液相： Z0： dF / dZ 0ilZ1：F ( 1) F ( 1)il ilD 1 d F ( 1)t g ill dZ(2) 进气温度为 Tg，液体为 Tl，气体热容 cg，液体热容 cp，液体的质量流量 Fl，气体的质量流量 Fg，单位时间反应热为 Q，C 为常数。鼓泡塔内温度均匀

32、，塔内温度和排出的气体液体温度均记为 T，根据能量守恒：Fl cp T Tl Fg cg T Tg Q CdTdt11填料塔广泛用于气体吸收，气液两相采用逆流操作，液体从塔顶均布后加入，沿填料表面成液膜下降，气体从塔底加入，沿塔上升并与液体实现逆流接触，气体中的活性组分被液体吸收后从塔底流出，净化后的气体从塔顶排出，如图所示。设从塔底加入的气体中含有待吸收组分 A 和惰性气体，惰气流量为 G（mol/s），从塔顶加入的液体惰性溶剂的流量为 L（mol/s），组分 A 在液相中以一级反应进行分解，给定塔的直径 D 和塔高 H、单位体积填料的液体持液量 L（m3/m3）和气液传质系数 kLa，以及

33、化学反应速率常数 kA、气液相 Henry 系数 HA，试用微元分析法建立一数学模型，描述气相浓度 yA (mol/mol 惰气)和液相浓度 xA (mol/mol 溶剂 )的沿塔分布，然后从模型中消去 xA，得到 yA 的单一方程，并给出适当的边界条件。提示：可假设在气液界面上满足 Henry 定律，则两相传质速率为 kLa (yAHAxA)。解：习题 11：填料吸收塔分析此填料塔用于气体的化学吸收，塔内物料平衡涉及两相：气相和液相。可分别对气相及液相中的待吸收组分作质量守恒， 守恒方程中将涉及到的未知量包括： 待吸收组分 A 在气相中的浓度 yA（molA/mol 惰性气体）和待吸收组分

34、A 在液相中的浓度 xA(molA/mol 溶剂 )。守恒方程数与未知量数均为 2，在给定边界条件下可以得到微分方程的特解。假设1 设待吸收组分 A 在液相中的反应速率为：只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除dcAr kxA Adt该反应为一级反应，反应速率常数为 k（单位 mol/m 3s）2 设吸收塔在连续操作过程中处于稳态，进而在质量衡算方程中涉及到的积累项均为 0解答步骤在填料塔任一高度 h 处取一厚度为 dh 的体积微元（见习题 11 图），分别考虑微元中气相和液相的待吸收组分 A 的质量守恒：气相输入项= y +dyA AG输出项 =y G+气液相间的质量传递

35、量a式中 气液相间的传质方向是由气相到液相，因而传递量属于输出项，其传质速率为：k a y H x ，其中 kL 是以气相摩尔分率差 yA H A xA 为总传质推动力的总传质系数L A A A（单位 kmol/(m 2sy)），a 是单位体积填料层所提供的有效传质面积（单位 m2 传质面积 /m3填料体积）。因而：气液相间的质量传递量2Dk a y H x dhL A A A42D式中 D 是填料塔直径， dh 为所取微元内的填料体积。42D输出项 =y G+k a y H x dh a L A A A4由于气相中不发生待吸收组分 A 的分解反应，因而：生成项 0吸收塔处于稳态操作，因而：积

36、累项 =0根据质量守恒：输入项 输出项 生成累 积累项即：2Dy dy G y G k a y H x dhA A A L A A A40 0整理得：2DGdy k a y H x dhA L A A A40只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除dy DAG k a y H xL A A Adh 420液相：输入项 气液相间传质量x LA式中气液相间的传质方向是从气相到液相， 因而传弟量属于输入项， 它的大小与气相输出项中的气液传质量相同，为：气液相间的质量传递量2Dk a y H x dhL A A A4因而：输入项2Dx L k a y H x dhA L A A A4

37、输出项x dx LA A生成项2Dr dhA L4式中-rA 为组分 A 在液相中的应反应速率，单位： mol/m3s。而2D4dh 是在所取填L料微元中的持液量。dxAr kxA Adt生成项2Dkx dhA L4k DL2生成项 dh4积累项 0根据质量守恒：输入项 输出项 生成累 积累项即：2 2D k DLx L k a y H x dh x dx L dhA L A A A A A4 40整理得：L2 2dx D k DA Lk a y H xL A A Adh 4 40只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除相应边界条件：在填料塔底部 h=0 yA=yb 式中 y

38、b 为进料气体中组分 A 的浓度顶部 h=H xA=xa 式中 xa为进料液体中组分 A 的浓度至此，联立关于气相与液相质量守恒的两个常微分方程，并加入边界条件：2dy DAG k a y H xL A A Adh 40 1 11 1L2dx D k DA Lk a y H xL A A Adh 4 420 1 11 2h 0 y yA ah Hx xA a求解此一阶常微分方程组即可得到气液相中组分 A 的浓度沿着塔高方向上的分布。下面对此常微分方程组进一步化简，消去 xA，以得到一个只含 yA 二阶常微分方程。对式（ 1-11-1）变形，得：xAyA4G dyA2k a D dhLHA代入式

39、（ 1-11-2）得到：4G dyAy2 A 2 2 2L dy 4G d y k a D dh D k DA A L Lk a y H2 2 L A AH dh k a D dh H 4 4A L A02 24GL d y L dy k DA A LG2 2k a D H dh H dh 4L A A0相应边界条件：h 0 y yA a4G dyAyA 2k a D dhLh H x xA aHA化简得：h 0 y yA a4G dyAh H y x H A a A2k a D dhL从而得到关于 yA 的一个二阶常微分方程及其边界条件：只供学习与交流此文档仅供收集于网络，如有侵权请联系网站

40、删除2 24G L d y L dy k DA A LG2 2k a D H dh H dh 4L A A0h 0 y yA a4G dyAh H y x H A a A2k a D dhL这是一个常系数二阶常微分方程， 将实际生产中的数据代入后可以很容易解出关于组分 A 在气相中的浓度 yA 沿塔高方向上的分布， 进而得到组分 A 在液相中的浓度 xA 沿塔高方向上的分布。12. 在一鼓泡容器内， 初始时刻装有体积为 V 的氨盐水 (NH 3NaCl H2O)，随后以流量 F通入 CO2 气体，在液相中即发生以下碳酸化反应CO2NH3k1 NH COO H (1)2NH COO2HO2k2

41、HCO NH33(2)HCO Na3k3 NaHCO3(3)这是制取纯碱 (Na2HCO3)的基本反应，是一个连串反应过程。其中反应 (1)为快速反应，发生在气液界面附近的液膜之内，而反应 (2) 、(3)为慢反应，发生在液相本体之中。设各步反应均为拟一级反应，其中 NH3 和 Na+ 大大过量，其浓度可近似考虑为常数。容器中的气含率 、气液比表面积 aV 和液膜厚度 均为已知量，液膜体积 aV1，试分别对液膜和液相本体导出该反应过程的数学模型，给出分批式和连续式两种操作情况下的反应器模型。提示： 液膜中的反应 -扩散过程可视为拟稳态过程， CO2 的气液界面浓度和本体浓度可设为已知， NH2

42、COO-的本体衡算方程需考虑液膜向液相本体的传递速率。解： 分析： 在此串连反应中 ,反应（ 1）为快反应，其反应速率远大于反应物 CO2 在 液相中的扩散速率，因而该反应只发生在气液界面附近的液膜之内，在液相本体中， CO2的浓度已为零。反应（ 2）、（3）为慢反应，因而液相中反应物的混合速率远大于反应速率，即由液膜中通过反应（ 1）生成的 NH 2COO- 在还未发生显著的化学反应之前就已混合均匀，可以认为NH2COO-以及其反应产物 HCO 3-在整个反应器中均匀分布。 而在反应器中， 液膜体积 av1，即液相本体的体积远大于液膜体积，所以反应（ 2）、（3）绝大部分都在液相本体中完成。

43、由于鼓泡反应器中液膜与液相本体中分别进行着不同的反应（对于快反应（ 1），主要发生在液膜中，而对于慢反应（ 2）、（3）则主要发生在液相本体中） ，因而需要分开讨论。对于液膜，可以采用拟稳态假设，即不论是分批式还是连续式操作，液膜内的物质浓度都不随时间变化，积累量均为零。对于液相本体，则需将分批式与连续式操作区分对待。此题中鼓泡反应器所涉及到的物质包括有： CO2、NH 3、NH2COO-、H+、H2O、HCO3-、+及 NaHCO 3 共八种，并且其浓度分布包括在液膜内浓度分布及液相本体中的浓度分布，共 Na涉及有八个变量在分批式及连续式两种操作情况下在两个不同相 （液膜及液相本体） 中的浓

44、度分布。但是，通过分析，可以先排除许多量：1 CO2：在本反应器中， CO2 作为反应物参加快反应（ 1），由前面的分析可知，CO2在液相本体中的浓度为零， 因而我们只需要求出在液膜中 CO2浓度分布即可，根据拟稳态假设， 所得到的浓度分布既可用于分批式操作， 也可用于连续式操作。2 NH 3、H2O、Na+：在本反应器中，这三种反应物都处于大大过量，其浓度可近似考虑为常数，因而可以在反应器模型中将其作为已知量处理，并视为常数。3 NH 2COO -、HCO3-：由前面的分析可知，这两种物质在液相本体中均处于均匀只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除分布，并且它们参与的反应

45、（ 2）、（3）都是在液相本体中完成，在液膜中的反应可以忽略，即液膜中 NH2COO-、HCO 3-的浓度分布对于反应器模型没有任影响。因而我们只需求得 NH 2COO-和 HCO3-在液相本体中的浓度随时间的分布即可 （在空间上均匀分布） ，而二者在液膜中的浓度分布无需涉及。+、NaHCO3：在本反应器中，这两种物质均只出现在反应产物中，并没有以 4 H反应物的形式参加反应， 并且在反应速率方程中也没有涉及，因而在建立的反应器模型中不涉及两者的浓度。综上所述，本反应器模型共涉及： CO2在液膜中的浓度分布， NH2COO-和 HCO 3-在液相本体中浓度随时间的分布共三个变量。 其中 CO2

46、在液膜中的浓度分布即可用于分批式操作，又用于连续式操作。而 NH 2COO -和 HCO 3-在液相本体中浓度随时间的分布则必须分为分批式操作和连续式操作进行讨论。假设1 在本题中，用符号 A 指代组分 CO2，符号 B 指代 NH2COO ，符号 C 指代 HCO3-。2 反应（ 1）、（2）、（3）均为拟一级反应，因而其反应速率分别为：dccor k c c21 1 co NHdt2 3dcNH COOr k c c22 2 H ONH COOdt22dcHCOr k c c33 3 HCO Nadt33 反应器体积为 V，而反应器中气含率为，易知反应器中的液体量为：（1-）V 。再由反应

47、器中的气液比表面积为 av（单位：有效气液接触面积 m2/液体体积 m3），得到反应器中的液膜面积（即气液接触面积）为 av（1- ）V ，液膜体积为： av（1- ）V。4 假设在反应过程中反应器内的液体体积不变，即忽略反应过程中的体积变化，因此在质量衡算方程中反应器内的液体休积（ 1-）V 为常数。解答步骤1 液膜中 CO2 的浓度分布：液膜中 CO2 的反应 -扩散过程可视为拟稳态过程，因而不论是分批式还是连续式均可以作为稳态处理。由于液膜厚度远小于气泡直径， 因而可以单位面积对液膜内的 CO2（用符号 A 代表）作质量守恒：输入项 =DAdCAdx1式中 1 表示单位面积， DA 为 CO2 在液膜中的传质系数， 由于液膜内始终处于稳态，浓度 CA 不随时间变化，因而 CA 仅是坐标 x 的函数， CO2 的质量守恒方程为一常微分方程。输出项2dC d dc dC d cA A A AD D dx 1 D D dx 1A A A A2dx dx dx dx dx只供学习与交流此文档仅供收集于网络，如有侵权请联系网站删除生成项 =k1cAcNH dx 13积累项 =0根据质量守恒：输入项 输出项 生成累 积累项2 d cAD dx k c cA 2 1 A NHdx3式中 CNH3 为液膜中 NH3 的浓度，