化工数学答案(全)

1、化工数学各章习题选解（仅供参考）第一章习题1.（V）在一个有效容积为V的半连续式搅拌反应器中，由原料A生产物质B,若浓度为5流量为Q的A溶液加入空反应器，反应遵循以下连串可逆步骤且所有的反应均为一级，证明在反应器中B的克分子数是以下微分方程的解dr+尊dt=C式中P=k2+k5R=血c=Qcoh证明：对A、B分别作质量衡算，有A：c.Q-k.Nk.N-（1）B：KNA-k2NB-k.NB=-由（2）得到：（IN&N严弦+十寻（3）（3）代入（2）,得:k&Q-灯儿=&+&)牛+i-iLP=&+kyR=kRyC=cqQd2NdN得令+P牛+RNb=C(5)

2、clrdt证毕。2冬天的池塘水面上结了一层厚度为/的冰层，冰层上方与温度为几的空气接触，下方与温度为0°C的池水接触。当7；.V0°C时，水的热量将通过冰层向空气中散发，散发的热量转化为冰层增加的厚度。已知水结冰的相变潜热为厶，冰的密度为P,导热系数为k,导温系数为Q,求：1) 当气温7；不随时间变化时，给出冰层厚度随时间变化的关系，若Lf=3.35X105J/kg,P=913kg/mk=2.22W/m°K,7；V=-1O°C，问冰冻三尺，需几口之寒？2) 当气温随时间变化时，设G=7(t)已知，导出冰层厚度变化的完整数学模型。解：(1)冰层的温度为O&

3、#39;C,水通过冰层向空气散发热量，记为0,该热量用于水结成冰。假设冰层面积为s，厚度为/根据导热方程，可得：Q=k(Jw-dt=Lfpsdl代入数值，Lf=3.35X105J/kg,P=913kg/m(V)在一个半分批式搅拌反应器中进行着一级放热化学反应，反应速率常数由Arrhenius关系式给出，反应热由釜内的冷却盘管移出，请自行设定有关的参数，导出该反应器的数学模型。解：设物料以恒定的体积流量F加入，则反应器中反应物浓度G与温度T由以下物料衡算与热量衡算方程给出物料衡算方程,k=2.22W/m°K,G=10°C,/=lm,求解积分上式得：fdt=f1dtJo913x

4、3.35xl05J。t=79.7天=80天若冰冻三尺，在几=一10°C时，需要约80天。若几一几，冰层厚度为/根据热量守恒：Q=k(Jw-dt=LfpsdlLtpldl=kTwdt两边枳分：=kTwdt厚度变化与几的关系为：”«zd（vcA）FCq+VrA=rA=kCAV=V0+Ft能屋衡算方程FpC/o-K(T-Te)A+kC(-AHr)=(2)合并，得数学模型为FPC血-K(TT)A+kCAV(-AH>.)=叮"<V=VFt(3)Ek*exp(-)KI5(0)=o,r(o)=tq式中K(T-Tc)A为冷却移热，kCAV(-Hr)为反应热。4. (V

5、)采用微元分析法推导出柱坐标系中的不定常热传导方程。解：考虑柱坐标系中热传导方程的形式。柱坐标系下的三个空间变量：向径儿经度角8,高度乙。在这三个方向上，与自变量的微分变化所对应的线段微元长度分别是dr=(dr,rdO.dz)(4)由偏导数的定义，温度梯度FT在三个方向的分量即温度在每个方向上的微元增量除以相应的线元长度，即竺竺竺drrdOdZ于是Fourier热传导定律在柱坐标系中的分量形式为接着考虑各方向输入和输出的微元通量，首先考虑r方向输入项:-k(rdO)(dz)(7)dr输出项:一R(rdO)(dz)-k(r)drd6dz(8)drdrdr于是方向的净输入通量为：k沁昨对0方向作同

6、样的分析，输入项:-k（d厂）（衣）rd6(10)prra输出项:一R-(dr)(dz)-k一-(-)drdOdzrdOrdOdO(11)净输入通量：kdCr）drdOdzrdOdO(12)AT输入项:_k一dz(13)z方向的分析，输出项:-k(rd0)(dr)-k(r)drdOdzdzdzdz(14)净输入通量：kr)drdOdz&dz(15)微元体内的积累项：pCpdri-dO）dz(16)将三个方向输入微元的热流净增量加和并令其等于枳累项，就得到孔卅2(止)+A咚+咚dtrdrdrrd0dv5. 风吹过皮肤表面时，人会有干燥凉爽的感觉，这是因为风的吹拂强化皮肤表面的对流传热与传

7、质，形成一个速度，温度，浓度（含水量）的边界层，设流动为层流（微风），考虑出汗的蒸发潜热，求：1）列出皮肤表面的三传问题的边界层方程，根据实际情况适当简化并给出问题的边界条件；2）将上述问题无量纲化，并解释所得到的各无量纲参数的物理意义；3）试分析速度分布，温度分布，含水量分布分别与哪些无量纲参数有关，并用简单的函数关系示意；4）根据所得结果定性的解释一些经验常识：为什么风越人越感觉到冷？为什么出汗后擦了汗感觉更凉快？当空气中湿度变化时，对表面散热会带来哪些影响？在冬天和夏天，人体对空气湿度的增加会有什么样的感觉？解：1）同时考虑流动传热传质时的边界层传递方程是0（"字+U竺）=-丁

