CN107589667A - 微界面强化反应器能量耗散构效调控模型建模方法 - Google Patents
微界面强化反应器能量耗散构效调控模型建模方法 Download PDFInfo
- Publication number
- CN107589667A CN107589667A CN201710766690.5A CN201710766690A CN107589667A CN 107589667 A CN107589667 A CN 107589667A CN 201710766690 A CN201710766690 A CN 201710766690A CN 107589667 A CN107589667 A CN 107589667A
- Authority
- CN
- China
- Prior art keywords
- gas
- liquid
- reactor
- area
- mix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000021715 photosynthesis, light harvesting Effects 0.000 title claims abstract description 18
- 239000007788 liquid Substances 0.000 claims abstract description 156
- 230000000694 effects Effects 0.000 claims abstract description 19
- 230000005587 bubbling Effects 0.000 claims abstract description 14
- 238000013461 design Methods 0.000 claims abstract description 7
- 238000005192 partition Methods 0.000 claims abstract description 4
- 239000000203 mixture Substances 0.000 claims description 88
- 239000012530 fluid Substances 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 14
- 230000003068 static effect Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000009423 ventilation Methods 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 4
- 230000035939 shock Effects 0.000 claims description 3
- 239000007921 spray Substances 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims 1
- 230000001186 cumulative effect Effects 0.000 claims 1
- 238000009795 derivation Methods 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 45
- 238000006243 chemical reaction Methods 0.000 description 34
- 238000012546 transfer Methods 0.000 description 18
- 239000012071 phase Substances 0.000 description 9
- 238000011160 research Methods 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 8
- 239000000463 material Substances 0.000 description 8
- 239000002245 particle Substances 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 7
- 238000005265 energy consumption Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000004581 coalescence Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000001965 increasing effect Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000007791 liquid phase Substances 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 241000790917 Dioxys <bee> Species 0.000 description 1
- 241001515806 Stictis Species 0.000 description 1
- 238000005273 aeration Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 238000007664 blowing Methods 0.000 description 1
- 210000000746 body region Anatomy 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 239000006143 cell culture medium Substances 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 238000005660 chlorination reaction Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000005984 hydrogenation reaction Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229920001983 poloxamer Polymers 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Physical Or Chemical Processes And Apparatus (AREA)
Abstract
本发明涉及一种微界面强化反应器能量耗散构效调控模型建模方法,将微界面强化反应器总的能量耗散率ε的计算划分为微界面强化反应器内三个不同区域能量耗散率的总和,包括反应器主体区鼓泡区的能量耗散率εR,气液破碎区的εmix以及气液出口区的εpl;确定εR、εmix和εpl各自的数值大小;获取ε与反应器结构参数相关的表达式。本发明的方法通过严谨的推导过程实现了微界面强化反应器能量耗散构效调控的数学模型的构建,为指导新型反应器的设计奠定理论基础。
Description
技术领域
本发明属于化工制造、反应器、建模技术领域,具体涉及一种微界面强化反应器能量耗散构效调控模型建模方法。
背景技术
氧化、加氢、氯化等多相反应在化工生产过程中广泛存在,其宏观反应速率一般受制于传质过程。气液反应的传质速率主要受液侧(或气侧)传质系数及气液相界面积a共同影响。已有研究表明,a对体积传质系数的影响程度更大,且容易调控。因此,增大a被视为提高受传质控制的气液反应体系反应效率的特别有效的途径。
气泡Sauter平均直径d32是决定a大小的关键参数之一,主要受气泡间及气液两相间相互作用力影响。气泡聚并和分裂分别是上述两种作用力的结果,并影响气泡直径的大小。因此,气泡聚并和破裂作为气泡的介观尺度行为,是决定a大小的深层次原因。关于气泡聚并和分裂行为的研究由来已久,普遍认为能量耗散率和d32是重要的影响因素。事实上,d32能够影响a及体积传质系数大小,是决定气液宏观反应速率的核心因素[1]。研究显示,当d32逐渐减小时,体积传质速率逐渐增大;特别是当d32小于1mm时,体积传质速率随d32的减小以类似于指数形式较快增大。因此尽可能减小d32能够强化气液传质并最终增大宏观反应速率。
鼓泡反应器和搅拌-鼓泡反应器是工业上最传统和常用的气液反应器。如PX氧化制TA的塔式鼓泡反应器,气泡直径通常大于10mm,乃至几厘米级,其传质界面面积十分有限,因此必须将反应器做得很大,以提高宏观反应速率,同时必须通过增加鼓气量来促进液体湍流,使气含率提高,进而增大界面面积,但此举必然降低空气中氧的利用率,增大压缩机功率和尾气排放,导致能耗过渡和物料损失及环境污染。从湍流动力学角度看,传统上用得最广泛的搅拌-鼓泡式气液反应器内大多形成对气泡宏观运动有影响但对气泡破碎作用甚微的大涡,气泡不能有效破碎,故气泡直径偏大,传质面积受限,以致反应效率偏低。为强化气液传质,塔式鼓泡反应器一般增设气体分布板、静态混合器等内件以加强混合,而搅拌釜则需安装不同结构的搅拌桨或内筒等结构,以增加液层的含气量。尽管如此,这两种反应器内的气泡直径通常为5~20mm,所提供的单位体积中的相界面积均十分有限,一般小于100m2/m3,故反应效率不可能获得突破性提高。因此,工业上经常通过高温高压和加大气量来提高气含率和相界面积,但这对反应过程的能耗、物耗及反应选择性都有重大的负面影响。
由于研发气泡的微破碎技术十分重要,故最近10年来,英、美、德、日等国的大学与研究机构开始关注和研发超细气泡技术[2-11],但其研究成果具有下列共同缺陷:
(1)采用机械破碎、流体撞击、超声等手段虽可得到一定量微米级尺度的气泡,但气液比(气体体积与液体体积之比)太低,一般低于1%,上限不超过5%。此外,产生微气泡的设备能耗和制造成本也太高。
(2)国内外尚未基于液相为连续相且高度湍流的微气泡体系特性,提出过系统化的微界面传质强化理论、微气泡测试与表征方法、微界面强化反应器构效调控理论及相关数学模型。
基于上述原因,虽有零星的应用实验结果发表,但尚未有规模化的工业应用报道,尤其在化工制造领域的应用基本上还处于空白状态。
当今时代的化工生产,基于创新、绿色、环保总体考虑,其生存和发展依赖于对材料与过程技术大幅创新。提高反应与分离过程的原子经济性对降低能耗、物耗、增强竞争力至关重要。基于此,我们提出“微界面传质强化反应-精细分离集成系统”新技术,试图从最基本的超细气液颗粒特性研究出发,解决在高度湍流状态下超细颗粒体系中,微界面化学反应器内流体流动、传质、反应、能量转换等构效调控全过程所涉及的理论、技术与应用问题。
本发明涉及的超细气液颗粒指的是超细气泡(或超细液滴),是颗粒当量直径处在1μm≤d0<1mm的微米级气液颗粒。在反应体系中,超细气液颗粒形成了超细界面(或微界面),超细界面的形成大幅强化了传质与反应速率,尤其是受传质控制的反应体系。
需要强调,经典的气液混合理论一般基于毫米-厘米级气液颗粒特征,目前最为合理的方法是多尺度能量最小化原理(EMMS)[12]。当前大多研究工作均是针对传统气液反应器内的毫米级上气泡[13,14],很少涉及超细颗粒体系。针对超细颗粒体系的混合、传质与反应特性,必须建立新的计算模型、测试与表征方法,以及构效调控模型,为此必须研究新的设备结构、能量输入方式与转换模式,从而形成全新的适合于超细颗粒反应体系的计算软件与硬件平台,为我国的过程工业生产技术上一个新台阶提供技术与装备支持。
已有研究指出[15,16],ε是决定气泡破裂机制的关键参数。一般认为,输入到体系中的能量以不同尺度的涡旋串级形式存在,最终在Kolmogorov尺度的耗散涡内耗散;大尺度涡主要提供体系宏观运动所需能量,尺寸小于气泡的含能涡通过撞击气泡表面而将能量传递给气泡,使气泡表面能增加,并导致气泡变形甚至破裂。ε增大,意味着气泡表面获得能量的速率增加,气泡更容易破裂。
关于能量耗散率较为经典的计算模型源于Cunningham基于热力学第一和第二定律对气液射流泵能量衡算的求解结果[17]。其研究的是射流断裂后与其携带的气体形成均匀气液混合物进行的能量传递,使气泡破碎。是针对传统反应器中产生气泡直径在1mm以上的气泡的。而对于微界面强化反应器不适用,后者是由于高速射流所携带的气体相互撞击进行能量交换而使气泡破碎的,两者机理不同。所以Cunningham基于热力学定律对εmin的表达式不适用。
所谓构效调控数学模型,是指将超细气液颗粒反应体系之反应效率(能效和物效)与体系理化特性、微界面特性、传质特性和反应器结构用数学方法关联起来,从而实现可通过调整结构参数和操作参数以获得反应过程能效物效的最大化目标,或者在给定反应目标(任务)和能耗物耗下,设计出高效的反应器结构。而对于微界面强化反应器而言,此方面工作在国际上仍为空白。
参考文献
[1]Levenspiel O.Chemical Reaction Engineering[M].Wiley New York etc.,1972.
[2]Xu JH,Li SW,Chen GG,LuoG..Formation of monodispersemicrobubbles ina microfluidic device[J].AIChE Journal,2006,52(6):2254-2259.
[3]Li P and Tsuge H.Ozone transfer in a new gas-induced contactorwith microbubbles[J].Journal of Chemical Engineering of Japan,2006,39(11):1213-1220.
[4]Muroyama K,Imai K,Oka Y,Hayashi J,Mass transfer properties in abubble column associated with micro-bubble dispersions[J].ChemicalEngineering Science,201,100:464-473.
[5]Maeda Y,Hosokawa S,Baba Y,Tomiyama Akio.Generation mechanism ofmicro-bubbles in a pressurized dissolution method[J].Experimental Thermal andFluid Science,2015,60:201-207.
[6]Hasegawa H,Nagasaka Y,Kataoka H.Electrical potential ofmicrobubble generated by shear flow in pipe with slits.Fluid DynamicsResearch,2008,40(7-8):554-564.
[7]Weber J and Agblevor F.Microbubble fermentation ofTrichodermareesei for cellulase production[J].Process Biochemistry,2005,40(2):669-676.
[8]Rehman F,Medley GJ,Bandulasena H,Zimmerman WB.Fluidic oscillator-mediated microbubble generation to provide cost effective mass transfer andmixing efficiency to the wastewater treatment plants[J].Environmentalresearch,2015,137:32-39.
[9]Stride E and Edirisinghe M.Novel microbubble preparationtechnologies[J].Soft Matter,2008,4(12):2350.
[10]Druzinec D,Salzig De,Kraume M,Czermak P.Micro-bubble aeration inturbulent stirred bioreactors:Coalescence behavior in Pluronic F68containingcell culture media[J].Chemical Engineering Science,2015,126:160-168.
[11]李宝璋,尚龙安,姜信真.倒喷射式环流反应器的研究[J].西北大学学报(自然科学版).1989,04:65-69.
[12]Chen JH,Yang N,Ge W,Li JH.Stability-driven structure evolution:exploring the intrinsic similarity between gas-solid and gas-liquid systems[J].Chinese Journal of Chemical Engineering.2012,20(1):167-177.
[13]Hinze JO.Fundamentals of the hydrodynamic mechanism of splittingin dispersion processes[J].AIChE Journal.1955,1(3):289-295.
[14]Zhong S,Zou X,Zhang ZB,Tian HZ.A flexible image analysis methodfor measuring bubble parameters[J].Chemical Engineering Science,2016,141(17):143-153.
[15]Baldy,S.,A Generation-Dispersion Model of Ambient and TransientBubbles in the Close Vicinity of Breaking Waves.Journal of GeophysicalResearch-Oceans,1993.98(C10):p.18277-18293.
[16]Prince,M.J.and H.W.Blanch,Bubble Coalescence and Break-up in Air-Sparged Bubble-Columns.Aiche Journal,1990.36(10):p.1485-1499.
[17]Cunningham,R.,Gas compression with the liquid jet pump.Journal offluids engineering,1974.96(3):203-215.
发明内容
本发明的目的在于,克服现有技术的缺陷,提供一种微界面强化反应器能量耗散构效调控模型建模方法。
针对以上背景技术,本发明的内容是构建微界面强化反应器(MTIR)构效调控数学模型。所谓构效调控数学模型,是指将超细气液颗粒反应体系之反应效率(能效和物效)与体系理化特性、微界面特性、传质特性和反应器结构用数学方法关联起来,从而实现可通过调整结构参数和操作参数以获得反应过程能效物效的最大化目标,或者在给定反应目标(任务)和能耗物耗下,设计出高效的反应器结构。而对于微界面强化反应器而言,此方面工作在国际上仍为空白。
为实现上述技术目的,本发明采用如下技术方案:
一种微界面强化反应器能量耗散构效调控模型建模方法,包括:
步骤100:将微界面强化反应器总的能量耗散率ε的计算划分为微界面强化反应器内三个不同区域能量耗散率的总和,包括反应器主体区鼓泡区的能量耗散率εR,气液破碎区的εmix以及气液出口区的εpl;
步骤110:其中,反应器主体区鼓泡区的能量耗散率εR采用如下方式计算:
反应器气体鼓泡过程中,依据气泡对体系做功,εR表示为:
其中,QG为反应器内通气体积流量,m3/s;S0为反应器横截面积,m2;
步骤120:计算气液破碎区的εmix:
基于εmix经典计算模型,假设气液混合为绝热过程并忽略液体势能变化,忽略气体质量流量,并使能量耗散率的单位统一为W/Kg,获取计算公式如下:
其中,Lmix为气泡破碎区长度,m;P0、P1分别为气泡破碎区入口液体静压及出口气液混合物压力,Pa;λ1为气液体积流量之比:K1为喷嘴直径与气泡破碎区直径的比值,K1=DN/D1;S1为气泡破碎器横截面积,m2;ρL为液体密度,kg/m3;QL为反应器内液体循环体积流量,m3/s;
λ1=QG/QL (3)
步骤121:计算气泡破碎区入口液体静压P0及出口气液混合物压力P1:
忽略气泡破碎区管壁摩擦损失,则:
其中,φmix为气泡破碎区气含率,按下式计算:
忽略管道摩擦及喷嘴处能量损失,依据能量守恒原理,体系实际获得的能量E0为:
即:
由式(6)(9)得:
步骤122:计算气泡破碎区长度Lmix:
Lmix通过测量气液破碎区管内壁压力突变来确定,或通过如下方式确定:
其中:PH为气液破碎区上方气压,Pa;ρMZ为气液破碎区内气液混合物密度,Kg/m3;vN为射流口的射流速度,m/s;Ue,max为气液破碎区涡旋最大返回速度,m/s;
PH由伯努利方程推得:
PH≈PG0 (13)
式中,PG0为供气压力,Pa;
ρMZ通过下式计算:
ρMZ=ρGφmix+ρL(1-φmix)≈ρL(1-φmix) (14)
式中,ρG为气体密度,g/m3;
考虑气液破碎区气液混合物流速的影响,Ue,max为射流口射流流速与气液破碎区气液混合物流速的矢量合成的结果,采用下式计算Ue,max:
将式(13)(16)代入式(12),并化简后可得:
获取反应器气泡破碎器长度Lb,并根据式(17)计算Lmix;
①当Lmix<Lb时,以式(17)的计算结果为Lmix的实际数值;
②当Lmix≥Lb时,说明射流能量近似完全在气泡破碎器区域耗,则:
Lmix=Lb (20)
步骤130:计算气液出口区的εpl;
假设气泡在气液出口区呈均匀分布状态,此区域的能量耗散速率εpl由下式计算:
反应器结构设计时保证λ1可调范围足够大,通过实验确定反应器基本结构参数间的关系:K1=0.5,Lb=13D1;代入前述相应表达式并化简可得:
步骤200:确定εR、εmix和εpl各自的数值大小;
步骤210:依据进入反应器的气液流量等于气泡破碎区出口的气液流量平衡原理,得到:
式中,CL为基于反应器内有效体积πD0 2H0/4的液体循环倍数,即每小时液体循环总体积与反应器有效体积的比值;u1为气泡破碎器出口气液混合物线速度,m/s;λ1取值0.1~0.5;
由式(25)可知:则u1增大时,反应器横截面积S0也增大;结合式(1)可知,此时εR减小;为对反应器不同区域的能量耗散率进行比较,假设:u1=3.0m/s;CL=20;H0=1.5m;由式(25)可得,当λ1=0.1~0.5时:
D0≈19D1 (26)
选定D1数值,计算并比较反应器不同区域在不同喷嘴液速下的能量耗散率,确定与气液破碎区的能量耗散率εmix相比,反应器主体区、气液出口区的能量耗散率可忽略不计,即εmix≈ε;则整个反应器的能量耗散率ε与反应器结构参数之间的数学关系,可由式(2)计算确定,即:
本发明的另一目的在于提供上述方法构建的数学模型。
本发明的又一目的在于提供上述方法在反应器设计中的应用。
具体包括:
设计反应器使其结构参数满足λ1=0.1~0.5、K1=0.5、Lb=13D1,以及下述方程:
其中,Lb为反应器气泡破碎器长度;D1为气泡破碎区直径;K1为喷嘴直径与气泡破碎区直径的比值,K1=DN/D1;λ1为气液体积流量之比;QL为反应器内液体循环体积流量;Lmix为气泡破碎区长度;D1为气泡破碎管直径;λ1为气液体积流量之比,λ1=QG/QL;QG为反应器内通气体积流量;P0为气泡破碎器入口处液体的静压;P1为气泡破碎区出口气液混合物压力;ε为能量耗散率;μL为液体动力粘度;σL为液体表面张力;ρL为液体密度。
本发明的方法适用于微界面强化反应器,其核心在于气泡破碎器。气泡破碎器的原理是高速射流所携带的气体相互撞击进行能量传递,使气泡破碎,其结构参数有Lb、D1,详细结构见附图1,除此之外该反应器的其它结构参数有D0、H0,具体反应器结构相关内容已经公布于发明人在先申请的专利CN10618766A中,本发明中不再赘述。
因有研究指出,能量耗散率ε是决定气泡破裂机制的关键参数。而ε的增大,意味着气泡更容易破裂,即可以得到尺度更小的气泡,但是针对微米级尺度的能量耗散率的表达式,目前相关研究还是空白,本发明的方法通过严谨的推导过程实现了微界面强化反应器能量耗散构效调控的数学模型的构建,为指导新型反应器的设计奠定理论基础。
本发明的方法通过构建微界面强化反应器的能量耗散构效调控模型,填补了气液体系微米级尺度能量耗散率在国际上的空白。而且构效调控的数学模型,能将超细气液颗粒反应体系之反应效率(能效和物效)与体系理化特性、微界面特性、传质特性和反应器结构用数学方法关联起来,从而实现可通过调整结构参数和操作参数以获得反应过程能效物效的最大化目标,或者在给定反应目标(任务)和能耗物耗下,设计出高效的反应器结构。
附图说明
图1是一种微界面强化反应器结构示意图,用于说明本发明建模方法在反应器装置中的应用;其中1-反应器,2-泵前阀,3-循环泵,4-泵后阀,5-液体流量计,6-换热器,7-气泡破碎器,8-测温仪,9-下降管,10-气体阀,11-气体流量计,12-气相入口,13-压力表,14-液位计;D0-反应器直径,H0-反应器内初始液位高度,D1-气泡破碎管直径,Lb-气泡破碎区长度;
图2是反应器内不同区域能量耗散率大小比较;图2(a)是反应器主体区鼓泡区的εR,图2(b)是气液出口区的εpl;
图3是3种代表性体系中d0与ε的关系。
具体实施方式
实施例1
本实施例具体说明本发明的微界面强化反应器能量耗散构效调控模型建模方法。
其步骤包括:
步骤100:对于基于高度湍动的射流能量破碎气体的微界面强化反应器,将其总的能量耗散率ε的计算划分为微界面强化反应器内三个不同区域能量耗散率的总和,包括反应器主体区鼓泡区的能量耗散率εR,气液破碎区的εmix以及气液出口区的εpl;
步骤110:其中,反应器主体区鼓泡区的能量耗散率εR采用如下方式计算:
反应器气体鼓泡过程中,依据气泡对体系做功,εR表示为:
其中,QG为反应器内通气体积流量,m3/s;S0为反应器横截面积,m2;
在实际鼓泡过程中,表观气速vG由于受液体静压的限制而不可能很大,因此,鼓泡过程εR值也较小,这是导致鼓泡塔内气泡较大的重要原因。因此,想要获得较小的气泡,尤其是微气泡,必须改变自然鼓泡方式而采用其它方式。
步骤120:计算气液破碎区的εmix:
气泡破碎区能量耗散率εmix经典计算模型源于Cunningham基于热力学第一和第二定律对气液射流泵能量衡算的求解结果。但是本发明应用的气泡破碎器与其有所不同:前者是射流断裂后与其携带的气体形成均匀气液混合物进行能量传递,使气泡破碎,而后者是高速射流所携带的气体相互撞击进行能量传递,使气泡破碎。
基于εmix经典计算模型,假设气液混合为绝热过程并忽略液体势能变化,忽略气体质量流量,并使能量耗散率的单位统一为W/Kg,获取计算公式如下:
其中,Lmix为气泡破碎区长度,m;P0、P1分别为气泡破碎区入口液体静压及出口气液混合物压力,Pa;λ1为气液体积流量之比:K1为喷嘴直径与气泡破碎区直径的比值,K1=DN/D1;S1为气泡破碎器横截面积,m2;ρL为液体密度,kg/m3;QL为反应器内液体循环体积流量,m3/s;
λ1=QG/QL (3)
步骤121:计算气泡破碎区入口液体静压P0及出口气液混合物压力P1:
Evans假设气泡破碎区压力各处相等,即:
P1≈P0 (4)
并对式(2)进行了修正,得到如下简化形式:
理论和实验结果均表明,严格意义上式(5)不合理。更为重要的是,体系总输入能量对εmix的影响未得到体现。若忽略气泡破碎区管壁摩擦损失,则:
其中,φmix为气泡破碎区气含率,按下式计算:
忽略管道摩擦及喷嘴处能量损失,依据能量守恒原理,体系实际获得的能量E0为:
即:
当液体饱和蒸汽压为PL *(Pa),为使喷嘴处液体不汽化,应满足:
当喷嘴及气泡破碎器结构尺寸一定情况下,上式决定了体系所需最小能量。换言之,若系统已选用确定型号的循环泵,则实际循环流量由不等式(10)所限。
由式(6)(9)得:
步骤122:计算气泡破碎区长度Lmix:
Lmix通过测量气液破碎区管内壁压力突变来确定,或通过如下方式确定:
其中:PH为气液破碎区上方气压,Pa;ρMZ为气液破碎区内气液混合物密度,Kg/m3;vN为射流口的射流速度,m/s;Ue,max为气液破碎区涡旋最大返回速度,m/s;
PH由伯努利方程推得:
PH≈PG0 (13)
式中,PG0为供气压力,Pa;
ρMZ通过下式计算:
ρMZ=ρGφmix+ρL(1-φmix)≈ρL(1-φmix) (14)
式中,ρG为气体密度,g/m3;
Ue,max一般可由下式计算:
Ue,max=0.085vN (15)
但式(15)仅考虑Ue,max与射流口液体流速有关,而并未考虑与气液破碎区气液混合物流速的影响。本发明研究认为,Ue,max实际应为射流口射流流速与气液破碎区气液混合物流速的矢量合成的结果。因此,考虑气液破碎区气液混合物流速的影响,采用下式计算Ue,max:
将式(13)(16)代入式(12),并化简后可得:
获取反应器气泡破碎器长度Lb,并根据式(17)计算Lmix,当气泡破碎器长度Lb(m)确定时,Lmix随操作工况的不同而发生改变,有可能出现如下两种情况,因此分情况讨论;
①当Lmix<Lb时,以式(17)的计算结果为Lmix的实际数值;
此时Lmix>0,K1还应同时满足下面不等式组:
K1 2(1+λ1)≤0.915 (19)
②当Lmix≥Lb时,说明射流能量近似完全在气泡破碎器区域耗,则:
Lmix=Lb (20)
步骤130:计算气液出口区的εpl;
假设气泡在气液出口区呈均匀分布状态,此区域的能量耗散速率εpl由下式计算:
上式是经化简后的结果。其中,νL为液体运动粘度,m2/s。不难看出,εpl主要与反应器操作参数(QL、QG)及气泡破碎器的气液出口管直径D1有关。
由以上分析可知,反应体系多因素(操作参数、物性参数和结构参数)的联合作用导致εR、εmix和εpl的变化较为复杂。而且进一步分析可知,气液体积流量比λ1是调控a的关键参数,可根据实际生产需要对λ1进行调节。为此,反应器结构设计时应保证λ1可调范围足够大,通过实验确定反应器基本结构参数间的关系:K1=0.5,Lb=13D1;代入前述相应表达式并化简可得:
由于本发明的超细气泡是在气泡破碎器中形成,因此应建立气泡破碎器的能量耗散率的计算模型。前文已指出,反应器内存在三个不同大小的区域能量耗散率,究竟采用何种形式,须首先对它们进行比较。
步骤200:确定εR、εmix和εpl各自的数值大小;
步骤210:要确定εR、εmix和εpl各自的数值大小,须确定反应器直径D0与气泡破碎器直径D1的关系,理论上D0与D1并不存在定量关系。为此,依据进入反应器的气液流量等于气泡破碎区出口的气液流量平衡,可得到
式中,CL为基于反应器内有效体积πD0 2H0/4的液体循环倍数,即每小时液体循环总体积与反应器有效体积的比值;u1为气泡破碎器出口气液混合物线速度,m/s;λ1取值0.1~0.5;
由式(25)可知:则u1增大时,反应器横截面积S0也增大;结合式(1)可知,此时εR减小;为对反应器不同区域的能量耗散率进行比较,假设:u1=3.0m/s;CL=20;H0=1.5m;由式(25)可得,当λ1=0.1~0.5时:
D0≈19D1 (26)
式(26)基本可以反映大多数情况下反应器直径与气泡破碎管直径的关系。在特殊情形时,可以适当调整D1前面的数值。
以D1=0.012m为例,计算并比较反应器不同区域在不同喷嘴液速下的能量耗散率,如图2所示。其中vN为喷嘴液速,m/s。由图2可知,与气液破碎区的能量耗散率εmix相比,反应器主体区、气液出口区的能量耗散率可忽略不计,即εmix≈ε;则整个反应器的能量耗散率ε与反应器结构参数之间的数学关系,可由式(2)计算确定,即:
实施例2
本实施例以图1所示的反应器为例,说明实施例1所述建模方法构建的模型在二氧化碳和水体系反应器中的应用。图1的反应器结构可为现有微界面强化反应器的结构,仅采用本发明的方法进行参数设计,本发明中对反应器的结构不再赘述。
根据实施例1构建的能量耗散构效调控模型如下:
式中,QL为反应器内液体循环体积流量;Lmix为气泡破碎区长度;D1为气泡破碎管直径;λ1为气液体积流量之比,λ1=QG/QL;QG为反应器内通气体积流量;P0为气泡破碎器入口处液体的静压;P1为气泡破碎区出口气液混合物压力;ε为能量耗散率;μL为液体动力粘度;σL为液体表面张力;ρL为液体密度。
实施例2选用的模型主要考虑Lmix小于Lb的情况,因为相反的情况不常见,比较极端。
根据实施例1,反应器结构参数还应满足:λ1=0.1~0.5、K1=0.5、Lb=13D1;
对于二氧化碳和水体系,当操作条件为:QL=2000L/h(5.56×10-4m3/s),气体流量QG=0.2QL,T=298K,PG0=1atm;而此体系中液相的物性参数为:ρL=1000kg/m3,μL=8.9×10-4Pa·s,σL=7.197×10-4N/m;反应器气泡破碎管直径D1=0.02m;E0表示体系输入的能量,即循环泵铭牌上的额定功率,取E0=1000W。根据操作条件和上述模型应用MIR计算得到产生气泡Sauter平均直径d32=0.1mm时的能量耗散率ε=1001.03W/kg。
Levenspiel认为,多相体系的宏观反应速率可由下式表示:
经化简后的气液反应宏观速率方程可简化为:
表1和表2是在MIR中相同体系的各个参数的计算情况:
表1在MIR反应器中模型公式计算的参数
表2在MIR反应器中模型公式计算的三种阻力(气膜、液膜、本征)
实施例3:
本实施例以图1所示的反应器为例,说明实施例1所述建模方法构建的模型在空气-水体系反应器中的应用,与现有装置产生的气泡Sauter平均直径d32=1mm时的能量耗散率的结果对比如下。
表3和表4是相同体系不同粒径下的各个参数的对比情况:
表3不同粒径下模型公式计算的参数
表4不同粒径下模型公式计算的三种阻力(气膜、液膜、本征)
如表3、表4所示,在相同条件下,产生的直径越小,所需的能量耗散率则越大。微米级气泡的产生,是由于高速射流所携带的气体相互碰撞进行能量传递,使能量耗散率最大的气液破碎区会产生很多小的旋涡及涡流,剧烈的湍动,较大的气泡会在相互碰撞和剪切力的作用下变成更小的气泡。因此会减小气泡直径,增大相界面积,从而减小了传质阻力和强化了气液传质速率。
能量耗散率ε是决定单一气泡直径d0的关键变量,且由反应器调控参数决定。对于鼓泡塔,有如下计算公式:
由于实际鼓泡过程vG受液体静压的限制而不可能很大,故鼓泡过程ε值一般较小,这决定了鼓泡塔内d0一般处在毫米-厘米级水平。但若采用强制气液喷射搅拌方式,ε值就会通过专门的结构设计变得很大。由喷射器理论可知,体系初始气泡的大小由εmix决定,而初始气泡的大小将直接决定反应器中气泡大小及分布。
当不考虑射流破裂情形时,从理论上考察4种代表性反应体系中ε与d0的关系,见图3。由图可知,当ε较小时,体系物性对d0的影响较显著;而随着ε的逐渐增大,体系物性的影响逐渐减小。
对于普通鼓泡反应器而言,由于其ε值一般小于1W/kg,故其d0一般都大于1mm。而微界面强化反应器可以实现ε=10~100W/kg之间,因此其气泡尺寸可以达到微米级。
Claims (4)
1.一种微界面强化反应器能量耗散构效调控模型建模方法,其特征在于,包括:
步骤100:将微界面强化反应器总的能量耗散率ε的计算划分为微界面强化反应器内三个不同区域能量耗散率的总和,包括反应器主体区鼓泡区的能量耗散率εR,气液破碎区的εmix以及气液出口区的εpl;
步骤110:其中,反应器主体区鼓泡区的能量耗散率εR采用如下方式计算:
反应器气体鼓泡过程中,依据气泡对体系做功,εR表示为:
其中,QG为反应器内通气体积流量,m3/s;S0为反应器横截面积,m2;
步骤120:计算气液破碎区的εmix:
基于εmix经典计算模型,假设气液混合为绝热过程并忽略液体势能变化,忽略气体质量流量,并使能量耗散率的单位统一为W/Kg,获取计算公式如下:
其中,Lmix为气泡破碎区长度,m;P0、P1分别为气泡破碎区入口液体静压及出口气液混合物压力,Pa;λ1为气液体积流量之比:K1为喷嘴直径与气泡破碎区直径的比值,K1=DN/D1;S1为气泡破碎器横截面积,m2;ρL为液体密度,kg/m3;QL为反应器内液体循环体积流量,m3/s;
λ1=QG/QL (3)
步骤121:计算气泡破碎区入口液体静压P0及出口气液混合物压力P1:
忽略气泡破碎区管壁摩擦损失,则:
其中,φmix为气泡破碎区气含率,按下式计算:
忽略管道摩擦及喷嘴处能量损失,依据能量守恒原理,体系实际获得的能量E0为:
即:
由式(6)(9)得:
步骤122:计算气泡破碎区长度Lmix:
Lmix通过测量气液破碎区管内壁压力突变来确定,或通过如下方式确定:
其中:PH为气液破碎区上方气压,Pa;ρMZ为气液破碎区内气液混合物密度,Kg/m3;vN为射流口的射流速度,m/s;Ue,max为气液破碎区涡旋最大返回速度,m/s;
PH由伯努利方程推得:
PH≈PG0 (13)
式中,PG0为供气压力,Pa;
ρMZ通过下式计算:
ρMZ=ρGφmix+ρL(1-φmix)≈ρL(1-φmix) (14)
式中,ρG为气体密度,g/m3;
考虑气液破碎区气液混合物流速的影响,Ue,max为射流口射流流速与气液破碎区气液混合物流速的矢量合成的结果,采用下式计算Ue,max:
将式(13)(16)代入式(12),并化简后可得:
获取反应器气泡破碎器长度Lb,并根据式(17)计算Lmix;
①当Lmix<Lb时,以式(17)的计算结果为Lmix的实际数值;
②当Lmix≥Lb时,说明射流能量近似完全在气泡破碎器区域耗,则:
Lmix=Lb (20)
步骤130:计算气液出口区的εpl;
假设气泡在气液出口区呈均匀分布状态,此区域的能量耗散速率εpl由下式计算:
反应器结构设计时保证λ1可调范围足够大,通过实验确定反应器基本结构参数间的关系为K1=0.5,Lb=13D1;代入前述相应表达式并化简可得:
步骤200:确定εR、εmix和εpl各自的数值大小;
步骤210:依据进入反应器的气液流量等于气泡破碎区出口的气液流量平衡原理,得到:
式中,CL为基于反应器内有效体积πD0 2H0/4的液体循环倍数,即每小时液体循环总体积与反应器有效体积的比值;u1为气泡破碎器出口气液混合物线速度,m/s;λ1取值0.1~0.5;
由式(25)可知:则u1增大时,反应器横截面积S0也增大;结合式(1)可知,此时εR减小;为对反应器不同区域的能量耗散率进行比较,假设:u1=3.0m/s;CL=20;H0=1.5m;由式(25)可得,当λ1=0.1~0.5时:
D0≈19D1 (26)
选定D1数值,计算并比较反应器不同区域在不同喷嘴液速下的能量耗散率,确定与气液破碎区的能量耗散率εmix相比,反应器主体区、气液出口区的能量耗散率可忽略不计,即εmix≈ε;则整个反应器的能量耗散率ε与反应器结构参数之间的数学关系,可由式(2)计算确定,即:
。
2.权利要求1所述方法构建的能量耗散构效调控模型,其特征在于,构建的能量耗散构效调控模型如下:
式中,QL为反应器内液体循环体积流量;Lmix为气泡破碎区长度;D1为气泡破碎管直径;λ1为气液体积流量之比,λ1=QG/QL;QG为反应器内通气体积流量;P0为气泡破碎器入口处液体的静压;P1为气泡破碎区出口气液混合物压力;ε为能量耗散率;μL为液体动力粘度;σL为液体表面张力;ρL为液体密度。
3.权利要求1所述方法在反应器设计中的应用。
4.根据权利要求3所述的应用,其特征在于,设计反应器使其结构参数满足λ1=0.1~0.5、K1=0.5、Lb=13D1,以及下述方程:
其中,Lb为反应器气泡破碎器长度;D1为气泡破碎区直径;K1为喷嘴直径与气泡破碎区直径的比值,K1=DN/D1;λ1为气液体积流量之比;QL为反应器内液体循环体积流量;Lmix为气泡破碎区长度;D1为气泡破碎管直径;λ1为气液体积流量之比,λ1=QG/QL;QG为反应器内通气体积流量;P0为气泡破碎器入口处液体的静压;P1为气泡破碎区出口气液混合物压力;ε为能量耗散率;μL为液体动力粘度;σL为液体表面张力;ρL为液体密度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710766690.5A CN107589667B (zh) | 2017-08-30 | 2017-08-30 | 微界面强化反应器能量耗散构效调控模型建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710766690.5A CN107589667B (zh) | 2017-08-30 | 2017-08-30 | 微界面强化反应器能量耗散构效调控模型建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107589667A true CN107589667A (zh) | 2018-01-16 |
CN107589667B CN107589667B (zh) | 2020-09-04 |
Family
ID=61050488
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710766690.5A Active CN107589667B (zh) | 2017-08-30 | 2017-08-30 | 微界面强化反应器能量耗散构效调控模型建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107589667B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109887550A (zh) * | 2019-01-09 | 2019-06-14 | 南京大学 | Miha纯气动操作条件下能量耗散调控模型建模方法 |
CN113035285A (zh) * | 2021-04-01 | 2021-06-25 | 南京延长反应技术研究院有限公司 | 微气泡尺寸对油品加氢脱硫效果影响的计算方法 |
CN113075091A (zh) * | 2021-02-24 | 2021-07-06 | 南京延长反应技术研究院有限公司 | 一种确定px氧化制ta体系中气泡大小的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7442360B2 (en) * | 2004-04-27 | 2008-10-28 | Velocys, Inc. | Hydrogen peroxide production in microchannel reactors |
CN103678904A (zh) * | 2013-12-06 | 2014-03-26 | 杭州电子科技大学 | 一种动态Gibbs反应器建模方法 |
CN106933262A (zh) * | 2015-12-29 | 2017-07-07 | 天津城建大学 | 一种在线压力信号解析的多相流场量测及气-液传质增效控制方法 |
-
2017
- 2017-08-30 CN CN201710766690.5A patent/CN107589667B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7442360B2 (en) * | 2004-04-27 | 2008-10-28 | Velocys, Inc. | Hydrogen peroxide production in microchannel reactors |
CN103678904A (zh) * | 2013-12-06 | 2014-03-26 | 杭州电子科技大学 | 一种动态Gibbs反应器建模方法 |
CN106933262A (zh) * | 2015-12-29 | 2017-07-07 | 天津城建大学 | 一种在线压力信号解析的多相流场量测及气-液传质增效控制方法 |
Non-Patent Citations (2)
Title |
---|
张娇: "微通道内微细颗粒增强气液传质的研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
谢灵丹: "微通道内微细颗粒对气液传质的影响研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109887550A (zh) * | 2019-01-09 | 2019-06-14 | 南京大学 | Miha纯气动操作条件下能量耗散调控模型建模方法 |
CN109887550B (zh) * | 2019-01-09 | 2022-06-21 | 南京大学 | Miha纯气动操作条件下能量耗散调控模型建模方法 |
CN113075091A (zh) * | 2021-02-24 | 2021-07-06 | 南京延长反应技术研究院有限公司 | 一种确定px氧化制ta体系中气泡大小的方法 |
CN113075091B (zh) * | 2021-02-24 | 2023-01-10 | 南京延长反应技术研究院有限公司 | 一种确定px氧化制ta体系中气泡大小的方法 |
CN113035285A (zh) * | 2021-04-01 | 2021-06-25 | 南京延长反应技术研究院有限公司 | 微气泡尺寸对油品加氢脱硫效果影响的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107589667B (zh) | 2020-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107563051B (zh) | 微界面强化反应器气泡尺度构效调控模型建模方法 | |
CN107561938A (zh) | 微界面强化反应器反应速率构效调控模型建模方法 | |
CN107335390A (zh) | 微界面强化反应器相界面积构效调控模型建模方法 | |
CN107346378A (zh) | 微界面强化反应器传质速率构效调控模型建模方法 | |
JP3231870U (ja) | マイクロ界面強化水素化反応システム | |
Cramers et al. | Hydrodynamics and mass transfer characteristics of a loop-venturi reactor with a downflow liquid jet ejector | |
CN107589667A (zh) | 微界面强化反应器能量耗散构效调控模型建模方法 | |
Poncin et al. | Hydrodynamics and volumetric gas–liquid mass transfer coefficient of a stirred vessel equipped with a gas-inducing impeller | |
Moraveji et al. | Gas‐Liquid Hydrodynamics and Mass Transfer in Aqueous Alcohol Solutions in a Split‐Cylinder Airlift Reactor | |
Bao et al. | Influence of impeller diameter on overall gas dispersion properties in a sparged multi-impeller stirred tank | |
Dong et al. | Effect of liquid phase rheology and gas–liquid interface property on mass transfer characteristics in bubble columns | |
CN113257367B (zh) | 一种基于苯酚废水的湿式氧化构效调控模型的建模方法 | |
Krichnavaruk et al. | Analysis of gas–liquid mass transfer in an airlift contactor with perforated plates | |
Gaddis et al. | The impinging—stream reactor: a high performance loop reactor for mass transfer controlled chemical reactions | |
Kleingeld et al. | The development and modelling of high-intensity impinging stream jet reactors for effective mass transfer in heterogeneous systems | |
Gourich et al. | Influence of operating conditions and design parameters on hydrodynamics and mass transfer in an emulsion loop–venturi reactor | |
Lu et al. | Operating and hydrodynamic characteristics of a cocurrent downflow bubble column reactor | |
Sivaiah et al. | Gas entrainment and holdup characteristics in a modified gas–liquid–solid down flow three-phase contactor | |
Ide et al. | Mass transfer characteristics in gas bubble dispersed phase generated by plunging jet containing small solute bubbles | |
Gavrilescu et al. | Concentric-tube airlift bioreactors: Part III: Effects of geometry on mass transfer | |
Yamagiwa et al. | Effects of liquid property on air entrainment and oxygen transfer rates of a plunging jet reactor | |
Meng et al. | Pierced cylindrical gas inlet device for sulfur dioxide removal from waste gas streams | |
Liu et al. | Region-dependent mass transfer behavior in a forced circulation airlift loop reactor | |
Majeed et al. | Design correlations of an air-lift tube reactor | |
Rong et al. | Numerical simulation of micro-nano bubble central intake aeration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |