CN112100944A - 基于cfd模拟与piv测量的多尺度条件下厌氧消化流场可视化方法及应用 - Google Patents
基于cfd模拟与piv测量的多尺度条件下厌氧消化流场可视化方法及应用 Download PDFInfo
- Publication number
- CN112100944A CN112100944A CN202011014677.2A CN202011014677A CN112100944A CN 112100944 A CN112100944 A CN 112100944A CN 202011014677 A CN202011014677 A CN 202011014677A CN 112100944 A CN112100944 A CN 112100944A
- Authority
- CN
- China
- Prior art keywords
- anaerobic digestion
- cfd
- flow field
- simulation
- piv
- 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
Images
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/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- 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
Abstract
本发明属于环境工程技术领域,涉及不同含固率厌氧消化系统的可视化,尤其涉及一种基于CFD模拟与PIV测量的多尺度条件下厌氧消化流场可视化方法及应用,利用CFD对厌氧消化的流场进行模拟,然后选取与厌氧消化沼液流变特性相似的透明流体模拟工作流体并利用PIV进行测量,所得测量结果与CFD的模拟结果进行验证;分析验证结果,对CFD模拟中相应的网格尺寸、湍流模型、数值算法进行调整,直至两者达到吻合。本发明操作步骤简洁,有助于指导反应器的设计及运行,减少多尺度条件下厌氧消化反应系统内的死区、提高搅拌效率并进一步增加反应器的处理能力。本发明的应用可提高厌氧消化流场的运行性能、运行稳定性以及经济性,并在实际工程中起到指导作用。
Description
技术领域
本发明属于环境工程技术领域,涉及不同含固率厌氧消化系统的可视化,尤其涉及一种基于CFD模拟与PIV测量的多尺度条件下厌氧消化流场可视化方法及应用。
背景技术
厌氧消化不仅能有效处理有机废弃物,还可以回收利用绿色清洁能源。然而,高固厌氧所消化(10%~20%)的沼液是粘稠的非牛顿流体,具有剪切稀化、触变性、粘弹性等复杂的流变特征;导致高固厌氧消化系统难以得到充分搅拌,流场中存在较大范围的短流、死区、停滞区域等,从而引起抑制性中间产物(如挥发性有机酸和氨氮)的累积,大大降低了高固厌氧消化的运行效率。此外,对于低含固率条件下(0.1%~10%)厌氧消化系统中传质、传热差也是厌氧消化反应效率低的原因之一。因此,厌氧消化流场的可视化有助于减少反应器的死区、强化系统的流动性,进而提升厌氧消化系统的运行性能和运行稳定性。
目前,计算流体力学(Computational Fluid Dynamics, CFD)可以用来预测流体的流动、传热、传质等现象,实现目标流场的可视化。CFD技术结合流体力学、数值计算和计算机科学技术,已运用于反应器内部流场的模拟之中。有研究基于污泥的流变学特性,在三维的 CFD 模型中利用层流模型模拟分析了厌氧消化反应器的搅拌性能,提出了用于评价反应器搅拌性能的均匀度指数相关方法。然而,现有研究中CFD对厌氧消化系统流场的模拟尚缺乏验证方法,即对应模拟结果的准确性难以得到有效的确认,模拟准确性难以信赖。
近年来发展的粒子图像测速法(Particle Image Velocimetry,PIV)能够实现空间速度场中全流速场无接触的实时测量,展现丰富的瞬时流速场信息。PIV技术曾用于废水在内循环厌氧消化(Internal Circulation Anaerobic Digestion, ICAD)的流场特性,以及沼气气泡在不同工况下废水厌氧消化中的特性。PIV多用于废水厌氧消化,但是对高固厌氧消化流场的检测极其受限,主要原因为PIV要求检测对象为透明流体,而高固厌氧消化沼液为灰黑色不透明流体。此外,PIV系统由于其技术特性,目前只限于实验室规模,无法应用于实际工程流场中的检测。
发明内容
鉴于CFD技术缺乏验证手段模拟结果准确性难以信赖以及PIV技术仅能运用于实验室规模的缺陷,本发明公开一种基于CFD模拟与PIV测量的厌氧消化流场的多尺度条件下厌氧消化流场,可实现不同含固率不同尺度条件下厌氧消化流场的准确、快速可视化。
技术方案
一种基于CFD模拟与PIV测量的厌氧消化流场的多尺度条件下厌氧消化流场可视化方法,利用CFD对厌氧消化的流场进行模拟,然后选取与不同含固率厌氧消化沼液流变特性相似的透明流体作为模拟工作流体并利用PIV进行测量,所得测量结果与CFD的模拟结果进行验证;分析验证结果,对CFD模拟中相应的网格尺寸、湍流模型、数值算法进行调整,直至两者达到吻合。
本发明较优公开例中,所述厌氧消化的流场含固率为0.1~20%。
本发明较优公开例中,所述利用CFD对厌氧消化的流场进行模拟,包括如下步骤:
步骤1:揭示有机废弃物(包括不限于污泥、畜禽粪便、餐厨垃圾)厌氧消化沼渣沼液流变学特性:利用动态振荡试验测量其粘弹性,在高剪切速率下,记录沼渣沼液的剪切稀释特性,以幂率模型模拟对应的流变学特性,所述剪切速率为0.01~100s-1;
步骤2:用扫描电镜和激光粒度仪分析厌氧消化沼渣沼液的物质结构,记录之。
根据上述步骤测得的厌氧消化沼渣沼液流变学特性及物质结构特征,确定与厌氧消化沼渣沼液流变学特性和物质结构特征相似的透明流体作为模拟工作流体。
进一步的,对于低固厌氧消化中含固率< 3%的沼液,采用去离子水作为工作模拟流体,具有高灵敏、快速、经济的效益;对于低固厌氧消化含固率在3%~10%,采用不同浓度的聚丙烯酰胺溶液作为工作模拟流体;对于含固率在10%~20%的高固厌氧消化可采用不同浓度锂藻土悬浮液进行模拟。
本发明较优公开例中,所述利用PIV进行测量,包括如下步骤:
步骤A. 按几何相似建立目标厌氧消化系统的冷模反应器,即将目标厌氧消化系统按照比例缩小至实验室尺寸并构建其对应模型,将模拟工作流体置入冷模反应器中,选取多个关键面以观察厌氧消化冷模反应器内的速度、剪切力、涡流和死区的分布规律,其中,所述关键面,除了轴对称中心平面,根据需要还可以关注其他不同切面;
步骤B. 利用PIV系统对选取的关键面进行测量,光束由柱面透镜发射,聚焦于反应器的对应观测面,CMOS相机放置于垂直于激光光束平面的位置;
步骤C. 对PIV系统所捕获的照片进行处理,通过DynamicStudio软件(DantecDynamics,丹麦)实现所述流场的可视化,测得厌氧消化全流场的速度分布、剪切力分布、涡流分布以及死区分布数据。
本发明较优公开例中,所述CFD的模拟结果,即利用CFD对厌氧消化流场进行数值模拟的结果;具体的,基于Gambit软件建立厌氧消化反应器三维模型,进行网格划分,输出计算文件,导入Fluent软件进行计算,计算收敛后导出结果,以CFD-Post进行后处理,获得多个关键面的速度、剪切力、涡流和死区分布规律的模拟结果。
将PIV测得的厌氧消化流场数据与CFD模拟结果进行对比,若与模拟结果匹配度不佳,则对CFD模拟中相应的网格尺寸、湍流模型、数值算法进行调整优化,直至两者结果达到吻合。
本发明还公开了所述方法的应用,即将吻合后的CFD模拟中相应的网格尺寸、湍流模型、数值算法对实际工程中厌氧消化进行CFD计算,采用CFD-Post对厌氧消化实际工程CFD计算结果进行后处理,分析并输出得到实际工程厌氧消化流场的水力学特性。
本发明通过PIV测量结果和相应的CFD数值计算进行对比,用PIV测量结果来约束数值模拟方法,如此使得CFD数值模拟的结果更准确有效;本发明的实施可优化和改进厌氧消化反应器构型、搅拌设备以及搅拌策略,能够增强厌氧消化的混合程度,从而强化厌氧消化的运行性能与运行稳定性。本发明突破了PIV测量技术仅限于实验室规模的局限,有效解决了CFD数值模拟技术缺乏验证方法的问题,可在实际工程中起到指导作用。
有益效果
本发明克服了CFD技术缺乏验证方法,模拟准确性不佳以及PIV技术仅局限于实验室规模的缺点,获得多尺度条件下厌氧消化流程水力学特征,特别是实际工程中厌氧消化流程的相关水力学特征。本发明操作步骤简洁,有助于指导反应器的设计及运行,减少厌氧消化反应系统内的死区、提高搅拌效率并进一步增加反应器的处理能力。本发明的应用可提高厌氧消化流程的运行性能、运行稳定性以及经济性,并在实际工程中起到指导作用。
附图说明
图1. 基于CFD模拟与PIV测量的多尺度条件下厌氧消化流场可视化流程图;
图2. 冷模试验反应器;
图3. 含固率为0.5 %的污泥低固厌氧消化流场;
图4.PIV技术原理示意图;
图5. (a)PIV测量结果,(b)CFD模拟结果;
图6. 叶轮A、叶轮B和叶轮C的示意图;
图7. 叶轮A不同转速条件下的流场图;
图8. 叶轮B不同转速条件下的流场图;
图9. 叶轮C不同转速条件下的流场图;
图10. 中尺度条件下高固厌氧消化反应器示意图;
图11. 中尺度条件下高固厌氧消化流场图;
图12. 大尺度条件下(实际工程)高固厌氧消化反应器示意图;
图13. 中尺度条件下(实际工程)高固厌氧消化流场图。
具体实施方式
下面结合实施例对本发明进行详细说明,以使本领域技术人员更好地理解本发明,但本发明并不局限于以下实施例。
实施例1
一种基于CFD模拟与PIV测量的多尺度条件下厌氧消化流场可视化方法,包括如下步骤:
步骤1、以含固率为0.5 %的污泥厌氧消化系统为对象,采用去离子水作为其模拟工作流体;
步骤2、按几何相似建立厌氧消化冷模模型,无生物化学反应,如图2,该冷模反应器为透明圆筒,直径为100mm,高度为200mm,试验液体填充高度为120mm,为了减小光折射指数的影响,将冷模反应器置于装满水的透明方形容器中;
将去离子水装入冷模反应器中,选取轴对称中心平面用以观察厌氧消化反应器内的速度、剪切力、涡流和死区的分布规律;
步骤3、采用PIV系统对厌氧消化流场进行测量,对PIV测得照片进行处理,通过DynamicStudio软件(Dantec Dynamics,丹麦)实现其流场的可视化;
步骤4、对厌氧消化方案进行数值模拟计算,通过CFD-Post进行后处理,对多个关键面的速度、剪切力、涡流和死区分布规律进行采集;
步骤5、将CFD模拟结果与采用PIV技术测量的厌氧消化流场进行对比,验证了CFD结果的可靠性,准确性达到89%以上。
输出相关网格尺寸、湍流模型、数值算法应用于实际工况条件下的流场可视化,可视化结果见图3。
实施例2
一种基于CFD模拟与PIV测量的多尺度条件下厌氧消化流场可视化方法,包括如下步骤:
步骤1、以含固率为14.3%的猪粪高固厌氧消化系统为对象,测量其沼渣沼液的流变特性,并采用幂率模型模拟对应的流变特性;采用扫描电镜和激光粒度仪对猪粪高固厌氧消化的物质结构特性进行测量,获得高固厌氧消化沼渣沼液的物质结构特征,据此,选用3.00wt %的透明锂藻土悬浮液作为模拟工作流体;
步骤2、按几何相似建立高固厌氧消化冷模模型,无生物化学反应,如图2,该冷模反应器为透明圆筒,直径为100mm,高度为200mm,试验液体填充高度为120mm,为了减小光折射指数的影响,将冷模反应器置于装满水的透明方形容器中;
将3.00 wt %的透明锂藻土悬浮液装入冷模反应器中,选取轴对称中心平面用以观察高固厌氧消化反应器内的速度、剪切力、涡流和死区的分布规律;
步骤3、如图4采用PIV系统对高固厌氧消化流场进行测量,对PIV测得照片进行处理,通过DynamicStudio软件(Dantec Dynamics,丹麦)实现其流场的可视化;
步骤4、对高固厌氧消化方案进行数值模拟计算,通过CFD-Post进行后处理,对多个关键面的速度、剪切力、涡流和死区分布规律进行采集;
步骤5、将CFD模拟结果与采用PIV技术测量的高固厌氧消化流场进行对比,对比结果见图5;在100rpm工况下叶轮A模拟HSAD混合流场与PIV测量结果具有相似性,即其混合区仅在叶轮周围且叶轮周围的流速在0.1m/s左右,这验证了CFD结果的可靠性。
输出相关网格尺寸、湍流模型、数值算法应用于实际工况条件下的流场可视化。
本实施例公开了不同搅拌桨不同转速条件下的搅拌特性,不同叶轮的示意图如图6,分别为叶轮A、叶轮B和叶轮C,三个叶轮的详情如下:
(1)搅拌桨A为双层叶轮,叶轮直径为50mm;
(2)搅拌桨B设计为直径90mm的三层叶轮;
(3)搅拌桨C为螺带叶轮,叶轮直径为90mm;
叶轮A、叶轮B和叶轮C不同转速条件下的模拟结果如图7、图8、图9,以下将三个不同叶轮的流场描述如下:
(1)叶轮A
不同转速下叶轮A的流场如图7所示。从图中可以看出:在50rpm、100rpm和150rpm的低转速下,混合区仅仅位于叶轮周围。随着转速进一步提高,两个叶轮之间的流动性不断增强。此外,即使转速高达300rpm,搅拌轴附近的流体也是处于停滞状态。
此外,提高转速可以增加流场的流速和混合区域,但是对于提高高固厌氧消化的混合区域效果有限。例如,当转速为50rpm时,x=-25 mm-25 mm,y=10 mm-30 mm,60 mm-80mm范围内混合效果较好,当转速增至300rpm时,x=-40 mm-40 mm,y=0 mm-90 mm范围内混合效果较好。然而,值得注意的是,即使在300rpm的转速下,区域(y=90 mm-150 mm)也不能混合。由此可见,轴向混合和径向混合都随着转速的提高而增强,但增强的程度有限。
(2)叶轮B
在转速为50rpm时,仅仅叶轮周围发生流体的流动,而叶轮之间的混合不明显。然而,随着转速的进一步提高,三个叶轮之间的区域混合逐渐增强。随着转速的增加,流场流速和混合区域都有明显的增加。在各个转速下都能实现径向混合,这表明径向混合在很大程度上取决于叶轮的直径。在转速为50rpm时,y=20 mm-140 mm范围内实现了轴向混合,但当转速增加到300rpm时,轴向混合增加到10 mm-150 mm,轴向搅拌效果较好。
(3)叶轮C
叶轮C搅拌下高固厌氧消化流场的混合相当充分,随着转速的提高混合逐渐由边缘向中心靠拢,其反应器边缘搅拌更为强烈。这与叶轮A和叶轮B搅拌条件下的流场存在很明显的差异,在叶轮A和叶轮B搅拌条件下,反应器中心的混合比边缘的混合更充分。当转速为50rpm时,x=-50 mm-30 mm和x=30 mm-50 mm范围可实现径向混合,而当转速为300rpm时,x=-50 mm-50 mm范围均可搅拌充分。此外,轴向混合即使在低转速下也足够充分。结果表明,螺带叶轮比桨叶轮更适合实高固厌氧消化的充分混合。
实施例3
一种基于CFD模拟与PIV测量的多尺度条件下厌氧消化流场可视化方法,包括如下步骤:
步骤1、以含固率为14.3%的猪粪高固厌氧消化系统为对象,测量其沼渣沼液的流变特性,并采用幂率模型模拟对应的流变特性;采用扫描电镜和激光粒度仪对猪粪高固厌氧消化的物质结构特性进行测量,获得高固厌氧消化沼渣沼液的物质结构特征,据此,选用3.00wt %的透明锂藻土悬浮液作为模拟工作流体;
步骤2、按几何相似建立高固厌氧消化冷模模型,无生物化学反应,如图2,该冷模反应器为透明圆筒,直径为100mm,高度为200mm,试验液体填充高度为120mm,为了减小光折射指数的影响,将冷模反应器置于装满水的透明方形容器中;
将3.00 wt %的透明锂藻土悬浮液装入冷模反应器中,选取轴对称中心平面用以观察高固厌氧消化反应器内的速度、剪切力、涡流和死区的分布规律;
步骤3、如图4采用PIV系统对高固厌氧消化流场进行测量,对PIV测得照片进行处理,通过DynamicStudio软件(Dantec Dynamics,丹麦)实现其流场的可视化;
步骤4、对高固厌氧消化方案进行数值模拟计算,通过CFD-Post进行后处理,对多个关键面的速度、剪切力、涡流和死区分布规律进行采集;
步骤5、将CFD模拟结果与采用PIV技术测量的高固厌氧消化流场进行对比,对比结果见图5;在100rpm工况下叶轮A模拟HSAD混合流场与PIV测量结果具有相似性,即其混合区仅在叶轮周围且叶轮周围的流速在0.1m/s左右,这验证了CFD结果的可靠性。
输出相关网格尺寸、湍流模型、数值算法应用于实际工况条件下的流场可视化。
以如图10的中型尺度反应器为研究对象,该反应器长×宽×高 = 3000 mm×2000mm×24000mm,其中有效水深为2000mm。搅拌器采用折板搅拌桨,桨叶角度为45°,搅拌桨厚度为100mm。研究其在0.25rpm、5rpm、13.5rpm 以及 50rpm条件下的搅拌特性,如图10所示。
如图11可见,随着搅拌速度的增大,该反应器流体的流速以及搅拌范围均有显著增加。在0.25 rpm条件下,流动主要集中于搅拌桨附近,且最大流速为0.03m/s; 当转速为5rpm时,流动有所增强,最大流速为0.3 m/s,然而搅拌的区域依旧仅限于两个搅拌桨之间;当转速为13.5 rpm时,搅拌桨之间的区域相对转速为5 rpm时有所增强,且最大流速为0.9m/s。转速为50 rpm时,反应器内的搅拌明显增强,最大流速为3 m/s。
实施例4
一种基于CFD模拟与PIV测量的多尺度条件下厌氧消化流场可视化方法,包括如下步骤:
步骤1、以含固率为14.3%的猪粪高固厌氧消化系统为对象,测量其沼渣沼液的流变特性,并采用幂率模型模拟对应的流变特性;采用扫描电镜和激光粒度仪对猪粪高固厌氧消化的物质结构特性进行测量,获得高固厌氧消化沼渣沼液的物质结构特征,据此,选用3.00wt %的透明锂藻土悬浮液作为模拟工作流体;
步骤2、按几何相似建立高固厌氧消化冷模模型,无生物化学反应,如图2,该冷模反应器为透明圆筒,直径为100mm,高度为200mm,试验液体填充高度为120mm,为了减小光折射指数的影响,将冷模反应器置于装满水的透明方形容器中;
将3.00 wt %的透明锂藻土悬浮液装入冷模反应器中,选取轴对称中心平面用以观察高固厌氧消化反应器内的速度、剪切力、涡流和死区的分布规律;
步骤3、如图4采用PIV系统对高固厌氧消化流场进行测量,对PIV测得照片进行处理,通过DynamicStudio软件(Dantec Dynamics,丹麦)实现其流场的可视化;
步骤4、对高固厌氧消化方案进行数值模拟计算,通过CFD-Post进行后处理,对多个关键面的速度、剪切力、涡流和死区分布规律进行采集;
步骤5、将CFD模拟结果与采用PIV技术测量的高固厌氧消化流场进行对比,对比结果见图5;在100rpm工况下叶轮A模拟HSAD混合流场与PIV测量结果具有相似性,即其混合区仅在叶轮周围且叶轮周围的流速在0.1m/s左右,这验证了CFD结果的可靠性。
输出相关网格尺寸、湍流模型、数值算法应用于实际工况条件下的流场可视化。
以如图12的实际尺度反应器为研究对象,该反应器长×宽×高= 36000mm×9000mm× 12000mm,其中有效水深为9000mm。搅拌器采用折板搅拌桨,桨叶角度为30°,搅拌桨厚度为900mm。研究其在0.25rpm、5rpm、13.5rpm以及50rpm条件下的搅拌特性,如图13所示。
如图13可见,随着搅拌速度的增大,该反应器流体的流速以及搅拌范围均有显著增加。在0.25 rpm条件下,流动主要集中于搅拌桨附近,且最大流速为0.2m/s; 当转速为5rpm时,流动有所增强,最大流速为2 m/s,然而反应器两端的搅拌有限;当转速为13.5 rpm时,流动相对转速为5 rpm时有所增强,且最大流速为5 m/s。转速为50 rpm时,反应器内的搅拌明显增强,最大流速为20 m/s。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (9)
1.一种基于CFD模拟与PIV测量的厌氧消化流场的多尺度条件下厌氧消化流场可视化方法,其特征在于:利用CFD对厌氧消化的流场进行模拟,然后选取与不同含固率厌氧消化沼液流变特性相似的透明流体作为模拟工作流体并利用PIV进行测量,所得测量结果与CFD的模拟结果进行验证;分析验证结果,对CFD模拟中相应的网格尺寸、湍流模型、数值算法进行调整,直至两者达到吻合。
2.根据权利要求1所述基于CFD模拟与PIV测量的厌氧消化流场的多尺度条件下厌氧消化流场可视化方法,其特征在于:厌氧消化的流场含固率为0.1~20%。
3.根据权利要求1所述基于CFD模拟与PIV测量的厌氧消化流场的多尺度条件下厌氧消化流场可视化方法,其特征在于,所述利用CFD对厌氧消化的流场进行模拟,包括如下步骤:
步骤1.利用动态振荡试验测量有机废弃物厌氧消化沼渣沼液粘弹性,在高剪切速率下,记录沼渣沼液的剪切稀释特性,以幂率模型模拟对应的流变学特性,所述剪切速率为0.01~100s-1,其中所述有机废弃物包括但不限于污泥、畜禽粪便、餐厨垃圾;
步骤2. 以扫描电镜和激光粒度仪分析厌氧消化沼渣沼液的物质结构,记录之。
4.根据权利要求3所述基于CFD模拟与PIV测量的厌氧消化流场的多尺度条件下厌氧消化流场可视化方法,其特征在于:依据所测得的厌氧消化沼渣沼液流变学特性及物质结构特征,确定与厌氧消化沼渣沼液流变学特性和物质结构特征相似的透明流体作为模拟工作流体。
5.根据权利要求4所述基于CFD模拟与PIV测量的厌氧消化流场的多尺度条件下厌氧消化流场可视化方法,其特征在于:含固率< 3%的沼液,采用去离子水作为工作模拟流体;含固率在3%~10%的沼液沼渣,采用不同浓度的聚丙烯酰胺溶液作为工作模拟流体;含固率在10%~20%的沼液沼渣,采用不同浓度锂藻土悬浮液进行模拟。
6.根据权利要求1所述基于CFD模拟与PIV测量的厌氧消化流场的多尺度条件下厌氧消化流场可视化方法,其特征在于,所述利用PIV进行测量,包括如下步骤:
步骤A. 按几何相似建立目标厌氧消化系统的冷模反应器,即将目标厌氧消化系统按照比例缩小至实验室尺寸并构建其对应模型,将模拟工作流体置入冷模反应器中,选取多个关键面以观察厌氧消化冷模反应器内的速度、剪切力、涡流和死区的分布规律,其中,所述关键面,除了轴对称中心平面,根据需要还应关注不同切面;
步骤B. 利用PIV系统对选取的关键面进行测量,光束由柱面透镜发射,聚焦于反应器的对应观测面,CMOS相机放置于垂直于激光光束平面的位置;
步骤C. 对PIV系统所捕获的照片进行处理,通过DynamicStudio软件(DantecDynamics,丹麦)实现所述流场的可视化,测得厌氧消化全流场的速度分布、剪切力分布、涡流分布以及死区分布数据。
7.根据权利要求1所述基于CFD模拟与PIV测量的厌氧消化流场的多尺度条件下厌氧消化流场可视化方法,其特征在于:所述CFD的模拟结果,即利用CFD对厌氧消化流场进行数值模拟的结果。
8.根据权利要求7所述基于CFD模拟与PIV测量的厌氧消化流场的多尺度条件下厌氧消化流场可视化方法,其特征在于:基于Gambit软件建立厌氧消化反应器三维模型,进行网格划分,输出计算文件,导入Fluent软件进行计算,计算收敛后导出结果,以CFD-Post进行后处理,获得多个关键面的速度、剪切力、涡流和死区分布规律的模拟结果。
9.一种如权利要求1-8任一所述基于CFD模拟与PIV测量的厌氧消化流场的多尺度条件下厌氧消化流场可视化方法的应用,其特征在于:将吻合后的CFD模拟中相应的网格尺寸、湍流模型、数值算法对实际工程中厌氧消化进行CFD计算,采用CFD-Post对厌氧消化实际工程CFD计算结果进行后处理,分析并输出得到实际工程厌氧消化流场的水力学特性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011014677.2A CN112100944B (zh) | 2020-09-24 | 2020-09-24 | 基于cfd模拟与piv测量的多尺度条件下厌氧消化流场可视化方法及应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011014677.2A CN112100944B (zh) | 2020-09-24 | 2020-09-24 | 基于cfd模拟与piv测量的多尺度条件下厌氧消化流场可视化方法及应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112100944A true CN112100944A (zh) | 2020-12-18 |
CN112100944B CN112100944B (zh) | 2022-05-24 |
Family
ID=73755279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011014677.2A Active CN112100944B (zh) | 2020-09-24 | 2020-09-24 | 基于cfd模拟与piv测量的多尺度条件下厌氧消化流场可视化方法及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112100944B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114029021A (zh) * | 2021-12-01 | 2022-02-11 | 兰州理工大学 | 一种聚偏氟乙烯微米级颗粒制备过程中物料运动情况的分析方法 |
CN115403218A (zh) * | 2022-08-24 | 2022-11-29 | 中冶赛迪工程技术股份有限公司 | 煤化工废水处理装置以及处理工艺 |
CN117213548A (zh) * | 2023-11-09 | 2023-12-12 | 达斯玛环境科技(北京)有限公司 | 一种搅拌效果评估方法、系统、设备及存储介质 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101919711A (zh) * | 2010-08-25 | 2010-12-22 | 四川省医学科学院(四川省人民医院) | 基于多普勒图像信息的心脏流场速度矢量场可视化描述方法 |
CN102156085A (zh) * | 2010-12-06 | 2011-08-17 | 昆明理工大学 | 一种化工领域多相混合流场均匀性的检测方法 |
CN102798511A (zh) * | 2012-02-06 | 2012-11-28 | 哈尔滨电机厂有限责任公司 | 水力机械内部流场精细化测试试验装置 |
CN102855342A (zh) * | 2012-05-29 | 2013-01-02 | 同济大学 | 一种厌氧连续流搅拌槽式生物制氢反应器的优化设计方法 |
CN103048114A (zh) * | 2012-12-19 | 2013-04-17 | 浙江大学 | 液力减速器内流场三维piv测试系统试验装置及方法 |
KR101402645B1 (ko) * | 2013-11-21 | 2014-06-03 | 대한민국 | 뱀장어 초기 자어 사육수조 청소시스템 |
CN103969053A (zh) * | 2014-04-02 | 2014-08-06 | 杭州汽轮动力集团有限公司 | 1.5级跨音速轴流压气机试验装置 |
CN104050330A (zh) * | 2014-06-26 | 2014-09-17 | 中国科学院生态环境研究中心 | 升流式厌氧发酵生物制氢反应器的优化设计方法及其应用 |
CN104122068A (zh) * | 2014-06-29 | 2014-10-29 | 吉林大学 | 一种液力变矩器内部流场分离式测试装置 |
CN105842040A (zh) * | 2016-04-25 | 2016-08-10 | 北京建筑大学 | 一种污泥相似溶液的检验方法及装置 |
CN106777589A (zh) * | 2016-11-29 | 2017-05-31 | 中国水利水电科学研究院 | 一种基于流量需求的锯齿型灌水器逆向设计方法 |
US9964495B1 (en) * | 2012-11-02 | 2018-05-08 | University Of Maryland | Method and system for spatially-resolved 3-dimensional characterization of near-field sprays |
CN108504548A (zh) * | 2018-04-03 | 2018-09-07 | 农业部规划设计研究院 | 一种全混式厌氧发酵反应器的设计优化方法 |
-
2020
- 2020-09-24 CN CN202011014677.2A patent/CN112100944B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101919711A (zh) * | 2010-08-25 | 2010-12-22 | 四川省医学科学院(四川省人民医院) | 基于多普勒图像信息的心脏流场速度矢量场可视化描述方法 |
CN102156085A (zh) * | 2010-12-06 | 2011-08-17 | 昆明理工大学 | 一种化工领域多相混合流场均匀性的检测方法 |
CN102798511A (zh) * | 2012-02-06 | 2012-11-28 | 哈尔滨电机厂有限责任公司 | 水力机械内部流场精细化测试试验装置 |
CN102855342A (zh) * | 2012-05-29 | 2013-01-02 | 同济大学 | 一种厌氧连续流搅拌槽式生物制氢反应器的优化设计方法 |
US9964495B1 (en) * | 2012-11-02 | 2018-05-08 | University Of Maryland | Method and system for spatially-resolved 3-dimensional characterization of near-field sprays |
CN103048114A (zh) * | 2012-12-19 | 2013-04-17 | 浙江大学 | 液力减速器内流场三维piv测试系统试验装置及方法 |
KR101402645B1 (ko) * | 2013-11-21 | 2014-06-03 | 대한민국 | 뱀장어 초기 자어 사육수조 청소시스템 |
CN103969053A (zh) * | 2014-04-02 | 2014-08-06 | 杭州汽轮动力集团有限公司 | 1.5级跨音速轴流压气机试验装置 |
CN104050330A (zh) * | 2014-06-26 | 2014-09-17 | 中国科学院生态环境研究中心 | 升流式厌氧发酵生物制氢反应器的优化设计方法及其应用 |
CN104122068A (zh) * | 2014-06-29 | 2014-10-29 | 吉林大学 | 一种液力变矩器内部流场分离式测试装置 |
CN105842040A (zh) * | 2016-04-25 | 2016-08-10 | 北京建筑大学 | 一种污泥相似溶液的检验方法及装置 |
CN106777589A (zh) * | 2016-11-29 | 2017-05-31 | 中国水利水电科学研究院 | 一种基于流量需求的锯齿型灌水器逆向设计方法 |
CN108504548A (zh) * | 2018-04-03 | 2018-09-07 | 农业部规划设计研究院 | 一种全混式厌氧发酵反应器的设计优化方法 |
Non-Patent Citations (2)
Title |
---|
王乐等: "气体搅拌下厌氧消化反应器CFD数值模拟及模型研究", 《农业机械学报》 * |
马吉龙: "基于CFD的厌氧发酵搅拌技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)农业科技辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114029021A (zh) * | 2021-12-01 | 2022-02-11 | 兰州理工大学 | 一种聚偏氟乙烯微米级颗粒制备过程中物料运动情况的分析方法 |
CN115403218A (zh) * | 2022-08-24 | 2022-11-29 | 中冶赛迪工程技术股份有限公司 | 煤化工废水处理装置以及处理工艺 |
CN117213548A (zh) * | 2023-11-09 | 2023-12-12 | 达斯玛环境科技(北京)有限公司 | 一种搅拌效果评估方法、系统、设备及存储介质 |
CN117213548B (zh) * | 2023-11-09 | 2024-01-19 | 达斯玛环境科技(北京)有限公司 | 一种搅拌效果评估方法、系统、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112100944B (zh) | 2022-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112100944B (zh) | 基于cfd模拟与piv测量的多尺度条件下厌氧消化流场可视化方法及应用 | |
Wu | CFD investigation of turbulence models for mechanical agitation of non-Newtonian fluids in anaerobic digesters | |
Cortada-Garcia et al. | Computational fluid dynamic studies of mixers for highly viscous shear thinning fluids and PIV validation | |
Wang et al. | The effect of high speed mixing and polymer dosing rates on the geometric and rheological characteristics of conditioned anaerobic digested sludge (ADS) | |
CN101813641A (zh) | 用于验证多相搅拌混合均匀状态及程度的方法 | |
Wiedemann et al. | Modeling mixing in anaerobic digesters with computational fluid dynamics validated by experiments | |
Wang et al. | Effect of impeller type and scale-up on spatial distribution of shear rate in a stirred tank | |
Montante et al. | Fluid dynamics characterization of a stirred model bio-methanation digester | |
Chung et al. | Reconstruction of 3-D flow field inside miniature stirred vessels using a 2-D PIV technique | |
Wang et al. | The effect of polymer doses and extended mixing intensity on the geometric and rheological characteristics of conditioned anaerobic digested sludge (ADS) | |
Yang et al. | Research on the homogenization evaluation of cemented paste backfill in the preparation process based on image texture features | |
CN103279433B (zh) | 一种表征颗粒团聚体的方法及装置 | |
Hu et al. | Visualization of mass transfer in mixing processes in high solid anaerobic digestion using Laser Induced Fluorescence (LIF) technique | |
Simmons et al. | Mixing in a model bioreactor using agitators with a high solidity ratio and deep blades | |
Arnau et al. | Evaluation of hydraulic mixing performance in a full-scale anaerobic digester with an external liquid recirculation system using CFD and experimental validation | |
Panckow et al. | Characterisation of particle stress in turbulent impeller flows utilising photo-optical measurements of a flocculation system–PART 1 | |
CN105842040B (zh) | 一种污泥相似溶液的检验方法及装置 | |
Maluta et al. | Hydrodynamics, power consumption and bubble size distribution in gas-liquid stirred tanks | |
Wu | CFD prediction of mixing time in anaerobic digesters | |
Vlaev et al. | Shear and skin friction on particles in power-law fluids agitated by flat-blade and fluid foil impellers | |
Zhang et al. | Experimental investigation on the uniformity optimization and chaos characterization of gas-liquid two-phase mixing process using statistical image analysis | |
Huber et al. | Positioning a measurement system for determining the mixing quality in biogas digesters | |
Li et al. | Rheology-dependent mass diffusion characteristics of high-solids anaerobic co-digestate of food waste and sewage sludge | |
Ramírez-Muñoz et al. | Hydrodynamics performance of Newtonian and shear-thinning fluids in a tank stirred with a high shear impeller: Effect of liquid height and off-bottom clearance | |
Park et al. | Examining the effect of hydraulic turbulence in a rapid mixer on turbidity removal with CFD simulation and PIV analysis |
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 |