化工数学答案全

1、化工数学各章习题选解（仅供参考）第一章习题1. （）在一个有效容积为V的半连续式搅拌反应器中，由原料生产物质，若浓度为c0流量为Q的溶液加入空反应器，反应遵循以下连串-可逆步骤且所有的反应均为一级，证明在反应器中的克分子数是以下微分方程的解式中证明：对A、B分别作质量衡算，有A：B：由（2）得到：（3）代入（2），得：令得证毕。2. 冬天的池塘水面上结了一层厚度为l的冰层，冰层上方与温度为的空气接触，下方与温度为0的池水接触。当0时，水的热量将通过冰层向空气中散发，散发的热量转化为冰层增加的厚度。已知水结冰的相变潜热为，冰的密度为，导热系数为k，导温系数为，求：1) 当气温不随时间变化时，给出

2、冰层厚度随时间变化的关系，若3.35×105，9133，k2.22 °K，10，问冰冻三尺，需几日之寒？2）当气温随时间变化时，设 (t)已知，导出冰层厚度变化的完整数学模型。解：(1) 冰层的温度为0，水通过冰层向空气散发热量，记为Q，该热量用于水结成冰。假设冰层面积为s，厚度为l 根据导热方程，可得：代入数值，3.35×105，9133，k2.22 °K，10，l1m，求解积分上式得：t79.7天80天若冰冻三尺，在10时，需要约80天。(2) 若 (t)，冰层厚度为l根据热量守恒：两边积分：厚度变化与的关系为：3. （）在一个半分批式搅拌反应器中进

3、行着一级放热化学反应，反应速率常数由关系式给出，反应热由釜内的冷却盘管移出，请自行设定有关的参数，导出该反应器的数学模型。解：设物料以恒定的体积流量F加入，则反应器中反应物浓度与温度T由以下物料衡算与热量衡算方程给出物料衡算方程能量衡算方程合并，得数学模型为式中K(T)A为冷却移热，(-)为反应热。4（）采用微元分析法推导出柱坐标系中的不定常热传导方程。解：考虑柱坐标系中热传导方程的形式。柱坐标系下的三个空间变量：向径r，经度角q，高度z。在这三个方向上，与自变量的微分变化所对应的线段微元长度分别是 由偏导数的定义，温度梯度ÑT在三个方向的分量即温度在每个方向上的微元增量除以相应的线

4、元长度，即于是热传导定律在柱坐标系中的分量形式为接着考虑各方向输入与输出的微元通量，首先考虑r方向 于是r方向的净输入通量为： （9）对q方向作同样的分析，z方向的分析，微元体内的积累项： （16）将三个方向输入微元的热流净增量加与并令其等于积累项，就得到 （17）5.风吹过皮肤表面时，人会有干燥凉爽的感觉，这是因为风的吹拂强化皮肤表面的对流传热与传质，形成一个速度，温度，浓度（含水量）的边界层，设流动为层流（微风），考虑出汗的蒸发潜热，求：1）列出皮肤表面的三传问题的边界层方程，根据实际情况适当简化并给出问题的边界条件；2）将上述问题无量纲化，并解释所得到的各无量纲参数的物理意义；3）试分析

5、速度分布，温度分布，含水量分布分别与哪些无量纲参数有关，并用简单的函数关系示意；4）根据所得结果定性的解释一些经验常识：为什么风越大越感觉到冷？为什么出汗后擦了汗感觉更凉快？当空气中湿度变化时，对表面散热会带来哪些影响？在冬天与夏天，人体对空气湿度的增加会有什么样的感觉？解：1）同时考虑流动传热传质时的边界层传递方程是表示重力在x方向的分量，为热膨胀系数，为密度变化系数 水汽化潜热 水蒸发速度由于可忽略，可忽略， 化简后边界条件y0， u0（皮肤表面气流速度） TT0（皮肤表面温度） cc0（皮肤表面的含水量）y1 uu（速度边界层外气流速度）y2 TT（温度边界层外气流温度）y3 cc（浓度

6、边界层外气流中含水量浓度）1，2，3分别为速度边界层，温度边界层，浓度边界层的厚度。2）无量纲化无量纲物理性质的比值无量纲化后边界条件在对于较大的或，热传导与扩散效应与黏性比较相对较弱，热边界层与扩散边界层位于速度边界层内部，反之，对于较小的或，热传导与扩散速率大于黏性传递速率，热与扩散边界层就有可能扩展到速度边界层之外。3）速度分布，温度分布，含水量分布的简单函数关系式4）风越大，皮肤表面的气体更新速度越快，水的蒸发速度变快，传热越快，感觉到冷出汗后感觉更凉快，是因为减小了汗水层的厚度，蒸发速度加快当空气中湿度变大时，皮肤表面水的蒸发速度变慢，不利于传热夏天空气湿度增加，汗水蒸发困难，人感觉

