1、 结晶器的选择与设计 在现代结晶器的设计中 考虑结晶动力学 设计产品的尺寸,尺寸 产品的质量。结晶器工艺设计 1.确定操作模型 (1).连续型 (2).间歇型 2.产生过饱和度的方法 冷却、蒸发、反应、其它 3。结晶过程质量衡算 质量、热量、粒数衡算(略)颗粒衡算与质量、热量衡算的关系 3.设计流程操作模型的选择 连续:大型生产一般产生的晶体的尺寸分布较宽可用一些控制手段来完成大晶体窄分布的产品。间歇:相对小的生产能力(50T/天)一般可产生分布较窄的晶体也需要适当的控制方法。操作过程的选择 冷却:溶解度随温度变化较大的体系,适宜用冷却的方法。蒸发:如果溶解度随温度变化不大的体系可用蒸发法。在
2、溶解度高时,为提高回收率用蒸发的方法。反应:能产生沉淀的反应体系。其它:溶剂萃取,高压结晶等。结晶器类型的选择 DTB结晶器(draft tube and baffle)图5.7(导流筒、挡板式)特征:固体悬浮较好,可带有结晶排除系统和产品颗粒分级系统,实现晶体粒度的控制。FC(forced circulation)强制循环 图5.11 5.17 5.18 强制循环结晶器可用于很多过程,可根据不同的过程来调节,其也可设置粒度控制系统。结晶器的计算 生产能力和产品质量,和粒度要求 计算要求的基础数据与模型 溶解度数据 溶液的基本性质和颗粒的基本性质 如密度、形状系数等 结晶成长速率模型与成核动力
3、学模型 溶液的初始浓度C0 计算:MT悬浮密度 C1溶液的终止浓度,根据最后排出体系的温度而定,即C1=f(Tf)VW蒸发水量 V0溶液的进料体积流率,此参数对多组分系统,可根据母液中杂质浓度而定,对单组分系统MT是一个操作参数,根据结晶器的流动状态而定。因此,可用上式计算体积、流率。当结晶器操作于MSMPR模型 结晶器体积:进料流率*停留时间 这样设计的结晶器比较简单,但因为是MSMPR的假设下,而以其结果很难适用于大型结晶器。11CTCM0VWV03GLd3GLd 假设在稳定连续操作的结晶过程中,其产品的颗粒的累积尺寸分布可用Rosion-Rammler 方程表示:*mm(L/L)xree
4、-=设计图及设计程序设计图及设计程序r相对产品颗粒个数的累积分布,当所有的晶体考虑在内 r=1 m改良的均一系数,不同的m值表示不同的分布函数,假设在一个过程中的产品其m值不变,也不随颗粒的尺寸而变。L晶体的尺寸。L*参考尺寸,其定义为在此尺寸下,r值为0.3679与m无关。x无因次晶体尺寸,x=L/L*在用方程(1)表示积累产品晶体分布时,颗粒的粒数密度(即在一定的颗粒尺寸范围内的晶体个数)为 (2)mm 1xdrdxmxe-=-因此,在此尺寸范围内的晶体质量为 (3)产品晶体的质量分布 kg/mh 晶体密度 kg/m3 晶体的体积形状系数 以个数为基准的晶体生产速率,它等于晶体的成核速率#
5、/h3*3drCvdxw(x)Fk x L()=rm*m2xCvFk Lmxe+-=rw(x)CrvkF 如果以产品中最大质量的晶体作为晶体尺寸来表示,这种晶体尺寸可以从方程(3)得出即 (4)注意:当m=1时,xd=3,其分布式和MSMPR的结果一致。当m=0时,产品晶体为同一尺寸,表示一个连续操作的带有理想产品分级的操作过程,对方程(3)进行积分,得到以质量为基础的生产速率。(5)m*3mm 1xdwCvdxFk L(m2mx)xe+-=r+-0=1m2dmx(1)=+d*3LdLx=3*3drvvdxLP=Fk L x()dxrm*3m 2-xv v0=F k L mxe dm+rF3#
6、/m hgFBV=mv*3m2-xPvV0L(Bk mxedx)+r=当 m=1mm 2x3x00mxedxx e dx6是全部结晶器内成核速率,如果用V表示结晶器体积,成核速率表示为B,即(6)1md2dmLx(1)3L*=+=33dLL*27方程(6)变成v*33PvvdV26Bk LBk L9r=这一结果和MSMPR的结果相同以上分析是从产品晶体部分来考虑,下面再从结晶器内部的晶体变化来考虑。在稳定操作状态下,细小的晶体或者从外部加入(加晶种)或者来源于系统内的成核,假设在操作过程中成长的晶体没有破碎或聚并。因此,在结晶器内一定尺寸的晶体的晶数密度与产品中大于此晶体尺寸的积累尺寸分布成正
7、比,也就是说,在这种假设下,在结晶器内的晶体密度分布正比于方程(1)的r。如果把晶种(成核)在dx范围内的个数表示为n0dx(注意:n0为核的粒数密度 n0,在无因次晶体尺寸x的颗粒数为 m-x00(n dx)r=(n dx)e4#/mhg 因此,在结晶器内全部晶体重量为 从另一方面讲,在结晶器内晶体的全部个数也可以表示为 (8)晶体的平均停留时间 从晶体的平均停留时间的定义 平均停留时间 一定尺寸的晶体个数 一定尺寸晶体的停留时间m3-xtvv00wk L n edx=rm*33-xvv00k L nx edx=rm-x00N=F=n edxQdrdrdxdx00dxdxdxdr(7)等号两
8、边同时等于结晶器内晶体的总个数(相对值,因为r为相对分布),从晶体生长的角度我们有 或 (9)由方程(8)和(9)粒数密度 (10)LGavddL)(Gav)d/dL(Ldxemxdxmavx1m0)d/dL(L0dxdrdxexmmav*x0m)d/dL(LdxexmFdxenNmav*mx0m)d/dL(Lx000mxav0mxm*dxe)d/dL(dxexmLF0n把方程(10)代入晶体总质量方程(7)由生产速率方程(方程(5)和方程(11)我们有 0mxav00mx3mxm4*Vcdxe)d/dL(dxexdxexmLkFtw0mx2m0mxav0mx30mxm4*dxexdxe)d/
9、dL(dxexdxexmPLtw(11)(12)从结晶器内的悬浮密度考虑,如果我们用表示溶液的含量分数,结晶器内晶体的质量可表示为:最后得到mm*mx3x00mmxm 2xC00PLx edxx edxdLavdVedxxedx(1)()tw(1)V(13)(14)分析方程(14)此方程给出生产强度(P/CV)悬浮密度(1-),晶体尺寸L(或X),成长速率,成核速率 方程(6),从产品角度分析而得到的粒数,质量衡算。方程(14)是从结晶器内的结晶过程而得到的粒数,质量衡算方程,这两个方程必须建立一定的联系,从此而得到一个结晶器的设计图。设 (15)(16)ddLdxexmImx2m100mx2
10、mmx00mx3mxmdxexdxedxexdxex2Imv*3m2-xPvV0L(Bk mxedx)+r=(6)用Ld代替L*方程(6)可表示为 方程(14)可表示为 3v3dVPIBkLCCdLPd 4VdavL I1m3)1(IIm213m1)1(IIm224(17)(18)作图把方程(17)取对数在一定的Ld下,与 为一条直线,斜率为1 这条直线取决于要求的产品粒度Ld和相应m值为特征。同一直线下,表示在一定的m值下的同一粒径,这条直线方程做特征线,就是说,这些平行线表示不同的粒径。另一方面,在一定的生产能力下,由方程(17),Ld和Bkv立方根成正比,因此,Bkv的轴可以和Ld轴平行
11、而置。cPv 3dvloglog(Bk I)3logLcPvlogv 3log(Bk I)作一个表示Ld,P/CV,Bkv的关系图,m=1作出特征线并确定Ld在Ld轴上的值,当m=1时,I3可被计算为0.225 当P/CV=0.01,Bkv=1108,Ld=7.63110-4 P/CV=0.01,Bkv=1109,Ld=3.54210-4 图8 用纵坐标表示P/CV(对数值)用一横坐标表示Bkv(对数值)另一横坐标表示Ld轴(对数值)在Bkv轴上找到1108和P/CV=0.01 作一直线通过此点,但斜率等于1标此直线与Ld轴交点为7.63110-4与上相同做出,P/CV=0.01,Bkv=11
12、09不同的m值做图。按道理讲,不同的m值,其P/CV、Bkv和Ld应遵循不同的规律,可以完全做成另一个图,然而,要是用一张图表示即存在的的信息总量较为简单。在特征方程中,如果我们把P/CV考虑成一个常数,Bkv=k(m)1/Ld3 这个方程表明,当P/CV不变时,Bkv与Ld3成反比,如不同的m值,它们的关系只改变其系数但比例关系不变,也就是说这条直线的位置有些变化,但变化趋势是一样的。这样我们就可以用一个Ld与Bkv之间的关系。3dL1V)m(kBk13d23d13d1323d232V1V)(LL1LILIBk(Bk(m)k)例如:m=1,P/CV=0.01,Bkv=1108,Ld=7.63
13、110-4 m=2,P/CV=0.01,Bkv=1108,Ld=6.01310-4 在例1图上作出m=2的Bkv轴和Ld轴 做图:选画一条辅助线 画一条横线为当m=2时的Bkv线从m=1的Bkv轴上找到Bkv=1108 通过此点作垂直线交m=2于H点作平行于辅助线的直线通过H点交P/CV=0.01的横线于G 通过G作斜率等于1的特征线,这条特征线的Ld必须等于6.01310-4从m=1的轴上找到6.01310-4,做垂直线交于此条直线位F 通过F点作出直线即为m=2得Ld轴这里强调一下,在这样的图中,如果m值不同BkV和Ld得值是不同的。到此,我们在一张图上把Ld,BkV,m,P/V的关系表示
14、出来。结晶器内的变化规律图 由方程 从产品方程可得 代入 如果假设 为常数 则 即从结晶器内的粒数与质量衡算avddL4dVP)(1(ILC313vC)(LIBk1VPdavddLIIBk1VP)(1()()(331431v34CavddLI1dVP)(1(L4C)m(C)()(31v34CBk1VPCP1VV4log()logC(m)logBk 与 为直线斜率斜率为 ,这条线与表示 与G乘积相同,所以做坐标()轴和G轴时,这个方程也一定要满足,要等于方程左面的乘积。即()或 可以任选一个值而另一个值必须满足上面的方程,把这些关系标在同一张图上,成为图7)(lgVPCVlgBk4111avdd
15、L4VPL)(1(ICd1avddL)(成核速率对晶体尺寸的影响分析 在生产强度一定的条件下,即,P/CV为常数 控制结晶器中的成核速率 成核速率增加,操作点向右移动,晶体尺寸减小。反之增加。另一方面,从晶体的成长与质量横算考虑,在图中有要满足操作线,即悬浮密度与成长速率的关系。生产速率增加 生产速率的增加,在设计图上操作点上移 产品的晶体尺寸可能增加,已可能降低。要看其造作是怎样提高生产能力的 在原操作点左方,晶体尺寸增加。在原操作点右方,晶体尺寸减小。保持产品粒度,而增加(或降低)生产能力 保持产品粒度,在设计图上是早同一条特征线,但其操作线 在变化。同一条特征线,在不同的生产能力下,要求
16、有不同的操作线。其悬浮密度必须保持恒定,以确保产品的粒度不变,因此,成长速率必须成比例增加。设计图的应用-用试验数据估计结晶动力学数据 试验可测得 粒度分布,生产速率,结晶器有效体积,悬浮密度 根据粒度分布,可计算粒数分布 获得平均尺寸Lw(即积累尺寸分布为36.8%的晶体尺寸)和在参数n,n为(Rosin-Rammlardiagram的斜率)由n可以得到m,Lw可转化为Ld 由测得的晶体尺寸分布可得到Ld,和m 从图中根据悬浮密度Ld,m和生产能力确定成长速率和成核速率(如果实验中通过改变不同的参数和测量起过饱和度,可得成核速率模型)m与n de 关系,和LD与LW,m 的关系 改善晶体尺寸
17、 已知:m,Lw,或m,Ld(产品要求)生产能力p,生产速率 计算:悬浮密度(质量衡算)设计方法:由悬浮密度和成长速率确定操作线 由Ld确定特征线 由两条线的交点即为系统点 有特征线上读出P/V得到V同时得到BKV值,即如果要达到产品的粒度要求成核速率必须控制,在此成核速率不能达到时,系统操作必须调节悬浮密度或成长速率。结晶器体积的设计结晶过程的模拟(基于平仓贤理论)缺欠 动力学数据(成核速率在生产上的要求与其规律没有联系)图的准确程度 用完全计算的方法实现结晶器的设计,考虑到体系的成核与成长速率 设计方程可以改写为 可用试差求解此两个方程 d 4TPL IV=GM3vd3PB=KVL I计算方法已知:p,c,m,Ld,MT(其中可给范围值)给一个成长速率G 计算V从方程a 计算成核速率从方程b 用成核速率模型计算成核速率 比较如不匹配,重设成长速率进行系统模拟。