8、+“马+pgg-TJ+Qg,g（C-cx）oxbyoxdypCp("M+U竺)=k马+“(字)，+aHi;oxovcv"dyzde.de.、*少c：表示重力在x方向的分量，0为热膨胀系数，f为密度变化系数/水汽化潜热£水蒸发速度由于歿可忽略，r=0,（字）2可忽略，dxdy化简后zdudu、d2uP（M+（?）=/ox莎dy応p（"責+公）=*石严比zde.de.、=52c-。（忖+u才m=oxcycv边界条件y=0,u=0（皮肤表面气流速度）T=T0（皮肤表面温度）c=c0（皮肤表面的含水量）y=6】u=uoo（速度边界层外气流速度）y=2T=T8（温

9、度边界层外气流温度）y=$3c=c8（浓度边界层外气流中含水量浓度）6】，§2，心分别为速度边界层，温度边界层，浓度边界层的厚度。2）无量纲化无量纲物理性质的比值=-=iAT=-=?ra=Sc无屋纲化后ann/v1ft,)vdx卜丄丄ddyvronTnrw1-+htu=7dxdyanrHuc+n.t>c-cdxdyc“52nv46厂k52nTaHi-边界条件y=o.n=o在y=oc,口=1对于较人的P或Sc,热传导与扩散效应与黏性比较相对较弱，热边界层和扩散边界层位于速度边界层内部，反之，对于较小的Pr或Sc,热传导与扩散速率大于黏性传递速率，热和扩散边界层就有可能扩展到速度边

10、界层之外。3）速度分布，温度分布，含水量分布的简单函数关系式幺=i（nv>i）-i）n7（n7.<i）（31）-=u（-i）nt.（nc<i）-=i（口i）4）风越犬，皮肤表面的气体更新速度越快，水的蒸发速度变快，传热越快，感觉到冷出汗后感觉更凉快，是因为减小了汗水层的厚度，蒸发速度加快当空气中湿度变人时，皮肤表面水的蒸发速度变慢，不利于传热夏天空气湿度增加，汗水蒸发困难，人感觉闷热冬天空气湿度增加，少量的汗水在皮肤表面使人感觉温暖。6. （V）在管式反应器模型（1.4.15）中，当PeTO时，相当于完全返混的情况。试从方程（4.15）出发，通过适当的体积积分和取极限PeTO

11、,导出均相釜式反应器模型。解：当PeTO时，由原方程（4.15）及边界条件可知，c=const,说明在完全返混的情况卞，反应器内具有均匀的浓度。对于任意的Peclet数，对方程4.15进行体积积分得到-Dacdz+西=_亦:+丄於七dOPedz010式中忘为反应器内的平均浓度。将边界条件（4.15）代入（31）,得到de<clO=-Dac-cz=i+l(32)。二0,c(0)=0上式对任意Peclet数均成立，仅当P0时，反应器内浓度均匀，c=c“上式成为无量纲的理想混合釜式反应器数学模型。7. (V)烯坯在Zieglar-Natta催化剂颗粒上的气相聚合过程可用最简单的固体核模型来描述

12、，如附图所示。气相中的烯坯单体在催化剂颗粒(图中阴影部分)表面聚合后生成一多孔的固体聚合物壳层并将催化剂包裹在内部，外部的气相烯坯单体只有扩散穿过此固体聚合物壳层后才能到达催化剂表面参与反应。试求：(i) 证明单体在壳层中的扩散及聚合物粒子的生长由以卞方程描述3M_1Q,r3M、=L)(r-)&t厂&rdR_MWDQMdT=oTI式中M为单体浓度(mol/m3),j为聚合物壳层的密度(kg/m3),D为单体在壳层中的扩散系数(m2/s),M缈为单体的分子量，R为聚合物颗粒的半径。(U)设催化剂核半径为比，单体在外部气相本体中的浓度为以上述参量为r和M的特征尺度，并引入适当的时间

