化工问题的建模与数学分析方法练习题带答案测试题题库各章节练习题.doc

各章习题选解
（仅供参考）
第一章习题
1. （√） 在一个有效容积为V的半连续式搅拌反应器中，由原料Ａ生产物质Ｂ，若浓度为c0流量为Q的Ａ溶液加入空反应器，反应遵循以下连串-可逆步骤

且所有的反应均为一级，证明在反应器中Ｂ的克分子数NB是以下微分方程的解

式中

证明：对A、B分别作质量衡算，有
A：
B：
由（2）得到：

（3）代入（2），得：

令
得
证毕。

2. 冬天的池塘水面上结了一层厚度为l的冰层，冰层上方与温度为Tw的空气接触，下方与温度为0℃的池水接触。当Tw＜0℃时，水的热量将通过冰层向空气中散发，散发的热量转化为冰层增加的厚度。已知水结冰的相变潜热为Lf，冰的密度为ρ，导热系数为k，导温系数为α，求：
1) 当气温Tw不随时间变化时，给出冰层厚度随时间变化的关系，若Lf＝3.35×105J/kg，ρ＝913kg/m3，k＝2.22W/m °K，Tw＝－10℃，问冰冻三尺，需几日之寒？
2）当气温随时间变化时，设Tw＝Tw (t)已知，导出冰层厚度变化的完整数学模型。
解：
(1) 冰层的温度为0℃，水通过冰层向空气散发热量，记为Q，该热量用于水结成冰。假设冰层面积为s，厚度为l 根据导热方程，可得：

代入数值，Lf＝3.35×105J/kg，ρ＝913kg/m3，k＝2.22W/m °K，Tw＝－10℃，l＝1m，
求解积分上式得：

t＝79.7天≈80天
若冰冻三尺，在Tw＝－10℃时，需要约80天。

(2) 若Tw＝Tw (t)，冰层厚度为l
根据热量守恒：

两边积分：

厚度变化与Tw 的关系为：

3. （√） 在一个半分批式搅拌反应器中进行着一级放热化学反应，反应速率常数由 Arrhenius关系式给出，反应热由釜内的冷却盘管移出，请自行设定有关的参数，导出该反应器的数学模型。
解：设物料以恒定的体积流量F加入，则反应器中反应物浓度CA与温度T由以下物料衡算与热量衡算方程给出
物料衡算方程

能量衡算方程

合并，得数学模型为

式中K(T－Tc)A为冷却移热，kCAV(-ΔHr)为反应热。
4．（√）采用微元分析法推导出柱坐标系中的不定常热传导方程。
解：考虑柱坐标系中热传导方程的形式。柱坐标系下的三个空间变量：向径r，经度角q，高度z。在这三个方向上，与自变量的微分变化所对应的线段微元长度分别是

由偏导数的定义，温度梯度ÑT在三个方向的分量即温度在每个方向上的微元增量除以相应的线元长度，即

于是Fourier热传导定律在柱坐标系中的分量形式为

接着考虑各方向输入和输出的微元通量，首先考虑r方向

于是r方向的净输入通量为： （9）
对q方向作同样的分析，

z方向的分析，
微元体内的积累项： （16）
将三个方向输入微元的热流净增量加和并令其等于积累项，就得到
（17）
5. 风吹过皮肤表面时，人会有干燥凉爽的感觉，这是因为风的吹拂强化皮肤表面的对流传热与传质，形成一个速度，温度，浓度（含水量）的边界层，设流动为层流（微风），考虑出汗的蒸发潜热，求：
1）列出皮肤表面的三传问题的边界层方程，根据实际情况适当简化并给出问题的边界条件；
2）将上述问题无量纲化，并解释所得到的各无量纲参数的物理意义；
3）试分析速度分布，温度分布，含水量分布分别与哪些无量纲参数有关，并用简单的函数关系示意；
4）根据所得结果定性的解释一些经验常识：为什么风越大越感觉到冷？为什么出汗后擦了汗感觉更凉快？当空气中湿度变化时，对表面散热会带来哪些影响？在冬天和夏天，人体对空气湿度的增加会有什么样的感觉？
解：
1）同时考虑流动传热传质时的边界层传递方程是

表示重力在x方向的分量，为热膨胀系数，为密度变化系数
水汽化潜热 水蒸发速度
由于可忽略，，可忽略，
化简后