7、闷热冬天空气湿度增加，少量的汗水在皮肤表面使人感觉温暖。6（）在管式反应器模型（1.4.15）中，当®0时，相当于完全返混的情况。试从方程（4.15）出发，通过适当的体积积分与取极限®0，导出均相釜式反应器模型。解：当®0时，由原方程（4.15）及边界条件可知，说明在完全返混的情况下，反应器内具有均匀的浓度。对于任意的数，对方程4.15进行体积积分得到 （31）式中为反应器内的平均浓度。将边界条件（4.15）代入（31），得到（32）上式对任意数均成立，仅当®0时，反应器内浓度均匀，上式成为无量纲的理想混合釜式反应器数学模型。7. （）烯烃在催化剂颗粒上

8、的气相聚合过程可用最简单的固体核模型来描述，如附图所示。气相中的烯烃单体在催化剂颗粒（图中阴影部分）表面聚合后生成一多孔的固体聚合物壳层并将催化剂包裹在内部，外部的气相烯烃单体只有扩散穿过此固体聚合物壳层后才能到达催化剂表面参与反应。试求： （i）证明单体在壳层中的扩散及聚合物粒子的生长由以下方程描述式中为单体浓度(3)，s为聚合物壳层的密度(3),为单体在壳层中的扩散系数(m2), M为单体的分子量，R为聚合物颗粒的半径。 （）设催化剂核半径为，单体在外部气相本体中的浓度为，以上述参量为r与M的特征尺度，并引入适当的时间尺度，将上述方程无量纲化。然后根据气相单体与固体聚合物密度之间的巨大差别

9、(g 103)将问题进一步简化。 （）设单体在催化剂核表面的浓度恒为0（瞬时反应），R的初始值为R0(R0> )，求解上述简化后的模型并给出聚合物粒子半径R随时间的变化关系。提示：对单体的浓度分布可采用拟稳态假定。解：为简化计算，令单体分子量的单位是（1） 问题建模如图1所示，对微元作物料衡算得 如图2.对微元作物料衡算得图1图2（2） 无量纲化与简化分析：本问题存在着两个特征时间尺度，一个是单体组分内扩散通过聚合物壳层的时间尺度，该尺度可以从内扩散方程（19）中得出，为此，近似取代替，代替，2代替，就可估算出；另一个是聚合物颗粒生长的特征时间尺度t2，可以从方程（20）中用类似的比值代

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

11、小参数，粒径将随时间变化，表示我们关注的是颗粒的生长。而单体内扩散方程（25）中的时间导数项含小参数e，可以略去，说明在慢时间尺度t2上考虑粒子生长时，单体的内扩散过程可以忽略时间变化项，内扩散可以作为拟稳态过程来考虑。从与中得到的不同简化模型说明时间尺度的选择需要根据建模目的来考虑，使简化后的模型能够代表所关注的过程的主要特征。（3）对单体的浓度分布作拟稳态（时间导数项为零）假设，即（25）中的，得解得代入（26），得解得8. 在缺乏数学模型的某些情况下，仅仅根据量纲分析或尺度比较也可以获得一些很有价值的结果，考虑以下例子：1） 对于固体颗粒在黏性流体中的流动问题，颗粒受到的阻力f仅仅与颗粒

12、尺度d，动力学粘度与速度有关，即ff（d，）根据量纲齐次化的要求，物理方程等式两边的量纲应该相同，而有参数d，组成的具有离地量纲的参量只可能是d，因此上述函数关系只可能取一下形式，f式中A是一个只与颗粒形状有关的常数，上式即为定律。现根据上述量纲分析方法分析湍流的消磁度运动。湍流中存在一系列大小不同的涡旋，能量从大尺度涡旋顺序传递给消磁度涡旋，同时将机械能耗散为热能，其中最小的涡旋尺度称为尺度，在这个尺度上，黏性与能量耗散占优，因此只有运动学粘度v (m2)与能量耗散速率(）两个产量起作用，其他物理量都可以用这两个量表示。试根据量纲齐次化原理推导出尺度及局部速度与v，的关系（可相差一个常数）习

13、题8：湍动射流2）流体在自由空间中的射流形成一个夹角为的圆锥型区域，如图所示，设UU（x）为距喷口x处的平均流速，RR（x）为x处的射流半径，试根据总动量沿x方向守恒的要求确定速度U与射流区总流量沿x的变化关系（可相差一个常数）3）对于放热反应，当反应器尺寸增大时，其体积按长度的三次方增长，而表面积却按平方增长，因此体积增大有利于热量的增加，而体积减小有利于冷却散热。这是化学工程中说明“放大效应”的一个典型例子。根据类似的道理解释为什么生活在寒冷地区的动物一般体型较大（例如北方人就比南方人高大），而且形状趋于圆滑，而热带地区的动物体型较小且趋于瘦长（例如南方人比北方人相对较瘦，且身体凸出部分的

14、轮廓更为明显）。解：1）分析：题中出现的符号意义如下尺度 m能量耗散速率运动学粘度v m2局部速度从f可知f量纲与d相同，又的量纲与的量纲相同（f力N，速度，密度3，V体积m3）可知V与3量纲相同，f与v量纲相同，与量纲相同，带入中2）分析：总动量沿x方向守恒，设C（C为常数）。又，将R带入总动量表达式中，3）分析：将人体看成一个反应器，食物在体内消化放出热量，除供人体正常活动所需的能量与储存在体内外，以热量的形式通过体表释放到体外。人体散热与人体表面积与外界温度与人体温度的差值成正比，表面积与人体尺寸的平方成正比。可记为散热A*温度差*表面积。人体放热与人体体积成正比，体积与人体尺寸的立方成

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

16、粘性管流的分布，求：1) 对于两端开口的毛细管，证明液位高度H随时间t的变化满足以下方程式中为液体的动力学粘度，g为重力。 2）对于上端封闭的毛细管，设总管长为l，管内气体满足理想气体状态方程，试推导相应的液位高度H的变化方程。 3）从上述方程中求出最大液位高度与时间变化关系H(t),据此讨论H变化的趋势。解：(1) 弯曲液面形成的毛细压强可以用以下 方程计算，同时又受到重力的作用产生压强，总的P为：式中为气液表面张力，为气液界面与固壁之间的接触角，R为毛细管半径。分布，体积流率的表达式为：式中为液体的动力学粘度，g为重力。(2) 若毛细管上端是封闭的，则P由三部分组成，还有一部分是液面上端产

17、生的压强。 式中为大气压强，L为毛细管长。(3) 当达到最高液位时，则：0积分得： 从上式可以看出，H先增高，到最大值后开始下降。10.气液两相的传质过程与色谱过程有许多类似之处，例如，气相通过反应器（鼓泡塔、板式塔、填料塔等）的流动可以看成是溶质通过固定相的运动，气液传质阻力可类比于气固传质阻力，气液两相的逆流操作模式也与移动床相似。此外，气液两相在界面上处于平衡状态，由定律表述，与§7.2节考虑的微孔分子筛的内扩散过程类似。与色谱问题不同的是，许多气液反应器（鼓泡塔与搅拌釜）中的液相或液固两相一般都处于全混流状态，而色谱柱中固定相是静止的，移动床中固体接近平推流。试根据与移动床的

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

19、数，给出温度T满足的方程。解：数学模型一般包括物料衡算、热量衡算与动量衡算，对于鼓泡塔的气液反应体系，气液两相的温度场计算是不必要，因为塔内混合良好，温差很小。动量衡算也非必要。物料衡算是要考虑到对流、相间传质、轴向分散及化学反应等影响因素。(1) 假设塔内等温，作气相及液相的物料衡算：气相： 流入项 传质项 返混项 积蓄项液项：流入项 传质项 返混项 反应项 积蓄项式中：假设 气相平推流 液相全混流 边界条件： 气相：Z0： Z1： 液相：Z0： Z1： (2) 进气温度为，液体为，气体热容，液体热容，液体的质量流量，气体的质量流量，单位时间反应热为Q，C为常数。鼓泡塔内温度均匀，塔内温度与

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

21、沿塔分布，然后从模型中消去，得到的单一方程，并给出适当的边界条件。提示：可假设在气液界面上满足定律，则两相传质速率为()。解：分析 此填料塔用于气体的化学吸收，塔内物料平衡涉及两相：气相与液相。可分别对气相及液相中的待吸收组分作质量守恒，守恒方程中将涉及到的未知量包括：待吸收组分A在气相中的浓度（惰性气体）与待吸收组分A在液相中的浓度(溶剂)。守恒方程数与未知量数均为2，在给定边界条件下可以得到微分方程的特解。假设 1 设待吸收组分A在液相中的反应速率为：该反应为一级反应，反应速率常数为k（单位 3·s）2 设吸收塔在连续操作过程中处于稳态，进而在质量衡算方程中涉及到的积累项均为0解

