CN115563905A - 一种泡状跨流型的全耦合群体平衡方法及系统 - Google Patents
一种泡状跨流型的全耦合群体平衡方法及系统 Download PDFInfo
- Publication number
- CN115563905A CN115563905A CN202211397846.4A CN202211397846A CN115563905A CN 115563905 A CN115563905 A CN 115563905A CN 202211397846 A CN202211397846 A CN 202211397846A CN 115563905 A CN115563905 A CN 115563905A
- Authority
- CN
- China
- Prior art keywords
- bubble
- phase
- flow
- turbulence
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000010168 coupling process Methods 0.000 title description 9
- 238000005859 coupling reaction Methods 0.000 title description 9
- 239000007788 liquid Substances 0.000 claims abstract description 107
- 230000005514 two-phase flow Effects 0.000 claims abstract description 72
- 238000003756 stirring Methods 0.000 claims abstract description 60
- 230000016507 interphase Effects 0.000 claims abstract description 51
- 238000012546 transfer Methods 0.000 claims abstract description 44
- 239000012530 fluid Substances 0.000 claims abstract description 26
- 230000000704 physical effect Effects 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims abstract description 16
- 239000012071 phase Substances 0.000 claims description 112
- 238000004581 coalescence Methods 0.000 claims description 83
- 230000007246 mechanism Effects 0.000 claims description 39
- 239000007791 liquid phase Substances 0.000 claims description 33
- 238000013467 fragmentation Methods 0.000 claims description 23
- 238000006062 fragmentation reaction Methods 0.000 claims description 23
- 230000001939 inductive effect Effects 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000005461 lubrication Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000001050 lubricating effect Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 78
- 238000012937 correction Methods 0.000 description 15
- 238000010586 diagram Methods 0.000 description 11
- 239000011800 void material Substances 0.000 description 7
- 239000000126 substance Substances 0.000 description 6
- 150000001875 compounds Chemical class 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000012821 model calculation Methods 0.000 description 4
- 239000003380 propellant Substances 0.000 description 3
- 238000013019 agitation Methods 0.000 description 2
- 230000009172 bursting Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 description 1
- 235000015842 Hesperis Nutrition 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 235000012633 Iberis amara Nutrition 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 210000000746 body region Anatomy 0.000 description 1
- 238000009835 boiling Methods 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000005501 phase interface Effects 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004781 supercooling Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明涉及一种泡状跨流型的全耦合群体平衡方法及系统。该方法包括:获取两相流场的含气率、气泡群数密度以及流体物性参数;确定两相流场中的气泡体积,根据流体物性参数确定变形临界体积;确定气液两相流的流型;根据不同流型的气泡形状相关参数以及两相速度确定不同流型的相间力系数,构建相间动量传递模型,更新含气率以及两相速度,选取气泡聚并破裂模型;更新气泡群数密度以更新相间动量传递模型,根据更新后的相间动量传递模型以及气泡聚并破裂模型确定气泡诱发液体湍流模型,求解气液两相湍流方程,更新两相湍流参数预测泡状流至所述湍流搅拌流转变过程中的湍流参数场、速度场以及相场。本发明能够准确描述泡状流至湍流搅拌流的演化过程。
Description
技术领域
本发明涉及气液两相流动过程模拟领域,特别是涉及一种泡状跨流型的全耦合群体平衡方法及系统。
背景技术
在石油、化工、核、航空、医药等工业领域中广泛出现含有泡状流、湍流搅拌流、弹状流等气液两相流动过程,准确预测由泡状流至湍流搅拌流的流型演变过程对工业系统的过程优化及工艺改造具有重要意义。
例如:液氧、液氢等低温液体被用于液体火箭推进剂,为运载火箭提供动力源。当推进剂流经火箭动力系统加注管路时,局部受热或压降易于诱发推进剂产生过冷沸腾或空化现象,形成气-液泡状流。在湍流输运条件下,气泡在流场中聚并、破裂,演化成多尺寸泡状流、甚至弹状流将造成动力系统内压力波动甚至阻塞输运管路的根源,威胁火箭的安全运载。
实验观测表明,泡状流流型转变气泡群尺寸非均匀性增强的结果;随着含气率增加,当流场中气泡体积增至某一数值,气泡形状由球形转变为帽形,流场中呈现球形、帽形(变形)两类气泡,复杂泡状流转变为湍流搅拌流。现有的气泡聚并、破裂机制模型均基于球形气泡假设推导获得,且未充分考虑含气率、气泡尺寸的影响,对于准确描述泡状流至湍流搅拌流的演化过程将不再完全适用。
发明内容
本发明的目的是提供一种泡状跨流型的全耦合群体平衡方法及系统,以准确描述泡状流至湍流搅拌流的演化过程。
为实现上述目的,本发明提供了如下方案:
一种泡状跨流型的全耦合群体平衡方法,包括:
基于两相流场的边界条件,获取所述两相流场的含气率、气泡群数密度以及流体物性参数;所述两相流场包括液相以及气相;所述流体物性参数包括密度以及粘度;
根据所述含气率以及所述气泡群数密度确定所述两相流场中的气泡体积,并根据所述流体物性参数确定变形临界体积;
根据所述气泡体积以及所述变形临界体积确定气液两相流的流型;所述气液两相流的流型包括泡状流流型以及湍流搅拌流流型;
根据不同的所述流型确定不同流型的气泡形状相关参数;所述气泡形状相关参数包括气泡形状直径以及界面面积密度;在所述泡状流流型中的气泡形状为球形气泡,在所述湍流搅拌流流型中的气泡形状为变形气泡;
根据不同流型的所述气泡形状相关参数以及所述两相流场的两相速度确定不同流型的相间力系数,并根据所述相间力系数构建相间动量传递模型;所述相间力系数包括曳力系数、虚拟质量力系数、侧升力系数、湍流耗散系数以及壁面润滑系数;
根据所述相间动量传递模型求解气液两相质量及动量守恒方程,更新所述含气率以及所述两相流场的两相速度,并根据不同流型选取气泡聚并破裂模型;所述气泡聚并破裂模型包括球形气泡聚并破裂模型以及变形气泡聚并破裂模型;
求解气泡群数密度传递方程,更新气泡群数密度,根据更新的气泡群数密度更新所述相间动量传递模型,并根据更新后的相间动量传递模型以及所述气泡聚并破裂模型确定气泡诱发液体湍流模型;
根据所述气泡诱发液体湍流模型求解气液两相湍流方程,更新两相湍流参数并预测所述泡状流至所述湍流搅拌流转变过程中的湍流参数场、速度场以及相场;所述两相湍流参数包括气相和液相的湍动能、湍流耗散率、湍流强度、湍流长度、湍流时间和湍流粘度。
可选的,所述根据所述气泡体积以及所述变形临界体积确定气液两相流的流型,具体包括:
判断所述气泡体积是否小于所述变形临界体积,得到第一判断结果;
若所述第一判断结果表示为所述气泡体积小于所述变形临界体积,确定气液两相流的流型为泡状流流型;
若所述第一判断结果表示为所述气泡体积不小于所述变形临界体积,确定气液两相流的流型为湍流搅拌流流型。
可选的,在所述泡状流流型中,所述气泡形状直径为:
所述界面面积密度为:
αi1=(36πnαi 2)1/3;其中,αi1为所述泡状流流型中液相或气相的界面面积密度;αi为液相或气相的含气率;当下角标i为l时,表示液相;当下角标i为g时,表示气相;
在所述湍流搅拌流流型中,所述气泡形状直径为:
所述界面面积密度为:
αi2=(36πnαi 2)1/3;其中,αi2为所述湍流搅拌流流型中液相或气相的界面面积密度。
可选的,所述相间动量传递模型为:
可选的,所述气液两相质量及动量守恒方程为:
其中,t为时间;αi为液相或气相的含气率;ρi为液相或气相的密度;ui为液相或气相的流度;g为重力加速度矢量;P'为气液两相共同的压力场;τi为应力张量。
可选的,所述球形气泡聚并破裂模型为:
其中,为所述泡状流流型中总气泡群数源项,N为气泡数量;为所述泡状流流型中湍流运动中气泡随机碰撞导致的气泡聚并的源项,为所述泡状流流型中气相速度梯度诱发气泡聚并的源项,为所述泡状流流型中尾涡夹带机制诱发气泡聚并的源项,为所述泡状流流型中湍流诱发的气泡破裂的源项;
所述变形气泡聚并破裂模型为:
其中,为所述湍流搅拌流流型中总气泡群数源项;为所述湍流搅拌流流型中湍流运动中气泡随机碰撞导致的气泡聚并的源项,为所述湍流搅拌流流型中气相速度梯度诱发气泡聚并的源项,为所述湍流搅拌流流型中尾涡夹带机制诱发气泡聚并的源项,为所述湍流搅拌流流型中湍流诱发的气泡破裂的源项。
可选的,所述气泡群数密度传递方程为:
可选的,所述根据所述气泡诱发液体湍流模型求解气液两相湍流方程,更新两相湍流参数并预测所述泡状流至所述湍流搅拌流转变过程中的湍流参数场、速度场以及相场,之后还包括:
判断所述气泡诱发液体湍流模型是否收敛,得到第二判断结果;
若所述第二判断结果表示为所述气泡诱发液体湍流模型收敛,输出所有计算参数;所述所有计算参数包括含气率、气泡群数密度、两相流速、两相湍流参数和气泡形状相关参数;
若所述第二判断结果表示为所述气泡诱发液体湍流模型未收敛,更新所述流体物性,并继续迭代更新。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供了一种泡状跨流型的全耦合群体平衡方法及系统,引入两相流场的含气率、气泡群数密度以及流体物性参数确定气液两相流的不同流型,并根据不同流型选取气泡聚并破裂模型;气泡聚并破裂模型包括球形气泡聚并破裂模型以及变形气泡聚并破裂模型,不仅依靠于球形气泡假设推导获得,考虑了含气率、气泡尺寸的影响,能够准确描述泡状流至湍流搅拌流的演化过程。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明所提供的泡状跨流型的全耦合群体平衡方法流程图;
图2为本发明所提供的基于流型判断的气泡内部与外部运动模型示意图;
图3为本发明所提供的考虑的帽型气泡形状参数示意图;
图4为本发明所提供的复杂泡状流主导机制分析图;图4(a)为径向气泡尺寸图;图4(b)为径向含气率图;
图5为本发明所提供的高液速工况不同入射含气率下径向气泡尺寸预测图;图5(a)为高液速工况下入射段图;图5(b)为高液速工况下充分发展段图;
图6为本发明所提供的低液速工况不同入射含气率下径向含气率预测图;图6(a)为低液速工况下入射段图;图6(b)为低液速工况下充分发展段图;
图7为本发明所提供的湍流搅拌流工况不同源项的主导机制分析图;图7(a)为湍流搅拌流工况下轴向平均气泡尺寸图;图7(b)为湍流搅拌流工况下径向气泡尺寸图;
图8为本发明所提供的不同的湍流搅拌流工况下径向气泡尺寸发展过程图;图8(a)为不同的湍流搅拌流工况下入射段图;图8(b)为不同的湍流搅拌流工况下过渡段图;图8(c)为不同的湍流搅拌流工况下充分发展段图;
图9为本发明所提供的另一种泡状跨流型的全耦合群体平衡方法流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种泡状跨流型的全耦合群体平衡方法及系统,能够准确描述泡状流至湍流搅拌流的演化过程。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例一
图1为本发明所提供的泡状跨流型的全耦合群体平衡方法流程图,如图1所示,一种泡状跨流型的全耦合群体平衡方法,包括:
步骤101:基于两相流场的边界条件,获取所述两相流场的含气率、气泡群数密度以及流体物性参数;所述两相流场包括液相以及气相;所述流体物性参数包括密度以及粘度。
步骤102:根据所述含气率以及所述气泡群数密度确定所述两相流场中的气泡体积,并根据所述流体物性参数确定变形临界体积。
在实际应用中,在迭代计算中,首先获取上一迭代步长中含气率及气泡群数密度,根据Vb=α/n计算控制容积内气泡体积,并根据更新的物性参数获得变形临界体积Vcr判断流型,如果流场中气泡体积Vb小于变形临界体积Vcr,则气液两相流为泡状流,如果超过即为湍流搅拌流,具体流程图见图2;在泡状流流型中,气泡形状参数与界面面积密度分别为:
αi1=(36πnαi 2)1/3 (2)
其中,db为气泡直径;Vb为气泡体积;αi1为所述泡状流流型中液相或气相的界面面积密度;αi为液相或气相的含气率;当下角标i为l时,表示液相;当下角标i为g时,表示气相。
在湍流搅拌流流型中,气泡形状参数与界面面积密度分别为:
αi2=(36πnαi 2)1/3 (4)
其中,Db为气泡直径;αi2为所述湍流搅拌流流型中液相或气相的界面面积密度。
步骤103:根据所述气泡体积以及所述变形临界体积确定气液两相流的流型;所述气液两相流的流型包括泡状流流型以及湍流搅拌流流型。
步骤104:根据不同的所述流型确定不同流型的气泡形状相关参数;所述气泡形状相关参数包括气泡形状直径以及界面面积密度;在所述泡状流流型中的气泡形状为球形气泡,在所述湍流搅拌流流型中的气泡形状为变形气泡。
步骤105:根据不同流型的所述气泡形状相关参数以及所述两相流场的两相速度确定不同流型的相间力系数,并根据所述相间力系数构建相间动量传递模型;所述相间力系数包括曳力系数、虚拟质量力系数、侧升力系数、湍流耗散系数以及壁面润滑系数。
在实际应用中,根据流型中气泡形状相关参数计算所有相间力系数,更新相间动量传递模型,模型计算中涉及的含气率、两相流速、湍动参数以及气泡形状相关参数均采用上一迭代步长内计算结果;该步骤模型包含以下部分:
其中各相间力系数具体如下:
(a)曳力系数
其中,Eo为厄特沃什数;E为平均纵横比;F(E)为平均纵横比相关系数。
(b)虚拟质量力系数
其中,CVM1为虚拟质量力系数。
(c)侧升力系数
Dh=ds(1+0.163Eo0.757)1/3 (12)
其中,Reb为气泡雷诺数;f(Eod)为修正厄特沃什数相关系数;Eod为修正厄特沃什数;Dh为水力直径,ds为径向平均气泡直径。
(d)湍流耗散系数
St=τg/τl (14)
τg=4db/(3CD|ug-ul|) (15)
其中,CTD为湍流耗散系数;Cμ为湍流模型常数;St为湍流斯托克斯数;CD为曳力系数;k为湍流能;ε为湍流耗散率;τg为气相应力张量;τl为液相应力张量。
(e)壁面润滑系数
CWC=10;CWD=6.8;p=1.7
其中,CW(Eo)为厄特沃什数相关系数;CWD为阻尼系数;CWC为截止系数;D为通道直径;yw为壁面间距。
步骤106:根据所述相间动量传递模型求解气液两相质量及动量守恒方程,更新所述含气率以及所述两相流场的两相速度,并根据不同流型选取气泡聚并破裂模型;所述气泡聚并破裂模型包括球形气泡聚并破裂模型以及变形气泡聚并破裂模型。
在实际应用中,求解气液两相质量、动量守恒方程,更新含气率、两相流速;方程主要包含以下几个部分:
下角标i=l或g表示液相或气相,α、ρ与u分别为液相或气相的含气率、密度和速度,g为重力加速度矢量,t为时间;而气液两相共同的压力场为P'、应力张量τ可分别表示为:
所述球形气泡聚并破裂模型为:
其中,为所述泡状流流型中总气泡群数源项,N为气泡数量;为所述泡状流流型中湍流运动中气泡随机碰撞导致的气泡聚并的源项,为所述泡状流流型中气相速度梯度诱发气泡聚并的源项,为所述泡状流流型中尾涡夹带机制诱发气泡聚并的源项,为所述泡状流流型中湍流诱发的气泡破裂的源项;
所述变形气泡聚并破裂模型为:
其中,为所述湍流搅拌流流型中总气泡群数源项;为所述湍流搅拌流流型中湍流运动中气泡随机碰撞导致的气泡聚并的源项,为所述湍流搅拌流流型中气相速度梯度诱发气泡聚并的源项,为所述湍流搅拌流流型中尾涡夹带机制诱发气泡聚并的源项,为所述湍流搅拌流流型中湍流诱发的气泡破裂的源项。
步骤107:求解气泡群数密度传递方程,更新气泡群数密度,根据更新的气泡群数密度更新所述相间动量传递模型,并根据更新后的相间动量传递模型以及所述气泡聚并破裂模型确定气泡诱发液体湍流模型。
在实际应用中,所述气泡群数密度传递方程为:
步骤108:根据所述气泡诱发液体湍流模型求解气液两相湍流方程,更新两相湍流参数并预测所述泡状流至所述湍流搅拌流转变过程中的湍流参数场、速度场以及相场;所述两相湍流参数包括气相和液相的湍动能、湍流耗散率、湍流强度、湍流长度、湍流时间和湍流粘度。
求解气液两相湍流方程,更新两相湍流参数;气液两相湍流方程具体如下:
(a)液相标准k-e湍流模型:
k-ε模型中剪切诱发湍流粘性系数由湍动能、湍流耗散率表示:
湍动能k和湍流耗散率ε方程如下:
式中,Gk代表平均速度梯度引起的湍流动能,可由Gk=τl:ul表示;Gb代表气泡诱发湍动能;YM为可压缩湍流中脉动体胀率对整个耗散率的贡献;Cε1、Cε2和Cε3代表湍流模型常数,σε和σk分别代表k和ε的湍流Prandtl数,分别为:Cε1=1.44,Cε2=1.92,Cμ=0.09,σK=1,σε=1.3。
(b)液相湍流方程修正
1)附加粘度
2)液体湍动能源项
Gk,bit=Ck3(Gk,bit1+Gk,bit2) (31)
式中,σ为表面张力,上标B和C分别代表破裂和聚并,Gk,bit为液体湍动能源项,Ck3为湍动能源项系数,Gk,bit2为聚并、破裂导致气泡释放的界面自由能主导的湍动能源项。Gk,bit1为由曳力做功主导的湍动能源项, 为基于平均气泡数数量的气泡破裂的源项,为基于平均气泡数数量的气泡聚并的源项,α为含气率,ai为单位体积的界面面积
3)液体湍流耗散频率源项
式中,为气泡驻留时间,为气泡碰撞时间,ds为径向气泡平均直径,C3ε1、C3ε2为修正系数,Gε,bit为液体湍流耗散频率源项,CVM为虚拟质量力系数,ds为径向平均气泡直径,αmax为含气率最大值。
(c)气相零方程模型
μt=ρfμUtlt (40)
式中,μt为分散相湍流黏度,lt为湍流长度标度,Ut为湍流速度标度,取流体域中的最大速度;lt为湍流长度标度,fμ为比例常数,VD为流体域体积。
式中,σg为气相的湍流普朗特数,值为1。
在实际应用中,所述步骤108之后还包括:判断所述气泡诱发液体湍流模型是否收敛,若是,输出所有计算参数;所述所有计算参数包括含气率、气泡群数密度、两相流速、两相湍流参数、气泡形状相关参数;若否,更新所述流体物性,并继续迭代更新;其中,收敛标准为:物理模型中所有求解变量的场叠加和变化均小于10-5。
在实际应用中,对于两种气泡的聚并破裂模型进行了修正,具体包括以下步骤:
改进的球形气泡聚并破裂模型:
本发明通过比较其他聚并破裂模型的构成与不足,采用下列统一表达式改进聚并、破裂模型:
考虑到气泡聚并现象是不同聚并机制累积的结果且不同的聚并机制具有不同的聚并时间与聚并效率,气泡聚并率为:
(a)湍流随机运动诱发气泡聚并:
本发明采用下列形式表示湍流随机运动诱发气泡聚并的模型:
1)湍流涡尺寸相关的聚并时间:
以气泡体积表示气泡间的相对脉动速度:
单气泡的平均滑移时间为:
其中,ub为气泡脉动速度。
2)物性相关的聚并效率:
采用Luo and Svendsen能量守恒关系得到的气泡碰撞时间:
引入韦伯数的概念:
将上述两式代入液膜拉伸模型得到气泡聚并效率为:
3)修正系数与最大含气率:
采用修正因子γRC校正气泡滑移时间:
其中,αmax为最大含气率。
其中两气泡之间的平均距离hb为:
气泡平均湍动距离与平均速度、滑移时间有关,可由下式表示:
此外,根据球形气泡假设得到理论最大含气率:
根据上述提出的聚并破裂模型的统一表达形式,湍流随机碰撞诱发气泡聚并率为:
(b)气泡速度差诱发气泡聚并:
由气相速度差造成的气泡聚并机制,聚并模型与湍流随机碰撞机制相同:
在气相速度差诱发气泡聚并机制中,气泡的滑移时间为:
与湍流随机碰撞机制不同,相邻气泡的接近速度与湍动速度无关,可采用气相速度梯度表示:
该聚并模型中,考虑到气泡体积对气泡滑移时间的影响,修正系数γVG=γRC;在相对运动速度下,气泡到达相邻气泡的可能性为:
式中,气泡滑移距离采用气泡相对速度与滑移时间表示:
(c)湍流涡冲击下气泡破裂:
根据统一表达式,采用破裂时间与破裂效率表示的破裂模型为:
1)基于最小湍流涡的破裂时间
修正系数γTI校正滑移时间:
气泡滑移时间表示为:
2)最小湍流涡:
根据力平衡与能量平衡关系推导得到诱发气泡破裂的最小湍流涡尺寸分别为:
3)不同控制区的破裂效率:
当Ve,min1和Ve,min2为诱发气泡破裂的最小湍流涡尺寸,气泡分别处于湍流冲击力控制区和湍流冲击能控制区,气泡破裂效率分别为:
改进的变形气泡聚并破裂模型:
考虑了尾涡夹带机制诱发的帽型气泡聚并现象,以及气泡变形因素对气泡与气泡间碰撞频率、湍流涡与气泡间碰撞频率、聚并效率、破裂效率等因素造成的影响,推导了帽型气泡聚并破裂模型。帽型气泡形状参数示意图如图3所示。
(a)尾涡夹带机制诱发气泡聚并:
其中,γWE为修正因子。
单气泡滑移时间可表示为:
得到帽型气泡尾涡内夹带气泡的平均滑移时间为:
(b)湍流随机碰撞诱发气泡聚并:
帽型气泡碰撞频率与聚并效率表达式如下:
其中,Sb为气泡群相互作用机制源项。
根据统一特征时间表达形式,湍流随机碰撞诱发气泡聚并率为:
式中帽型气泡的滑移时间与碰撞时间分别为:
其中,气泡平均湍动距离与平均湍动速度、滑移时间有关,可由下式表示:
(c)气泡速度差诱发气泡聚并:
将上述模型中湍动速度采用气泡速度差代替,得到气相速度梯度诱发气泡聚并率:
在以上模型中帽型气泡速度差造成气泡滑移时间、碰撞时间与聚并效率分别为:
式中,urg2为两帽型气泡在碰撞过程中速度差,可由下式表示:
修正气泡滑移时间的气泡体积影响因子γVG形式不变。与气泡运动距离有关的可能性函数为:
(d)湍流涡冲击下气泡破裂:
气泡滑移时间可表示为:
根据力平衡与能量平衡关系推导得到诱发气泡破裂的最小湍流涡尺寸分别为:
当Ve,min1和Ve,min2为诱发气泡破裂的最小湍流涡尺寸,气泡分别处于湍流冲击力控制区和湍流冲击能控制区,气泡破裂效率分别为:
如图4-图8所示,本发明将流场气泡界定为球形、帽形两类气泡,经过研究获得了两类气泡的聚并、破裂机制模型,球形气泡聚并、破裂机制模型与变型气泡聚并、破裂机制模型均由湍流随机碰撞、气泡速度差、尾涡夹带机制诱发气泡聚并机制以及湍流涡碰撞诱发气泡破裂机制构成。将上述气相质量传递机制模型与相间动量传递机制模型,分别嵌入双流体模型、湍流模型、群数密度传递方程构建的全耦合群体平衡模型框架,形成了一种适用于涵盖泡状流-湍流搅拌流的全耦合群体平衡模型及数值模拟装置。相对于传统的泡状流数学模型,泡状跨流型的全耦合群体平衡模型基于双流体模型和湍流模型及群数密度传递方程框架,修正了球形与变形两类气泡聚并、破裂机制模型;同时,该模型将预测范围由泡状流拓展至高含气率泡状流、以及泡状流至湍流搅拌流流型,改进后模型对流场中气泡尺寸、气液流速及含气率等参数预测精度提升了50%。
相对于CFD-PBM全耦合群体平衡模型,该模型采用群数密度传递方程代替多组气泡群平衡方程,对于气泡尺寸范围较宽的流型单位控制容积内控制方程数量大大减少、模型计算效率高。
图9为本发明所提供的另一种泡状跨流型的全耦合群体平衡方法流程图,如图9所示,具体包括如下步骤:
步骤2,在迭代计算中,首先获取上一迭代步长中含气率及气泡群数密度,根据Vb=α/n计算控制容积内气泡体积,并根据更新的物性参数获得Vcr判断流型。
步骤3,根据流型中气泡形状相关参数计算所有相间力系数,更新相间动量传递模型,模型计算中涉及的含气率、两相流速、湍动参数以及气泡形状相关参数均采用上一迭代步长内计算结果。
步骤4,求解气液两相质量、动量守恒方程,更新含气率、两相流速。
步骤5,根据更新的流动参数重新计算步骤2与步骤3。
步骤6,并根据流型选取气泡聚并、破裂机制模型;模型计算中涉及的含气率、两相流速为更新参数,湍动参数以及气泡形状相关参数均采用上一迭代步长内计算结果。
步骤7,求解气泡群数密度传递方程,更新气泡群数密度。
步骤8,根据更新的流动参数以及气泡形状参数重新计算步骤2、步骤3和步骤6;并根据更新的相间动量传递模型与气泡聚并破裂模型计算气泡诱发液体湍动模型。
步骤9,求解气液两相湍流方程,更新两相湍流参数。
步骤10,判断模型是否收敛,收敛标准为:物理模型中所有求解变量的场叠加和变化均小于10-5;若模型收敛,输出所有计算参数,否则更新流体物性继续迭代。
实施例二
为了执行上述实施例一对应的方法,以实现相应的功能和技术效果,下面提供一种泡状跨流型的全耦合群体平衡系统。
一种泡状跨流型的全耦合群体平衡系统,包括:
参数获取模块,用于基于两相流场的边界条件,获取所述两相流场的含气率、气泡群数密度以及流体物性参数;所述两相流场包括液相以及气相;所述流体物性参数包括密度以及粘度。
气泡体积以及变形临界体积确定模块,用于根据所述含气率以及所述气泡群数密度确定所述两相流场中的气泡体积,并根据所述流体物性参数确定变形临界体积。
气液两相流的流型确定模块,用于根据所述气泡体积以及所述变形临界体积确定气液两相流的流型;所述气液两相流的流型包括泡状流流型以及湍流搅拌流流型。
在实际应用中,所述气液两相流的流型确定模块,具体包括:第一判断单元,用于判断所述气泡体积是否小于所述变形临界体积,得到第一判断结果;泡状流流型确定单元,用于若所述第一判断结果表示为所述气泡体积小于所述变形临界体积,确定气液两相流的流型为泡状流流型;湍流搅拌流流型确定单元,用于若所述第一判断结果表示为所述气泡体积不小于所述变形临界体积,确定气液两相流的流型为湍流搅拌流流型。
气泡形状相关参数确定模块,用于根据不同的所述流型确定不同流型的气泡形状相关参数;所述气泡形状相关参数包括气泡形状直径以及界面面积密度;在所述泡状流流型中的气泡形状为球形气泡,在所述湍流搅拌流流型中的气泡形状为变形气泡。
相间动量传递模型构建模块,用于根据不同流型的所述气泡形状相关参数以及所述两相流场的两相速度确定不同流型的相间力系数,并根据所述相间力系数构建相间动量传递模型;所述相间力系数包括曳力系数、虚拟质量力系数、侧升力系数、湍流耗散系数以及壁面润滑系数。
气泡聚并破裂模型选取模块,用于根据所述相间动量传递模型求解气液两相质量及动量守恒方程,更新所述含气率以及所述两相流场的两相速度,并根据不同流型选取气泡聚并破裂模型;所述气泡聚并破裂模型包括球形气泡聚并破裂模型以及变形气泡聚并破裂模型。
气泡诱发液体湍流模型确定模块,用于求解气泡群数密度传递方程,更新气泡群数密度,根据更新的气泡群数密度更新所述相间动量传递模型,并根据更新后的相间动量传递模型以及所述气泡聚并破裂模型确定气泡诱发液体湍流模型。
预测模块,用于根据所述气泡诱发液体湍流模型求解气液两相湍流方程,更新两相湍流参数并预测所述泡状流至所述湍流搅拌流转变过程中的湍流参数场、速度场以及相场;所述两相湍流参数包括气相和液相的湍动能、湍流耗散率、湍流强度、湍流长度、湍流时间和湍流粘度。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
Claims (10)
1.一种泡状跨流型的全耦合群体平衡方法,其特征在于,包括:
基于两相流场的边界条件,获取所述两相流场的含气率、气泡群数密度以及流体物性参数;所述两相流场包括液相以及气相;所述流体物性参数包括密度以及粘度;
根据所述含气率以及所述气泡群数密度确定所述两相流场中的气泡体积,并根据所述流体物性参数确定变形临界体积;
根据所述气泡体积以及所述变形临界体积确定气液两相流的流型;所述气液两相流的流型包括泡状流流型以及湍流搅拌流流型;
根据不同的所述流型确定不同流型的气泡形状相关参数;所述气泡形状相关参数包括气泡形状直径以及界面面积密度;在所述泡状流流型中的气泡形状为球形气泡,在所述湍流搅拌流流型中的气泡形状为变形气泡;
根据不同流型的所述气泡形状相关参数以及所述两相流场的两相速度确定不同流型的相间力系数,并根据所述相间力系数构建相间动量传递模型;所述相间力系数包括曳力系数、虚拟质量力系数、侧升力系数、湍流耗散系数以及壁面润滑系数;
根据所述相间动量传递模型求解气液两相质量及动量守恒方程,更新所述含气率以及所述两相流场的两相速度,并根据不同流型选取气泡聚并破裂模型;所述气泡聚并破裂模型包括球形气泡聚并破裂模型以及变形气泡聚并破裂模型;
求解气泡群数密度传递方程,更新气泡群数密度,根据更新的气泡群数密度更新所述相间动量传递模型,并根据更新后的相间动量传递模型以及所述气泡聚并破裂模型确定气泡诱发液体湍流模型;
根据所述气泡诱发液体湍流模型求解气液两相湍流方程,更新两相湍流参数并预测所述泡状流至所述湍流搅拌流转变过程中的湍流参数场、速度场以及相场;所述两相湍流参数包括气相和液相的湍动能、湍流耗散率、湍流强度、湍流长度、湍流时间和湍流粘度。
2.根据权利要求1所述的泡状跨流型的全耦合群体平衡方法,其特征在于,所述根据所述气泡体积以及所述变形临界体积确定气液两相流的流型,具体包括:
判断所述气泡体积是否小于所述变形临界体积,得到第一判断结果;
若所述第一判断结果表示为所述气泡体积小于所述变形临界体积,确定气液两相流的流型为泡状流流型;
若所述第一判断结果表示为所述气泡体积不小于所述变形临界体积,确定气液两相流的流型为湍流搅拌流流型。
6.根据权利要求5所述的泡状跨流型的全耦合群体平衡方法,其特征在于,所述球形气泡聚并破裂模型为:
其中,为所述泡状流流型中总气泡群数源项,N为气泡数量;为所述泡状流流型中湍流运动中气泡随机碰撞导致的气泡聚并的源项,为所述泡状流流型中气相速度梯度诱发气泡聚并的源项,为所述泡状流流型中尾涡夹带机制诱发气泡聚并的源项,为所述泡状流流型中湍流诱发的气泡破裂的源项;
所述变形气泡聚并破裂模型为:
8.根据权利要求7所述的泡状跨流型的全耦合群体平衡方法,其特征在于,所述根据所述气泡诱发液体湍流模型求解气液两相湍流方程,更新两相湍流参数并预测所述泡状流至所述湍流搅拌流转变过程中的湍流参数场、速度场以及相场,之后还包括:
判断所述气泡诱发液体湍流模型是否收敛,得到第二判断结果;
若所述第二判断结果表示为所述气泡诱发液体湍流模型收敛,输出所有计算参数;所述所有计算参数包括含气率、气泡群数密度、两相流速、两相湍流参数和气泡形状相关参数;
若所述第二判断结果表示为所述气泡诱发液体湍流模型未收敛,更新所述流体物性,并继续迭代更新。
9.一种泡状跨流型的全耦合群体平衡系统,其特征在于,包括:
参数获取模块,用于基于两相流场的边界条件,获取所述两相流场的含气率、气泡群数密度以及流体物性参数;所述两相流场包括液相以及气相;所述流体物性参数包括密度以及粘度;
气泡体积以及变形临界体积确定模块,用于根据所述含气率以及所述气泡群数密度确定所述两相流场中的气泡体积,并根据所述流体物性参数确定变形临界体积;
气液两相流的流型确定模块,用于根据所述气泡体积以及所述变形临界体积确定气液两相流的流型;所述气液两相流的流型包括泡状流流型以及湍流搅拌流流型;
气泡形状相关参数确定模块,用于根据不同的所述流型确定不同流型的气泡形状相关参数;所述气泡形状相关参数包括气泡形状直径以及界面面积密度;在所述泡状流流型中的气泡形状为球形气泡,在所述湍流搅拌流流型中的气泡形状为变形气泡;
相间动量传递模型构建模块,用于根据不同流型的所述气泡形状相关参数以及所述两相流场的两相速度确定不同流型的相间力系数,并根据所述相间力系数构建相间动量传递模型;所述相间力系数包括曳力系数、虚拟质量力系数、侧升力系数、湍流耗散系数以及壁面润滑系数;
气泡聚并破裂模型选取模块,用于根据所述相间动量传递模型求解气液两相质量及动量守恒方程,更新所述含气率以及所述两相流场的两相速度,并根据不同流型选取气泡聚并破裂模型;所述气泡聚并破裂模型包括球形气泡聚并破裂模型以及变形气泡聚并破裂模型;
气泡诱发液体湍流模型确定模块,用于求解气泡群数密度传递方程,更新气泡群数密度,根据更新的气泡群数密度更新所述相间动量传递模型,并根据更新后的相间动量传递模型以及所述气泡聚并破裂模型确定气泡诱发液体湍流模型;
预测模块,用于根据所述气泡诱发液体湍流模型求解气液两相湍流方程,更新两相湍流参数并预测所述泡状流至所述湍流搅拌流转变过程中的湍流参数场、速度场以及相场;所述两相湍流参数包括气相和液相的湍动能、湍流耗散率、湍流强度、湍流长度、湍流时间和湍流粘度。
10.根据权利要求9所述的泡状跨流型的全耦合群体平衡系统,其特征在于,所述气液两相流的流型确定模块,具体包括:
第一判断单元,用于判断所述气泡体积是否小于所述变形临界体积,得到第一判断结果;
泡状流流型确定单元,用于若所述第一判断结果表示为所述气泡体积小于所述变形临界体积,确定气液两相流的流型为泡状流流型;
湍流搅拌流流型确定单元,用于若所述第一判断结果表示为所述气泡体积不小于所述变形临界体积,确定气液两相流的流型为湍流搅拌流流型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211397846.4A CN115563905A (zh) | 2022-11-09 | 2022-11-09 | 一种泡状跨流型的全耦合群体平衡方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211397846.4A CN115563905A (zh) | 2022-11-09 | 2022-11-09 | 一种泡状跨流型的全耦合群体平衡方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115563905A true CN115563905A (zh) | 2023-01-03 |
Family
ID=84768898
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211397846.4A Pending CN115563905A (zh) | 2022-11-09 | 2022-11-09 | 一种泡状跨流型的全耦合群体平衡方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115563905A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170220050A1 (en) * | 2014-10-22 | 2017-08-03 | Landmark Graphics Corporation | Flow regime identification apparatus, methods, and systems |
CN107335390A (zh) * | 2017-08-30 | 2017-11-10 | 南京大学 | 微界面强化反应器相界面积构效调控模型建模方法 |
CN111191394A (zh) * | 2019-12-23 | 2020-05-22 | 国家海洋环境预报中心 | 一种气泡羽流湍流闭合的方法 |
-
2022
- 2022-11-09 CN CN202211397846.4A patent/CN115563905A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170220050A1 (en) * | 2014-10-22 | 2017-08-03 | Landmark Graphics Corporation | Flow regime identification apparatus, methods, and systems |
CN107335390A (zh) * | 2017-08-30 | 2017-11-10 | 南京大学 | 微界面强化反应器相界面积构效调控模型建模方法 |
CN111191394A (zh) * | 2019-12-23 | 2020-05-22 | 国家海洋环境预报中心 | 一种气泡羽流湍流闭合的方法 |
Non-Patent Citations (2)
Title |
---|
XI-BAO ZHANG 等: "Effects of bubble coalescence and breakup models on the simulation of bubble columns", 《CHEMICAL ENGINEERING SCIENCE》, 6 June 2020 (2020-06-06), pages 1 - 14 * |
段欣悦 等: "两种群体平衡模型在大规模多粒径泡状流中的应用研究", 《西安交通大学学报》, vol. 45, no. 12, 31 December 2011 (2011-12-31), pages 92 - 97 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kassab et al. | Air-lift pumps characteristics under two-phase flow conditions | |
Taha et al. | CFD modelling of slug flow in vertical tubes | |
Gutiérrez et al. | Numerical study of Taylor bubbles rising in a stagnant liquid using a level-set/moving-mesh method | |
Kaiping et al. | A contribution to study on the lift of ventilated supercavitating vehicle with low Froude number | |
Abubakar et al. | Linear stability analysis of Taylor bubble motion in downward flowing liquids in vertical tubes | |
Fatjo | New dimensionless number to predict cavitation in accelerated fluid | |
CN115563905A (zh) | 一种泡状跨流型的全耦合群体平衡方法及系统 | |
Zhou et al. | Evolution of free surface in the formation of thermo-solutocapillary convection within an open cavity | |
CN115659692A (zh) | 一种湍流流动状态下泡状流数值预测方法及系统 | |
Pinto et al. | Experimental study on the transition in the velocity of individual Taylor bubbles in vertical upward co-current liquid flow | |
Maranzoni | Galilean-invariant expression for Bernoulli’s equation | |
Pakhomov et al. | Simulation of the turbulent structure of a flow and heat transfer in an ascending polydisperse bubble flow | |
Li et al. | MUSIG modeling and evaluation of nitrogen bubble coalescence in a bottom-closed vertical tube | |
CN111723451B (zh) | 一种牛顿-拉弗森算法的子通道求解方法 | |
Wang et al. | Soft hydraulics in channels with thick walls: The finite-Reynolds-number base state and its stability | |
Pakhomov et al. | Modeling the effect of bubbles on a pattern and heat transfer in a turbulent polydisperse upward two-phase flow after sudden enlargement in a tube | |
Chegini et al. | Comparison of various turbulence models for violent geysers in vertical pipes | |
Asiagbe et al. | Large eddy simulation of microbubble transport in vertical channel flows | |
Ashrafian et al. | Wall boundary conditions for rough walls | |
Goncharova et al. | Numerical modeling of the tangential stress effects on convective fluid flows in an open cavity | |
Isaev et al. | Modeling a decrease in hydraulic losses during turbulent flow in a U-bend channel with a circular cavern with a large opening angle | |
Ramdin et al. | CFD for Multiphase Flow Transport in Pipelines | |
Darzi et al. | Gravitational effects on hydrodynamic characteristics of horizontal and vertical upward two-phase flows | |
Goyal | Fluid Mechanics and Hydraulic Machines | |
Yoon | Dimensional 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 |