边界条件
y＝0， u＝0（皮肤表面气流速度） T＝T0（皮肤表面温度） c＝c0（皮肤表面的含水量）
y＝δ1 u＝u∞（速度边界层外气流速度）
y＝δ2 T＝T∞（温度边界层外气流温度）
y＝δ3 c＝c∞（浓度边界层外气流中含水量浓度）
δ1，δ2，δ3分别为速度边界层，温度边界层，浓度边界层的厚度。

2）无量纲化

无量纲物理性质的比值

无量纲化后

边界条件
在
对于较大的Pr或Sc，热传导与扩散效应与黏性比较相对较弱，热边界层和扩散边界层位于速度边界层内部，反之，对于较小的Pr或Sc，热传导与扩散速率大于黏性传递速率，热和扩散边界层就有可能扩展到速度边界层之外。

3）速度分布，温度分布，含水量分布的简单函数关系式

4）风越大，皮肤表面的气体更新速度越快，水的蒸发速度变快，传热越快，感觉到冷
出汗后感觉更凉快，是因为减小了汗水层的厚度，蒸发速度加快
当空气中湿度变大时，皮肤表面水的蒸发速度变慢，不利于传热
夏天空气湿度增加，汗水蒸发困难，人感觉闷热
冬天空气湿度增加，少量的汗水在皮肤表面使人感觉温暖。
6．（√）在管式反应器模型（1.4.15）中，当Pe®0时，相当于完全返混的情况。试从方程（4.15）出发，通过适当的体积积分和取极限Pe®0，导出均相釜式反应器模型。
解：当Pe®0时，由原方程（4.15）及边界条件可知，c=const，说明在完全返混的情况下，反应器内具有均匀的浓度。对于任意的Peclet数，对方程4.15进行体积积分
得到 （31）
式中为反应器内的平均浓度。
将边界条件（4.15）代入（31），得到
（32）
上式对任意Peclet数均成立，仅当Pe®0时，反应器内浓度均匀，，上式成为无量纲的理想混合釜式反应器数学模型。

7. （√） 烯烃在Zieglar－Natta催化剂颗粒上的气相聚合过程可用最简单的固体核模型来描述，如附图所示。气相中的烯烃单体在催化剂颗粒（图中阴影部分）表面聚合后生成一多孔的固体聚合物壳层并将催化剂包裹在内部，外部的气相烯烃单体只有扩散穿过此固体聚合物壳层后才能到达催化剂表面参与反应。试求：
（i）证明单体在壳层中的扩散及聚合物粒子的生长由以下方程描述

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

得
如图2.对微元dR作物料衡算

得

图1
图2

（2） 无量纲化与简化
分析：本问题存在着两个特征时间尺度，一个是单体组分内扩散通过聚合物壳层的时间尺度，该尺度可以从内扩散方程（19）中得出，为此，近似取代替，M/r代替，M/r2代替，就可估算出；另一个是聚合物颗粒生长的特征时间尺度t2，可以从方程（20）中用类似的比值代替微分的办法估算出。在对问题进行无量纲化时，不同时间尺度的选择代表着所关注的不同过程。
如果选取为时间尺度，式（20）和（21）可分别无量纲化为（仍然用当前变量表示无量纲变量）：

此时式（23）中出现一个小参数e 。时间尺度t1称为快时间尺度，选择这一尺度所得到的方程（22）中不含有小参数e，表示我们关注的是单体M通过聚合物壳层的不定常扩散而不是粒子的生长。略去（23）中的小参数项后得到，说明在考虑单体内扩散时，由于时间较短，可以将粒子半径作为常数考虑。
因此，选择t1为时间尺度显然不妥，得到的不是我们希望关注的问题。
如果选取为时间尺度，（20）、（21）式可无量纲化为：

此时粒子生长方程（26）不含小参数，粒径将随时间变化，表示我们关注的是颗粒的生长。而单体内扩散方程（25）中的时间导数项含小参数e，可以略去，说明在慢时间尺度t2上考虑粒子生长时，单体的内扩散过程可以忽略时间变化项，内扩散可以作为拟稳态过程来考虑。
从和中得到的不同简化模型说明时间尺度的选择需要根据建模目的来考虑，使简化后的模型能够代表所关注的过程的主要特征。
（3）对单体的浓度分布作拟稳态（时间导数项为零）假设，即（25）中的，得