22、答步骤在填料塔任一高度h处取一厚度为的体积微元（见习题11图），分别考虑微元中气相与液相的待吸收组分A的质量守恒：气相式中 气液相间的传质方向是由气相到液相，因而传递量属于输出项，其传质速率为：，其中是以气相摩尔分率差为总传质推动力的总传质系数（单位 (m2·s·y)），a是单位体积填料层所提供的有效传质面积（单位 m2传质面积3填料体积）。因而：式中 D是填料塔直径，为所取微元内的填料体积。由于气相中不发生待吸收组分A的分解反应，因而：吸收塔处于稳态操作，因而：根据质量守恒：即： 整理得：液相：式中气液相间的传质方向是从气相到液相，因而传弟量属于输入项，它的大小与气相输出

23、项中的气液传质量相同，为：因而：式中为组分A在液相中的应反应速率，单位：3·s。而是在所取填料微元中的持液量。根据质量守恒：即：整理得：相应边界条件：在填料塔底部 0 式中为进料气体中组分A的浓度 顶部 式中为进料液体中组分A的浓度至此，联立关于气相与液相质量守恒的两个常微分方程，并加入边界条件：求解此一阶常微分方程组即可得到气液相中组分A的浓度沿着塔高方向上的分布。下面对此常微分方程组进一步化简，消去，以得到一个只含二阶常微分方程。对式（1-11-1）变形，得：代入式（1-11-2）得到：相应边界条件：化简得：从而得到关于的一个二阶常微分方程及其边界条件：这是一个常系数二阶常微分方

24、程，将实际生产中的数据代入后可以很容易解出关于组分A在气相中的浓度沿塔高方向上的分布，进而得到组分A在液相中的浓度沿塔高方向上的分布。12.在一鼓泡容器内，初始时刻装有体积为V的氨盐水(3H2O)，随后以流量F通入2气体，在液相中即发生以下碳酸化反应 (1) (2)(3)这是制取纯碱(23)的基本反应，是一个连串反应过程。其中反应(1)为快速反应，发生在气液界面附近的液膜之内，而反应(2)、(3)为慢反应，发生在液相本体之中。设各步反应均为拟一级反应，其中3 与大大过量，其浓度可近似考虑为常数。容器中的气含率、气液比表面积与液膜厚度均为已知量，液膜体积<<1，试分别对液膜与液相本体

25、导出该反应过程的数学模型，给出分批式与连续式两种操作情况下的反应器模型。提示：液膜中的反应-扩散过程可视为拟稳态过程，2的气液界面浓度与本体浓度可设为已知，2的本体衡算方程需考虑液膜向液相本体的传递速率。解： 分析： 在此串连反应中,反应（1）为快反应，其反应速率远大于反应物2在液相中的扩散速率，因而该反应只发生在气液界面附近的液膜之内，在液相本体中，2的浓度已为零。反应（2）、（3）为慢反应，因而液相中反应物的混合速率远大于反应速率，即由液膜中通过反应（1）生成的2在还未发生显著的化学反应之前就已混合均匀，可以认为2以及其反应产物3-在整个反应器中均匀分布。而在反应器中，液膜体积<&l

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

27、膜及液相本体）中的浓度分布。但是，通过分析，可以先排除许多量：1 2：在本反应器中，2作为反应物参加快反应（1），由前面的分析可知，2在液相本体中的浓度为零，因而我们只需要求出在液膜中2浓度分布即可，根据拟稳态假设，所得到的浓度分布既可用于分批式操作，也可用于连续式操作。2 3、H2O、：在本反应器中，这三种反应物都处于大大过量，其浓度可近似考虑为常数，因而可以在反应器模型中将其作为已知量处理，并视为常数。3 2、3-：由前面的分析可知，这两种物质在液相本体中均处于均匀分布，并且它们参与的反应（2）、（3）都是在液相本体中完成，在液膜中的反应可以忽略，即液膜中2、3-的浓度分布对于反应器模型没

28、有任影响。因而我们只需求得2与3-在液相本体中的浓度随时间的分布即可（在空间上均匀分布），而二者在液膜中的浓度分布无需涉及。4 、3：在本反应器中，这两种物质均只出现在反应产物中，并没有以反应物的形式参加反应，并且在反应速率方程中也没有涉及，因而在建立的反应器模型中不涉及两者的浓度。综上所述，本反应器模型共涉及：2在液膜中的浓度分布，2与3-在液相本体中浓度随时间的分布共三个变量。其中2在液膜中的浓度分布即可用于分批式操作，又用于连续式操作。而2与3-在液相本体中浓度随时间的分布则必须分为分批式操作与连续式操作进行讨论。假设 1 在本题中，用符号A指代组分2，符号B指代2，符号C指代3-。2