13、尺度，将上述方程无量纲化。然后根据气相单体与固体聚合物密度之间的巨大差别(Q./Qgl(P)将问题进一步简化。(皿)设单体在催化剂核表面的浓度恒为0(瞬时反应)，/?的初始值为R°(Ro>0，求解上述简化后的模型并给出聚合物粒子半径R随时间的变化关系。提示：对单体的浓度分布可采用拟稳态假定。解：为简化计算，令单体分子量的单位是kg(1)问题建模如图1所示，对微元也作物料衡算AA(18)(MA)dr+(AJ)dr=0dtdrJ=-D,4=4龙厂dr得如图2对微元dR作物料衡算dtrdrdr(19)47rR'p、dR11drf=R(20)(21)dR_M、D0Mdtpsdr

14、（2）无量纲化与简化分析：本问题存在着两个特征时间尺度，一个是单体组分内扩散通过聚合物壳层的时间尺度r2McMcMq亠,该尺度可以从内扩散方程（19）中得出，为此，近似取代替,M/r代替,Drxdtdrd2Mr2M/卢代替就可估算出=另一个是聚合物颗粒生长的特征时间尺度忌可以从方6厂1D程（20）中用类似的比值代替微分的办法估算出J=在对问题进行无量纲化时,2不同时间尺度的选择代表着所关注的不同过程。如果选取"=¥为时间尺度，式（20）和（21）可分别无量纲化为（仍然用当前变量表示无量纲变量）：6M1d.dMr八0q）（22）dtrordrdRdMI”、=（23）dtdrr

15、cp=mhmw,=<1Ps此时式（23）中出现一个小参数时间尺度珀称为快时间尺度，选择这一尺度所得到的方程（22）中不含有小参数&表示我们关注的是单体M通过聚合物壳层的不定常扩散而不是粒子的生长。略去（23）中的小参数项后得到R=const，说明在考虑单体内扩散时，由于时间较短，可以将粒子半径作为常数考虑。因此，选择6为时间尺度显然不妥，得到的不是我们希望关注的问题。如果选取，务为吋间尺度，3、皿式可无量纲化为:dM£dtdMdr)(25)dRdMdldr(26)此时粒子生长方程（26）不含小参数，粒径将随时间变化.表示我们关注的是颗粒的生长。而单体内扩散方程（25）中

16、的时间导数项含小参数&可以略去，说明在慢时间尺度邑上考虑粒子生长时，单体的内扩散过程可以忽略时间变化项，内扩散可以作为拟稳态过程来考虑。从（D和中得到的不同简化模型说明时间尺度的选择需要根据建模目的来考虑，使简化后的模型能够代表所关注的过程的主要特征。对单体的浓度分布作拟稳态（时间导数项为零）假设，即（25）中的葺"得1dr=tM=0(27)解得代入（26）,得(28)dR_rrdtRR-rc)R(0)=£(29)解得&在缺乏数学模型的某些情况下,仅仅根据量纲分析或尺度比较也可以获得一些很有价值的结果，考虑以下例子：1）对于固体颗粒在黏性流体中的Stock流

17、动问题，颗粒受到的阻力f仅仅与颗粒尺度d,动力学粘度卩和速度U有关，即f=f（d,U,U）根据量纲齐次化的要求，物理方程等式两边的量纲应该相同，而有参数d,口，u组成的具有离地量纲的参量只可能是duu,因此上述函数关系只可能取一下形式，f=AdUu式中A是一个只与颗粒形状有关的常数，上式即为Stock定律。现根据上述量纲分析方法分析湍流的消磁度运动。湍流中存在一系列人小不同的涡旋，能屋从人尺度涡旋顺序传递给消磁度涡旋，同时将机械能耗散为热能，其中最小的涡旋尺度称为Kolmogorov尺度，在这个尺度上，黏性和能量耗散占优，因此只有运动学粘度v（nr/s）和能量耗散速率e（W/Kg）两个产量起作

18、用，其他物理量都可以用这两个量表示。试根据量纲齐次化原理推导出Kolmogorov尺度X及局部速度u'与v,e的关系（可相差一个常数）U.习题8：湍动射流2）流体在自由空间中的射流形成一个夹角为a的圆锥型区域，如图所示，设U=U（x）为距喷IIx处的平均流速，R=R（x）为x处的射流半径，试根据总动量pU2R2沿x方向守恒的要求确定速度U和射流区总流量g=UR2沿x的变化关系（可相差一个常数）3）对于放热反应，当反应器尺寸增人时，其体积按长度的三次方增长，而表面积却按平方增长，因此体积增人有利于热量的增加，而体积减小有利于冷却散热。这是化学工程中说明“放人效应”的一个典型例子。根据类似

19、的道理解释为什么生活在寒冷地区的动物一般体型较X（例如北方人就比南方人高人），而且形状趋于圆滑，而热带地区的动物体型较小且趋于瘦长（例如南方人比北方人相对较瘦，且身体凸出部分的轮廓更为明显）。解：1）分析：题中出现的符号意义如下Kolmogorov尺度入m能量耗散速率£W/Kg运动学粘度vnr/s局部速度qni/s从f=AdUu可知f量纲与dUu相同，又£的量纲W/Kg与马的量纲相同Pv（f力一N,u速度一m/s,P密度一Kg/nP,V体枳一）可知V与入3量纲相同，f与入vq量纲相同，u与八量纲相同,v2）分析：总动量pU2R2沿x方向守恒，设pU2R2=C（C为常数）。又

20、R=fgax,将R带入总动量表达式中,pU'Qgax)2=Ccp(/gax)Q=UR2=Cz宀小fgax式中Q为气液表面张力,o为气液界面与固壁之间的接触角，R为毛细管半径。3）分析：将人体看成一个反应器，食物在体内消化放出热量，除供人体正常活动所需的能量和储存在体内外，以热量的形式通过体表释放到体外。人体散热与人体表面枳和外界温度和人体温度的差值成正比，表面积与人体尺寸的平方成正比。可记为散热=A*温度差*表面积。人体放热与人体体积成正比，体积与人体尺寸的立方成正比，可记为放热=E*体积（A,B为关于人体散热放热的常数）。为保持人体温度一定，放热与散热需要平衡，散热=放热。寒冷地区的

21、外界温度和人体温度的差值比热带地区人，所以需要的体枳与表面枳比也人，即表现为人体尺寸人，也就是人的体型较人，而形状圆滑是为了减小表面枳所致。热带地区的人体型瘦长，身体凸出部分轮廓明显，可以增人表面积，方便散热。9.在水平液面上垂直插入一个半径为R的毛细管，此时液体将在表面张力的拉动下沿着管中上升。弯曲液面形成的毛细压强可以用以下Yong-Laplace方程计算an2gCosOP=R式中C为气液表面张力，0为气液界面与固壁之间的接触角，管中流体一方面受到毛细压强的驱动而上升，一方面又受到重力和粘性阻力的作用，设流动速度遵从粘性管流的Poiseuille分布，求：1）对于两端开II的毛细管，证明液

22、位高度H随时间f的变化满足以下方程式中u为液体的动力学粘度，pg为重力。2）对于上端封闭的毛细管，设总管长为/,管内气体满足理想气体状态方程，试推导相应的液位高度H的变化方程。3）从上述方程中求出最犬液位高度和时间变化关系H（/）,据此讨论H变化的趋势。解：（1）弯曲液面形成的毛细压强可以用以FYong-Laplace方程计算，同时又受到重力的作用产生压强，总的4P为：IgCosO_R一如Poiseuille分布，体枳流率的表达式为:dV_R、7idH万厂dt_式中为液体的动力学粘度，pg为重力。pdH(JdtR2APR22aCos&(2)若毛细管上端是封闭的，则AP由三部分组成，还有

23、一部分是液面上端产生的压强。RL-H式中几为人气压强，厶为毛细管长。/dH<jdtr2apr2EHbIgCosO8/“RL-HPOH)12R8PR(L_H0(3) 当达到最高液位凤时,pdHadtR'P_SHaR2(IgCosO8HoR2aCos0PgRdH_adt8/2oRCos&-pgR'H8/z2aRCose-pgR'H积分得：(一pgR'H°+2aRCos0(l-n(2aRCos&-pgRHJ)+IcRCosOloRCosO-1)PUR('从上式可以看出，H先增高，到最大值竺竺？后开始下降。PgRio气液两相的传质

24、过程与色谱过程有许多类似之处，例如，气相通过反应器(鼓泡塔.板式塔、填料塔等）的流动可以看成是溶质通过固定相的运动，气液传质阻力可类比于气固传质阻力，气液两相的逆流操作模式也与移动床相似。此外，气液两相在界面上处于平衡状态，由Henry定律表述，与§7.2节考虑的微孔分子筛的内扩散过程类似。与色谱问题不同的是，许多气液反应器（鼓泡塔与搅拌釜）中的液相或液固两相一般都处于全混流状态，而色谱柱中固定相是静止的，移动床中固体接近平推流。试根据与移动床的类比建立如图所示的鼓泡塔反应器的稳态数学模型，图中气体从塔底加入，经分布器之后形成分散的气泡并在液体中浮升，最后从容器的上部输出；液体则从塔

25、顶加入，从底部流出。气相中的组分A被液体吸收后在液相中发生一级化学反应。鼓泡塔中气相的流动可考虑为平推流，液相考虑为全混流。其它已知的参数为：鼓泡塔液位高度I，气含率6，空塔气速U,加入液体的质量流率F,单位体积气液传质系数km，一级反应动力学常数kA,Herny系数Ha。所建立的数学模型要求包括以下内容：G习题10：鼓泡塔反应器1）设“和c分别为反应组分在气相和液相中的浓度，给出其方程和边界条件；2）如果是强放热反应，反应热通过溶剂蒸发和气液相的连续流动移出，请自行设定有关物性参数，给出温度丁满足的方程。解：数学模型一般包括物料衡算、热屋衡算和动量衡算，对于鼓泡塔的气液反应体系，气液两相的温

26、度场计算是不必要，因为塔内混合良好，温差很小。动量衡算也非必要。物料衡算是要考虑到对流、相间传质、轴向分散及化学反应等影响因素。（1）假设塔内等温，作气相及液相的物料衡算:流入项A传质项返混项积蓄项气相:液项:需+WW齐护+怜（7见传质项流入项返混项反应项积蓄项式中:假设气相平推流=0液相全混流E.a=边界条件：气相：Z=0：%=FZ=l：dFJdZ=0液相：Z=0：dFjJdZ=0Z=l:吊(+1)=吊(-1)+Dg；')切(2)进气温度为A，液体为71，气体热容“，液体热容Cp,液体的质量流量Fi，气体的质量流量件,单位时间反应热为0，C为常数。鼓泡塔内温度均匀，塔内温度和排出的气

27、体液体温度均记为T，根据能量守恒：11.填料塔广泛用于气体吸收，气液两相采用逆流操作，液体从塔顶均布后加入，沿填料表面成液膜卞降，气体从塔底加入，沿塔上升并与液体实现逆流接触，气体中的活性组分被液体吸收后从塔底流出，净化后的气体从塔顶排出，如图所示。设从塔底加入的气体中含有待吸收组分A和惰性气体，惰气流量为G(mol/s),从塔顶加入的液体惰性溶剂的流量为L(mol/s),组分A在液相中以一级反应进行分解，给定塔的直径D和塔高H、单位体积填料的液体持液量2(m3/m3)和气液传质系数kLa,以及化学反应速率常数灯、气液相Heniy系数HA,试用微元分析法建立一数学模型，描述气相浓度%(moVm

28、ol惰气)和液相浓度心(moVmol溶剂)的沿塔分布，然后从模型中消去心，得到”的单一方程，并给出适当的边界条件。提示：可假设在气液界面上满足Henry定律，则两相传质速率为kLaSHaxa).解：分析液体进厶&气体G%EH如十g必.X4液体出习题11:填料吸收塔此填料塔用于气体的化学吸收，塔内物料平衡涉及两相：气相和液相。可分别对气相及液相中的待吸收组分作质量守恒，守恒方程中将涉及到的未知量包括：待吸收组分A在气相中的浓度*(molA/mo1惰性气体)和待吸收组分A在液相中的浓度XA(moLVmol溶剂)。守恒方程数与未知量数均为2,在给定边界条件下可以得到微分方程的特解。假设1设待

29、吸收组分A在液相中的反应速率为：该反应为一级反应，反应速率常数为k(单位moLin3-s)2设吸收塔在连续操作过程中处于稳态，进而在质量衡算方程中涉及到的积累项均为0解答步骤在填料塔任一高度h处取一厚度为dh的体枳微元(见习题11图)，分别考虑微元中气相和液相的待吸收组分A的质量守恒：气相输入项二(ya+dyA)G输出项二汕+气液相间的质量传递量式中气液相间的传质方向是由气相到液相，因而传递量属于输出项，其传质速率为：kLa(yA-HAxA),其中灶是以气相摩尔分率差(-HAxA)为总传质推动力的总传质系数(单位kiiiol/(m2-sAv),a是单位体积填料层所提供的有效传质面积(单位nF传

30、质面积/nP填料体枳)。因而：气液相间的质量传递量=kLa(yA-Haxa)手dh式中D是填料塔直径，也d/?为所取微元内的填料体积。4输出项=yfi+kLci(儿-Haxa)耳-dh由于气相中不发生待吸收组分A的分解反应，因而：生成项=0吸收塔处于稳态操作，因而：积累项二0根据质量守恒：输入项-输出项+生成累=积累项即：(%+如)。-yAG+kLa(yA-HAxA)dh+0=0整理得：Gd*-kLa(yA-HAxA)-dh=04=°液相：输入项=xaL+气液相间传质量式中气液相间的传质方向是从气相到液相，因而传弟量属于输入项，它的人小与气相输出项中的气液传质量相同，为：气液相间的质

31、量传递量-kLa（yAHxxadh因而:输入项-x上+灯（比-/內）4汕输出项=（耳+如）厶生成项=-rAxdhxeL4式中心为组分A在液相中的应反应速率，单位：mol/m3-so而dhx£L是在所取填4料微元中的持液量。生成项=-kgx7TP2"T"生成项=一區竺dh4积累项=0根据质量守恒:输入项-输出项+生成累=积累项即：X/+S伉H"dh(xA+dxA)LdhQ44整理得:q闷仇-）44=0相应边界条件：在填料塔底部11=0VA=Vb式中yb为进料气体中组分A的浓度顶部h=HXA=Xa式中Xa为进料液体中组分A的浓度至此，联立关于气相与液相质量守

32、恒的两个常微分方程，并加入边界条件:=0(1-11-1)(1-11-2)求解此一阶常微分方程组即可得到气液相中组分A的浓度沿着塔高方向上的分布。卞面对此常微分方程组进一步化简，消去XA，以得到一个只含二阶常微分方程。对式(1-11-1)变形,得：J飞Ha4GclyA4GdyALdyA4Gd2yA)八kLaD2dhdhkLa;rD2dlr相应边界条件:r/?=oyA=儿kLci7rDHAd/厂(=0yA=儿4Gd*kLa7rD2dh4Gd%)4kLa7rD2dh化简得：/?=0<h=H从而得到关于VA的一个二阶常微分方程及其边界条件:kLci7rD-HA4G厶yA=ya这是一个常系数二阶常

33、微分方程，将实际生产中的数据代入后可以很容易解出关于组分A在气相中的浓度yA沿塔高方向上的分布，进而得到组分A在液相中的浓度XA沿塔高方向上的分布。12.在一鼓泡容器内，初始时刻装有体积为V的氨盐水（NH3-NaCl-H2O）,随后以流量F通入CO?气体，在液相中即发生以下碳酸化反应CO、+NH.NH,COO+HNH2COO+H.OHCO'i+NHHCO；+Na+NaHCO,这是制取纯碱（NgHCCh）的基本反应，是一个连串反应过程。其中反应（1）为快速反应，发生在气液界面附近的液膜之内，而反应（2）、（3）为慢反应，发生在液相本体之中。设各步反应均为拟一级反应，其中NHs和N犷人大过

34、量，其浓度可近似考虑为常数。容器中的气含率£、气液比表面积"和液膜厚度<5均为已知量，液膜体枳a«l，试分别对液脸和液相本体导出该反应过程的数学模型，给出分批式和连续式两种操作情况下的反应器模型。提示：液膜中的反应-扩散过程可视为拟稳态过程，CO2的气液界面浓度和本体浓度可设为已知，NH.COO-的本体衡算方程需考虑液膜向液相本体的传递速率。解：分析：在此串连反应中，反应（1）为快反应，其反应速率远人于反应物CO：在液相中的扩散速率，因而该反应只发生在气液界面附近的液膜之内，在液相本体中，co：的浓度已为零。反应（2）、（3）为慢反应，因而液相中反应物的混合

35、速率远人于反应速率，即由液脸中通过反应（1）生成的NHzCOO-在还未发生显著的化学反应之前就已混合均匀，可以认为NH2COO-以及其反应产物HCOs-在整个反应器中均匀分布。而在反应器中，液膜体积av5«1,即液相本体的体积远人于液膜体积，所以反应（2）、（3）绝人部分都在液相本体中完成。由于鼓泡反应器中液膜与液相本体中分别进行着不同的反应（对于快反应（1）,主要发生在液膜中，而对于慢反应（2）、（3）则主要发生在液相本体中），因而需要分开讨论。对于液脸，可以采用拟稳态假设，即不论是分批式还是连续式操作，液膜内的物质浓度都不随时间变化，积累量均为零。对于液相本体，则需将分批式与连续

36、式操作区分对待。此题中鼓泡反应器所涉及到的物质包扌舌有：CO,NH3、NH2COO>IT、H2O.HCO.、Na+及NaHCO3共八种，并且其浓度分布包括在液膜内浓度分布及液相本体中的浓度分布，共涉及有八个变量在分批式及连续式两种操作情况下在两个不同相（液膜及液相本体）中的浓度分布。但是，通过分析，可以先排除许多量：1CO2：在本反应器中，CO2作为反应物参加快反应（1）,由前面的分析可知，CO2在液相本体中的浓度为零，因而我们只需要求出在液脸中CO2浓度分布即可,根据拟稳态假设，所得到的浓度分布既町用于分批式操作，也町用于连续式操作。2NH3、HQ、Na+：在本反应器中，这三种反应物都

37、处于人大过量，其浓度可近似考虑为常数，因而可以在反应器模型中将其作为已知量处理，并视为常数。3NH2COO>HCOc由前面的分析可知，这两种物质在液相本体中均处于均匀分布，并且它们参与的反应（2）、（3）都是在液相本体中完成，在液脸中的反应可以忽略，即液脸中NHCOO-、HCO3-的浓度分布对于反应器模型没有任影响。因而我们只需求得NHMOO-和HCCh-在液相本体中的浓度随时间的分布即口1（在空间上均匀分布），而二者在液膜中的浓度分布无需涉及。4HJNaHCOj：在本反应器中，这两种物质均只出现在反应产物中，并没有以反应物的形式参加反应，并且在反应速率方程中也没有涉及，因而在建立的反应

38、器模型中不涉及两者的浓度。综上所述，本反应器模型共涉及：CO?在液膜中的浓度分布，NH.COO-和HCOy在液相本体中浓度随时间的分布共三个变量。其中CO2在液膜中的浓度分布即可用于分批式操作，又用于连续式操作。而NH2COOfflHCO亍在液相本体中浓度随时间的分布则必须分为分批式操作和连续式操作进行讨论。假设1在本题中,用符号A指代组分CO2,符号B指代NH2COO,符号C指代HCO3o2反应（1）、（2）、（3）均为拟一级反应，因而其反应速率分别为：%切-fi=也伽dcNH.coo-,一乙=一匕=k2CcOOCHzO.-HC还_3反应器体积为V,而反应器中气含率为£,易知反应器

39、中的液体量为：（1-£）V。再由反应器中的气液比表面积为av（单位：有效气液接触面积m-/液体体积mJ,得到反应器中的液膜面枳（即气液接触面枳）为av（l-e）V,液膜体积为：a、（1-£）V6。4假设在反应过程中反应器内的液体体枳不变，即忽略反应过程中的体枳变化，因此在质量衡算方程中反应器内的液体休枳（1-£）V为常数。解答步骤1液膜中CO：的浓度分布：液膜中CO2的反应-扩散过程可视为拟稳态过程，因而不论是分批式还是连续式均可以作为稳态处理。由于液膜厚度远小于气泡直径，因而可以单位面积对液膜内的CO2（用符号A代表）作质量守恒：输入项=_D4xldx式中1表示

40、单位面积，Da为CO?在液脸中的传质系数，由于液膜内始终处于稳态，浓度6不随时间变化，因而6仅是坐标x的函数，CO）的质量守恒方程为一常微分方程。输出项J"孕_?仏纠伸字祭亦1dxdx£L丿Idxeix）生成项二-比dxx1积累项二0根据质量守恒：输入项-输出项+生成累=积累项式中63为液膜中NHs的浓度，可视为常数。与之相应的边界条件为：fx=0cA=co2f.co为气液界面处co：的浓度L=/cA=0在液相本体中cq的浓度为o因而可以得到关于CO2(A)在液脸中浓度分布的二阶常微分方程:DAJdx=klCACNHi<x=0cA=co2tx=65=0可以解得:此式即

41、是CO三在液膜中的浓度分布，而液膜内COj总的反应速率为：-R严J一心加=J(一*心5企)小'企个液肢V企个液肢V式中-R1是在鼓泡反应器所有液膜中的反应(1)的总反应速率，单位moVso将前面得到的CO：浓度沿x方向上的分布代入体枳枳分中，很容易就可以把体积积分变为关于x的线积分，求得-Ri,但这里我们利用稳态模型，即CO?在液膜内总的反应速率-R】应等于CO：通过气液相界面，自气相向液相的传递量：J(-匕圧册J加=-。XSX字|俺个液聊农I*式中S是反应器中总的气液接触面积，即液膜面枳:s=av(1-s)V利用Ca沿x方向上的分布函数，易得到:式中负号表示co：在液膜中的反应是作为

42、反应物不断消耗的。2分批式操作中NH：COO-和HCO3-在液相本体中浓度随时间的分布：在液相本体中对NHzCOO-作质量守恒方程：输入项=尺式中-R】是在所有液膜中反应(1)的总反应速率(单位mol/s),它也就等于在所有液膜中nh2coo-的生成速率，而液膜内始终处于稳态，即所有生成的NH2COO-都从液膜传递到液相本体中，因而在液相本体中NH2COO-的输入项就是-Ri。在分批式操作中鼓泡反应器内液体没产出，因而输出项=0生成项二“一训=-k2cBcHO(1-£)V式中(1-e)V为反应器中的液相体枳，由于NH2COO-在反应(2)中是反应物，在反应过程中不断消耗，因而在生成项

43、中为负。积累项=牟dt根据质屋守恒：输入项-输出项+生成累=积累项整理得:”,匝de1+Y5才(1一0)V=K-kscy(1-£)V=KcNHDAav(1-£)vcoJ-心“(1-s)V消去总体枳(l-s)V,得到:dt卜皿yllCNH3aV1CBCHZO这是一个常微分方程，加入相应的初始条件:即在反应器的初始时刻氨盐水中不含NH.COO-dcB1+£Y5dt得到:J&Cnh、Daa、co2.Sbchzo/=0cB=0对此常微分方程求解，即可得到在分批式操作条件下nh2coo-浓度随时间的分布。同样，在液相本体中对HCOs作质量守恒：输入项=0输出项=0生

44、成项=(k“c.cHo-k.c.c十)(1一$)卩2NH.COOH)O3HCO、N&八丿积累项=牟(1Y)Udt根据质量守恒：输入项-输出项+生成累=积累项消去总体积，得到:这也是一个常微分方程，其中C可以通前面NHzCOO-的质量恒算得到,NH(K)C址。、Cw+均为常数，可由实验得到。加入相应的初始条件:t=0cc=0即在反应器的初始时刻氨盐水中不含HCO,得到：t=0cc=0对此常微分方程求解，即可得到在分批式操作条件卞HCOH农度随时间的分布。3连续式操作中NHMOO-和HCO3-在液相本体中浓度随时间的分布：假设在连续式操作过程中反应液的体枳流量为q°同样，在液相本

45、体中对NHzCOO-作质量守恒方程：输入项=尺由于进料的氨盐水中不含NH2COO-,因而输入项只包含NH.COO-自液脸扩散到液相本体中的量出。输出项=go生成项=-k2cBcHzO由于是连续操作，反应器处于稳态：积累项=0根据质量守恒：输入项-输出项+生成累=积累项得到：&-徑-他唧址。（1-级=0化简得：c严"彳+心5,。（1-叩在连续操作过程中，在液相本体里NHzCOO-的浓度不随时间变化，为一常数。同样，在液相本体中对HCO3-作质量守恒方程：由于进料的氨盐水中不含HCOy：输入项=0输出项=%生成项=（g叱。5。-显怙）（1-叩积累项=0根据质量守恒：输入项-输出项

46、+生成累=积累项得到：0_%+（k2cNHCOO.cHzO-臥怙）0化简得：在连续式操作中，在液相本体里HCOs-的浓度也不随时间变化，为一常数。综上所述：在分批式操作中，本反应器的模型为：液膜中:"q=液相本体中：NH2COOSl<de+0YnJ&Cnh、qavCO2i"2CBCH.O檢勺=0HCO；满足dtt=0df2d1-ef=0T=2CnhzCOO-CHzO*30%+cc=0在连续式操作中，本反应器的模型为:液膜中:Cg液相本体中:CNHZCOO'&9+%。(1_£)卩其中&=1一£ylCNH3AaV(1-C

47、O2z人97/<00-切：0(1一£)"q+3h(i-wW化工数学第二章习题1. 求以卞微分方程的解1 +2Vtnnx=4sinxdx2 胃+2少+2g=厂cosxd©dx3 x2+2力购+2g=©+Ind©dx解：+2ytanx=4sinxdx一I2tanxdxy=eJ解:j4sinxe2ulxdx+c=cos2xJ4sinxeCQZxV1dx+c=cos2x(4+cJCOS"X=cos2cosx+c二4cosx+Ceos2xc为任意常数(2)d£ydy-rv+2+2y=ecosxdxdx解：先求通解：特征方程为：V

48、+2X+2=0=-l+i=1z故方程通解为：ynQGcosx+C/Smx)由于a=-l+i为方程一根设特解为：=代入方程得：A=-i2)=一*ixel+i)x=-iixex(cosx+isinx)=扌xex(/cosx-suix)取其实部：7*=丄xesHix2方程的解为：y=ex(Gcosx+C2sinx)+丄xexsmxx2+2x+2y=x+lnxdxdx解：设兀=&原式可化为:孕+£+2y*+fdrdt设特解为:)'；=加特征方程为：22+2+2=0y；=At+B代入方程得:所以，方程的解为：y二acosr+C2sinr)+-z+-t-2x2z424即y=x2(