解得

代入（26），得

解得

8. 在缺乏数学模型的某些情况下，仅仅根据量纲分析或尺度比较也可以获得一些很有价值的结果，考虑以下例子：
1） 对于固体颗粒在黏性流体中的Stock流动问题，颗粒受到的阻力f仅仅与颗粒尺度d，动力学粘度μ和速度υ有关，即
f＝f（d，μ，υ）
根据量纲齐次化的要求，物理方程等式两边的量纲应该相同，而有参数d，μ，υ组成的具有离地量纲的参量只可能是dμυ，因此上述函数关系只可能取一下形式，
f＝Adμυ
式中A是一个只与颗粒形状有关的常数，上式即为Stock定律。
现根据上述量纲分析方法分析湍流的消磁度运动。湍流中存在一系列大小不同的涡旋，能量从大尺度涡旋顺序传递给消磁度涡旋，同时将机械能耗散为热能，其中最小的涡旋尺度称为Kolmogorov尺度，在这个尺度上，黏性和能量耗散占优，因此只有运动学粘度v (m2/s)和能量耗散速率ε(W/Kg）两个产量起作用，其他物理量都可以用这两个量表示。试根据量纲齐次化原理推导出Kolmogorov尺度λ及局部速度υλ与v，ε的关系（可相差一个常数）
习题8：湍动射流
2）流体在自由空间中的射流形成一个夹角为α的圆锥型区域，如图所示，设U＝U（x）为距喷口x处的平均流速，R＝R（x）为x处的射流半径，试根据总动量沿x方向守恒的要求确定速度U和射流区总流量沿x的变化关系（可相差一个常数）
3）对于放热反应，当反应器尺寸增大时，其体积按长度的三次方增长，而表面积却按平方增长，因此体积增大有利于热量的增加，而体积减小有利于冷却散热。这是化学工程中说明“放大效应”的一个典型例子。根据类似的道理解释为什么生活在寒冷地区的动物一般体型较大（例如北方人就比南方人高大），而且形状趋于圆滑，而热带地区的动物体型较小且趋于瘦长（例如南方人比北方人相对较瘦，且身体凸出部分的轮廓更为明显）。
解：
1）分析：题中出现的符号意义如下
Kolmogorov尺度λ m
能量耗散速率ε W/Kg
运动学粘度v m2/s
局部速度 m/s

从f＝Adμυ可知f量纲与dμυ相同，
又ε的量纲W/Kg与的量纲相同
（f力－N，υ速度－m/s，ρ密度－Kg/m3，V体积－m3）
可知V与λ3量纲相同，f与λv 量纲相同，υ与υλ量纲相同，
带入ε＝中

2）分析：
总动量沿x方向守恒，设＝C（C为常数）。
又，将R带入总动量表达式中，

3）分析：
将人体看成一个反应器，食物在体内消化放出热量，除供人体正常活动所需的能量和储存在体内外，以热量的形式通过体表释放到体外。
人体散热与人体表面积和外界温度和人体温度的差值成正比，表面积与人体尺寸的平方成正比。可记为散热＝A*温度差*表面积。
人体放热与人体体积成正比，体积与人体尺寸的立方成正比，可记为放热＝B*体积（A，B为关于人体散热放热的常数）。
为保持人体温度一定，放热与散热需要平衡，散热＝放热。
寒冷地区的外界温度和人体温度的差值比热带地区大，所以需要的体积与表面积比也大，即表现为人体尺寸大，也就是人的体型较大，而形状圆滑是为了减小表面积所致。
热带地区的人体型瘦长，身体凸出部分轮廓明显，可以增大表面积，方便散热。

9. 在水平液面上垂直插入一个半径为R的毛细管，此时液体将在表面张力的拉动下沿着管中上升。弯曲液面形成的毛细压强可以用以下Yong－Laplace 方程计算

式中σ为气液表面张力，θ为气液界面与固壁之间的接触角，管中流体一方面受到毛细压强的驱动而上升，一方面又受到重力和粘性阻力的作用，设流动速度遵从粘性管流的Poiseuille分布，求：
1) 对于两端开口的毛细管，证明液位高度H随时间t的变化满足以下方程