29、反应（1）、（2）、（3）均为拟一级反应，因而其反应速率分别为：3 反应器体积为V，而反应器中气含率为，易知反应器中的液体量为：（1-）V。再由反应器中的气液比表面积为（单位：有效气液接触面积m2/液体体积m3），得到反应器中的液膜面积（即气液接触面积）为（1-）V，液膜体积为：（1-）V。4 假设在反应过程中反应器内的液体体积不变，即忽略反应过程中的体积变化，因此在质量衡算方程中反应器内的液体休积（1-）V为常数。解答步骤1 液膜中2的浓度分布：液膜中2的反应-扩散过程可视为拟稳态过程，因而不论是分批式还是连续式均可以作为稳态处理。由于液膜厚度远小于气泡直径，因而可以单位面积对液膜内的2（用

30、符号A代表）作质量守恒：式中1表示单位面积，为2在液膜中的传质系数，由于液膜内始终处于稳态，浓度不随时间变化，因而仅是坐标x的函数，2的质量守恒方程为一常微分方程。根据质量守恒：式中3为液膜中3的浓度，可视为常数。与之相应的边界条件为：因而可以得到关于2（A）在液膜中浓度分布的二阶常微分方程：可以解得：此式即是2在液膜中的浓度分布，而液膜内2总的反应速率为：式中1是在鼓泡反应器所有液膜中的反应（1）的总反应速率，单位。将前面得到的2浓度沿x方向上的分布代入体积积分中，很容易就可以把体积积分变为关于x的线积分，求得1，但这里我们利用稳态模型，即2在液膜内总的反应速率1应等于2通过气液相界面，自气

31、相向液相的传递量：式中S是反应器中总的气液接触面积，即液膜面积：利用沿x方向上的分布函数，易得到：式中负号表示2在液膜中的反应是作为反应物不断消耗的。2 分批式操作中2与3-在液相本体中浓度随时间的分布：在液相本体中对2作质量守恒方程：式中1是在所有液膜中反应（1）的总反应速率（单位），它也就等于在所有液膜中2的生成速率，而液膜内始终处于稳态，即所有生成的2都从液膜传递到液相本体中，因而在液相本体中2的输入项就是1。在分批式操作中鼓泡反应器内液体没产出，因而式中（1-）V为反应器中的液相体积，由于2在反应（2）中是反应物，在反应过程中不断消耗，因而在生成项中为负。根据质量守恒：整理得：消去总体

32、积，得到：这是一个常微分方程，加入相应的初始条件：得到：对此常微分方程求解，即可得到在分批式操作条件下2浓度随时间的分布。同样，在液相本体中对3-作质量守恒：根据质量守恒：得到：消去总体积，得到：这也是一个常微分方程，其中可以通前面2的质量恒算得到，、均为常数，可由实验得到。加入相应的初始条件：得到：对此常微分方程求解，即可得到在分批式操作条件下3-浓度随时间的分布。3 连续式操作中2与3-在液相本体中浓度随时间的分布：假设在连续式操作过程中反应液的体积流量为q。同样，在液相本体中对2作质量守恒方程：由于进料的氨盐水中不含2，因而输入项只包含2自液膜扩散到液相本体中的量R1。由于是连续操作，反

33、应器处于稳态：根据质量守恒：得到：化简得：在连续操作过程中，在液相本体里2的浓度不随时间变化，为一常数。同样，在液相本体中对3-作质量守恒方程：由于进料的氨盐水中不含3-：根据质量守恒：得到：化简得：在连续式操作中，在液相本体里3-的浓度也不随时间变化，为一常数。综上所述：在分批式操作中，本反应器的模型为：在连续式操作中，本反应器的模型为：化工数学第二章习题1. 求以下微分方程的解解：(1)（2）（3）2. 求解第一章给出的连续结晶器的稳态数学模型式中，成核速率B，生长速率G，流量F均可考虑为常数，加入流体的粒数分布 为l的任意函数(l)。 解： 求解以下方程：的解。分离变量得：积分得：（C为

34、常数）设，将n的表达式代入原微分方程，得：则，C1为常数。代入n表达式得原微分方程解：l0时，则原方程解为：3. 电极加热炉中石墨电极棒的传热问题可用以下方程描述式中D，U，A，T0均为常数，但导热系数k为温度的线性函数，= k0-aT，试求出上述方程的通解（建议：采用改进的p变换，可使求解更为简捷）。解：设,则有 原方程变为:, 其中为任意常数;将带入,变形得:积分得关系式为: 其中为任意常数。4. 设反应物A，B在液膜中发生以下瞬时反应v为化学计量系数，该反应受扩散限制，试导出相应的变换并求出反应物A，B与产物P的浓度分布。解：此反应受扩散限制 （1） （2） (1)-(2)得（3）令（4

