,,,天津工业大学,,中空纤维膜材料与膜过程教育部重点实验室, ,,环境与化学工程学院,天津,,,,,,,,, 中国科学院生态环境研究中心污染控制研究室,北京,,,,,,,,,北京市交通运输环境保护中心,北京,,,,,,,,,天津膜天膜科技股份有限公司,天津,,,,,,,
摘 要,采用,,,,,,,,多相流模型对膜生物反应器进行气液两相流数值模拟,对“对齐”和“对齐导流”两种反应器构型内部流场气含率、速度场和膜面液体流速进行了分析比较,并就布气方式进行了分析和优化,同时借助缩小实验模型对模拟结果进行了实验验证。 结果表明,两种反应器构型流场内的气泡呈现出汇集于膜组件中心位置,然后在膜组件顶端散开的流动状态,导流作用对反应器内气含率分布的影响不大,但对速度场分布特性影响显著,体积缩小,,,倍的对齐导流模拟装置中的气液流动状态与,,,模拟结果基本一致,通过对布气方式的优化模拟发现,不均匀布气方式可以改善气体分布状况,提高反应器内气含率和流场的湍流强度。
为了减缓膜面污染及维持高浓度的污泥生化降解能力,膜生物反应器,,,,,,,,, ,,,,,,,,,,,,,,,运行过程中需要高强度的曝气作用。 因此与其它污水处理工艺相比,,,,运行中需要很高的能耗,其中膜擦洗曝气、生化工艺曝气能耗占系统运行能耗的,,,,,, 。 因此了解各种布气方式下,,,内部流场特性,对于合理设计,,,构型及布气方式以达到节能降耗的目的至关重要。 由于反应器内部流场复杂多变,因此仅通过实验手段获得整体流场分布特性受到一定的局限性。近年来,随着计算流体力学,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,技术及其软件的迅速发展,已在反应器设计、优化中得到了广泛的应用,,,,, ,它的优势在于应用数值软件求解描述过程,可以实现过程设计、优化以及放大。如今已成为解决工程问题的一种新的手段,,,,, 。
,,,内部流场数值模拟属于气液两相流范畴,,,,,, ,气液两相流是常见的一种物质流动,广泛应用于工程和环境领域中,,,,,,, 。在,,,运行过程中,曝气是重要的环节,曝气的作用一方面是向反应器内充空气,以保证好氧微生物所需的溶解氧,另一方面是吹扫膜组件,气泡对膜面的冲刷剪切以减缓膜污染,,,,,,, 。对于,,,构型及布气方式的优化,由于工程实际装置的复杂性及实验条件的局限性,有些问题在实验室难以完成。而计算流体力学模拟可以实现及预见工程放大可能出现的一些问题。 ,,,,,等,,,,运用,,,对应用最广泛的,,,反应池构型进行了二维数值模拟,对其进行布气方式优化,考察了特定曝气强度和布气方式对,,,运行特性的影响。
在上述研究的基础上,运用,,,数值模拟方法对曝气引发的,,,内流场的状态进行了分析,同时比较了曝气管正对齐膜组件下方,简称“对齐”布气,及“对齐”布气时在膜单元两侧加上曝气导流板,简称“对齐”导流布气,构型的改进给,,,内部流场带来的改变,由此提出了布气方式的优化方案,以便为大型,,,工程设计提供相应的理论基础。
由于本研究的重点是曝气引发的反应器内流场状态,且膜组件过滤引发的液体流动速度远小于曝气引发的气、液混合流动速度,因此将膜简化为不透水的无滑移壁面,模拟中不考虑膜过滤。
借助,,,,,, ,, ,, ,,划分反应器网格结构,网格决定着模型的准确度。生成网格不但要考虑其细密程度来提高模拟过程的准确程度,又要避免由于网格过于复杂带来的计算复杂性。该研究采用四面体网格,分别画了网格数为,,,万和,,,万,种,模拟结果相同,证明了网格的无关性,网格整体与局部图如图,。反应器模拟活性污泥混合液,其混合液粘度及密度均为实测值,分别为,, ,,, ,, ,,·,,, · ,,,和, ,,,,,·,,, 。从流体实际流动情况出发,采用非定常流动模式,简单且计算精度较高的标准,,ε两方程湍流模型,该模型在一定程度上考虑了流动场中各点的湍动能传递和流动的历史作用, 。求解器选用求解不可压流场且需要内存较小的压力速度耦和,,,,,,算法,稳态,,分离显式解算器,对流项选取一阶迎风离散格式。设计液面高度为, ,,,,,,曝气孔进气速度为,, ,,, ,,· ,,, 。设置入口边界条件为速度入口边界条件,出口边界条件为压力出口边界条件,过滤器壁面及膜组件壁面为无滑移壁面边界条件,反应器为流场边界条件。
研究以,,,中反应器内的流体为分析对象,以天津膜天膜科技股份有限公司的,,,,Ⅱ,,型帘式膜组件为基础模型。反应器及膜组件结构参数如下,反应器尺寸为, ,,,,,×, ,,,,,×, ,,,,,,膜组件尺寸为,,,,,×,,,,,×, ,,,,,,单个组件膜面积,,,, ,组件间距,,,,,膜组件距反应池边距,,,,,,底部采用穿孔管曝气,曝气孔,,,,每个膜组件下均匀排布十个曝气孔,气孔间距,,,,,排间距,,,,,。 图,为两种布气方式下的三维反应器简化模型。对“对齐导流”布气如图,,,,及“对齐”布气如图,,,,进行三维数值模拟,分析其内部流场特性。 为模拟求解方便,假设反应器内均匀。
模拟过程中将污泥的固相和水的液相视作“假想的均一混合相”,则可将三相流问题简化为“气液两相流”问题。通过对反应器不同布气方式的模拟,结果得到反应器内多个截面的气液两相速度场及速度流线等信息,以此来了解反应器内水力情况。通
过实验室缩小模型与模拟结果的对比,验证模拟结果,并对改进的布气方式进行了模拟计算。
,, ,, , 反应器内气含率分布 图,为模型的简化示意图,实验选取了,,,,,,,,,,,几个截面从横向和纵向分板反应器内流场的变化。 图,和,为反应器不同截面及轴向上的气含率分布情况,从图,可以看出反应器曝气板底部的气含率很低,从曝气板到膜组件底部位置气含率开始显著增加。在曝气器上方的位置气体沿膜组件中间位置比较集中的上升,然后在膜组件顶端散开,这一模拟结果与,,,,,,,模拟结果一致。气体在流场中的这种分布可能是由于从曝气孔出来的小气泡形态比较稳定,对周围气泡的形态及速度变化没有明显的影响,组件之间的阻挡作用使大量的小气泡在膜组件中心位置比较集中的缓慢上升,随着气泡的逐渐上升,到达组件顶部后,气体没有组件的阻挡作用,气泡尺寸越来越大,较大的气泡容易变形破裂或与周围小气泡融合,最终大量大气泡在膜组件顶部破裂散开,从膜组件顶端散开的气体在水流的作用下沿着壁面被带到反应器下部,在反应器内部形成循环流动。从图,也可以看出,气体主要集中在反应器中心位置,反应器底部和反应器两侧气含率较低。
,, ,, , 反应器内液速分布 为进一步分析流体在膜面的流动特性,研究考察了曝气槽中平行膜面以及垂直膜面几个典型截面上流体的速度等值线及矢量场分布情况。 图,为各截面液体速度云图,图,为各截面速度矢量图。速度云图及流线图中的颜色深浅代表速度绝对值的大小,箭头指向代表流体的运动方向。从图,中可以看出,膜单元中间,及膜单元顶部存在大量的高流速区域,而反应器内底部、靠近反应器壁面及膜单元外围的膜组件存在一定范围的水力死区。曝气带动的混合液流动呈现出了不均匀的状态,这会致使膜面污染的不均匀性,外侧膜组件受到曝气清洗的作用较弱,导致污染较严重。
分析图,可知,流体在反应器内呈现为循环运动,漩涡的形成有利于气液两相的充分混合。气泡的轨迹倾向于流向膜组件中心位置,然后在膜组件顶端散开。从平行膜面截面速度矢量图,,,,可以
看出,气泡带动液体在膜组件之间形成的一束高流速流体比较集中的向上流动,到达组件上部时由于大气泡的相互作用形成大量的涡旋,大量涡旋的扰动使流体在组件之间形成内循环流动,这有利于提高膜面的冲刷作用。但膜组件下部两侧存在一定范围的低流速区域,这些区域涡旋不明显,流体扰动不强烈,流体对膜面的冲刷作用较弱。
结合图,和,可以看出,流体在反应器内的流动呈现一定的不均匀性,组件中间的升流区流速较大而两侧的流速较小,这使得膜污染的分布也呈现一种不均匀性。 同时流体在膜单元中间位置形成流速较大的升流区,到达膜组件顶部扩散,流体往反应器两侧向下流动,再与膜单元外围膜组件间隔区域形成循环流动。,, , 导流作用对反应器内部流场的影响
,, ,, , 导流板对气含率分布的影响 图,、图,为导流结构的反应器不同截面及轴向上的气含率分布情况。对比图,和,可以看出导流板的增加,使气体分布有所改变,但变化不大。分析图,数据也可以看出,反应器轴向上气含率分布基本一致,导流布气时气含率略高,但总的来说导流作用对反应器内气含率分布的影响不是很大。
图,, “对齐”导流布气下反应器各截面流体速度分布云图,截面选取同图,,
由图,,可知,加导流板之后,流体在反应器内仍为循环运动,与“对齐”布气不同的是,导流板的设置使流场内的循环流动明显增强。流场内液相质点形成的旋涡非常明显,导流板的设置增强了反应器内流体的循环流动,形成了非常明显的升流区和降流区,这有利于流体浓度及速度的分布均匀。通过图,,中这两种反应器构型不同截面及膜表面平均速度的对比也可以发现,导流结构能提高反应器内部的液体流速,改变流场分布情况,使得膜表面的流体速度显著增加,且分布更均匀,有利于更均匀的对膜组件冲刷与水力清洗。
,, ,, , 导流板对流场液速分布的影响 从图,,可以看出,在导流板的作用下,反应器内速度场分布趋于均匀,高流速区域明显增加,底部流体流速有所提升,没有出现大范围的低流速区。对比图,和,,的反应器内流速分布可以看出,导流板在一定程度上提升了膜单元外侧组件周围的混合液流速,能够使外侧的组件受到更为均匀的水力冲刷作用,对保障膜单元中各组件的污染均匀性具有一定作用。
通过对“对齐”与“对齐导流”这两种反应器构型的模拟发现,均匀排布气孔时,膜单元流场中气泡轨迹流向膜单元中心位置,然后在膜单元顶端散开。集中上升的气泡不利于氧的传质和膜组件的清洗。为了改善反应器内部的气含率分布情况,研究对曝气孔的排布进行了改进。对布气孔位置进行了调整,将,,个曝气孔分成,组,,组,个曝气孔均匀排布,,组之间间距,,,,,。为了求解方便,实验对单个膜组件的小型反应器进行了模拟,实验模型如图,,。其网格划分、计算模型选择与全尺寸模型模拟相同。
匀布气时,气体分成两束上升,分布较广,组件上部区域的速度要小于均匀布气方式,但循环流非常明显。通过对比图,,,,, 、,,,,,可以发现,均匀布气时,流体速度主要集中在膜组件中心及组件上部,没有在反应器两侧形成循环流动,使得反应器底部也几乎没有循环流动,流体在这一区域速度较低,不均匀布气时,在,束气流的作用下,液体在膜组件两侧形成了,个方向相反的涡流。涡流提高了反应器内流场的湍流强度,使得反应器壁面附近流体速度变大,如图,,,,,所示。 由于漩涡将流体带至反应器底部,使底部壁面附近流体产生了一定的流速,这能够有效地改善反应器底部低流速区。
改进曝气形式的内部流场气含率分布如图,,。对比图,可以看出,这一布气方式使反应器内形成,束气流,这,束气流在膜组件上部汇合并散开。不均匀布气方式改善了气体分布状况,使气体分布更均匀,提高了反应器内气含率。
从图,,、,,中对比“对齐”均匀布气和改进布气两种布气方式近膜面处的流体速度及速度矢量分布可知,均匀布气时,流场内高流速区域集中在膜组件中心位置与膜组件上部,膜组件两侧、反应器壁面及底部有大范围的低流速区域,不均匀布气时,流场内速度分布均匀,反应器内膜组件及壁面附近均有较大的流速。反应器底部也有较大的流速。通过对图,,中各特征区域局部放大的流线轨迹,如图,,所示,进行分析,可以对上述,种速度分布做出解释。分析图,,,,, 、,,,,,发现,均匀布气时组件上端有较大的流速,循环流动不明显,这是由于均匀曝气时小气泡集中上升,流体速度主要集中在这一区域,不均
影响反应器内部流场的因素很多,譬如,反应器构型,有无导流板, ,曝气孔的排布方式,均匀和非均匀布置曝气孔, ,曝气孔的大小等。而本研究的重点是从宏观上考察一定曝气强度下反应器不同构型及不同布气方式对反应器内部流场的影响,并且结合小型实验验证数值模拟结果准确性。研究对体积缩小,,,倍的对齐导流模拟装置进行了流动实验,图,,,,,为对齐导流曝气小型实验模型,模型尺寸为,, ,,,,×,, ,,,,×,, ,,,,,曝气孔直径为,, ,,,,,实验各项参数的获得均通过佛洛德数相
似性准则计算得到。为了便于观察反应器内部流场运动特性,选取密度与流体相近的压实海绵颗粒进行示踪,并通过高清相机,,,,,,,,,,,,,,,,,,,连续采集反应器内流体流动状态的图片与数值模拟结果进行对比分析。
从图,,,,,中可看出压实海绵示踪颗粒在反应器内部做旋转运动,这与模拟结果图,,中的的速度矢量图中液体运动方向一致。 图,,,,, 、 ,,,中的观测到从曝气孔出来的大量小气泡在组件及导流板的作用下在膜组件中心位置比较集中的缓慢上升。随着气泡的逐渐上升,气泡尺寸越来越大,大气泡上升速度相对较快并且容易变形破裂或与周围小气泡融合,最终大量大气泡在膜组件顶部破裂散开,从膜组件顶端散开的气体在水流的作用下沿着壁面被带到反应器下部,在反应器内部形成循环流动。从图,,,,,可看出延曝气孔往上流体紊动程度明显增大。这可能是由于组件及导流板对气泡的束缚作用减弱,且较大的气泡对周围气泡的形态及速度变化产生的影响较大,对于其上方的气泡有顶托作用,对位于的气泡其尾涡会加速气泡的上升,而对位于其侧面稍远的气泡则有抑制作用甚至会使其侧面的气泡小幅度下降。而组件底端距曝气口较近的地方,刚从曝气孔出来的大量小气泡受组件及导流板的阻滞作用,上升及扩散速度相对于大气泡较慢,形态比大气泡稳定,但小气泡对周围气泡的形态及速度变化没有明显的影响。气泡轨迹与模拟结果相同,由此说明整个模拟过程中参数、计算模型及计算公式的选择都比较合理,能使模拟结果和实际情况很好的吻合。
,,通过模拟气孔均匀排布的布气方式发现反应器内部气泡轨迹倾向于流向膜组件中心位置,然后在膜组件顶端散开,从膜组件顶端散开的气体在水流的作用下,沿着壁面被带到反应器下部,在整个反应器内形成循环。
,,对比“对齐”与“对齐导流”反应器内流体速度分布场发现,导流板的导流作用,对流场特性影响显著,能很好的改善反应器内部速度分布不均的情况。导流板的设置增加了流场局部区域内的紊流强度,强化了反应器内流体的循环流动,形成了非常明显的升流区和降流区,有利于气体在反应器内的扩散和流体浓度分布均匀。
,,通过对不均匀布气的改进结构的模拟发现,气泡轨迹仍为沿着膜组件中心集中,但改进结构改善了气体分布状况,提高了反应器内气含率和流场的湍流强度,液体在膜组件两侧形成了,个方向相反的涡流,使反应器内流体速度分布更均匀,能改善膜面的不均匀污染。
,,对体积缩小,,,倍的对齐导流小型模型装置进行了流动实验,模拟结果与实验结果基本一致,即模拟过程中参数、计算模型及计算公式的选择都比较合理。888集团官网入口