式中μ为液体的动力学粘度，ρg为重力。
2）对于上端封闭的毛细管，设总管长为l，管内气体满足理想气体状态方程，试推导相应的液位高度H的变化方程。
3）从上述方程中求出最大液位高度Ho和时间变化关系H(t),据此讨论H变化的趋势。

解：
(1) 弯曲液面形成的毛细压强可以用以下Yong－Laplace 方程计算，同时又受到重力的作用产生压强，总的ΔP为：

式中σ为气液表面张力，θ为气液界面与固壁之间的接触角，R为毛细管半径。
Poiseuille分布，体积流率的表达式为：

式中μ为液体的动力学粘度，ρg为重力。

(2) 若毛细管上端是封闭的，则ΔP由三部分组成，还有一部分是液面上端产生的压强。
式中Po为大气压强，L为毛细管长。

(3) 当达到最高液位Ho时，，则：
＝0

积分得：
从上式可以看出，H先增高，到最大值后开始下降。

10. 气液两相的传质过程与色谱过程有许多类似之处，例如，气相通过反应器（鼓泡塔、板式塔、填料塔等）的流动可以看成是溶质通过固定相的运动，气液传质阻力可类比于气固传质阻力，气液两相的逆流操作模式也与移动床相似。此外，气液两相在界面上处于平衡状态，由Henry定律表述，与§7.2节考虑的微孔分子筛的内扩散过程类似。与色谱问题不同的是，许多气液反应器（鼓泡塔与搅拌釜）中的液相或液固两相一般都处于全混流状态，而色谱柱中固定相是静止的，移动床中固体接近平推流。试根据与移动床的类比建立如图所示的鼓泡塔反应器的稳态数学模型，图中气体从塔底加入，经分布器之后形成分散的气泡并在液体中浮升，最后从容器的上部输出；液体则从塔顶加入，从底部流出。气相中的组分A被液体吸收后在液相中发生一级化学反应。鼓泡塔中气相的流动可考虑为平推流，液相考虑为全混流。其它已知的参数为：鼓泡塔液位高度l，气含率εg，空塔气速U，加入液体的质量流率F，单位体积气液传质系数kLa，一级反应动力学常数kA，Henry系数HA。所建立的数学模型要求包括以下内容：
习题10：鼓泡塔反 应器
1）设cg 和cL 分别为反应组分在气相和液相中的浓度，给出其方程和边界条件；
2）如果是强放热反应，反应热通过溶剂蒸发和气液相的连续流动移出，请自行设定有关物性参数，给出温度T满足的方程。
解：数学模型一般包括物料衡算、热量衡算和动量衡算，对于鼓泡塔的气液反应体系，气液两相的温度场计算是不必要，因为塔内混合良好，温差很小。动量衡算也非必要。物料衡算是要考虑到对流、相间传质、轴向分散及化学反应等影响因素。

(1) 假设塔内等温，作气相及液相的物料衡算：

气相：
流入项 传质项 返混项 积蓄项
液项：流入项 传质项 返混项 反应项 积蓄项

式中：

假设 气相平推流 液相全混流
边界条件：
气相：Z＝0：
Z＝1：
液相：Z＝0：
Z＝1：

(2) 进气温度为Tg，液体为Tl，气体热容cg，液体热容cp，液体的质量流量Fl，气体的质量流量Fg，单位时间反应热为Q，C为常数。
鼓泡塔内温度均匀，塔内温度和排出的气体液体温度均记为T，根据能量守恒：

习题11：填料吸收塔
11．填料塔广泛用于气体吸收，气液两相采用逆流操作，液体从塔顶均布后加入，沿填料表面成液膜下降，气体从塔底加入，沿塔上升并与液体实现逆流接触，气体中的活性组分被液体吸收后从塔底流出，净化后的气体从塔顶排出，如图所示。设从塔底加入的气体中含有待吸收组分A和惰性气体，惰气流量为G（mol/s），从塔顶加入的液体惰性溶剂的流量为L（mol/s），组分A在液相中以一级反应进行分解，给定塔的直径D和塔高H、单位体积填料的液体持液量eL（m3/m3）和气液传质系数kLa，以及化学反应速率常数kA、气液相Henry系数HA，试用微元分析法建立一数学模型，描述气相浓度yA (mol/mol惰气)和液相浓度xA (mol/mol溶剂)的沿塔分布，然后从模型中消去xA，得到yA的单一方程，并给出适当的边界条件。
提示：可假设在气液界面上满足Henry定律，则两相传质速率为kLa (yA－HAxA)。
解：
分析
此填料塔用于气体的化学吸收，塔内物料平衡涉及两相：气相和液相。可分别对气相及液相中的待吸收组分作质量守恒，守恒方程中将涉及到的未知量包括：待吸收组分A在气相中的浓度yA（molA/mol惰性气体）和待吸收组分A在液相中的浓度xA(molA/mol溶剂)。守恒方程数与未知量数均为2，在给定边界条件下可以得到微分方程的特解。
假设
1 设待吸收组分A在液相中的反应速率为：