49、qcos-hix+C2sinInx)+x+Inx-2. 求解第一章给出的连续结晶器的稳态数学模型式中，成核速率B,生长速率G,流量F均可考虑为常数，加入流体的粒数分布为/的任意函数"加二m(/)。解：求解以卞方程：的解。nG分离变量得:积分得：则心启d/+GC/为常数。JG代入"表达式得原微分方程解：/=0时,"（°）=J仏曲+C严菩则原方程解为:/?=吨一用/=073. 电极加热炉中石墨电极棒的传热问题可用以卜方程描述普（DA式中D,U,4,几均为常数，但导热系数心为温度的线性函数，kT=kQ-的通解（建议：采用改进的"变换p=kT,可使求解

50、更为简捷）。dxdx，试求出上述方程解:原方程务唸2呼（T）变为:务警1117?=1+CGn=Ce（C为常数）设C=C（/）,将兀的表达式代入原微分方程，得:ttDUA(T-Tq).畔=¥门aTA=籌=+(九)+5-k仏-尹+&兮心厂_人乙打+q,其中q为任意常数;赛土2ttDUA-F3+込鱼严2将kt=k.-aT带入，变形得:dX=±dT2兀DUAk°-aT32+q积分得XT关系式为:X=±J2ttDUA伙o-aT)_a_r+alQ+kQTi_kjj+$clT+c.其中qc,为任意常数。4. 设反应物A,3在液膜中发生以下瞬时反应A+vB=P”

