CN107346378A - 微界面强化反应器传质速率构效调控模型建模方法 - Google Patents
微界面强化反应器传质速率构效调控模型建模方法 Download PDFInfo
- Publication number
- CN107346378A CN107346378A CN201710766695.8A CN201710766695A CN107346378A CN 107346378 A CN107346378 A CN 107346378A CN 201710766695 A CN201710766695 A CN 201710766695A CN 107346378 A CN107346378 A CN 107346378A
- Authority
- CN
- China
- Prior art keywords
- bubble
- reactor
- formula
- liquid
- gas
- 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
- 238000012546 transfer Methods 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000007788 liquid Substances 0.000 claims abstract description 210
- 238000006243 chemical reaction Methods 0.000 claims abstract description 43
- 230000000694 effects Effects 0.000 claims abstract description 19
- 238000013461 design Methods 0.000 claims abstract description 7
- 239000000203 mixture Substances 0.000 claims description 77
- 239000012530 fluid Substances 0.000 claims description 28
- 230000003068 static effect Effects 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 15
- 230000033001 locomotion Effects 0.000 claims description 15
- 238000009792 diffusion process Methods 0.000 claims description 14
- 239000007791 liquid phase Substances 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 13
- 230000005587 bubbling Effects 0.000 claims description 11
- 239000002245 particle Substances 0.000 claims description 11
- 230000004044 response Effects 0.000 claims description 10
- 238000009423 ventilation Methods 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 7
- 230000010354 integration Effects 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 5
- 238000009776 industrial production Methods 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 241000197192 Bulla gouldiana Species 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000005192 partition Methods 0.000 claims description 3
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 239000002904 solvent Substances 0.000 claims description 3
- 239000007921 spray Substances 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 11
- 239000000463 material Substances 0.000 abstract description 9
- 239000007789 gas Substances 0.000 description 71
- 239000010410 layer Substances 0.000 description 10
- 239000012071 phase Substances 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 4
- 238000005265 energy consumption Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000004581 coalescence Methods 0.000 description 3
- 230000001965 increasing effect 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
- 238000009825 accumulation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007664 blowing Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 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
- 238000001311 chemical methods and process Methods 0.000 description 1
- 238000005660 chlorination reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram 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
- 239000007792 gaseous phase Substances 0.000 description 1
- 239000003292 glue Substances 0.000 description 1
- 238000005984 hydrogenation reaction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000002609 medium Substances 0.000 description 1
- VIKNJXKGJWUCNN-XGXHKTLJSA-N norethisterone Chemical compound O=C1CC[C@@H]2[C@H]3CC[C@](C)([C@](CC4)(O)C#C)[C@@H]4[C@@H]3CCC2=C1 VIKNJXKGJWUCNN-XGXHKTLJSA-N 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
- 230000009467 reduction Effects 0.000 description 1
- 238000000926 separation method 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
- 230000007723 transport mechanism Effects 0.000 description 1
- 239000005436 troposphere Substances 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
Abstract
本发明公开了一种微界面强化反应器传质速率构效调控模型建模方法，通过严谨的理论推导分别建立了建立气侧传质系数计算模型和液侧传质系数计算模型。采用本发明的建模方法构建的传质速率构效调控模型可以直观的看出传质速率和气泡大小的关系，为研究微界面体系奠定了理论基础，从而也可通过调整结构参数和操作参数以实现获得反应过程能效物效的最大化目标，或者在给定反应目标(任务)和能耗物耗下，设计出高效的反应器结构。
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]，很少涉及超细颗粒体系。针对超细颗粒体系的混合、传质与反应特性，必须建立新的计算模型、测试与表征方法，以及构效调控模型，为此必须研究新的设备结构、能量输入方式与转换模式，从而形成全新的适合于超细颗粒反应体系的计算软件与硬件平台，为我国的过程工业生产技术上一个新台阶提供技术与装备支持。
传统反应器的传质建模是基于实验数据拟合得到体积传质参数(如液相体积传质系数kLa、气相体积传质系数kGa等)经验或半经验的关联式。此类模型较多，其缺陷在于：所得关联式是在特定条件下(包括特定反应器结构、搅拌方式、气液物性、操作工况、测量方法等)得到的，难以普遍适用。
基于连续介质假设的Navier-Stokes方程(N-S方程)，理论上可以得到设计参数。但由于工业气液反应体系中气液界面存在复杂的时空变化，即使采用上述方法能获得数值解，其结果的适用性仍令人生疑。主要原因在于求解方程所依赖的边界条件受反应器结构(外壳和内件、搅拌器结构等)和流体运动方式的影响较大，且很难确定。
气液界面传质阻力由kLa和kGa大小决定。但由于多数气液反应，如加氢和氧化，气体在液相中的溶解度不大，气侧传质阻力较液侧一般要小几个数量级[15,16]，故多数情况下可以忽略气侧传质阻力。因而气液界面传质建模一般是针对反应器结构、气液物性、操作工况和kLa之间的关系展开[17-23]。但就传质过程而言，kLa是一个宏观参数，它并没有充分反映气液传质界面及其两侧的详细信息。因此，只有将kLa中的两项(kL和a)分别加以研究才能获取更多微观信息。为建立普遍适用的气液界面传质模型，还须考虑气侧传质阻力，一概忽略该阻力是不科学的，应依据不同体系的计算或测定结果决定是否可忽略。因此，为建立精确普适性的反应体系气液传质模型，需要对kL、kG和气液相界面积a分别进行研究[24-26]。
经典kG数学模型一般膜理论。据此，决定其大小的参数有温度T、亨利系数HA、气相扩散系数DG及有效气膜厚度δG。由于影响δG大小的因素较多，且不易测量，故尚未见其理论计算模型。除此之外，实际可采用Gedde静止球模型[27]：其中，t32为直径为d32的气泡在体系中的停留时间，s。本发明也采用的此表达式。但因所采用的d32的表达式不同，所以，修改后的气侧传质速率也将不同。
相对于kG而言，kL的理论研究较多。自上世纪20年代发展至今，已有包括双膜理论[28]、渗透理论[29]及表面更新理论[30]三种经典理论及在它们基础上的一些改进型模型[31-35]，如速度滑移模型[36]，界面湍流理论[37,38]，单涡模型[39]，表面更新伸展模型[40]等。而基于膜理论的kL模型中，有效液膜厚度δL的准确测量存在与δG相同的测量难度的问题，因此用于实际理论计算的意义不大[41,42]。而基于表面更新理论的kL渗透理论模型，则由于液相扩散系数DG和液体微元暴露时间tc(s)可以分别建立相应理论计算模型，故被广泛采用。
所谓构效调控数学模型，是指将超细气液颗粒反应体系之反应效率(能效和物效)与体系理化特性、微界面特性、传质特性和反应器结构用数学方法关联起来，从而实现可通过调整结构参数和操作参数以获得反应过程能效物效的最大化目标，或者在给定反应目标(任务)和能耗物耗下，设计出高效的反应器结构。而对于微界面强化反应器而言，此方面工作在国际上仍为空白。
参考文献
[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]Akita K,Yoshida F.Bubble size,interfacial area,and liquid-phasemass-transfer coefficient in bubble columns[J].Industrial&EngineeringChemistry Process Design and Dvelopment.1974,13(1):84-91.
[16]Schumpe A,Deckwer W D.Organic liquids in a bubble column:holdupsand mass transfer coefficients[J].AIChE journal.1987,33(9):1473-1480.
[17]Jordan U,Schumpe A.The gas density effect on mass transfer inbubble columns with organic liquids[J].Chemical Engineering Science.2001,56(21):6267-6272.
[18]Shah Y T,Kelkar B G,Godbole S P,et al.Design parametersestimations for bubble column reactors[J].AIChE Journal.1982,28(3):353-379.
[19]Sada E,Kumazawa H,Lee C,et al.Gas holdup and mass-transfercharacteristics in a three-phase bubble column[J].Industrial&EngineeringChemistry Process Design and Development.1986,25(2):472-476.
[20]Mandal A,Kundu G,Mukherjee D.Interfacial Area and Liquid‐SideVolumetric Mass Transfer Coefficient in a Downflow Bubble Column[J].TheCanadian Journal of Chemical Engineering.2003,81(2):212-219.
[21]Kawase Y,Halard B,Moo-Young M.Theoretical prediction ofvolumetric mass transfer coefficients in bubble columns for Newtonian andnon-Newtonian fluids[J].Chemical Engineering Science.1987,42(7):1609-1617.
[22]Vázquez G,Cancela M A,Riverol C,et al.Application of theDanckwerts method in a bubble column:Effects of surfactants on mass transfercoefficient and interfacial area[J].Chemical Engineering Journal.2000,78(1):13-19.
[23]Velan M,Ramanujam T K.Gas-liquid mass transfer in a down flow jetloop reactor[J].Chemical engineering science.1992,47(9):2871-2876.
[24]Bouaifi M,Hebrard G,Bastoul D,et al.A comparative study of gashold-up,bubble size,interfacial area and mass transfer coefficients instirred gas-liquid reactors and bubble columns[J].Chemical Engineering andProcessing.2001,40(2SI):97-111.
[25]Kantarci N,Borak F,Ulgen K O.Bubble column reactors[J].ProcessBiochemistry.2005,40(7):2263-2283.
[26]Martín M,Montes F J,Galán M A.Mass transfer from oscillatingbubbles in bubble column reactors[J].Chemical Engineering Journal.2009,151(1-3):79-88.
[27]Charpentier,J.-C.,Mass Transfer Rates in Gas-LiquidAbsorbers.Advances in chemical engineering,1981.11.
[28]Lewis W K,Whitman W G.Principles of gas absorption[J].INDUSTRIALAND ENGINEERING CHEMISTRY.1924,16:1215-1220.
[29]Higbie R.The rate of absorption of a pure gas into still liquidduring short periods of exposure[M].1935.
[30]Danckwerts P V.SIGNIFICANCE OF LIQUID-FILM COEFFICIENTS IN GASABSORPTION[J].INDUSTRIAL AND ENGINEERING CHEMISTRY.1951,43(6):1460-1467.
[31]Kishinevsky M K.Two approaches to the theoretical analysis ofabsorption processes[J].J.Appl.Chem.USSR.1955,28:881-886.
[32]Kishinevsky M K,Serebryansky V T.The mechanism of mass transferat the gas-liquid interface with vigorous stirring[J].J.appl.Chem.USSR.1956,29:29-33.
[33]Dobbins W E.Mechanism of gas absorption by turbulent liquids[M].New York University,Department of Civil Engineering,1964.
[34]Perlmutter D D.Surface-renewal models in mass transfer[J].Chemical Engineering Science.1961,16(3):287-296.
[35]Kulkarni A A.Mass Transfer in Bubble Column Reactors:Effect ofBubble Size Distribution[J].Industrial&Engineering Chemistry Research.2007,46(7):2205-2211.
[36]Calderba.P H.REVIEW SERIES NO 3-GAS ABSORPTION FROM BUBBLES[J].TRANSACTIONS OF THE INSTITUTION OF CHEMICAL ENGINEERS AND THE CHEMICALENGINEER.1967,45(8):C209.
[37]Levich V G,Technica S.Physicochemical hydrodynamics[M].Prentice-hall Englewood Cliffs,NJ,1962.
[38]Davies J T.Turbulent phenomena[J].Academic,New York.1972.
[39]Luk S,Lee Y H.Mass transfer in eddies close to air‐waterinterface[J].AICHE journal.1986,32(9):1546-1554.
[40]Jajuee B,Margaritis A,Karamanev D,et al.Application of surface-renewal-stretch model for interface mass transfer[J].Chemical engineeringscience.2006,61(12):3917-3929.
[41]Vasconcelos J,Rodrigues J,Orvalho S,et al.Effect of contaminantson mass transfer coefficients in bubble column and airlift contactors[J].CHEMICAL ENGINEERING SCIENCE.2003,58(8):1431-1440.
[42]Colombet D,Legendre D,Cockx A,et al.Experimental study of masstransfer in a dense bubble swarm[J].Chemical Engineering Science.2011,66(14):3432-3440.
发明内容
本发明的目的在于，克服现有技术的缺陷，提供微界面强化反应器传质速率构效调控模型建模方法。
为实现上述技术目的，本发明采用如下技术方案：
微界面强化反应器微界面强化反应器一种微界面强化反应器传质速率构效调控模型建模方法，包括：
步骤10：建立气侧传质系数kG计算模型；
根据膜理论，气侧传质系数kG大小由气相扩散系数DG及有效气膜厚度δG决定；
基于Chapman-Enskog动理论，DG的一般形式为：
式中，T为气体温度；MA、MB分别为气体A及溶剂B的摩尔质量；PG为反应器内气泡内气体平均压力；
采用Gedde静止球模型建立kG理论计算模型：
式中，d32为反应器内气泡Sauter平均直径；t32为直径为d32的气泡在体系中的停留时间，其表达式为：
式中，H0为反应器内初始液位高度；v32为与d32对应的反应器内气泡平均上升速度；vG为反应器内的表观气速，QG为反应器内通气体积流量；
将式(3)代入上式(2)后得：
步骤20：建立液侧传质系数计算模型；
建立基于表面更新理论的渗透模型，具体如下：
式中tc为液体微元暴露时间；DL为液相扩散系数；
依据速度滑移模型的定义，将tc利用式(6)定义：
tc＝d32/vS (6)
式中vs为气泡滑移速度；
将式(6)代入式(5)得：
根据气泡滑移速度vs的定义：
vS＝v32-u (8)
式中u为气泡间流体速度；
当液体循环流量为QL，对于均匀上升气流，u的表达式为：
式中，vL为表观液速，S0为反应器横截面积，D0为反应器直径；φG为反应器内气含率，假定为φG＝vG/v32；
将式(9)代入式(8)得：
将式(10)代入式(7)得：
式中，液相扩散系数DL基于修正的Wilke–Chang方程获取，如式(12)所示：
基于式(4)，式(11)，式(12)计算获取微界面强化反应器传质速率。
其中，所述反应器内气泡Sauter平均直径d32采用如下方式获取：
(1)以微界面强化反应器最大气泡直径dmax和最小气泡直径dmin为自变量，气泡Sauter平均直径d32为因变量，建立dmax、dmin和d32间的关系；具体步骤如下：
设x,m,n分别为反应器气液体系中的气泡粒径、气泡粒径几何自然对数的均值和标准差，获取气泡粒径x的概率密度函数：
气泡粒径满足此分布时的气泡Sauter平均直径d32为：
d32＝exp(m+2.5n2) (14)
气泡粒径x呈对数正态分布，因此lnx的数学期望(算术平均值)为：根据气泡粒径x的概率密度函数绘制气泡粒径概率密度图，当时，概率密度最大；即此处的一阶导数为0：
将方程(15)代入(13)得到方程(16)：
由(15)、(16)可得：
由于：
将方程(13)代入(18)并化简后可得：
令：则上式简化为：
方程(20)左端为误差函数，与标准误差函数的差别在于积分限的不同，将式(17)分别代入上述积分上下限，并将方程(20)转化为标准误差函数后可得：
方程(21)中，erf(·)为误差函数；
对于如下形式的误差函数：
采用泰勒级数展开进行近似计算，泰勒级数展开依据误差函数自变量的取值范围不同而采用不同的形式，当z≤4时，erf(z)可展开为：
由于：
当dmax/dmin为1000时：
而根据方程(23)：
因此，当：
即：
时，方程(21)近似成立；
此外，方程(21)成立的条件与n及dmax/dmin的大小有关，且n受dmax/dmin的大小的制约；构建气泡粒径累积概率密度g(n)以考察n和dmax/dmin对方程(21)成立条件的影响，令气泡粒径累积概率密度g(n)为：
绘制g(n)～n关系曲线；获取确保方程(21)成立的n的可取值范围与dmax/dmin的关系；
取不等式(28)的等号条件，即：
由式(17)及(30)确定m和n，进而由方程(14)建立d32基本数学模型；其结果如下：
(2)基于Kolmogorov-Hinze理论，构建微界面强化反应器最大气泡直径dmax、最小气泡直径dmin与反应器参数间的关系；
能使气泡破裂的最小湍流涡尺度是Kolmogorov尺度的11.4～31.4倍，假设此倍率为11.4，由于湍流涡仅能破碎直径大于其尺度的气泡，因此，气泡直径最小值dmin与该湍流涡尺度一致，即：
dmin＝11.4(μL/ρL)0.75ε-0.25 (33)
基于Kolmogorov-Hinze理论，最大气泡直径dmax由下式(34)确定：
dmax＝ε-2/5(σLWecrit/2ρL)3/5 (34)
其中，ε为能量耗散率；μL为液体动力粘度；σL为液体表面张力；ρL为液体密度；Wecrit为气泡破碎临界weber数；
基于气泡破碎的共振理论确定Wecrit：
其中，α2为气泡体积模量，α2＝2,3,...；当α2越大，气泡高阶振动越激烈，气泡就越小，对于超细气泡颗粒选择α2＝2，即Wecrit＝1.24；
此时：
dmax＝0.75(σL/ρL)0.6ε-0.4 (36)
将式(33)、式(36)代入式(32)求解d32。
进一步的，所述能量耗散率ε采用如下方式获取：
步骤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 (39)
步骤121：计算气泡破碎区入口液体静压P0及出口气液混合物压力P1：
忽略气泡破碎区管壁摩擦损失，则：
其中，φmix为气泡破碎区气含率，按下式计算：
忽略管道摩擦及喷嘴处能量损失，依据能量守恒原理，体系实际获得的能量E0为：
即：
由式(40)、(43)得：
步骤122：计算气泡破碎区长度Lmix：
Lmix通过测量气液破碎区管内壁压力突变来确定，或通过如下方式确定：
其中：PH为气液破碎区上方气压，Pa；ρMZ为气液破碎区内气液混合物密度，Kg/m3；vN为射流口的射流速度，m/s；Ue,max为气液破碎区涡旋最大返回速度，m/s；
PH由伯努利方程推得：
PH≈PG0 (46)
式中，PG0为供气压力，Pa；
ρMZ通过下式计算：
ρMZ＝ρGφmix+ρL(1-φmix)≈ρL(1-φmix) (47)
式中，ρG为气体密度，g/m3；
考虑气液破碎区气液混合物流速的影响，Ue,max为射流口射流流速与气液破碎区气液混合物流速的矢量合成的结果，采用下式计算Ue,max：
将式(46)(48)代入式(45)，并化简后可得：
获取反应器气泡破碎器长度Lb，并根据式(49)计算Lmix；
①当Lmix＜Lb时，以式(49)的计算结果为Lmix的实际数值；
②当Lmix≥Lb时，说明射流能量近似完全在气泡破碎器区域耗，则：
Lmix＝Lb (50)
步骤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；
由式(55)可知：则u1增大时，反应器横截面积S0也增大；结合式(37)可知，此时εR减小；为对反应器不同区域的能量耗散率进行比较，假设：u1＝3.0m/s；CL＝20；H0＝1.5m；由式(55)可得，当λ1＝0.1～0.5时：
D0≈19D1 (56)
选定D1数值，计算并比较反应器不同区域在不同喷嘴液速下的能量耗散率，确定与气液破碎区的能量耗散率εmix相比，反应器主体区、气液出口区的能量耗散率可忽略不计，即εmix≈ε；则整个反应器的能量耗散率ε与反应器结构参数之间的数学关系，可由式(38)计算确定，即：
进一步的，所述d32对应的反应器内气泡平均上升速度v32采用如下方式获取：
对于微界面强化反应器中的微米级气泡，考虑升力、气泡周围压力梯度、气泡对其周围湍流场响应时间τb(s)对曳力的影响，微气泡运动的矢量方程可表示为式(58)：
式中，u为气泡间流体速度；y(t)表示t时刻气泡的位置；g为重力加速度；ω为气泡涡量；τb为气泡响应时间；
对于微界面强化反应器的微气泡体系，反应器内为无旋势流场，因此ω＝0；假设气泡进入体系后宏观上为匀速上升，并且稳态时，反应器内湍流各向同性，气泡稳定上升，则Du/Dt＝0；dv32/dt＝0；气泡响应时间τb表示为：
式中，v0为无限大静止液相中单个气泡上升速度；
仅考虑反应器内气泡一维运动情况并结合以上假设，将式(58)简化：
v32＝v0±u (60)
分析反应器内气泡周围液体的三种运动情况：
a.气泡周围液体宏观竖直向下运动；
此时：
v32＝v0-u (61)
当体系中气泡较小时，基于反应器工业生产需求，须满足：
v0≥u (62)
u采用式(9)表示；
结合式(9)(61)，有：
且：
对式(63)化简后得：
v32 2+(vL-vG-v0)v32+vGv0＝0 (65)
求解式(65)，获得唯一实解：
b.气泡周围液体宏观运动速度为零；
根据微界面强化反应器中微气泡主导的反应体系，将Fan通过实验拟合获取的v0的表达式化简得：
式中，σL为液体表面张力；μL为液体动力粘度；ρL为液体密度；
c.气泡周围液体宏观竖直向上运动：
此时：
v32＝v0+u (68)
并且：
v32 2-(vL+vG+v0)v32+vGv0＝0 (69)
整理后得：
根据所述a、b、c，气泡上升速度v32依据液体宏观流动方向表示为：
式中，vL前的“+”表示气泡周围液体宏观竖直向上运动的情况，“-”则表示气泡周围液体宏观竖直向下运动的情况；上式中v0依据式(67)计算。
本发明的另一目的在于提供上述方法构建的反应器传质速率构效调控模型，所述传质速率构效调控模型如下：
dmin＝11.4(μL/ρL)0.75ε-0.25 (33)
dmax＝0.75(σL/ρL)0.6ε-0.4 (36)
式中，d32为反应器内气泡Sauter平均直径；dmin、dmax分别为反应器中最大气泡和最小气泡直径；v32为与d32对应的反应器内气泡平均上升速度；vG为反应器内的表观气速；vL为反应器内的表观液速；v0：无限大静止液相中单个气泡上升速度；S0：反应器横截面积；D0：反应器直径；QL为反应器内液体循环体积流量；Lmix为气泡破碎区长度；D1为气泡破碎管直径；λ1为气液体积流量之比，λ1＝QG/QL；QG为反应器内通气体积流量；P0为气泡破碎器入口处液体的静压；P1为气泡破碎区出口气液混合物压力；ε为能量耗散率；μL为液体动力粘度；σL为液体表面张力；ρL为液体密度。
本发明的又一目的在于提供上述方法在反应器设计中的应用。
本发明的方法适用于微界面强化反应器，其核心在于气泡破碎器。气泡破碎器的原理是高速射流所携带的气体相互撞击进行能量传递，使气泡破碎，其结构参数有Lb、D1，详细结构见附图1，除此之外该反应器的其它结构参数有D0、H0，具体反应器结构相关内容已经公布于发明人在先申请的专利CN10618766A中，本发明中不再赘述。
本发明公开了微界面强化反应器(MIR)构效调控模型的构建方法，所谓构效调控数学模型，是指将超细气液颗粒反应体系之反应效率(能效和物效)与体系理化特性、微界面特性、传质特性和反应器结构用数学方法关联起来，从而实现可通过调整结构参数和操作参数以获得反应过程能效物效的最大化目标，或者在给定反应目标(任务)和能耗物耗下，设计出高效的反应器结构。而对于微界面强化反应器而言，此方面工作在国际上仍为空白。
采用本发明的方法构建的微界面强化反应器传质速率构效调控数学模型有如下有益效果：
(1)可以直观的看出传质速率和气泡大小的关系，不仅填补了国际上此方面的空白，为研究微界面体系奠定了理论基础，从而也可通过调整结构参数和操作参数以实现获得反应过程能效物效的最大化目标，或者在给定反应目标(任务)和能耗物耗下，设计出高效的反应器结构。
(2)本发明的方法构建了dmax、dmin和d32的直接计算关系，而不再采用实验拟合的方式获取d32的具体数值，极大减少了在反应器中应用时产生的误差；现有的方法构建的d32模型多针对的是鼓泡反应器(Bubble column,BC)和鼓泡搅拌釜反应器(Bubbling-stirredreactor,BSR)内的气液体系，或射流泵(gas-liquid jet bump,GLJB)内的空气-水体系。而对于工业微界面反应器(MIR)，则不一定适用，其原因在于：①MIR内气泡破碎机制有别于上述反应器；②工业气液反应体系中可能涉及到高粘度液体，而现有技术并没有考虑液体粘度对d32的影响；而本发明的方法适用于工业微界面反应器(MIR)，其通用性更好。
(3)本发明的方法还构建了能量耗散率ε的表达式，ε是决定气泡破裂机制的关键参数。而ε的增大，意味着气泡更容易破裂，即可以得到尺度更小的气泡，但是针对微米级尺度的能量耗散率的表达式，目前相关研究还是空白，本发明的方法通过严谨的推导过程实现了微界面强化反应器能量耗散率ε的数学模型的构建，为指导新型反应器的设计奠定理论基础。
附图说明
图1是一种反应器结构示意图，用于说明本发明建模方法在反应器装置中的应用。
其中1-反应器，2-泵前阀，3-循环泵，4-泵后阀，5-液体流量计，6-换热器，7-气泡破碎器，8-测温仪，9-下降管，10-气体阀，11-气体流量计，12-气相入口，13-压力表，14-液位计；D0-反应器直径，H0-反应器内初始液位高度，D1-气泡破碎管直径，Lb-气泡破碎区长度。
具体实施方式
实施例1
本发明的方法具体包括如下步骤：
步骤10：建立气侧传质系数kG计算模型；
根据膜理论，气侧传质系数kG大小由气相扩散系数DG及有效气膜厚度δG决定；
基于Chapman-Enskog动理论，DG的一般形式为：
式中，T为气体温度；MA、MB分别为气体A及溶剂B的摩尔质量，Kg/mol；PG为反应器内气泡内气体平均压力，Pa；
对于中高温气体，上式准确性较高，一般误差在5～10％以内。由式可知，当气体组成一定时，DG随温度升高而增大，随气体分压的增大而减小。低密度气体的DG几乎与气体组成无关，一般为5×10-6～1×10-5m2·s-1，较液相扩散系数DL(10-10～10-9m2·s-1)高几个数量级。
δG通常约定平行于界面的速度分量与相应外流速度相差1％的位置为其外边界。但影响δG的因素较多，且不易测量，故依据膜理论建立kG理论计算模型较困难。实际可采用Gedde静止球模型，其数学形式如下：
上式中，d32为反应器内气泡Sauter平均直径；t32为直径为d32的气泡在体系中的停留时间，s。其表达式为：
式中，H0为反应器内初始液位高度；v32为与d32对应的反应器内气泡平均上升速度；vG为反应器内的表观气速，QG为反应器内通气体积流量；
将式(3)代入上式(2)后得：
步骤20：建立液侧传质系数计算模型；
kL体现溶于液相中的气体分子在气液边界层内的传递速率大小，本质上，经典理论均基于边界层假设，不同之处在于气体分子的传递机理有差异。双膜理论认为在静止边界层内，气体分子是通过分子扩散传递，而渗透理论及表面更新理论则多考虑到边界层外流场的影响，认为是液体微元(或湍流涡)更新边界层的结果。Prandtl边界层理论将它们进行了统一。该理论认为，气液边界层实际由静止的粘滞层和对流层组成，质量传递过程由分子扩散和对流传质机理控制。Danckwerts也指出，kL由分子扩散速率和对流传质速率两部分组成。但通常由于气液体系中气泡易变形且存在尾流，其边界层流动实际十分复杂，故基于Prandtl边界层理论建立kL精确数学模型仍然十分困难，必须进行合理的简化。
本发明中建立基于表面更新理论的渗透模型，具体如下：
式中tc为液体微元暴露时间；DL为液相扩散系数；
依据速度滑移模型的定义，考虑到气泡周围液体的绕流，表面更新可近似认为是气泡表面原有液体被新鲜液体替换的过程。进一步而言，气泡(其运动方向上的投影直径为d32)相对于其周围液相滑移(滑移速度vS,m.s-1)是气泡表面更新的原因。因此将tc利用式(6)定义：
tc＝d32/vS (6)
将式(6)代入式(5)得：
根据气泡滑移速度vs的定义：
vS＝v32-u (8)
式中u为气泡间流体速度；
当液体循环流量为QL，对于均匀上升气流，u的表达式为：
式中，vL为表观液速，S0为反应器横截面积，D0为反应器直径；φG为反应器内气含率，假定为φG＝vG/v32；
将式(9)代入式(8)得：
将式(10)代入式(7)得：
式中，液相扩散系数DL基于修正的Wilke–Chang方程获取，如式(12)所示：
基于式(4)，式(11)，式(12)计算获取微界面强化反应器传质速率。
现有d32的具体数值采用实验拟合的方式获取，在反应器中应用时可能产生很大的误差，本发明设计的d32算法构建了dmax、dmin和d32的直接计算关系，极大地减小了在反应器中应用时的误差。
具体的，其计算步骤包括：
(1)以微界面强化反应器最大气泡直径dmax和最小气泡直径dmin为自变量，气泡Sauter平均直径d32为因变量，建立dmax、dmin和d32间的关系；具体步骤如下：
设x,m,n分别为反应器气液体系中的气泡粒径、气泡粒径几何自然对数的均值和标准差，获取气泡粒径x的概率密度函数：
气泡粒径满足此分布时的气泡Sauter平均直径d32为：
d32＝exp(m+2.5n2) (14)
气泡粒径x呈对数正态分布，因此lnx的数学期望(算术平均值)为：根据气泡粒径x的概率密度函数绘制气泡粒径概率密度图，当时，概率密度最大；即此处的一阶导数为0：
将方程(15)代入(13)得到方程(16)：
由(15)、(16)可得：
由于：
将方程(13)代入(18)并化简后可得：
令：则上式简化为：
方程(20)左端为误差函数，与标准误差函数的差别在于积分限的不同，将式(17)分别代入上述积分上下限，并将方程(20)转化为标准误差函数后可得：
方程(21)中，erf(·)为误差函数；
对于如下形式的误差函数：
采用泰勒级数展开进行近似计算，泰勒级数展开依据误差函数自变量的取值范围不同而采用不同的形式，当z≤4时，erf(z)可展开为：
由于：
当dmax/dmin为1000时：
而根据方程(23)：
因此，当：
即：
时，方程(21)近似成立；
此外，方程(21)成立的条件与n及dmax/dmin的大小有关，且n受dmax/dmin的大小的制约；构建气泡粒径累积概率密度g(n)以考察n和dmax/dmin对方程(21)成立条件的影响，令气泡粒径累积概率密度g(n)为：
绘制g(n)～n关系曲线；获取确保方程(21)成立的n的可取值范围与dmax/dmin的关系；
取不等式(28)的等号条件，即：
由式(17)及(30)确定m和n，进而由方程(14)建立d32基本数学模型；其结果如下：
(2)基于Kolmogorov-Hinze理论，构建微界面强化反应器最大气泡直径dmax、最小气泡直径dmin与反应器参数间的关系；
能使气泡破裂的最小湍流涡尺度是Kolmogorov尺度的11.4～31.4倍，假设此倍率为11.4，由于湍流涡仅能破碎直径大于其尺度的气泡，因此，气泡直径最小值dmin与该湍流涡尺度一致，即：
dmin＝11.4(μL/ρL)0.75ε-0.25 (33)
基于Kolmogorov-Hinze理论，最大气泡直径dmax由下式(34)确定：
dmax＝ε-2/5(σLWecrit/2ρL)3/5 (34)
其中，ε为能量耗散率；μL为液体动力粘度；σL为液体表面张力；ρL为液体密度；Wecrit为气泡破碎临界weber数；
基于气泡破碎的共振理论确定Wecrit：
其中，α2为气泡体积模量，α2＝2,3,...；当α2越大，气泡高阶振动越激烈，气泡就越小，对于超细气泡颗粒选择α2＝2，即Wecrit＝1.24；
此时：
dmax＝0.75(σL/ρL)0.6ε-0.4 (36)
将式(33)、式(36)代入式(32)求解d32。
其中，所述能量耗散率ε是决定气泡破裂机制的关键参数。本发明基于微界面强化反应器，设计了更适用于微气泡体系的ε计算方法。
所述能量耗散率ε采用如下方式获取：
步骤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 (39)
步骤121：计算气泡破碎区入口液体静压P0及出口气液混合物压力P1：
忽略气泡破碎区管壁摩擦损失，则：
其中，φmix为气泡破碎区气含率，按下式计算：
忽略管道摩擦及喷嘴处能量损失，依据能量守恒原理，体系实际获得的能量E0为：
即：
由式(40)、(43)得：
步骤122：计算气泡破碎区长度Lmix：
Lmix通过测量气液破碎区管内壁压力突变来确定，或通过如下方式确定：
其中：PH为气液破碎区上方气压，Pa；ρMZ为气液破碎区内气液混合物密度，Kg/m3；vN为射流口的射流速度，m/s；Ue,max为气液破碎区涡旋最大返回速度，m/s；
PH由伯努利方程推得：
PH≈PG0 (46)
式中，PG0为供气压力，Pa；
ρMZ通过下式计算：
ρMZ＝ρGφmix+ρL(1-φmix)≈ρL(1-φmix) (47)
式中，ρG为气体密度，g/m3；
考虑气液破碎区气液混合物流速的影响，Ue,max为射流口射流流速与气液破碎区气液混合物流速的矢量合成的结果，采用下式计算Ue,max：
将式(46)(48)代入式(45)，并化简后可得：
获取反应器气泡破碎器长度Lb，并根据式(49)计算Lmix；
①当Lmix＜Lb时，以式(49)的计算结果为Lmix的实际数值；
②当Lmix≥Lb时，说明射流能量近似完全在气泡破碎器区域耗，则：
Lmix＝Lb (50)
步骤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；
由式(55)可知：则u1增大时，反应器横截面积S0也增大；结合式(37)可知，此时εR减小；为对反应器不同区域的能量耗散率进行比较，假设：u1＝3.0m/s；CL＝20；H0＝1.5m；由式(55)可得，当λ1＝0.1～0.5时：
D0≈19D1 (56)
选定D1数值，计算并比较反应器不同区域在不同喷嘴液速下的能量耗散率，确定与气液破碎区的能量耗散率εmix相比，反应器主体区、气液出口区的能量耗散率可忽略不计，即εmix≈ε；则整个反应器的能量耗散率ε与反应器结构参数之间的数学关系，可由式(38)计算确定，即：
本实施例中，所述d32对应的反应器内气泡平均上升速度v32采用如下方式获取：
对于微界面强化反应器中的微米级气泡，考虑升力、气泡周围压力梯度、气泡对其周围湍流场响应时间τb(s)对曳力的影响，微气泡运动的矢量方程可表示为式(58)：
式中，u为气泡间流体速度；y(t)表示t时刻气泡的位置；g为重力加速度；ω为气泡涡量；τb为气泡响应时间；
对于微界面强化反应器的微气泡体系，反应器内为无旋势流场，因此ω＝0；假设气泡进入体系后宏观上为匀速上升，并且稳态时，反应器内湍流各向同性，气泡稳定上升，则Du/Dt＝0；dv32/dt＝0；气泡响应时间τb表示为：
式中，v0为无限大静止液相中单个气泡上升速度；
仅考虑反应器内气泡一维运动情况并结合以上假设，将式(58)简化：
v32＝v0±u (60)
分析反应器内气泡周围液体的三种运动情况：
a.气泡周围液体宏观竖直向下运动；
此时：
v32＝v0-u (61)
当体系中气泡较小时，基于反应器工业生产需求，须满足：
v0≥u (62)
u采用式(9)表示；
结合式(9)(61)，有：
且：
对式(63)化简后得：
v32 2+(vL-vG-v0)v32+vGv0＝0 (65)
求解式(65)，获得唯一实解：
b.气泡周围液体宏观运动速度为零；
根据微界面强化反应器中微气泡主导的反应体系，将Fan通过实验拟合获取的v0的表达式化简得：
式中，σL为液体表面张力；μL为液体动力粘度；ρL为液体密度；
c.气泡周围液体宏观竖直向上运动：
此时：
v32＝v0+u (68)
并且：
v32 2-(vL+vG+v0)v32+vGv0＝0 (69)
整理后得：
根据所述a、b、c，气泡上升速度v32依据液体宏观流动方向表示为：
式中，vL前的“+”表示气泡周围液体宏观竖直向上运动的情况，“-”则表示气泡周围液体宏观竖直向下运动的情况；上式中v0依据式(67)计算。
实施例2
本实施例以图1所示的反应器为例，说明实施例1所述建模方法构建的模型在二氧化碳和水体系反应器中的应用。图1的反应器结构可为现有反应器的结构，仅采用本发明的方法进行参数设计，本发明中对反应器的结构不再赘述。
根据实施例1构建的传质速率构效调控模型如下：
dmin＝11.4(μL/ρL)0.75ε-0.25 (33)
dmax＝0.75(σL/ρL)0.6ε-0.4 (36)
式中，d32为反应器内气泡Sauter平均直径，m；dmin、dmax分别为反应器中最大气泡和最小气泡直径，m；v32为与d32对应的反应器内气泡平均上升速度，m/s；vG为反应器内的表观气速,m/s；vL为反应器内的表观液速,m/s；v0：无限大静止液相中单个气泡上升速度，m/s；S0：反应器横截面积，m2；D0：反应器直径，m；QL为反应器内液体循环体积流量；Lmix为气泡破碎区长度；D1为气泡破碎管直径；λ1为气液体积流量之比，λ1＝QG/QL；QG为反应器内通气体积流量；P0为气泡破碎器入口处液体的静压；P1为气泡破碎区出口气液混合物压力；ε为能量耗散率；μL为液体动力粘度；σL为液体表面张力；ρL为液体密度。
实施例3选用的模型主要考虑Lmix小于Lb的情况，因为相反的情况不常见，比较极端。反应器结构参数还需满足：λ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时的气侧和液侧的传质速率分别为kG＝1.78×10-6mol/Pa·m3·s，kL＝5.75×10-4m/s。
Levenspiel认为，多相体系的宏观反应速率可由下式表示：
经化简后的气液反应宏观速率方程可简化为：
表1和表2是在MIR中相同体系的各个参数的计算情况：
表1在MIR反应器中模型公式计算的参数
表2在MIR反应器中模型公式计算的三种阻力(气膜、液膜、本征)
实施例3：
本实施例以图1所示的反应器为例，说明实施例1所述建模方法构建的模型在空气-水体系反应器中的应用，与现有装置产生的气液传质速率kL和kG相比的结果。
表3和表4是相同体系不同粒径下的各个参数的对比情况：
表3不同粒径下模型公式计算的参数
表4不同粒径下模型公式计算的三种阻力(气膜、液膜、本征)
如表3、表4所示，在相同条件下，通过模型计算，对于传统反应器产生的直径最小为1mm的气泡相比，MIR产生的气泡直径仅为0.1mm，为原来的1/10，通过改变了气液传质系数，最终影响了其反应的宏观反应速率，使其增大了将近4倍。此结果在工业生产中具有指导性的意义。
Claims (6)
1.一种微界面强化反应器传质速率构效调控模型建模方法，其特征在于，包括：
步骤10：建立气侧传质系数kG计算模型；
根据膜理论，气侧传质系数kG大小由气相扩散系数DG及有效气膜厚度δG决定；
基于Chapman-Enskog动理论，DG的一般形式为：
式中，T为气体温度；MA、MB分别为气体A及溶剂B的摩尔质量；PG为反应器内气泡内气体平均压力；
采用Gedde静止球模型建立kG理论计算模型：
式中，d32为反应器内气泡Sauter平均直径；t32为直径为d32的气泡在体系中的停留时间，其表达式为：
式中，H0为反应器内初始液位高度；v32为与d32对应的反应器内气泡平均上升速度；vG为反应器内的表观气速，QG为反应器内通气体积流量；
将式(3)代入上式(2)后得：
步骤20：建立液侧传质系数计算模型；
建立基于表面更新理论的渗透模型，具体如下：
式中tc为液体微元暴露时间；DL为液相扩散系数；
依据速度滑移模型的定义，将tc利用式(6)定义：
tc＝d32/vS (6)
式中vs为气泡滑移速度；
将式(6)代入式(5)得：
根据气泡滑移速度vs的定义：
vS＝v32-u (8)
式中u为气泡间流体速度；
当液体循环流量为QL，对于均匀上升气流，u的表达式为：
式中，vL为表观液速，S0为反应器横截面积，D0为反应器直径；φG为反应器内气含率，假定为φG＝vG/v32；
将式(9)代入式(8)得：
将式(10)代入式(7)得：
式中，液相扩散系数DL基于修正的Wilke–Chang方程获取，如式(12)所示：
基于式(4)，式(11)，式(12)计算获取微界面强化反应器传质速率。
2.根据权利要求1所述的方法，其特征在于，所述反应器内气泡Sauter平均直径d32采用如下方式获取：
(1)以微界面强化反应器最大气泡直径dmax和最小气泡直径dmin为自变量，气泡Sauter平均直径d32为因变量，建立dmax、dmin和d32间的关系；具体步骤如下：
设x,m,n分别为反应器气液体系中的气泡粒径、气泡粒径几何自然对数的均值和标准差，获取气泡粒径x的概率密度函数：
气泡粒径满足此分布时的气泡Sauter平均直径d32为：
d32＝exp(m+2.5n2) (14)
气泡粒径x呈对数正态分布，因此lnx的数学期望(算术平均值)为：根据气泡粒径x的概率密度函数绘制气泡粒径概率密度图，当时，概率密度最大；即此处的一阶导数为0：
将方程(15)代入(13)得到方程(16)：
由(15)、(16)可得：
由于：
将方程(13)代入(18)并化简后可得：
令：则上式简化为：
方程(20)左端为误差函数，与标准误差函数的差别在于积分限的不同，将式(17)分别代入上述积分上下限，并将方程(20)转化为标准误差函数后可得：
方程(21)中，erf(·)为误差函数；
对于如下形式的误差函数：
采用泰勒级数展开进行近似计算，泰勒级数展开依据误差函数自变量的取值范围不同而采用不同的形式，当z≤4时，erf(z)可展开为：
由于：
当dmax/dmin为1000时：
而根据方程(23)：
因此，当：
即：
时，方程(21)近似成立；
此外，方程(21)成立的条件与n及dmax/dmin的大小有关，且n受dmax/dmin的大小的制约；构建气泡粒径累积概率密度g(n)以考察n和dmax/dmin对方程(21)成立条件的影响，令气泡粒径累积概率密度g(n)为：
绘制g(n)～n关系曲线；获取确保方程(21)成立的n的可取值范围与dmax/dmin的关系；
取不等式(28)的等号条件，即：
由式(17)及(30)确定m和n，进而由方程(14)建立d32基本数学模型；其结果如下：
(2)基于Kolmogorov-Hinze理论，构建微界面强化反应器最大气泡直径dmax、最小气泡直径dmin与反应器参数间的关系；
能使气泡破裂的最小湍流涡尺度是Kolmogorov尺度的11.4～31.4倍，假设此倍率为11.4，由于湍流涡仅能破碎直径大于其尺度的气泡，因此，气泡直径最小值dmin与该湍流涡尺度一致，即：
dmin＝11.4(μL/ρL)0.75ε-0.25 (33)
基于Kolmogorov-Hinze理论，最大气泡直径dmax由下式(34)确定：
dmax＝ε-2/5(σLWecrit/2ρL)3/5 (34)
其中，ε为能量耗散率；μL为液体动力粘度；σL为液体表面张力；ρL为液体密度；Wecrit为气泡破碎临界weber数；
基于气泡破碎的共振理论确定Wecrit：
其中，α2为气泡体积模量，α2＝2,3,...；当α2越大，气泡高阶振动越激烈，气泡就越小，对于超细气泡颗粒选择α2＝2，即Wecrit＝1.24；
此时：
dmax＝0.75(σL/ρL)0.6ε-0.4 (36)
将式(33)、式(36)代入式(32)求解d32。
3.根据权利要求2所述的方法，其特征在于，所述能量耗散率ε采用如下方式获取：
步骤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 (39)
步骤121：计算气泡破碎区入口液体静压P0及出口气液混合物压力P1：
忽略气泡破碎区管壁摩擦损失，则：
其中，φmix为气泡破碎区气含率，按下式计算：
忽略管道摩擦及喷嘴处能量损失，依据能量守恒原理，体系实际获得的能量E0为：
即：
由式(40)、(43)得：
步骤122：计算气泡破碎区长度Lmix：
Lmix通过测量气液破碎区管内壁压力突变来确定，或通过如下方式确定：
其中：PH为气液破碎区上方气压，Pa；ρMZ为气液破碎区内气液混合物密度，Kg/m3；vN为射流口的射流速度，m/s；Ue,max为气液破碎区涡旋最大返回速度，m/s；
PH由伯努利方程推得：
PH≈PG0 (46)
式中，PG0为供气压力，Pa；
ρMZ通过下式计算：
ρMZ＝ρGφmix+ρL(1-φmix)≈ρL(1-φmix) (47)
式中，ρG为气体密度，g/m3；
考虑气液破碎区气液混合物流速的影响，Ue,max为射流口射流流速与气液破碎区气液混合物流速的矢量合成的结果，采用下式计算Ue,max：
将式(46)(48)代入式(45)，并化简后可得：
获取反应器气泡破碎器长度Lb，并根据式(49)计算Lmix；
①当Lmix＜Lb时，以式(49)的计算结果为Lmix的实际数值；
②当Lmix≥Lb时，说明射流能量近似完全在气泡破碎器区域耗，则：
Lmix＝Lb (50)
步骤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；
由式(55)可知：则u1增大时，反应器横截面积S0也增大；结合式(37)可知，此时εR减小；为对反应器不同区域的能量耗散率进行比较，假设：u1＝3.0m/s；CL＝20；H0＝1.5m；由式(55)可得，当λ1＝0.1～0.5时：
D0≈19D1 (56)
选定D1数值，计算并比较反应器不同区域在不同喷嘴液速下的能量耗散率，确定与气液破碎区的能量耗散率εmix相比，反应器主体区、气液出口区的能量耗散率可忽略不计，即εmix≈ε；则整个反应器的能量耗散率ε与反应器结构参数之间的数学关系，可由式(38)计算确定，即：
。
4.根据权利要求3所述的方法，其特征在于，所述d32对应的反应器内气泡平均上升速度v32采用如下方式获取：
对于微界面强化反应器中的微米级气泡，考虑升力、气泡周围压力梯度、气泡对其周围湍流场响应时间τb(s)对曳力的影响，微气泡运动的矢量方程可表示为式(58)：
式中，u为气泡间流体速度；y(t)表示t时刻气泡的位置；g为重力加速度；ω为气泡涡量；τb为气泡响应时间；
对于微界面强化反应器的微气泡体系，反应器内为无旋势流场，因此ω＝0；假设气泡进入体系后宏观上为匀速上升，并且稳态时，反应器内湍流各向同性，气泡稳定上升，则Du/Dt＝0；dv32/dt＝0；气泡响应时间τb表示为：
式中，v0为无限大静止液相中单个气泡上升速度；
仅考虑反应器内气泡一维运动情况并结合以上假设，将式(58)简化：
v32＝v0±u (60)
分析反应器内气泡周围液体的三种运动情况：
a.气泡周围液体宏观竖直向下运动；
此时：
v32＝v0-u (61)
当体系中气泡较小时，基于反应器工业生产需求，须满足：
v0≥u (62)
u采用式(9)表示；
结合式(9)(61)，有：
且：
对式(63)化简后得：
v32 2+(vL-vG-v0)v32+vGv0＝0 (65)
求解式(65)，获得唯一实解：
b.气泡周围液体宏观运动速度为零；
根据微界面强化反应器中微气泡主导的反应体系，将Fan通过实验拟合获取的v0的表达式化简得：
式中，σL为液体表面张力；μL为液体动力粘度；ρL为液体密度；
c.气泡周围液体宏观竖直向上运动：
此时：
v32＝v0+u (68)
并且：
v32 2-(vL+vG+v0)v32+vGv0＝0 (69)
整理后得：
根据所述a、b、c，气泡上升速度v32依据液体宏观流动方向表示为：
式中，vL前的“+”表示气泡周围液体宏观竖直向上运动的情况，“-”则表示气泡周围液体宏观竖直向下运动的情况；上式中v0依据式(67)计算。
5.根据权利要求4所述方法构建的反应器传质速率构效调控模型，其特征在于，所述传质速率构效调控模型如下：
dmin＝11.4(μL/ρL)0.75ε-0.25 (33)
dmax＝0.75(σL/ρL)0.6ε-0.4 (36)
式中，d32为反应器内气泡Sauter平均直径；dmin、dmax分别为反应器中最大气泡和最小气泡直径；v32为与d32对应的反应器内气泡平均上升速度；vG为反应器内的表观气速；vL为反应器内的表观液速；v0：无限大静止液相中单个气泡上升速度；S0：反应器横截面积；D0：反应器直径；QL为反应器内液体循环体积流量；Lmix为气泡破碎区长度；D1为气泡破碎管直径；λ1为气液体积流量之比，λ1＝QG/QL；QG为反应器内通气体积流量；P0为气泡破碎器入口处液体的静压；P1为气泡破碎区出口气液混合物压力；ε为能量耗散率；μL为液体动力粘度；σL为液体表面张力；ρL为液体密度。
6.权利要求1-4任一项所述方法在反应器设计中的应用。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710766695.8A CN107346378B (zh) | 2017-08-30 | 2017-08-30 | 微界面强化反应器传质速率构效调控模型建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710766695.8A CN107346378B (zh) | 2017-08-30 | 2017-08-30 | 微界面强化反应器传质速率构效调控模型建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107346378A true CN107346378A (zh) | 2017-11-14 |
| CN107346378B CN107346378B (zh) | 2020-09-04 |
Family
ID=60257624
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710766695.8A Active CN107346378B (zh) | 2017-08-30 | 2017-08-30 | 微界面强化反应器传质速率构效调控模型建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107346378B (zh) |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108563919A (zh) * | 2018-03-19 | 2018-09-21 | 中国石油化工股份有限公司 | 聚合物凝胶颗粒孔隙尺度运移的直接跟踪方法 |
| CN109806815A (zh) * | 2019-01-14 | 2019-05-28 | 浙江大学 | 基于液速波动的射流鼓泡反应器气液分散状态检测方法 |
| CN110245434A (zh) * | 2019-06-18 | 2019-09-17 | 南京延长反应技术研究院有限公司 | Miha纯气动操作条件下氢气传输调控模型建模方法 |
| CN111523194A (zh) * | 2019-02-04 | 2020-08-11 | 南京大学 | Miha纯气动操作条件下反应速率及转化率调控模型建模方法 |
| CN112175656A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接液化的悬浮床强化反应系统及方法 |
| CN112175649A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤与生物质的多级液化系统及工艺 |
| CN112175654A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接液化的沸腾床强化反应系统及工艺 |
| CN112175652A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接液化的乳化床强化反应系统及方法 |
| CN112175651A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接液化循环溶剂的制备系统及方法 |
| CN112175650A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种洗油和煤直接液化油混合加工系统及方法 |
| CN112175655A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接液化的强化反应系统及方法 |
| CN112175646A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接加氢液化富产芳烃的系统及方法 |
| CN112175653A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤液化的强化反应系统及方法 |
| US20210069666A1 (en) * | 2020-03-31 | 2021-03-11 | Nanjing Yanchang Reaction Technology Research Institute Co., Ltd. | Built-in micro interfacial enhanced reaction system and process for pta production with px |
| CN113035285A (zh) * | 2021-04-01 | 2021-06-25 | 南京延长反应技术研究院有限公司 | 微气泡尺寸对油品加氢脱硫效果影响的计算方法 |
| CN113066535A (zh) * | 2021-02-24 | 2021-07-02 | 南京延长反应技术研究院有限公司 | 一种确定微界面反应体系中气泡大小的方法 |
| CN113087255A (zh) * | 2019-12-23 | 2021-07-09 | 南京延长反应技术研究院有限公司 | 一种吡啶废水的处理系统及方法 |
| CN113223627A (zh) * | 2021-04-07 | 2021-08-06 | 南京延长反应技术研究院有限公司 | 一种微界面强化加氢反应的宏观反应速率计算方法 |
| CN113257367A (zh) * | 2021-05-24 | 2021-08-13 | 南京延长反应技术研究院有限公司 | 一种基于苯酚废水的湿式氧化构效调控模型的建模方法 |
| WO2021196385A1 (zh) * | 2020-03-31 | 2021-10-07 | 南京延长反应技术研究院有限公司 | Px生产pta的外置微界面机组强化氧化系统 |
| CN113971988A (zh) * | 2021-11-08 | 2022-01-25 | 南京延长反应技术研究院有限公司 | 一种丙烯氢甲酰化制丁醛微界面反应强化程度的评测方法 |
| WO2023077844A1 (zh) | 2021-11-08 | 2023-05-11 | 南京延长反应技术研究院有限公司 | 一种蒽醌法制双氧水体系中微气泡强化的评测方法 |
| CN116644688A (zh) * | 2023-06-19 | 2023-08-25 | 北京科技大学 | 一种基于破碎模式预测的雾化生产顺行控制方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1979111A (zh) * | 2005-11-30 | 2007-06-13 | 中国科学院过程工程研究所 | 一种测量颗粒流体两相流反应器内流动参数分布的方法 |
| CN101244868A (zh) * | 2008-03-06 | 2008-08-20 | 同济大学 | 一种射流曝气污水处理反应器的优化设计方法 |
| CN102034000A (zh) * | 2010-12-10 | 2011-04-27 | 华东理工大学 | 一种工业装置乙炔催化加氢反应过程优化运行方法 |
| CN106268544A (zh) * | 2016-08-05 | 2017-01-04 | 南京大学 | 塔式超细气泡反应器 |
| WO2017099687A1 (en) * | 2015-12-08 | 2017-06-15 | Turkiye Petrol Rafinerileri Anonim Sirketi Tupras | A control system for modeling and optimizing a diesel hydroprocessing unit |
| CN106933262A (zh) * | 2015-12-29 | 2017-07-07 | 天津城建大学 | 一种在线压力信号解析的多相流场量测及气-液传质增效控制方法 |
-
2017
- 2017-08-30 CN CN201710766695.8A patent/CN107346378B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1979111A (zh) * | 2005-11-30 | 2007-06-13 | 中国科学院过程工程研究所 | 一种测量颗粒流体两相流反应器内流动参数分布的方法 |
| CN101244868A (zh) * | 2008-03-06 | 2008-08-20 | 同济大学 | 一种射流曝气污水处理反应器的优化设计方法 |
| CN102034000A (zh) * | 2010-12-10 | 2011-04-27 | 华东理工大学 | 一种工业装置乙炔催化加氢反应过程优化运行方法 |
| WO2017099687A1 (en) * | 2015-12-08 | 2017-06-15 | Turkiye Petrol Rafinerileri Anonim Sirketi Tupras | A control system for modeling and optimizing a diesel hydroprocessing unit |
| CN106933262A (zh) * | 2015-12-29 | 2017-07-07 | 天津城建大学 | 一种在线压力信号解析的多相流场量测及气-液传质增效控制方法 |
| CN106268544A (zh) * | 2016-08-05 | 2017-01-04 | 南京大学 | 塔式超细气泡反应器 |
Non-Patent Citations (5)
| Title |
|---|
| CHARPENTIER ET AL: "Mass Transfer Rates in Gas-Liquid Absorbers", 《ADVANCES IN CHEMICAL ENGINEERING》 * |
| JAJUEE B.ET AL: "Application of surface-renewal-stretch model for interface mass transfer", 《CHEMICAL ENGINEERING SCIENCE》 * |
| 张志炳等: "气液反应体系相界面传质强化研究", 《化学工程》 * |
| 杨程等: "液相中气泡群粒径分布的检测与分析", 《南京大学学报（自然科学）》 * |
| 王云鹏等: "强化气液喷射搅拌反应器传质特性研究", 《南京大学学报（自然科学）》 * |
Cited By (35)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108563919A (zh) * | 2018-03-19 | 2018-09-21 | 中国石油化工股份有限公司 | 聚合物凝胶颗粒孔隙尺度运移的直接跟踪方法 |
| CN108563919B (zh) * | 2018-03-19 | 2022-04-19 | 中国石油化工股份有限公司 | 聚合物凝胶颗粒孔隙尺度运移的直接跟踪方法 |
| CN109806815A (zh) * | 2019-01-14 | 2019-05-28 | 浙江大学 | 基于液速波动的射流鼓泡反应器气液分散状态检测方法 |
| CN111523194A (zh) * | 2019-02-04 | 2020-08-11 | 南京大学 | Miha纯气动操作条件下反应速率及转化率调控模型建模方法 |
| CN111523194B (zh) * | 2019-02-04 | 2023-11-07 | 南京大学 | Miha纯气动操作条件下反应速率及转化率调控模型建模方法 |
| CN110245434A (zh) * | 2019-06-18 | 2019-09-17 | 南京延长反应技术研究院有限公司 | Miha纯气动操作条件下氢气传输调控模型建模方法 |
| CN110245434B (zh) * | 2019-06-18 | 2023-06-27 | 南京延长反应技术研究院有限公司 | Miha纯气动操作条件下氢气传输调控模型建模方法 |
| CN112175651A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接液化循环溶剂的制备系统及方法 |
| CN112175652A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接液化的乳化床强化反应系统及方法 |
| CN112175650A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种洗油和煤直接液化油混合加工系统及方法 |
| CN112175655A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接液化的强化反应系统及方法 |
| CN112175646A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接加氢液化富产芳烃的系统及方法 |
| CN112175653A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤液化的强化反应系统及方法 |
| CN112175654A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接液化的沸腾床强化反应系统及工艺 |
| CN112175649A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤与生物质的多级液化系统及工艺 |
| CN112175656A (zh) * | 2019-07-04 | 2021-01-05 | 南京延长反应技术研究院有限公司 | 一种煤直接液化的悬浮床强化反应系统及方法 |
| CN113087255A (zh) * | 2019-12-23 | 2021-07-09 | 南京延长反应技术研究院有限公司 | 一种吡啶废水的处理系统及方法 |
| CN113087255B (zh) * | 2019-12-23 | 2022-04-22 | 南京延长反应技术研究院有限公司 | 一种吡啶废水的处理系统及方法 |
| US20220395803A1 (en) * | 2020-03-31 | 2022-12-15 | Nanjing Yanchang Reaction Technology Research Institute Co., Ltd. | Strengthening oxidation system of external micro-interfacial unit for producing pta with px |
| US20210069666A1 (en) * | 2020-03-31 | 2021-03-11 | Nanjing Yanchang Reaction Technology Research Institute Co., Ltd. | Built-in micro interfacial enhanced reaction system and process for pta production with px |
| WO2021196385A1 (zh) * | 2020-03-31 | 2021-10-07 | 南京延长反应技术研究院有限公司 | Px生产pta的外置微界面机组强化氧化系统 |
| US11628415B2 (en) * | 2020-03-31 | 2023-04-18 | Nanjing Yanchang Reaction Technology Research Institute Co., Ltd. | Built-in micro interfacial enhanced reaction system and process for PTA production with PX |
| US11607663B2 (en) * | 2020-03-31 | 2023-03-21 | Nanjing Yanchang Reaction Technology Research Institute Co., Ltd. | Strengthening oxidation system of external micro-interfacial unit for producing PTA with PX |
| CN113066535A (zh) * | 2021-02-24 | 2021-07-02 | 南京延长反应技术研究院有限公司 | 一种确定微界面反应体系中气泡大小的方法 |
| CN113066535B (zh) * | 2021-02-24 | 2023-03-07 | 南京延长反应技术研究院有限公司 | 一种确定微界面反应体系中气泡大小的方法 |
| CN113035285A (zh) * | 2021-04-01 | 2021-06-25 | 南京延长反应技术研究院有限公司 | 微气泡尺寸对油品加氢脱硫效果影响的计算方法 |
| CN113223627A (zh) * | 2021-04-07 | 2021-08-06 | 南京延长反应技术研究院有限公司 | 一种微界面强化加氢反应的宏观反应速率计算方法 |
| CN113223627B (zh) * | 2021-04-07 | 2023-10-31 | 南京延长反应技术研究院有限公司 | 一种微界面强化加氢反应的宏观反应速率计算方法 |
| CN113257367B (zh) * | 2021-05-24 | 2023-10-27 | 南京延长反应技术研究院有限公司 | 一种基于苯酚废水的湿式氧化构效调控模型的建模方法 |
| CN113257367A (zh) * | 2021-05-24 | 2021-08-13 | 南京延长反应技术研究院有限公司 | 一种基于苯酚废水的湿式氧化构效调控模型的建模方法 |
| CN113971988A (zh) * | 2021-11-08 | 2022-01-25 | 南京延长反应技术研究院有限公司 | 一种丙烯氢甲酰化制丁醛微界面反应强化程度的评测方法 |
| CN113971988B (zh) * | 2021-11-08 | 2023-05-05 | 南京延长反应技术研究院有限公司 | 一种丙烯氢甲酰化制丁醛微界面反应强化程度的评测方法 |
| WO2023077844A1 (zh) | 2021-11-08 | 2023-05-11 | 南京延长反应技术研究院有限公司 | 一种蒽醌法制双氧水体系中微气泡强化的评测方法 |
| CN116644688A (zh) * | 2023-06-19 | 2023-08-25 | 北京科技大学 | 一种基于破碎模式预测的雾化生产顺行控制方法 |
| CN116644688B (zh) * | 2023-06-19 | 2023-12-15 | 北京科技大学 | 一种基于破碎模式预测的雾化生产顺行控制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107346378B (zh) | 2020-09-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107346378A (zh) | 微界面强化反应器传质速率构效调控模型建模方法 | |
| CN107335390B (zh) | 微界面强化反应器相界面积构效调控模型建模方法 | |
| CN107561938A (zh) | 微界面强化反应器反应速率构效调控模型建模方法 | |
| CN107563051B (zh) | 微界面强化反应器气泡尺度构效调控模型建模方法 | |
| Rebrov | Two-phase flow regimes in microchannels | |
| Forrester et al. | Bubble formation from cylindrical, flat and concave sections exposed to a strong liquid cross-flow | |
| Tan et al. | A non-spherical model for bubble formation with liquid cross-flow | |
| Liu et al. | CFD modelling of uneven flows behaviour in flat-sheet membrane bioreactors: From bubble generation to shear stress distribution | |
| Wutz et al. | Predictability of kLa in stirred tank reactors under multiple operating conditions using an Euler–Lagrange approach | |
| Bao et al. | Influence of impeller diameter on overall gas dispersion properties in a sparged multi-impeller stirred tank | |
| Chakravorty | Process intensification by pulsation and vibration in miscible and immiscible two component systems | |
| Li et al. | CFD analysis of CO2 absorption in a microporous tube-in-tube microchannel reactor with a novel gas-liquid mass transfer model | |
| Dong et al. | Effect of liquid phase rheology and gas–liquid interface property on mass transfer characteristics in bubble columns | |
| Roy et al. | CFD simulation of flow and axial dispersion in external loop airlift reactor | |
| Ammar et al. | Effect of the tank design on the flow pattern generated with a pitched blade turbine | |
| Li et al. | Dynamic characteristics of ventilated bubble moving in micro scale venturi | |
| CN107589667B (zh) | 微界面强化反应器能量耗散构效调控模型建模方法 | |
| Sun et al. | Liquid flow and breakage behaviors of two liquid jets impacting on the wire mesh with different impinging angles | |
| Alam et al. | Population balance and computational fluid dynamics modeling of swirl flow microbubble generator | |
| Gourich et al. | Comparison of hydrodynamic and mass transfer performances of an emulsion loop-venturi reactor in cocurrent downflow and upflow configurations | |
| Li et al. | Investigation into drag coefficient for sub-millimeter bubbles in gas–liquid bubbly flow: Experiments and CFD simulations | |
| Sato et al. | Effects of azimuthal angle of aeration hole on flows inside and outside an air diffuser pipe | |
| Alam et al. | Unsteady numerical simulation of gas-liquid flow in dual chamber microbubble generator | |
| Koutita et al. | A theoretical fluid dynamic model for estimation of the hold-up and liquid velocity in an external loop airlift bioreactor | |
| Li et al. | Experimental study of Taylor bubble flow in non-Newtonian liquid in a rectangular microchannel |
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 |