该反应为一级反应，反应速率常数为k（单位 mol/m3·s）
2 设吸收塔在连续操作过程中处于稳态，进而在质量衡算方程中涉及到的积累
项均为0
解答步骤
在填料塔任一高度h处取一厚度为dh的体积微元（见习题11图），分别考虑微元中气相和液相的待吸收组分A的质量守恒：
气相

式中 气液相间的传质方向是由气相到液相，因而传递量属于输出项，其传质速率为：，其中kL是以气相摩尔分率差为总传质推动力的总传质系数（单位 kmol/(m2·s·Δy)），a是单位体积填料层所提供的有效传质面积（单位 m2传质面积/m3填料体积）。因而：

式中 D是填料塔直径，为所取微元内的填料体积。

由于气相中不发生待吸收组分A的分解反应，因而：

吸收塔处于稳态操作，因而：

根据质量守恒：

即：

整理得：

液相：

式中气液相间的传质方向是从气相到液相，因而传弟量属于输入项，它的大小与气相输出项中的气液传质量相同，为：

因而：

式中-rA为组分A在液相中的应反应速率，单位：mol/m3·s。而是在所取填料微元中的持液量。

根据质量守恒：

即：

整理得：

相应边界条件：
在填料塔底部 h=0 yA=yb 式中yb为进料气体中组分A的浓度
顶部 h=H xA=xa 式中xa为进料液体中组分A的浓度
至此，联立关于气相与液相质量守恒的两个常微分方程，并加入边界条件：

求解此一阶常微分方程组即可得到气液相中组分A的浓度沿着塔高方向上的分布。
下面对此常微分方程组进一步化简，消去xA，以得到一个只含yA二阶常微分方程。
对式（1-11-1）变形，得：

代入式（1-11-2）得到：

相应边界条件：

化简得：

从而得到关于yA的一个二阶常微分方程及其边界条件：

这是一个常系数二阶常微分方程，将实际生产中的数据代入后可以很容易解出关于组分A在气相中的浓度yA沿塔高方向上的分布，进而得到组分A在液相中的浓度xA沿塔高方向上的分布。