51、为化学计量系数，该反应受扩散限制，试导出相应的W变换并求出反应物A,3和产物P的浓度分布。解：A+vB=P此反应受扩散限制=灼D樂"D舲"叫(2)V(1)-(2)得(yDAcA-Dbcb)=0令W=vDaca-Dbcb(4)易由(3)可知W=Ax+B由定义式(4)可得W的边界条件为x=O:W=vDacao:x=SW=-Dbcbq:其中5为液膜厚度由以上边界条件可得B=呜5。人_凶50+D/boSW=_心。5s。兀+叫甌(5)8在X=XR处，5=5=0,此时(5)中W=0,得Y一创go人RvDacao+D/bo将定义式(4)两边平方，利用瞬时反应4,疗不能共存的条件知，55=

52、0,于是有W，=(vDaca-Dbcb)=(vDaca+Drcb)2将上式两边开平方，得=g+%(6)根据(4),(6),可以将°用新变量W表出(VV>0)(7)1w5=万(网-旳一，(”<°)将(5)代入(7),(8)得X+5o(0<X<AjvDaca.+DrCboV3%_r.1+叽讥Dbb°IDrCro)5DrCbo(xR<x<)产物p的浓度分布为(1)5. Ch和CCh在生物组织中的传递过程对于呼吸作用和光合作用具有根本的重要性。在生物组织中，溶解在液体(血液、组织液)中的O,通过渗透与扩散两种机制输送到组织内部供细胞呼吸

53、，研究表明，在某一临界溶氧浓度/以上，单位体积的氧消耗速率为常数，设为因此，对于厚度为/的一片组织，代谢过程中氧的衡算方程为d'cdeD-U-q=Odkax式中U为液体渗透流率，D为溶氧扩散系数。设该组织外部的溶氧浓度为常数，则边界条件町表示为x=O,x=Z:c=cQ试求出溶氧浓度沿组织厚度方向的分布，据此判断氧浓度在何处达到最小值？渗透速率"需满足什么条件才能保证在组织内部不会出现缺氧的情况(c(Q>F)?解：dedxd2cUde_qdx2DdxD齐次方程的特征方程为A2A=0D解得所以齐次方程的通解为(Uc=cl+c2exp(万x利用比较系数法，求得非齐次方程的特解C=-X(2)U所以，非齐次方程的通解为U)ac(x)=q+c?exp(万xJX边界条件为所以溶氧浓度沿组织厚度方向的分布为C(x)=cQ+c2Xu对c(x)关于X求一次，及二次导数式中c?>0dedxUc、expD令半=0得ax所以,氧浓度在A0=lllexp(%)l达到