[发明专利]一种泡状跨流型的全耦合群体平衡方法及系统

说明书

本发明涉及一种泡状跨流型的全耦合群体平衡方法及系统。该方法包括：获取两相流场的含气率、气泡群数密度以及流体物性参数；确定两相流场中的气泡体积，根据流体物性参数确定变形临界体积；确定气液两相流的流型；根据不同流型的气泡形状相关参数以及两相速度确定不同流型的相间力系数，构建相间动量传递模型，更新含气率以及两相速度，选取气泡聚并破裂模型；更新气泡群数密度以更新相间动量传递模型，根据更新后的相间动量传递模型以及气泡聚并破裂模型确定气泡诱发液体湍流模型，求解气液两相湍流方程，更新两相湍流参数预测泡状流至所述湍流搅拌流转变过程中的湍流参数场、速度场以及相场。本发明能够准确描述泡状流至湍流搅拌流的演化过程。

技术领域

本发明涉及气液两相流动过程模拟领域，特别是涉及一种泡状跨流型的全耦合群体平衡方法及系统。

背景技术

在石油、化工、核、航空、医药等工业领域中广泛出现含有泡状流、湍流搅拌流、弹状流等气液两相流动过程，准确预测由泡状流至湍流搅拌流的流型演变过程对工业系统的过程优化及工艺改造具有重要意义。

例如：液氧、液氢等低温液体被用于液体火箭推进剂，为运载火箭提供动力源。当推进剂流经火箭动力系统加注管路时，局部受热或压降易于诱发推进剂产生过冷沸腾或空化现象，形成气－液泡状流。在湍流输运条件下，气泡在流场中聚并、破裂，演化成多尺寸泡状流、甚至弹状流将造成动力系统内压力波动甚至阻塞输运管路的根源，威胁火箭的安全运载。

实验观测表明，泡状流流型转变气泡群尺寸非均匀性增强的结果；随着含气率增加，当流场中气泡体积增至某一数值，气泡形状由球形转变为帽形，流场中呈现球形、帽形(变形)两类气泡，复杂泡状流转变为湍流搅拌流。现有的气泡聚并、破裂机制模型均基于球形气泡假设推导获得，且未充分考虑含气率、气泡尺寸的影响，对于准确描述泡状流至湍流搅拌流的演化过程将不再完全适用。

发明内容

本发明的目的是提供一种泡状跨流型的全耦合群体平衡方法及系统，以准确描述泡状流至湍流搅拌流的演化过程。

为实现上述目的，本发明提供了如下方案：

一种泡状跨流型的全耦合群体平衡方法，包括：

基于两相流场的边界条件，获取所述两相流场的含气率、气泡群数密度以及流体物性参数；所述两相流场包括液相以及气相；所述流体物性参数包括密度以及粘度；

根据所述含气率以及所述气泡群数密度确定所述两相流场中的气泡体积，并根据所述流体物性参数确定变形临界体积；

根据所述气泡体积以及所述变形临界体积确定气液两相流的流型；所述气液两相流的流型包括泡状流流型以及湍流搅拌流流型；

根据不同的所述流型确定不同流型的气泡形状相关参数；所述气泡形状相关参数包括气泡形状直径以及界面面积密度；在所述泡状流流型中的气泡形状为球形气泡，在所述湍流搅拌流流型中的气泡形状为变形气泡；

根据不同流型的所述气泡形状相关参数以及所述两相流场的两相速度确定不同流型的相间力系数，并根据所述相间力系数构建相间动量传递模型；所述相间力系数包括曳力系数、虚拟质量力系数、侧升力系数、湍流耗散系数以及壁面润滑系数；

根据所述相间动量传递模型求解气液两相质量及动量守恒方程，更新所述含气率以及所述两相流场的两相速度，并根据不同流型选取气泡聚并破裂模型；所述气泡聚并破裂模型包括球形气泡聚并破裂模型以及变形气泡聚并破裂模型；

求解气泡群数密度传递方程，更新气泡群数密度，根据更新的气泡群数密度更新所述相间动量传递模型，并根据更新后的相间动量传递模型以及所述气泡聚并破裂模型确定气泡诱发液体湍流模型；