12. 在一鼓泡容器内，初始时刻装有体积为V的氨盐水(NH3－NaCl－H2O)，随后以流量F通入CO2气体，在液相中即发生以下碳酸化反应
(1)
(2)
(3)
这是制取纯碱(Na2HCO3)的基本反应，是一个连串反应过程。其中反应(1)为快速反应，发生在气液界面附近的液膜之内，而反应(2)、(3)为慢反应，发生在液相本体之中。设各步反应均为拟一级反应，其中NH3 和Na+ 大大过量，其浓度可近似考虑为常数。容器中的气含率ε、气液比表面积aV和液膜厚度δ均为已知量，液膜体积aVδ<<1，试分别对液膜和液相本体导出该反应过程的数学模型，给出分批式和连续式两种操作情况下的反应器模型。
提示：液膜中的反应-扩散过程可视为拟稳态过程，CO2的气液界面浓度和本体浓度可设为已知，NH2COO-的本体衡算方程需考虑液膜向液相本体的传递速率。
解： 分析： 在此串连反应中,反应（1）为快反应，其反应速率远大于反应物CO2 在液相中的扩散速率，因而该反应只发生在气液界面附近的液膜之内，在液相本体中，CO2的浓度已为零。反应（2）、（3）为慢反应，因而液相中反应物的混合速率远大于反应速率，即由液膜中通过反应（1）生成的NH2COO-在还未发生显著的化学反应之前就已混合均匀，可以认为NH2COO-以及其反应产物HCO3-在整个反应器中均匀分布。而在反应器中，液膜体积avδ<<1，即液相本体的体积远大于液膜体积，所以反应（2）、（3）绝大部分都在液相本体中完成。
由于鼓泡反应器中液膜与液相本体中分别进行着不同的反应（对于快反应（1），主要发生在液膜中，而对于慢反应（2）、（3）则主要发生在液相本体中），因而需要分开讨论。对于液膜，可以采用拟稳态假设，即不论是分批式还是连续式操作，液膜内的物质浓度都不随时间变化，积累量均为零。对于液相本体，则需将分批式与连续式操作区分对待。
此题中鼓泡反应器所涉及到的物质包括有：CO2、NH3、NH2COO-、H+、H2O、HCO3-、Na+及NaHCO3共八种，并且其浓度分布包括在液膜内浓度分布及液相本体中的浓度分布，共涉及有八个变量在分批式及连续式两种操作情况下在两个不同相（液膜及液相本体）中的浓度分布。但是，通过分析，可以先排除许多量：
1 CO2：在本反应器中，CO2作为反应物参加快反应（1），由前面的分析可知，CO2在液相本体中的浓度为零，因而我们只需要求出在液膜中CO2浓度分布即可，根据拟稳态假设，所得到的浓度分布既可用于分批式操作，也可用于连续式操作。
2 NH3、H2O、Na+：在本反应器中，这三种反应物都处于大大过量，其浓度可近似考虑为常数，因而可以在反应器模型中将其作为已知量处理，并视为常数。
3 NH2COO-、HCO3-：由前面的分析可知，这两种物质在液相本体中均处于均匀分布，并且它们参与的反应（2）、（3）都是在液相本体中完成，在液膜中的反应可以忽略，即液膜中NH2COO-、HCO3-的浓度分布对于反应器模型没有任影响。因而我们只需求得NH2COO-和HCO3-在液相本体中的浓度随时间的分布即可（在空间上均匀分布），而二者在液膜中的浓度分布无需涉及。
4 H+、NaHCO3：在本反应器中，这两种物质均只出现在反应产物中，并没有以反应物的形式参加反应，并且在反应速率方程中也没有涉及，因而在建立的反应器模型中不涉及两者的浓度。
综上所述，本反应器模型共涉及：CO2在液膜中的浓度分布，NH2COO-和HCO3-在液相本体中浓度随时间的分布共三个变量。其中CO2在液膜中的浓度分布即可用于分批式操作，又用于连续式操作。而NH2COO-和HCO3-在液相本体中浓度随时间的分布则必须分为分批式操作和连续式操作进行讨论。
假设
1 在本题中，用符号A指代组分CO2，符号B指代NH2COO，符号C指代HCO3-。
2 反应（1）、（2）、（3）均为拟一级反应，因而其反应速率分别为：

3 反应器体积为V，而反应器中气含率为ε，易知反应器中的液体量为：
（1-ε）V。再由反应器中的气液比表面积为av（单位：有效气液接触面积m2/液体体积m3），得到反应器中的液膜面积（即气液接触面积）为av（1-ε）V，液膜体积为：av（1-ε）Vδ。
4 假设在反应过程中反应器内的液体体积不变，即忽略反应过程中的体积变化，因此在质量衡算方程中反应器内的液体休积（1-ε）V为常数。
解答步骤

1 液膜中CO2的浓度分布：
液膜中CO2的反应-扩散过程可视为拟稳态过程，因而不论是分批式还是连续式均可以作为稳态处理。
由于液膜厚度远小于气泡直径，因而可以单位面积
对液膜内的CO2（用符号A代表）作质量守恒：

式中1表示单位面积，DA为CO2在液膜中的传质系数，由于液膜内始终处于稳态，浓度CA不随时间变化，因而CA仅是坐标x的函数，CO2的质量守恒方程为一常微分方程。

根据质量守恒：

式中CNH3为液膜中NH3的浓度，可视为常数。
与之相应的边界条件为：

因而可以得到关于CO2（A）在液膜中浓度分布的二阶常微分方程：

可以解得：

此式即是CO2在液膜中的浓度分布，而液膜内CO2总的反应速率为：

式中-R1是在鼓泡反应器所有液膜中的反应（1）的总反应速率，单位mol/s。
将前面得到的CO2浓度沿x方向上的分布代入体积积分中，很容易就可以把
体积积分变为关于x的线积分，求得-R1，但这里我们利用稳态模型，即CO2
在液膜内总的反应速率-R1应等于CO2通过气液相界面，自气相向液相的传递
量：

