制备载体催化剂的浸渍过程数学模型及有限元正交配置解.pdfVIP

制备载体催化剂的浸渍过程数学模型 及有限元正交配置解 张柏青*赵修仁 大（连理工大学化工学院） 浸渍型催化剂在近代工业催化剂开发中的地位日益重要。根据反应的特点，可 以对 催化剂活性组分的分布提出特定的要求。目前在工艺上已有一些方法可制备活性组分符 合不同分布规律的催化剂，如所谓的蛋壳型，蛋白型及蛋黄型。 Maatman等[1]首先提出了设计和控制浸渍参数以制备具有预期活性组分分布的催 化剂的观点。 Weisz等[2]首先建立并求解了浸渍过程的数学模型。以后Lee等3[]又对浸渍过程各 个步骤的数学模型及其解法作了系统论述。但这些工作大多都是针对低含量贵金属盐催 化剂如（Pt,Pd等）而做的，对普通金属盐催化剂，只有Komiyama4[]研究过镍盐浸渍过程 的数学模型。 鉴于此过程的复杂性和求取模型参数的困难，本文仅研究了硝酸镍在球形A-O23l 载 体上的浸渍过程。 浸渍过程模型的数值解法，文献中多采用差分法，本文采用有限元正交配置法。 采用Lee等3[]的忽略扩散项的单组份浸渍模型来模拟硝酸镍在球形 -Al2O3上的浸 渍过程，经无因次化处理后此模型为： 初值： 边值： 式中=C/C0,=n/ns,=t/t,=1－r（/R),=ns/C0,K1=k1C0t,K2=k2t,Q()=V (t)t/4 R2t=8R2/kdP。其中u：浸渍组分无因次浓度，C：浸渍液浓度m（ol/m3或M ， θ龋吸附组分无因次表面浓度，n和ns分别为：表面浓度和表面饱和浓度(mol/m)2,：无 因次 时间，t和t分别鹞为时间和液体渗透的特征时间(s) ：径向无因次坐标，r：径向坐标(m)， R：球半径(m)，：相对吸附能力，K1和K2分别为无因次吸附和脱附的速率常数，k1和2 分别为吸附速率常数(m3/mol·s)和脱附速率常数(S)， 液禾迳傅奈抟虼翁寤俾剩 1987年10月13日收到。 现在中国化工进出口公司大连分公司。 V液体渗透的体积速率(m3/s)，：载体空隙率，:液体粘度 厘（泊），kd：载体的渗透度， P：大气压力。 如令 则(1)，(2)式可写成： 式中 可以证明(7)式的特征是一组垂直于 轴的直线。而由(6)式的特征可求得： 式中下标f为液体前沿，a为积分常数，0 a 1。其含义是不同时刻进入粒子的液体前沿 所留下的轨迹，如a=0表示首先进入粒子的液体前沿的轨迹。事实上，a就是在它所代表 的液体前沿以前进入粒子的无因次体积 （液体体积／粒子孔体积）。 利用链式法则及上述两个特征的性质可将(,)坐标转换 成(,a)坐标。从而将 (6)，( )式转换成： 初值： 边值 ： 这一组微分方程的求解前人都采用一阶差分法。本文采用有限元正交配置法[5]离散 (1)式，离散过程中采用由Lagrange多项式求得的非对称内配置点及其相应的离散矩 阵，将(1)式降阶为代数方程后用Gauss主元素消去法求解离散点上的u值。再将解得的 u值作为已知值用Gear法求解(2)式。 在计算时，根据所求曲线梯度的不同调节元的密度，但每个元内的内配置点数均取为 2，而每个元内配置点的位置均为0;0.2113;0.7888;1.0。 (1)式离散后的方程式为： 式中：IE一元的序号，hA(IE)－元IE的长度，N－内配置点数，IP－内配置点序号，m= 3(IE－1)+I离散点的序号。 相邻两元之间存在下列关系： 式中 由(14)和(16)式可见，对每一个元来说都可以得到(N+1)个独立的代数方程组。对 第一个元来说，u1=1。对其它元来说，每一个元的始点为前一个元的终点，其序号K相 同。有关此解法的细节及程序见另文[6。 结果与讨论如下： 1．为验证所用有限元正交配置方法的可靠性，对0.05MH2PtCl6溶液在A-O23l 上 的单组份浸渍曲线进行了计算，并将结果与Lee等用差分方法所得的结果一起列于图1 中。可见两法所得结果相近。但有限元正交配置法所需机时明显少于差分法。 2．计算了不同浓度硝酸镍溶液对A-O2l 3球进行单组分浸渍后，活性组分在粒内 的浓度分布