35、）易由（3）可知由定义式（4）可得的边界条件为；其中为液膜厚度由以上边界条件可得（5）在处，此时（5）中，得将定义式（4）两边平方，利用瞬时反应，不能共存的条件知，于是有将上式两边开平方，得（6）根据（4），（6），可以将，用新变量表出，（7），（8）将（5）代入（7），（8）得产物的浓度分布为5O2与2在生物组织中的传递过程对于呼吸作用与光合作用具有根本的重要性。在生物组织中，溶解在液体（血液、组织液）中的O2通过渗透与扩散两种机制输送到组织内部供细胞呼吸，研究表明，在某一临界溶氧浓度c* 以上，单位体积的氧消耗速率为常数，设为q，因此，对于厚度为l的一片组织，代谢过程中氧的衡算方程为式中U

36、为液体渗透流率，D为溶氧扩散系数。设该组织外部的溶氧浓度为常数，则边界条件可表示为试求出溶氧浓度沿组织厚度方向的分布，据此判断氧浓度在何处达到最小值？渗透速率U需满足什么条件才能保证在组织内部不会出现缺氧的情况（c(x)> c*）？解：齐次方程的特征方程为解得所以齐次方程的通解为（1）利用比较系数法，求得非齐次方程的特解（2）所以，非齐次方程的通解为（3）边界条件为所以溶氧浓度沿组织厚度方向的分布为（4）对关于求一次，及二次导数式中令得所以，氧浓度在达到最小值为保证组织内部不会出现缺氧的情况，要求即需要满足6（）求以下变系数方程的级数解（a）（b）(c) （a）解：将幂级数（4.5）代入

37、方程，逐项比较系数，令首项的系数为0，得到指标方程为指标方程的两个根为，属于第一种情况，可以将c代入递推公式确定各系数。令项的系数为0，得递推公式为首先将代入递推公式，有可得 则方程的第一解为 接着将代入递推公式，整理得则方程的第二解为最后得到方程的通解为（b）指标方程，重根c12=0递推公式为,将递推公式表示成对参数c的函数形式于是含有任意参数c的幂级数y()由下式给出当0时，上式给出方程的第一解第二解即对各项求导，并令c0得所以方程的通解为（c）将幂级数代入后比较系数得到指标方程 递推公式 指标方程的根为 两根相差一个整数m时，递推公式中的系数将成为0而使之无法确定。首先确定由大根得到的级

38、数，递推公式为然后将c10代入得到第一解为注：级数推导详见微积分下册第274页考虑以下含任意常数c的级数y()上述级数在c-2处有奇异性，第二解y2由下式给出方程的通解为 7. 环形法兰上的散热问题可用以下方程描述式中k与h分别为法兰的导热系数与向周边环境的传热系数，T0为环境温度。边界条件为在内圆边界 r = r1处:T = T1在外圆边界 r = r2处: T = T2试用有关的函数给出上述问题的通解并说明如何由边界条件确定通解中的任意常数。 （提示：作变换，化为标准形式）解：8.（）用矩阵解法求以下一阶线性微分方程组的通解，并将通解用实函数表示。（1）（2）（1）解：系数矩阵A的特征方程

39、为解得A的特征值当时，特征向量方程为式中的两个方程线性相关，取x1为独立变量，令x11，得到相应的特征向量为类似的得到对应的特征向量为因此，微分方程的通解为（2）解：系数矩阵A的特征方程为解得A的特征值当时，特征向量方程为式中的方程线性相关，取x1为独立变量，令x11，得到相应的特征向量为,相应的复数形式的解为上述实部与虚部为方程的两个线性无关特解，因此方程组的通解为9.（）对于上题（1）中的系数矩阵,请用下述方法求矩阵函数（1） 待定系数法（2） 插值法（1） 解：将特征根代入方程，得到解得 矩阵函数为 （2） 将特征根代入方程（7.24）得10. 如所示,两相互联接的搅拌釜中装有体积分别为

40、V1与V2的溶液,初始时刻釜中溶质浓度分别为y10与y20 ,从t = 0开始,两釜中的溶液以流量q通过管道泵送而相互交换,管道体积可以忽略,求两釜中的溶质浓度随时间的变化关系。 图6.1 两串联搅拌釜之间的质量交换qqV1y1V2y2.习题10：两互联搅拌釜的动态响应解：由条件得，对第一个釜：对第二个釜：以上方程组成如下方程组：设，有：其特征方程为：展开得：其特征值为：，为相异实根。与对应的特征向量方程为：取为独立变量，则。令，得特征向量：与对应的特征向量方程为：取为独立变量，则。令，得特征向量：综上，微分方程组的通解为：、为常数由初始条件，得：解得：代入得到原方程组的通解为：其中11设l=