式中S是反应器中总的气液接触面积，即液膜面积：

利用CA沿x方向上的分布函数，易得到：

式中负号表示CO2在液膜中的反应是作为反应物不断消耗的。
2 分批式操作中NH2COO-和HCO3-在液相本体中浓度随时间的分布：
在液相本体中对NH2COO-作质量守恒方程：

式中-R1是在所有液膜中反应（1）的总反应速率（单位mol/s），它也就等于在所有液膜中NH2COO-的生成速率，而液膜内始终处于稳态，即所有生成的NH2COO-都从液膜传递到液相本体中，因而在液相本体中NH2COO-的输入项就是-R1。
在分批式操作中鼓泡反应器内液体没产出，因而

式中（1-ε）V为反应器中的液相体积，由于NH2COO-在反应（2）中是反应物，在反应过程中不断消耗，因而在生成项中为负。

根据质量守恒：

整理得：

消去总体积，得到：

这是一个常微分方程，加入相应的初始条件：

得到：

对此常微分方程求解，即可得到在分批式操作条件下NH2COO-浓度随时间的分布。
同样，在液相本体中对HCO3-作质量守恒：

根据质量守恒：

得到：

消去总体积，得到：

这也是一个常微分方程，其中可以通前面NH2COO-的质量恒算得到，、均为常数，可由实验得到。加入相应的初始条件：

得到：

对此常微分方程求解，即可得到在分批式操作条件下HCO3-浓度随时间的分布。
3 连续式操作中NH2COO-和HCO3-在液相本体中浓度随时间的分布：
假设在连续式操作过程中反应液的体积流量为q。
同样，在液相本体中对NH2COO-作质量守恒方程：

由于进料的氨盐水中不含NH2COO-，因而输入项只包含NH2COO-自液膜扩散到液相本体中的量R1。

由于是连续操作，反应器处于稳态：

根据质量守恒：

得到：

化简得：

在连续操作过程中，在液相本体里NH2COO-的浓度不随时间变化，为一常数。
同样，在液相本体中对HCO3-作质量守恒方程：
由于进料的氨盐水中不含HCO3-：

根据质量守恒：

得到：

化简得：

在连续式操作中，在液相本体里HCO3-的浓度也不随时间变化，为一常数。
综上所述：
在分批式操作中，本反应器的模型为：

在连续式操作中，本反应器的模型为：

第二章习题
1. 求以下微分方程的解

解：
(1)

（2）

（3）

2. 求解第一章给出的连续结晶器的稳态数学模型

式中，成核速率B，生长速率G，流量F均可考虑为常数，加入流体的粒数分布 为l的任意函数nin=nin(l)。
解：
求解以下方程：

的解。
分离变量得：
积分得：
（C为常数）
设，将n的表达式代入原微分方程，得：

则，C1为常数。
代入n表达式得原微分方程解：

l＝0时，

则原方程解为：
3. 电极加热炉中石墨电极棒的传热问题可用以下方程描述

式中D，U，A，T0均为常数，但导热系数kＴ为温度的线性函数，kT = k0 -aT，试求出上述方程的通解（建议：采用改进的p变换，可使求解更为简捷）。
解：
设,则有
原方程
变为:

, 其中为任意常数;

将带入,变形得:

积分得X~T关系式为:

其中 为任意常数。
4. 设反应物A，B在液膜中发生以下瞬时反应

v为化学计量系数，该反应受扩散限制，试导出相应的W-变换并求出反应物A，B和产物P的浓度分布。
解：

此反应受扩散限制

（1）
（2）
(1)-(2)得
（3）
令
（4）
易由（3）可知

由定义式（4）可得的边界条件为
；；其中为液膜厚度
由以上边界条件可得

（5）
在处，，此时（5）中，得

将定义式（4）两边平方，利用瞬时反应，不能共存的条件知，，于是有

将上式两边开平方，得
（6）
根据（4），（6），可以将，用新变量表出
，（7）
，（8）
将（5）代入（7），（8）得
，
，
产物的浓度分布为

