CN112836444B - 基于cfd的干式厌氧沼气搅拌系统的优化方法 - Google Patents
基于cfd的干式厌氧沼气搅拌系统的优化方法 Download PDFInfo
- Publication number
- CN112836444B CN112836444B CN202110045514.9A CN202110045514A CN112836444B CN 112836444 B CN112836444 B CN 112836444B CN 202110045514 A CN202110045514 A CN 202110045514A CN 112836444 B CN112836444 B CN 112836444B
- Authority
- CN
- China
- Prior art keywords
- model
- stirring system
- phase
- gas
- calculation
- 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.)
- Active
Links
- 238000003756 stirring Methods 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000005457 optimization Methods 0.000 title claims abstract description 44
- 238000004364 calculation method Methods 0.000 claims abstract description 43
- 239000007789 gas Substances 0.000 claims description 48
- 239000012071 phase Substances 0.000 claims description 45
- 239000007788 liquid Substances 0.000 claims description 22
- 239000007791 liquid phase Substances 0.000 claims description 18
- 239000012530 fluid Substances 0.000 claims description 17
- 239000007787 solid Substances 0.000 claims description 13
- 230000005514 two-phase flow Effects 0.000 claims description 13
- 238000009826 distribution Methods 0.000 claims description 9
- 239000007790 solid phase Substances 0.000 claims description 9
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 8
- 238000011282 treatment Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 5
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 4
- 238000005094 computer simulation Methods 0.000 claims description 2
- 238000007872 degassing Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 20
- 239000000463 material Substances 0.000 abstract description 13
- 238000012423 maintenance Methods 0.000 abstract description 9
- 230000008021 deposition Effects 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 11
- 230000008569 process Effects 0.000 description 11
- 238000004088 simulation Methods 0.000 description 11
- 238000000855 fermentation Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 238000010438 heat treatment Methods 0.000 description 5
- 230000009467 reduction Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 239000007921 spray Substances 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- 238000012800 visualization Methods 0.000 description 4
- 238000004140 cleaning Methods 0.000 description 3
- 230000007774 longterm Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000010802 sludge Substances 0.000 description 3
- 239000002910 solid waste Substances 0.000 description 3
- 238000012546 transfer Methods 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
- 230000003321 amplification Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 239000008358 core component Substances 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 230000029087 digestion Effects 0.000 description 2
- 230000004151 fermentation Effects 0.000 description 2
- 239000010806 kitchen waste Substances 0.000 description 2
- 238000010907 mechanical stirring Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000004064 recycling Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 238000012916 structural analysis Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 150000001335 aliphatic alkanes Chemical class 0.000 description 1
- 230000003851 biochemical process Effects 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000000306 component Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000010791 domestic waste Substances 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000007792 gaseous phase Substances 0.000 description 1
- 244000144972 livestock Species 0.000 description 1
- 239000010871 livestock manure Substances 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000010813 municipal solid waste Substances 0.000 description 1
- 239000013618 particulate matter Substances 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000004381 surface treatment Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E50/00—Technologies for the production of fuel of non-fossil origin
- Y02E50/30—Fuel from waste, e.g. synthetic alcohol or diesel
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Mixers Of The Rotary Stirring Type (AREA)
Abstract
一种基于CFD的干式厌氧沼气搅拌系统的优化方法,包括对待优化干式厌氧沼气搅拌系统现状结构分析并获取核心特征结构,建立搅拌系统模型;将所述搅拌系统模型划分网格;对计算模型简化,设定边界条件后求解;若求解结果不满足优化目的,则修改搅拌系统模型的模型参数,直到求解结果满足优化目,此时的搅拌系统模型即为优化后的搅拌系统。本发明实现了干式厌氧沼气搅拌系统的快速优化设计,优化后的结构能减少沼气搅拌喷嘴周围的低速区,从而减少厌氧物料在此区域沉积,有利于有效反应容积的维持。
Description
技术领域
本发明属于环境工程固体废弃物处置和资源化技术领域,特别涉及一种基于CFD的干式厌氧沼气搅拌系统的优化方法。
背景技术
干式厌氧发酵由于相较于湿式厌氧消化工艺具有处理负荷强,容积产气率高、节约用水量、产沼液量少或不排放沼液、后处理简单以及运行费用低等优点,在欧洲等国外经济发达市场越来越多被用于畜禽废弃物、厨余垃圾等有机固体废弃物的处理,近年来在中国干式厌氧发酵工艺也更多地被用于缓解资源化政策需求及生活垃圾等有机固废存在处理能力缺口的矛盾。
厌氧过程尤其是固体含量较高的干式厌氧发酵过程,对物料的搅拌是厌氧产气效率的关键因素,搅拌直接决定了物料分布的均匀度进而影响厌氧微生物与物料的接触并对生化过程造成影响,同时也对物料的传热、传质等产生正、负反馈。工程实践中发现,对比传统的机械搅拌,干式厌氧发酵更适宜采用气体搅拌的方式,用回流沼气增压后再输入干式厌氧发酵罐,使厌氧物料得到搅拌、扰动,解决了传统机械搅拌高耗能、磨损大、故障率高等问题。然而,目前的沼气搅拌系统结构设计尚有不足,沼气输送管出口喷嘴深入厌氧罐内部、突出罐底表面一定距离的设计,导致突出的管道周围物料不能很好被增加沼气吹起,造成物料在此区域表现出易堆积,长期运行后将缩小罐体有效容积;对于堆积物料的清理维护,也因罐底密布突出表面的尖头喷嘴而不能采用机械车需采用人工进行清罐,大大延长了清罐维护时间,导致工程停机过长,不利于稳定运行和调试启动。因此,有必要针对干式厌氧的沼气搅拌系统进行优化。
计算流体力学(CFD)是将流体流动传热传质与计算机数值计算结合在一起形成的一门新兴学科,它的主要作用是对流场进行数值模拟计算,进而得到整个流场信息。应用CFD,不但基础成本低,而且获得的信息量大同时可以节约大量的试验经费,因而其被广泛的应用到各个领域。近年来,随着计算流体动力学的高速发展,CFD模拟仿真技术已经开始运用到反应器的流场分析及反应器的优化设计的工作之中,并有效地完成了对污泥等不透明流体的流场可视化分析,极大提高了反应器流场研究的效率以及数据分析的全面性。目前,国内外针对干式厌氧过程的反应器流场模拟研究较少,基于CFD对封闭黑箱反应流场可视化的强大优势,可采用合理简化,通过数值模拟及数据分析,实现对结构设计的快速优化。
发明内容
有鉴于此,本发明的主要目的之一在于提出一种基于CFD的干式厌氧沼气搅拌系统的优化方法,以期至少部分地解决上述技术问题中的至少之一。
为了实现上述目的,本发明提供了一种基于CFD的干式厌氧沼气搅拌系统的优化方法,包括:
(1)对待优化干式厌氧沼气搅拌系统现状结构分析并获取核心特征结构,建立搅拌系统模型;
(2)将所述搅拌系统模型划分网格;
(3)对计算模型简化,设定边界条件后求解;
(4)若求解结果不满足优化目的,则修改搅拌系统模型的模型参数,直到求解结果满足优化目,此时的搅拌系统模型即为优化后的搅拌系统。
基于上述技术方案可知,本发明的基于CFD的干式厌氧沼气搅拌系统的优化方法相对于现有技术至少具有以下优势之一:
(1)通过对结构特征的分析,采用对称分析法、降维简化法、突出局部特征法,大大减少几何建模和网格划分的工作,实现快速设计;
(2)采用欧拉单相流或欧拉气液两相流、气固两相流替换气、液、固三相流,大幅降低计算的复杂度,实现快速计算;
(3)对于工程放大和现有工程优化,采用计算流体力学可大幅降低用于实践的小试-中试-工程的建造成本,减少大量时间成本,和实验探索成本,实现优化结果快速可视化,有助于干式厌氧沼气工程的进一步优化;
(4)实现了干式厌氧沼气搅拌系统的快速优化设计,优化后的结构能减少沼气搅拌喷嘴周围的低速区,从而减少厌氧物料在此区域沉积,有利于有效反应容积的维持;
(5)优化后的沼气搅拌系统气体喷嘴管不在突出于罐体内部地面,在干式厌氧长期运行后,采用机械车清洁罐底淤泥代替原始人工清淤成为了可能,提高清淤效率,清淤时间从30多天减少为几天,大大缩短停机维护时间,有利于干式厌氧稳定工况的维持。
附图说明
图1为本发明快速优化方法流程图;
图2为干式厌氧沼气搅拌系统关键部件结构示意图;
图3为本发明快速优化方法过程图;
图4A为本发明实施例中原构型方案;
图4B为本发明实施例中优化构型方案示意图;
图5A至图5C为本发明实施例中三种构型单相流模拟速度分布云图;
图6A为本发明实施例中优化构型设计在喷嘴倾角50°的单相流模拟速度分布云图;
图6B为本发明实施例中优化构型设计在喷嘴倾角50°的单相流模拟速度矢量图;
图7为本发明实施例中优化构型设计的气-液两相流模拟速度矢量数据图。
附图标记说明:
1-倾斜喷管;2-数值喷管;3-反应罐底面;4-倾斜喷管出口截面;5-反应罐内部与喷管形成凹槽;6-喷嘴渐缩段;7-原构型喷管与罐底面形成低速区夹角;8-优化构型I凹槽与喷管形成低速区域。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明作进一步的详细说明。
本发明的干式厌氧沼气搅拌系统的快速优化方法,具体涉及城市生活垃圾、厨余、污泥、畜禽粪便等进行高含固和干式厌氧发酵资源化处理过程中采用气体搅拌设计的优化方法。所述快速优化方法是基于计算流体力学,通过对结构的特征剖析,在几何建模过程中采用降维简化和局部建模相结合,并在模型选择中考虑主相作用为重点,简化多相流模型,通过多重简化,大大节省计算模拟的时间,快速获取优化构型和原构型的流场信息,结果的可视化处理和综合比选,从而获得同时满足延长干式厌氧工程运行维护间隔时间和缩短维护时的清淤时间的沼气搅拌的气体喷嘴设计,为干式厌氧沼气搅拌系统的优化设计提供了一种经济、高效的方法。本发明技术方法成熟稳定,并节约了工艺设计从中试到工程放大的实体建造成本,具有快速、经济、高效的突出特点。
本发明公开了一种干式厌氧沼气搅拌系统的优化方法,包括:
(1)对待优化干式厌氧沼气搅拌系统现状结构分析并获取核心特征结构,建立搅拌系统模型;
(2)将所述搅拌系统模型划分网格;
(3)对计算模型简化,设定边界条件后求解;
(4)若求解结果不满足优化目的,则修改搅拌系统模型的模型参数,直到求解结果满足优化目,此时的搅拌系统模型即为优化后的搅拌系统。
在本发明的一些实施例中,步骤(1)中建立所述搅拌系统模型时将三维结构降维成二维模型。
在本发明的一些实施例中,步骤(1)中所述核心特征结构包括喷嘴相关构型。
在本发明的一些实施例中,步骤(2)中划分所述网格采用非结构化四面体网格生成方法。
在本发明的一些实施例中,步骤(2)中对所述核心特征结构进行局部网格加密处理;同时根据计算精度需求,划分多个网格尺寸不同的计算网格,进行网格无关化比较,选择流场稳定条件下的最少网格数的网格作为下一步的计算基础。
在本发明的一些实施例中,步骤(3)中所述简化方法包将气、液、固三相计算模型简化为两相计算模型或单相计算模型。
在本发明的一些实施例中,步骤(3)中所述计算模型包括用于模拟气、液两相流运动规律的欧拉-欧拉两相流模型和湍流模型。
在本发明的一些实施例中,所述欧拉-欧拉两相流模型在只考察速度场的分布信息时,简化为单相流模型。
在本发明的一些实施例中,步骤(3)中所述求解步骤中包括对求解结果可视化处理。
在本发明的一些实施例中,步骤(4)中所述求解结果包括搅拌系统模型结构的多种流场信息。
在一个示例性实施例中,本发明的干式厌氧沼气搅拌系统的快速优化方法,包括如下步骤:
(1)对待优化干式厌氧沼气搅拌系统现状结构分析,以实施效果或目的剖析结构、提取核心结构特征;
其中,干式厌氧沼气搅拌系统涉及气、液、固相三相,或气液两相、液固两相,涉及依据现有的几何结构,采用对称分析、降维简化方法,根据实施目的和优化效果,提取核心问题,进行单一化精准剖析。
(2)快速数字模型化与计算求解,按照核心结构特征简化的要求建模并划分网格,选用适用的简化模型求解;
(3)结果分析——比对流场信息,优选结构。
其中,快速数字模型化与计算求解包括如下步骤:快速前处理——局部特征几何建模及快速网格划分;快速计算求解——简化计算模型的选择、边界条件的设定与求解;后处理——流场信息获取,具体步骤如下:
1)快速前处理
(i)反应器局部特征几何建模
其中,步骤(1),单一化精准剖析解读反应器几何结构,采用建模软件(ANSYS18.0SpaceClaim、SolidWorks等)仅对反应器的所需优化的特征结构进行几何建模,中心/轴对称结构构建二维几何模型,降低模型复杂程度,节约后续网格划分、计算求解所耗费资源;
(ii)快速网格划分
采用ANSYS前处理软件将上述几何进行网格划分,采用非结构化网格和局部加密法实现网格的快速生成,同时确保网格质量和必要的计算精度,导出目标几何的网格mesh文件;
2)快速求解计算
(i)求解模型的选择与简化
所述的一种干式厌氧沼气搅拌系统,通常涉及沼气对液相和固相的扰动,模型的选择包含气、液、固三相,其中液相为主相,气、固归为第二相、第三相,所述简化模型是将气、液、固三相简化为气、液两相,液、固处理为拟流体,固相特征在液相中体现;
将mesh文件导入ANSYS FLUENT 18.0软件,定义计算模型:由采用欧拉多相流模型模拟干式厌氧沼气搅拌系统的气、液、固三相流简化为欧拉-欧拉两相流模型模拟气、液两相流运动规律,采用标准的湍动能耗散率(k-ε)模型描述湍动黏度μt,建立湍流模型;气、液两相的运动由对应的质量守恒方程和动量守恒方程控制;
所述控制方程如下:
顾量守恒方程即连续性方程:
气、液相动量守恒方程为:
假设气、液相间作用力较弱,不考虑相间作用力。
进一步,如无需观测气、液相交界面变化及作用力,可简化为单相流模型,液相为主相,气相的扰动特征由动能输入即入口速度来体现。
(ii)边界条件设定与计算求解
在模型数值计算中,设定相关参数,气液两相模型中将底部沼气入口作为气体入口,以固定速度为入口边界条件,顶部采用degassing(脱气)作为出口边界条件,其余作为wall(壁面),设置无滑移壁面边界条件,且所有壁面采用不可渗透壁面处理。
利用FLUENT求解器,动量、湍动能耗散采用二阶迎风格式,均方根残差收敛标准判据为1.0E-3,进行稳态迭代计算,选取迭代10000步后的稳定流场结果进行保存。
3)后处理
将上述计算结果保存后,利用后处理软件进行可视化处理,获取不同结构的流场信息。
其中,通过可视化流场信息,观察对比不同结构的流态数据,采用直观分析、选取特征点对比分析各项数值,通过比选确定优化设计方案。
其中,在干式厌氧沼气搅拌优化设计的应用。
其中,进一步,在干式厌氧沼气搅拌系统中通过加热搅拌的沼气实现冬季中温/高温厌氧发酵中的应用。
所述的一种干式厌氧沼气搅拌系统的快速优化设计方法在干式厌氧沼气搅拌系统中通过加热搅拌的沼气实现冬季中温/高温厌氧发酵中的应用,采用该快速优化设计方法,通过能量方程及输运方程,实现入口沼气温度的模拟以及内部温度场的数值模拟。
其中,沼气加热的温度范围在45℃~85℃。
本发明的目的之一在于提供一种干式厌氧沼气搅拌系统的快速优化设计方法,以经济、快速实现减少干式厌氧工程长期运行过程厌氧罐底部物料的堆积或底部堆积区快速清淤的目的,大大缩短工程停机维护的时间。
本发明的目的之二在于沼气搅拌系统中采用热源沼气的方式替代进料加热,并利用目的一的优化方法,实现温度分布的快速模拟。
以下通过具体实施例结合附图对本发明的技术方案做进一步阐述说明。需要注意的是,下述的具体实施例仅是作为举例说明,本发明的保护范围并不限于此。
请参阅图1-7,本发明提出的干式厌氧沼气搅拌系统的快速优化方法,是基于计算流体力学的数值模拟软件与结果,对比研究不同结构的沼气搅拌系统流场特征,优化后的设计同步实现干式厌氧罐沼气搅拌系统喷嘴周边物料沉积量的减少,以延长干式厌氧维护的周期间隔,同时缩短厌氧工程长期运行后的停机清淤维护时间。分析沼气搅拌系统结构构型,针对核心部件——沼气搅拌喷嘴的结构如图2所示,根据快速优化流程(图3)各步骤方法,顺序开展将几何建模时三维结构降维至二维,获得原构型和优化构型的二维建模图(图4A和图4B),并采用欧拉多相流模型或单相流模型获得主相在不同结构构型中的速度分布流场信息(图5A-图7),结合减少沉降物料量和缩短清淤时间的需求进行综合分析,从而获得优化的沼气搅拌系统。在对沼气进行能量(内能,包括加热、热辐射、热传导)输入的基础上,模拟温度场的变化。本发明提出的方法基于一套操作逻辑。快速优化的关键步骤的通过逻辑的执行与判断来实施。
1.现状结构分析,以实施效果或目的剖析结构、提取核心问题
针对干式厌氧沼气搅拌系统,采用对称法进行结构分析,主要包括储气罐、连通沼气气体管、沼气气体管与厌氧罐底部相连接、阀门等。所述沼气气体管为插入厌氧罐底部,并突出厌氧罐底部内部面一段长度,且气体管罐底内部面的夹角为60~90度,气体管横截面为圆形,夹角为90度的气体管出口处做渐缩处理,夹角60度的气体管喷嘴出口截面垂直于地平面,经分析,喷嘴直插入厌氧罐内部,直接影响气体出口的状态(速度、方向等),因此喷嘴相关构型为沼气搅拌的核心部件。
2.快速数字模型化与计算求解,按照特征几何简化的要求建模并划分网格,选用适用的简化模型求解
(1)几何建模和网格划分
通过结构剖析,利用计算流体力学的几何建模软件,快速选取喷嘴的中心截面构建二维几何模型,取气体管直径的50~100倍长度作为计算域边界,忽略整体厌氧罐的特征,重点突出局部特征,可大大节省计算资源和计算时间;采用非结构化四面体网格生成方法生成网格,对环喷嘴区域进行局部网格加密处理,确保计算精度和网格质量;同时可根据计算精度需求,划分多个网格尺寸不同的计算网格,进行网格无关化比较,选择流场稳定条件下的最少网格数的网格作为下一步的计算基础。
(2)模型选择
沼气搅拌系统涉及气体(沼气)、液体(厌氧消化液)、固体(颗粒物)三相,采用欧拉-欧拉多相模型可实现模拟计算;实际过程中固体液体相混合,因此可将固体和液体混合相视为一种流体,也就简化为气-液两相流;甚至,在只考察速度场的分布信息时,无需对气体和液体进行相区分,可进一步简化为单相流模型。
在温度场的模拟中涉及相间热传导,即气相到液相的热传导,由于相间传导系数不同,不可简化为单一流体。
在欧拉多相流模型中,二维计算域内求解得到各相控制方程。
所述控制方程,包括1)各相的质量守恒方程(连续性方程):
2)气、液相动量守恒方程为:
假设气相为不可压缩流体。气-液相间作用力——曳力较小,可忽略。
湍流模型的选择方面采用标准k-ε模型来描述湍动过程。由于湍流过程的精确模拟非常复杂,因此需考虑特征流态,尽可能减少次相流体对主相流体的影响。
(3)边界条件
在单相流或两相流计算过程中,入口流体设置为固定速度边界条件,边界为5%紊流强度。出口设置为上方,出口设定为大气压下的压力出口,两相流时为脱气出口;喷嘴和管内接触表面设置为内部面,其余管壁、计算域边界等均设置为壁面,且对液体为无滑移壁面,不渗漏壁面。
温度场计算所涉及不同材料间的导热系数,根据傅里叶定律获得;加压下混合气体(沼气)的导热系数采用加压下纯气体的导热系数公式计算,此时采用混合气体的虚拟临界常数,其他分子量、气体密度等均采用混合气体的数值。
所述混合气体(沼气)是指甲烷含量60%和二氧化碳含量40%的混合气体。
(4)计算求解
利用FLUENT求解器,动量、湍动能耗散采用二阶迎风格式,均方根残差收敛标准判据为1.0E-3,进行稳态迭代计算,选取迭代10000步后的稳定流场结果进行保存。
3.结果分析——比对流场信息,优选结构。
通过计算软件进行可视化处理分析,得到不同结构构型下的液相速度场、气相分率、气相速度场等详细流场信息,对比优化构型和原构型的流场数据,综合比选后获得优化的结构设计。根据温度场分布,获得入口气体的温度范围。最终实现干式厌氧沼气搅拌系统的优化设计。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种干式厌氧沼气搅拌系统的优化方法,包括:
(1)对待优化干式厌氧沼气搅拌系统现状结构分析并获取核心特征结构,建立搅拌系统模型;
(2)将所述搅拌系统模型划分网格;
(3)对计算模型简化,设定边界条件后求解;
(4)若求解结果不满足优化目的,则修改搅拌系统模型的模型参数,直到求解结果满足优化目的,此时的搅拌系统模型即为优化后的搅拌系统;
其中,步骤(2)中对所述核心特征结构进行局部网格加密处理;同时根据计算精度需求,划分多个网格尺寸不同的计算网格;
其中,步骤(3)中的所述对计算模型进行简化包括:将气、液、固三相简化为气、液两相,液、固处理为拟流体,固相特征在液相中体现,其中,在模型数值计算中,设定相关参数,气液两相模型中将底部沼气入口作为气体入口,以固定速度为入口边界条件,顶部采用脱气作为出口边界条件,其余作为壁面,设置无滑移壁面边界条件,且所有壁面采用不可渗透壁面处理,其中,动量、湍动能耗散采用二阶迎风格式。
2.根据权利要求1所述的优化方法,其特征在于,
步骤(1)中建立所述搅拌系统模型时将三维结构降维成二维模型。
3.根据权利要求1所述的优化方法,其特征在于,
步骤(1)中所述核心特征结构包括喷嘴相关构型。
4.根据权利要求1所述的优化方法,其特征在于,
步骤(2)中划分所述网格采用非结构化四面体网格生成方法。
5.根据权利要求1所述的优化方法,其特征在于,
步骤(2)中对所述计算网格进行网格无关化比较,选择流场稳定条件下的最少网格数的网格作为下一步的计算基础。
6.根据权利要求1所述的优化方法,其特征在于,
步骤(3)中所述简化方法包将气、液、固三相计算模型简化为单相计算模型。
7.根据权利要求1所述的优化方法,其特征在于,
步骤(3)中所述计算模型包括用于模拟气、液两相流运动规律的欧拉-欧拉两相流模型和湍流模型。
8.根据权利要求7所述的优化方法,其特征在于,
所述欧拉-欧拉两相流模型在只考察速度场的分布信息时,简化为单相流模型。
9.根据权利要求1所述的优化方法,其特征在于,
步骤(3)中所述求解步骤中包括对求解结果可视化处理。
10.根据权利要求1所述的优化方法,其特征在于,
步骤(4)中所述求解结果包括搅拌系统模型结构的多种流场信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110045514.9A CN112836444B (zh) | 2021-01-13 | 2021-01-13 | 基于cfd的干式厌氧沼气搅拌系统的优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110045514.9A CN112836444B (zh) | 2021-01-13 | 2021-01-13 | 基于cfd的干式厌氧沼气搅拌系统的优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112836444A CN112836444A (zh) | 2021-05-25 |
CN112836444B true CN112836444B (zh) | 2024-04-05 |
Family
ID=75928103
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110045514.9A Active CN112836444B (zh) | 2021-01-13 | 2021-01-13 | 基于cfd的干式厌氧沼气搅拌系统的优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112836444B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114029021B (zh) * | 2021-12-01 | 2023-04-14 | 兰州理工大学 | 一种聚偏氟乙烯微米级颗粒制备过程中物料运动情况的分析方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101164915A (zh) * | 2006-10-16 | 2008-04-23 | 林仁吉 | 两相一体化多层次内循环厌氧生物反应器 |
CN102855342A (zh) * | 2012-05-29 | 2013-01-02 | 同济大学 | 一种厌氧连续流搅拌槽式生物制氢反应器的优化设计方法 |
WO2017084106A1 (zh) * | 2015-11-20 | 2017-05-26 | 田川 | 一种数值模拟飞行器流场的系统及方法 |
CN108504548A (zh) * | 2018-04-03 | 2018-09-07 | 农业部规划设计研究院 | 一种全混式厌氧发酵反应器的设计优化方法 |
CN108795747A (zh) * | 2018-09-05 | 2018-11-13 | 扬州首创投资有限公司 | 一种立式厌氧发酵罐及厌氧发酵系统和厌氧发酵方法 |
-
2021
- 2021-01-13 CN CN202110045514.9A patent/CN112836444B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101164915A (zh) * | 2006-10-16 | 2008-04-23 | 林仁吉 | 两相一体化多层次内循环厌氧生物反应器 |
CN102855342A (zh) * | 2012-05-29 | 2013-01-02 | 同济大学 | 一种厌氧连续流搅拌槽式生物制氢反应器的优化设计方法 |
WO2017084106A1 (zh) * | 2015-11-20 | 2017-05-26 | 田川 | 一种数值模拟飞行器流场的系统及方法 |
CN108504548A (zh) * | 2018-04-03 | 2018-09-07 | 农业部规划设计研究院 | 一种全混式厌氧发酵反应器的设计优化方法 |
CN108795747A (zh) * | 2018-09-05 | 2018-11-13 | 扬州首创投资有限公司 | 一种立式厌氧发酵罐及厌氧发酵系统和厌氧发酵方法 |
Non-Patent Citations (4)
Title |
---|
Effect of Sulfate on the Migration and Transformation of Methylmercury in Advanced Anaerobic Digestion of Sludge;Wei yuandong等;《Journal 》;20200309;第41卷(第3期);1425-1431 * |
利用CFD方法优化沼气发酵罐内流场形态的研究综述;黄如一;赵鑫;李江;熊霞;郭亭;薛庆文;罗涛;龙恩深;梅自力;;中国沼气;20180820(第04期);25-30 * |
利用CFD模拟沼气发酵罐内流场形态的计算方法研究综述;王军;王子云;薛庆文;郭亭;梅自力;罗涛;黄如一;;中国沼气(第05期);62-68 * |
污泥厌氧消化反应器CFD数值模拟研究进展;曹秀芹;徐国庆;袁海光;江坤;仇付国;尹伟齐;付昆明;;环境工程学报;20181105(第11期);9-23 * |
Also Published As
Publication number | Publication date |
---|---|
CN112836444A (zh) | 2021-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shen et al. | Improving the mixing performances of rice straw anaerobic digestion for higher biogas production by computational fluid dynamics (CFD) simulation | |
CN109583104B (zh) | 一种基于modelica规范的核反应堆热工水力仿真架构方法 | |
Zhang et al. | Syngas production from microwave-assisted air gasification of biomass: Part 2 model validation | |
CN112836444B (zh) | 基于cfd的干式厌氧沼气搅拌系统的优化方法 | |
Koerich et al. | Optimization of bioreactor operating conditions using computational fluid dynamics techniques | |
Velásquez et al. | Optimization of the basin and inlet channel of a gravitational water vortex hydraulic turbine using the response surface methodology | |
Daryus et al. | CFD simulation of turbulent flows in proto X-3 bioenergy micro gas turbine combustor using std k-ε and rng k-ε model for green building application | |
Zeng et al. | A novel matrixed nonlinear grey Bernoulli model for interval prediction of power generation | |
Zhao et al. | MDO strategy for meridian plane design to improve energy conversion capability of LFR main coolant pump | |
GUO et al. | Study on bifurcation to chaos of surface tension gradient driven flow | |
Müller et al. | Optimization of an externally mixed biogas plant using a robust CFD method | |
Fei et al. | Review of numerical simulation research on submersible mixer for sewage | |
Du et al. | Optimization of the internal circulating fluidized bed using computational fluid dynamics technology | |
Cheng et al. | Study of the comprehensive evaluation of energy efficiency of an oilfield heating furnace based on combination weighting method | |
CN104268789A (zh) | 一种确定钢铁企业能源综合调度可行解的方法 | |
Chen et al. | Numerical Analysis of the Combustion in Micro Gas Turbine with Methane/Biogas Fuels | |
Wang et al. | Optimized Geometric Structure of a Rotational Hydrodynamic Cavitation Reactor | |
Kshirsagar et al. | Mixing performance improvement by passive modifications in an anaerobic digester design | |
Kamiya et al. | Development of JAEA Advanced Multi-Physics Analysis Platform for Nuclear Systems | |
CN112380769A (zh) | 基于稀疏检测和径向基函数插值的虚拟样本生成方法 | |
Sprague et al. | Outcomes from the DOE workshop on turbulent flow simulation at the exascale | |
Cadinu et al. | Relating system-to-CFD coupled code analyses to theoretical framework of a multiscale method | |
Kulbyakina et al. | System Analysis of Energy Efficiency at Oil and Gas Industry Enterprises | |
Kozaka et al. | The 3-D optimization of an axial mixer propeller through high fidelity simulations | |
Chen | Research on material big data processing based on optimization algorithm |
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 |