优秀学习资料欢迎下载化工数学第三章习题1 求下列函数的laplace 变换（1）2sin()ntT； （2）cos()t； （3）2sin t； （4）sin2tata； （ 5）11 exp()tt解：（1）22222sin()1 2111()22222()nnititTTnLtTL eL einnisisiTTnTnsT（2）22cos()cos( ()sLtLtses（3）22sin1 (1)cos2 21 1()24LtLLtsss（4）22222sin212()tLatadaa ds sassa（5）01 1exp()1 ()1ln01lnttLttLeetssss2 利用 Laplace 变换的微分性质求下列函数的Laplace 变换（1）3( )sin 2tf ttet； （ 2）30( )sin2tf tted精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 16 页优秀学习资料欢迎下载解： 、（1）322222sin2 2(3)422(3)(3)44(3)(3)4tL tetddssssss（2）30222222222sin212(3)41214(3)()(3)4(3)424(3)(3)4(3)4tL teddds sssssssssss s3 求下列函数的Laplace 逆变换（1）41s； （2）2125s； （3）4(1)(2)ss； （4）3(1)(3)sss解：（1）144 10311(41)!16tLstet（2）1211sin5255Lts（3）11244444(1)(2)12ttLLeessss（4）1131331322(1)(3)1322ttsLLeessss4证明二阶常系数线性常微分方程当其特征根为重根12ss时，方程的基本解为1stE tte当特征根为两共轭复根1,2si时，方程的基本解为1( )sintE tet精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 16 页优秀学习资料欢迎下载给出以上两种情况下非齐次方程的特解的一般表达式证；依题意，是求取以下问题的解答( )(0)(0)0aybycytyy对上式进行拉氏变换，仍记( ) ( )G sL y t，得到21211( )()()G sasbscssss当12ss时，211( )()G sss111211( )( )1s ts tE tLG stete特解的一般表达式为1()0( )( )()( )tstyE tf ttefd当特征根为两共轭复根1,2si时，2211( )()()()G ssisis11( )( )sintE tLG set特解的一般表达式为()01( )( )sin()( )ttyE tf tetfd5 求以下控制回路的传递函数F 1( )G s2( )Gs1( )G s+ Y 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 16 页优秀学习资料欢迎下载解：1321223()()()()()()1()()FGsYGsGsYGs GsYGsFGsGs6、在气固催化动力学实验研究中，催化剂表面的吸附与化学反应是两种较难区分的因素，因为这两种效应都发生在固体表面，一般的稳态实验难以鉴别，而采用脉冲动态实验就有可能区分和测定这两种参数，这对于了解催化过程机理和指导催化剂开发具有重要意义。设固体表面的化学反应为一级反应0,0Arrk ntn并设表面吸附为传质过程的控制步骤，忽略内、外扩散阻力，则颗粒模型由下式给出arnnkck ntK在不考虑固定床轴向扩散的情况下，试给出床层出口气体浓度的各阶时间矩与吸附、反应参数之间的关系。解：该气固催化模型可视为一个有反应发生的色谱过程，不考虑轴向扩散，则DZ=0，由于忽略内、外扩散阻力，表面吸附为传质过程的控制步骤，故颗粒上溶质的平均浓度=n(固相浓度 )。床层模型：10(0, )( )( ,0)0rccnk ctxcttc x(1) 颗粒模型：0,0arnnkck ntKtn(2) ka为吸附速率常数， n/K 为贴近颗粒表面处的流体相浓度，与固相浓度n 处于吸附平衡状态，令1对颗粒模型（ 2）进行 Laplace 变换，并仍记( )( ),( )( )L n tn s L c tc s得：( )( )( )( )arn ssn skc sk n sK(3) 则，( )( )()/arakn sc ssKkkK由 5.5.7 对颗粒传递函数M(s) 的定义得：精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 16 页优秀学习资料欢迎下载( )( )( )( )/1parammrcKkn sKM sc sc sKsKkkt stt式中,mrtt具有时间量纲，maKtk代表吸附过程的时间尺度；1rrtk代表反应过程的时间尺度。对床层方程（ 1）进行 Laplace 变换，得：0(0, )1rcs cnk cxcs(4) 将 M(s) 代入（ 4）中消去n(s)得：1( )0(0, )1rcsM sk cxcs(1( )( )( , )rsMskG sc l se，l出口气体浓度的各阶时间矩为：/01022002(0)(0)1/12(0)1/1(/1)rrktmrmmrmrGeeKGttKtKGtttt7、当采用上述脉冲实验方法来研究非均相催化反应时，由催化剂装填而成的反应柱往往难以将反应物 -产物的色谱峰完全分离。为解决这一问题，一般都在反应柱之后串连一色谱柱，将反应物 -产物分离后由检测器检测，如图所示。 这种实验技术称之为催化反应色谱方法（陈诵英，高荫本，彭少逸，石油学报（石油加工），1988,4(1):2937 ） 。由于检测器测得的信号是经过多个单元传递的结果，且包含各种物理因素的影响，因此在数据处理时必须注意剔除与反应无关的其它因素以得到真实的反映表面吸附与化学反应作用的信息。试根据串连过程的时间参数加和性质与方差加和性质说明如何设计相关实验获得吸附和反应动力学参数。解：由第 6 题结果可得反应柱传递函数：11(1( )1( )rsMskG se11111( )/1mmrKMststt111maKtk代表吸附过程的时间尺度；1rrtk代表反应过程的时间尺度。反应柱的各种时间矩为：精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页，共 16 页优秀学习资料欢迎下载11/0,11111,110,1112111112,110,110,11211(0)(0)1/12(0)1/1(/1)rrktmrmmrmrGeeKGttK tKGtttt忽略色谱过程的轴向扩散及内、外扩散阻力：床层模型：0(0, )( )( ,0)0ccntxcttc x颗粒模型：0,0annkctKtn对颗粒模型进行Laplace 变换得：( )( )( )( )( )/1aaamkn sn sKsn skc sKc sskKt s2222222( )( )( )1mnsKMscsts对床层模型进行Laplace 变换得：221( )20( )(0, )1Mscs cnGsexcs色谱柱的各阶时间矩为：0,221,222222222,222222222(0)1(0)1(0)1mGGKGKK t串联系统的总传递函数12( )( )( )G sGs Gs，令 s=0，得：0,0,10,2T根据时间和方差的加和性得：11121222111/1TmrKtttKtt精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 16 页优秀学习资料欢迎下载2222111111121222211121/1/1/1mTmmrmrmrtKKK ttttttt121212,KK均为可测量得到的常数，通过实验测定出1,mrtt值则可获得吸附和反应动力学的参数,arkk实验过程设计：首先以不与催化剂单独作用发生发应的反应物之一为研究对象进行脉冲实验， 获得各阶矩的数据后求出,1Tt和2,1T， 然后以能发生反应的混合物为对象进行脉冲实验，获得另一组数据,2Tt和2,2T。,111122211,2122211111/1TTmrtKKKtKtt1,1,211111/1mrTTmrmrttttKtttt的值22111111,1,2111111111112211/1/1/1mTTmmrmrmrtKKKtKtttttt从而得1mt和rt，代入111maKtk和1rrtk中得吸附速率常数ak和反应速率常数rk8停留时间分布或其它物理量分布的时间矩值除了可由脉冲实验来测定之外，也可由阶跃输入的方法来测定，特别在气体中溶质含量低、固体吸附量大的情况下，采用连续输入的阶跃实验技术上更容易实现。为了处理阶跃实验数据，理论上要求建立阶跃实验的动态输出曲线 yU (t)与脉冲输入时的各阶矩值k之间的关系，以便通过阶跃实验来测定脉冲输入的矩值。试按照以下步骤推导阶跃输入时的有关求矩公式：1）根据基本解的定义和性质证明，对于单位阶跃输入函数U(t)，相应的输出函数yU (t)与基本解 E(t)之间满足以下关系2）根据脉冲输入的各阶矩的定义将 E(t) 用 yU (t) 代换，证明式中 ys= yU()为阶跃实验的输出达到稳态以后的值。解：1) 单位阶跃函数dttEtkk)(0)()(tEdttdyU0sydttyytkUskk)(10精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页，共 16 页优秀学习资料欢迎下载00( )10tU tt用基本解构造单位阶跃输入函数的U(t) 的输出函数，可得：0( )()( )tuytE tUd设tx000( )( ) ()( )( )tuttytE x d txE x dxE x dx则可得( )( )udytE tdt2) 000( )( )()(0)uuuE t dtdytyy有齐次边界条件可知(0)0uy所以0()usyy0000( )|kkkkkuuut E t dtt dyt yy dt00()kkuydty dt100 ()( )( )kkusuyy t dtk tyyt dt9. 鼓泡塔是最简单的气体吸收器和气液反应器，动力学研究用的实验室鼓泡反应器一般都采用分批式操作，如图所示， 气体从底部通入装有液体的容器，经分布板之后形成分散的气泡并在液体中浮升，最终从容器的上部输出。鼓泡塔中气相的流动可考虑为平推流，液相则一般考虑为全混流。分别对气、 液两相组分建立物料衡算方程，并略去气相组分的时间导数项（该项较小，可以忽略），得到式中 cg和 cL分别为反应组分在气相和液相中的浓度，g为气含率， 即气泡所占的体积分率，U为空塔气速， kLa 为气液传质系数，kr为一级反应动力学常数，ci为气液界面处的液相浓0()()(1)gLiLlLLiLrLgcUk a cczdck accdzk cdtl精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页，共 16 页优秀学习资料欢迎下载度，由 Henry 定律给出当气相反应物采用脉冲输入时，初始与边界条件如下试采用矩量分析方法求出反应器顶部输出的气相浓度的0、1、2 阶时间矩， 并说明如何设计实验测定参数kLa, H, kr。解：对gc和Lc做关于时间t 的 Laplace 变换，设：( , )( , )ggL ct zcs z( )( )LLL c tcs( , ) /igccs zH边界条件：0:1gzc原方程可化为：()giLLdcUk a ccdz(1) 0()(1)lLLiLLrgk ascccdzk cl(2) (1)式是关于gc的一阶常微分方程，分离变量求解：(/)ggLLdcUk a cHcdzgLgLk adcdzHUcH c( )ln1gLLLk alclH cHUH cHccgi/0:0:0Lggccttcz精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页，共 16 页优秀学习资料欢迎下载设：( )1Lk algLHULclHceAH c得：( )(1)gLclAA H c(3) 将(1)代入 (2)得：( )0(0)()(1)(1)ggcllgLLgrgLgck aU dcUskcdzdclk a dzl()(1( )(1)LgrgUsk ccll设：(1)1( )()()ggrrLlclskB skUc(4) 由(3)，(4)两式可得：()(1)( )()(1)rgrAB skA HclB skA H由矩的定义可知：反应器顶部的的气相浓度的0 阶时间矩为：0(1)(1)(1)(0, )(1)1(1)(1)LLLk alk alHUHUgrrgk alrHUgrHUeelkABkA HclBkA HHUelk222122(1)(1)(1)(0, )(1)1(1)(1)LLk alHUgrgk alrHUgrHUelkABHclBkA HHUelk精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页，共 16 页优秀学习资料欢迎下载23222332(1)(1)2(1)(0, )(1)1(1)(1)LLk alHUgrgk alrHUgrHUelkB HAclBkA HHUelk由矩的物理意义可知，零阶矩表示区间内所有物理量的总和，一阶矩与零阶矩的比值，10t表示了分布的平均值，2220t或2200211表示了分布函数的方差。为了求出 0， 1，2 阶矩，需要得到反应器顶部输出的气相浓度对时间的变化曲线。在反应器入口打入原料气的一个脉冲输入，同时在反应器出口记录尾气的浓度，可以得到气相浓度对时间的变化曲线，利用微积分和概率统计的知识，可以求出浓度曲线下的面积，浓度的平均值以及方差，由定义得到各阶矩。然后代入0，1，2 阶矩的表达式中，就可以得到含有rk，H以及Lk a的三个方程。通过求解这个方程组可以得到这三个参数。10附图为一两级连续结晶器流程，由两个体积相同的结晶器串连组成。设第一个结晶器输入的料液中不含晶体，结晶器处于稳态操作，液体流量相同， 每个结晶器的成核与晶体生长速率已知（但不一定相同），试求最后一级结晶器输出的晶体重均粒径。nin= 0 解：由教材 P294 页连续结晶过程的矩量分析可知：第一个结晶器的各阶矩为：0,11B1,110,1G2,111,12G3,112,13G4,113,14G其中 G1和 B1分别是第一个结晶器的晶体的成长和成核速率。第二个结晶器的各阶矩为：0,20,12B1,21,120,2G精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页，共 16 页优秀学习资料欢迎下载2,22,121,22G3,23,122,23G4,24,123,24G其中 G2和 B2分别是第一个结晶器的晶体的成长和成核速率。所以第二个级结晶器输出的晶体重均粒径为：4,24,123,23,23,122,243WGlG将上面两个结晶器各阶矩代入，化简得：5432234411112112112122243223311112112122224()6()WBGB G GBG GBG GBGB GlBGBG GBG GBGB G4322344111121121121222322331111211212224 ()()BGBG GBG GBG GBGB GB GBG GBG GBGB G11已知聚合物的链长分布可由方程（5.6.50） 、 （5.6.51）描述，即1jjjdPPPdt（j1, 2, ）100,002jPIPj当 j1 时，可将j 考虑为连续变量并取如下近似21212jjjPPPPjj于是链长分布P(, j)由以下连续方程描述220012,0,00,0jPPPjjPPljPj（11-1）试对上述方程作关于变量j 的 Laplace 变换，然后求出链长分布的04 阶矩（注：上述近似仅在j1 时才能得到合理的结果）解：对式（ 11-1）做关于变量j 的 Laplace 变换，得20012(0)0dPsPIs PsIdP（11-2）精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 12 页，共 16 页优秀学习资料欢迎下载可得式（ 11-2）的解为01exp12IPAsss又(0)0P001exp12IIPssss（11-3）000(0)lim( )sPP sI式（ 11-3）可变换成如下形式：001exp12sPIIss（11-4）对式（ 11-4）求关于s 的二阶导数，得01exp112PsPIsss2200112exp1exp1122PsPIssIsss当 s 0时，有2101(0)2PI（11-5）同理有332300113 3exp11exp1122PsPIsssIsss32201(0)33PI（11-6）23423004401143exp16exp11221exp112PsPIssIsssIsss3432301(0)634PI（11-7）45303545001515exp1121110exp11exp1122PsPIsssIsssIsss精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 13 页，共 16 页优秀学习资料欢迎下载4543401(0)10155PI（11-8）12设无限长色谱反应柱中的溶质运动由以下方程描述221,0rcccKDk ctxxc xx（12-1）试按照本章 7 节所介绍的方法由空间矩和Hermite 多项式构造本问题的解答，讨论当吸附参数 K 变化时，色谱峰的移动速度、峰面积、峰宽呈何种趋势变化；如果输入的同一个脉冲中含有两个不同的独立组分，二组分具有相同的参数, D, kr, 但有不同的吸附平衡常数Ki（i1，2）试画一草图表示两组分色谱峰的运动图象，要求标出两峰在同一时刻的空间位置、峰面积、峰宽上的差异。解：对式（ 12-1）做关于x的 Fourier 变换，仍记F c(x, t) = c(, t)，21,01rdcKiciDck cdtc（12-2）可得式（ 12-2）的解为2,exp1rDiktctK可得其各阶矩值如下00,exp1rk tctK（12-3）01 0,1tictK（12-4）0122 0,11DttctKK（12-5）3321340,11tDtictKK（12-6）432460,11tDtctKK（12-7）一般有21210,11kkkkkDttctKK（12-8）该分布函数的形心位置x和第 04 阶中心矩 分别为的精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 14 页，共 16 页优秀学习资料欢迎下载101txK01（12-9）10（ 12-10）222021DtxK（12-11）332300320xx（12-12）2243424000463121DtxxxK（12-13）我们定义如下空间变量2221xxxxXDtK则浓度分布函数c（ x，t）关于 X 的 Hermite 多项式级数式（12-14）所示20,expnnnc x tXbt HX（12-14）其中222,2!2expnnnnnc x t HX dXc x t HXdxbtnHXXdX则其前五项系数分布为22001,22bc x t dx（12-15）2220111,202 222xxbc x tdx（ 12-16）222222020411,422024 2 228 2xxbc x tdx（12-17）精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 15 页，共 16 页优秀学习资料欢迎下载222223330121,8128 3! 2221,03!48xxxxbc x tdxc x txxdx（12-18）2222242404021,16481216 4!22241304!422xxxxbc x tdx（12-19）进一步计算可证明，更高阶的系数bn（n4）也恒为0 222002,expexp2exp1exp4211rc x tbXXk txxKDtDtKK（12-20）色谱峰的移动速度1dxvdtK（12-21）峰宽（标准偏差）102022211DKt（12-22）峰面积2exp1,exp4211exp1rrk txxKSc x t dxdxDtDtKKk tK（12-23）由式（ 12-21）至式（ 12-23）可知某一时刻色谱峰的移动速度v 随 K 的增大而变小，而峰面积 S和峰宽 0都随 K 的增大而增大。如果输入的同一脉冲中含有两个不同平衡常数Ki（i1， 2） 的独立组分， 不妨假设 K1K2，则某一时刻两峰的空间位置、峰面积、峰宽如图所示：精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 16 页，共 16 页