5．O2和CO2在生物组织中的传递过程对于呼吸作用和光合作用具有根本的重要性。在生物组织中，溶解在液体（血液、组织液）中的O2通过渗透与扩散两种机制输送到组织内部供细胞呼吸，研究表明，在某一临界溶氧浓度c* 以上，单位体积的氧消耗速率为常数，设为q，因此，对于厚度为l的一片组织，代谢过程中氧的衡算方程为

式中U为液体渗透流率，D为溶氧扩散系数。设该组织外部的溶氧浓度为常数，则边界条件可表示为

试求出溶氧浓度沿组织厚度方向的分布，据此判断氧浓度在何处达到最小值？渗透速率U需满足什么条件才能保证在组织内部不会出现缺氧的情况（c(x)> c*）？
解：

齐次方程的特征方程为

解得

所以齐次方程的通解为
（1）
利用比较系数法，求得非齐次方程的特解
（2）
所以，非齐次方程的通解为
（3）
边界条件为

所以溶氧浓度沿组织厚度方向的分布为

（4）
对关于求一次，及二次导数
式中

令得

所以，氧浓度在达到最小值
为保证组织内部不会出现缺氧的情况，要求
即需要满足

6．（√）求以下变系数方程的级数解
（a）
（b）
(c)
（a）解：将幂级数（4.5）代入方程，逐项比较系数，令首项的系数为0，得到指标方程为
指标方程的两个根为，属于第一种情况，可以将c代入递推公式确定各系数。
令项的系数为0，得递推公式为

首先将代入递推公式，有

可得
则方程的第一解为
接着将代入递推公式，整理得
则方程的第二解为
最后得到方程的通解为
（b）指标方程，重根c1=c2=0
递推公式为,将递推公式表示成an对参数c的函数形式

于是含有任意参数c的幂级数y(x,c)由下式给出

当c=0时，上式给出方程的第一解

第二解即

对各项求导，并令c＝0得

所以方程的通解为
（c）将幂级数代入后比较系数得到
指标方程
递推公式
指标方程的根为
两根相差一个整数m时，递推公式中am的系数将成为0而使之无法确定。首先确定由大根得到的级数，递推公式为
然后将c1＝0代入得到第一解为

注：级数推导详见微积分下册第274页

考虑以下含任意常数c的级数y(x,c)

上述级数在c＝-2处有奇异性，第二解y2由下式给出

方程的通解为

7. 环形法兰上的散热问题可用以下方程描述

式中k和h分别为法兰的导热系数和向周边环境的传热系数，T0为环境温度。边界条件为
在内圆边界 r = r1处: T = T1
在外圆边界 r = r2处: T = T2
试用有关的Bessel函数给出上述问题的通解并说明如何由边界条件确定通解中的任意常数。
（提示：作变换，化为标准形式）
解：

8.（√）用矩阵解法求以下一阶线性微分方程组的通解，并将通解用实函数表示。
（1）
（2）
（1）解：系数矩阵A的特征方程为

解得A的特征值
当时，特征向量方程为

式中的两个方程线性相关，取x1为独立变量，令x1＝1，得到相应的特征向量为

类似的得到对应的特征向量为
因此，微分方程的通解为

（2）解：系数矩阵A的特征方程为

解得A的特征值
当时，特征向量方程为

式中的方程线性相关，取x1为独立变量，令x1＝1，得到相应的特征向量为
,相应的复数形式的解为

上述实部和虚部为方程的两个线性无关特解，因此方程组的通解为

9.（√）对于上题（1）中的系数矩阵,请用下述方法求矩阵函数expAt
（1） 待定系数法
（2） lagrange插值法
（1） 解：将特征根代入方程，得到

解得
矩阵函数为
（2） 将特征根代入方程（7.24）得

10. 如所示, 两相互联接的搅拌釜中装有体积分别为V1和V2的溶液, 初始时刻釜中溶质浓度分别为y10与y20 , 从t = 0开始, 两釜中的溶液以流量q通过管道泵送而相互交换, 管道体积可以忽略, 求两釜中的溶质浓度随时间的变化关系。

图6.1 两串联搅拌釜之间的质量交换
q
q
V1
y1
V2
y2
.习题10：两互联搅拌釜的动态响应

解：
由条件得，对第一个釜：

对第二个釜：
以上方程组成如下方程组：设，有：

其特征方程为：
展开得：
其特征值为：，为相异实根。
与对应的特征向量方程为：

取为独立变量，则。令，得特征向量：

与对应的特征向量