41、 li为矩阵A的特征值，试根据特征向量方程(2.4.6)证明定理(2.5.24) 。解:对于求齐次方程的通解,有式2.4.6: 定义 方程2.4.6欲有非零解,其充分必要条件是 可表示为: 由定理 有 ,其中m>n; 则任何由矩阵级数定义的函数f(A)都可用不高于1次的A的多项式表示:考虑其参数形式, 若欲确定其中的(0,11),常用的方法是用函数上的n个点上的值(1,2)插值法来得到, 因是1次的多项式,故可通过插值的方法来准确表达,即 其中 将参数代换为A上式同样成立,有此式即为定理.12.不同形状的催化剂颗粒上的反应-扩散问题可用以下方程描述式中s为颗粒的形状指数，s=0为片形，s

42、=1为长圆柱形，s=2为球形，D为内扩散系数，k为一级反应速率常数，为外表面传质系数，为流体相本体浓度。(a) 选择适当的特征尺度将问题无量纲化；(b) 分别求取s=0，1，2时的粒内浓度分布；(c) 求催化剂有效系数与模数与数之间的关系，并讨论这两个参数对h的影响趋势。解：(1)取颗粒的半径为x的特征尺度。取为的特征尺度。可将问题无量纲化。选择R与分别为自变量x与因变量的特征尺度，作变换代入方程与边界条件，整理得：(2)0 1将上式对c求导，再令0,得第二解： 均为任意常数。2(3)所以与成正比，与成反比。13. 对于可逆反应，反应的转化率由于受到化学平衡的限制而难以提高，采用反应分离耦合操

43、作的方法就可以打破这一限制。设在一催化剂颗粒内部发生上述双组分可逆反应，同时也筛选出了某种吸附剂，使得A、B的吸附性能呈现较大差别，这样，我们就可以将催化剂与吸附剂掺混，然后采用逆流移动床来实现反应分离耦合操作。 设吸附等温线为线性，反应关于固相浓度为一级，忽略颗粒内外的传质阻力与床层返混，则逆流移动床数学模型由以下方程描述式中，c 与n分别表示流体与固体相的浓度，u, us 分别为流体与固体的运动速度，求1） 解上述方程，给出浓度A、B的沿塔分布。2）如果令，则二者分别表示床层任一截面上组分A、B的净流率。研究表明，只有当与的方向相反时，才能实现A、B的分离，此时边界上的间断条件表述为塔底（

44、z = 0）：塔顶（z ）： 式中上标“”与“”分别表示塔顶与塔底边界外部的值，根据上述条件确定流体相与固体相的出口浓度。3）如果取以下参数计算相关的浓度分布并作图考察各参数变化的影响。解：1)原方程组可化为如下形式：其特征方程为：展开得：其特征值为：，为相异实根。与对应的特征向量方程为：取为独立变量，得。令，得到特征向量：与对应的特征向量方程为：取为独立变量，得。令，得到特征向量：综上，微分方程组的通解为：、为常数即：且时，知：解得：则得到原方程组的解：时，知：解得：则得到原方程组的解：即A、B的沿塔分布。2)时，因，否则静流率为0，则由条件得：时，同上原因，由条件得：3)由条件得：则代入原

45、浓度分布为：14（）发生在催化剂颗粒上的放热反应总是与传热相互耦合，从而导致多重稳态的产生与相关的稳定性问题，在铂催化剂表面的氧化反应曾被作为一个典型的体系而得到广泛的研究。现考虑在铂金属丝表面的氧化，相应的质量与热量衡算方程可表示为试从上述方程出发，在稳态解附近作线性近似，导出相关的稳定性判据。解：方程右端为零时，所得的解是系统的稳态解，设为（）,考虑系统的小扰动并令代入原方程，即得扰动（x，y）满足的瞬态方程将k展为泰勒级数，并忽略二阶及其后的高阶项，得到线性近似方程相应的特征方程为特征根实部的正负由参数决定，分为以下情况：（1） 当<0,D>0时，两特征根的实部为负值，系统是

46、渐进稳定的；（2） 当>0,或者D<0时，至少有一个特征根的实部为正值，系统是不稳定的；（3） 当0,D>0时，两特征根为相异的虚根，x与y将围绕稳态进行振荡，系统是稳定的但不是渐进稳定的。系统渐进稳定的物理变量判据为第二种解法：方程右端为零时，所得的解是系统的稳态解，设为（）,考虑系统的小扰动原方程化为将k展成泰勒级数，并忽略二阶以上项方程化为略去方程的高阶项再按稳定性的常规方法判断。15 考虑催化剂颗粒的内扩散阻力时，一般情况下颗粒热稳定性问题的分析将变得十分复杂，但是，若对内扩散过程采用拟稳态假定，则问题就可以大大简化。试从第一章给出的催化剂颗粒简化模型(1.5.15）

47、、(1.5.16）出发，对热量衡算方程采用稳态附近的线性近似，然后针对薄片型催化剂颗粒（s0）导出相关的失稳条件（斜率条件）。解：当时，有：特征方程：得：则方程的解为：代入热量衡算方程为：对系统：设系统的稳态解为，考虑小扰动代入以上方程得到扰动满足的瞬态方程：将展成泰勒级数，并忽略二阶及其后的高阶项，得到线性近似方程：其特征方程为以下代数方程的解：其中当时系统渐近稳定，得到相应的斜率条件：16. 在釜式反应器的热稳定性问题中，如果采用控制装置来对移热过程进行自动调节，就有可能把一个本来不稳定的反应器状态变成稳定状态。设控制装置对反应器增加了一个正比于温度偏差的移热速率，使得总的移热速率成为试给

48、出相应条件下釜式反应器的热稳定性判据并结合图2.6分析在何种情况下系统只存在唯一且稳定的稳态解。解：反应器的质量与热量衡算方程：设系统的稳定解为，考虑系统的小扰动：代入衡算方程，得到扰动满足的瞬态方程：式中将与展开为泰勒级数，并忽略二阶及其后的高阶项，得到线性近似方程：上述方程的特征根为以下代数方程的解：式中特征根实部的正负由参数与决定，分为以下情况：(1) 当时，两特征根的实部均为负值，系统渐近稳定。(2) 当或时，至少有一个特征根的实部为正值，系统不稳定。(3) 当时，两特征根为相异的虚根，x与y将围绕稳态进行振荡，系统稳定但不是渐近稳定。因此渐近稳定时的判据为：当直线形的移热曲线与S形的

49、放热曲线只有一个交点时，表明系统只存在唯一且稳定的稳态解。化工数学第三章习题1.（） 设函数，是一阶拟线性偏微分方程（1.1）的解，是任意常数，求证：的任意函数也是方程（1.1）的解。证明：由分别对x，y求导，得由于g1，g2满足偏微分方程（1.1），因此有将方程分别对x，y求导 （7）将（7）代入偏微分方程（1.1）左端，再利用关系式（6），得 （8）说明偏微分方程（1.1）也得到满足，因此也是方程的解。2.求下列方程的通解（1）解： 由于与存在函数关系，我们得到方程的通解为，为任意函数（2）解： 方程通解表示为隐函数，为任意函数（3）解： 由（2）得： 代入（1）得：整理得： 两边同乘以得

50、：即： 由于与存在函数关系，我们得到方程的通解为，为任意函数3 （）求下列问题的解解：该问题的特征线方程为 （1）初始曲线表示为 （2）由（1）、（2）联立解出 （3） 消去s,x，得 （4）解：特征线方程为 （1） （2）由方程（1）解出 （3）由初值（2）得代回方程（3），消去s,x，得4（）在固定化酶反应器中，生物活性细胞被制备成颗粒状填充于固定床中，反应动力学由方程给出，设流动为平推流，初、边值由给出，求浓度分布。解：依题意，该固定床反应器的数学模型为 （1）上述问题的特征线方程为 （2）初始曲线I由方程（1）解出初值影响区的解消去x，得边值影响区的解消去h，得5. 变截面色谱柱中的线

51、性吸附曲线由以下方程表述式中，流量为常数，为截面积，若定义以下无量纲变量则可将方程化为以下形式现考虑环形平面区域的色谱过程，为溶质注入端的环形边界半径，设问题为圆对称，初始浓度为0，注入的溶质浓度为常数，请在平面上绘出边值影响区的特征线，并讨论截面积的变化对波在r方向传播的影响。解： 其特征线方程为:在边值影响区，特征线为用r、t表示后为边值影响区的特征线即如上图所示。由上图可看出，随r增大，A(r)增大，增大，减小，在r方向上的波速减小。此外，A(r)还跟该环形区域的厚度有关，即式中的值，当增大时，A(r)增大，特征曲线的曲率增大，增大，减小，在r方向上的波速减小。所以，总的来说，随A(r)

52、增大，特征曲线的曲率增大，增大，减小，在r方向上的波速减小。6 （）河口附近涨潮后形成水位h的梯级分布落潮时的水位遵循以下方程请描述水位h在落潮时的衰减过程。解：该问题的特征线方程为 （1）初始曲线 （2）由（1）解得 （3）代入（2）中消去初始变量x，就得到用隐函数分段表示的解h()。7. 移动床平衡色谱模型在吸附情况下由以下方程给出求：1）对于单组分色谱问题，根据波的追赶现象讨论移动床中激波生成的条件，并给出激波间断关系。2）对于双组分吸附分离问题，给出稳态条件下的解答。设进料流股中含有A、B两种组分，试确定合适的操作速度ur以便将A、B分离。解： 1）对于单组分色谱问题，模型简化为其中，令 得模型的特征曲线方程为由前两式可得波速的方程当 时，波速为正方向，此时