两相流模拟器 返回
流体工程·化学工程

两相流模拟器(气液流动)

利用Baker图判断气液两相流流型(气泡流、段塞流、环状流等)。用Lockhart-Martinelli关联式计算压降。

流体参数

计算结果

流型-
空泡率 α-
干度 x-
滑移比 S-
M-M参数 X_tt-
两相流乘数 φ²_L-
压降 [Pa/m]-
Lockhart-Martinelli关联式
$X_{tt}= \left(\frac{1-x}{x}\right)^{0.9}\left(\frac{\rho_G}{\rho_L}\right)^{0.5}\left(\frac{\mu_L}{\mu_G}\right)^{0.1}$
$\phi_L^2 = 1 + \frac{C}{X_{tt}}+ \frac{1}{X_{tt}^2}$
Baker图(流型图)
压降 vs 干度 x

什么是两相流流型与压降

🧑‍🎓
「两相流流型」是什么?听起来好复杂。
🎓
简单来说,就是气体和液体在管道里一起流动时,会“排成”不同的队形。比如在空调的制冷剂管路里,液态和气态的制冷剂混在一起流动,有时候是小气泡(气泡流),有时候会形成像香肠一样一节气体一节液体的样子(段塞流)。你试着在模拟器里把“气相质量分数”滑块往右拉高,看看流型图里的点怎么移动,就能直观感受到流量变化对流型的影响了。
🧑‍🎓
诶,真的吗?那不同的“队形”对工程师有什么影响呢?
🎓
影响可大了!不同的流型,管道受到的振动、阻力(也就是压降)和传热效率都完全不同。比如段塞流,液塞会像锤子一样周期性撞击管壁和阀门,容易造成疲劳损坏。工程现场常见的是,设计不当的蒸汽锅炉给水管路就可能出现段塞流,引发剧烈振动。你现在可以同时调高“液相密度”和“气相流速”,看看会不会更容易进入段塞流区域,压降曲线也会猛地跳上去哦。
🧑‍🎓
原来这么危险!那怎么提前算出这个阻力(压降)来避免问题呢?
🎓
这时候就要用到Lockhart-Martinelli关联式了。它是个经验公式,核心思想是先分别算出假如只有液体流动的压降和只有气体流动的压降,然后用一个叫Martinelli参数 $X_{tt}$ 的东西把两者关联起来,最后乘以一个“两相流乘数”就得到真实的两相流压降。你改变上面“液相粘度”和“气相粘度”的滑块,下面公式计算出的 $X_{tt}$ 和压降乘数 $\phi_L^2$ 就会实时变化,你能立刻看到哪个物性参数对压降影响最大。

物理模型与关键公式

首先,我们用一个无量纲参数——Martinelli参数 ($X_{tt}$) 来关联单相流压降。它考虑了气液两相的质量分数、密度和粘度比。

$$X_{tt}= \left(\frac{1-x}{x}\right)^{0.9}\left(\frac{\rho_G}{\rho_L}\right)^{0.5}\left(\frac{\mu_L}{\mu_G}\right)^{0.1}$$

其中,$x$ 是气相质量分数,$\rho_G$ 和 $\rho_L$ 分别是气相和液相密度,$\mu_G$ 和 $\mu_L$ 分别是气相和液相粘度。$X_{tt}$ 越小,表示气相主导作用越强。

得到 $X_{tt}$ 后,通过Lockhart-Martinelli关联式计算两相流压降乘数 $\phi_L^2$。这个乘数代表了实际两相流压降是纯液相流动压降的多少倍。

$$\phi_L^2 = 1 + \frac{C}{X_{tt}}+ \frac{1}{X_{tt}^2}$$

其中,$C$ 是一个经验常数,其取值取决于两相流是均为湍流(C=20)还是其他流态。最终的两相流压降 $\Delta P_{TP} = \phi_L^2 \cdot \Delta P_L$,$\Delta P_L$是假设只有液相流动时的压降。

现实世界中的应用

石油与天然气输送:在油气混输管道中,准确预测段塞流的发生和对应的巨大压降至关重要,这直接关系到泵站的设计间距和管道壁厚的选择,以避免管道振动破裂。

制冷与空调系统:制冷剂在蒸发器和冷凝器中经历剧烈的相变过程,流型从环状流到雾状流变化。通过模拟压降,可以优化管路设计,确保制冷效率和压缩机稳定运行。

核反应堆冷却系统:反应堆堆芯内的沸腾水流属于典型的两相流。精确计算不同流型下的压降和空泡份额,是评估冷却剂流量、确保堆芯安全冷却的关键。

化工精馏塔:塔板上的气液接触本质上也是两相流。了解流型有助于设计更高效的塔板结构,改善气液传质效率,从而提升产品分离纯度并降低能耗。

常见误解与注意事项

开始使用本模拟器时,有几个容易陷入的误区需要注意。首先是“物性参数直接使用默认值即可”的想法,这很危险。例如,若用常温常压下的空气-水物性计算高温高压的蒸汽-水两相流,密度比和粘度比将完全不同,导致流型判定和压降计算出现严重偏差。请尽量输入接近实际工艺条件的参数,最好通过物性数据库查询具体数值。

其次是“贝克流型图的边界线是绝对标准”的误解。该图仅是基于典型实验数据的“参考基准”。实际管道倾斜度、内壁粗糙度、入口形状都会影响流型转变。即使模拟器显示“段塞流”,稍调整管道布局就可能转变为“气泡流”。建议将工具输出视为“一种可能性参考”,尤其当判定为危险流型时,应遵循安全侧设计原则。

最后是关于压力损失计算中“C常数”的选择。洛克哈特-马蒂内利关联式中的 $C$ 值虽随流型变化,但实际流动在过渡区域往往存在模糊性。例如,若计算位置处于“环状流”与“波状流”边界附近,工程实践中的铁律是分别用两种 $C$ 值计算并取更保守值(压降较大者)。工具虽可自动判定,但了解背后的原理能让你在决策时更有把握。

相关工程领域

两相流模拟技术支撑的领域范围可能超乎你的想象。首先它与“热流体工程”整体密不可分——尤其在涉及沸腾或冷凝的换热器设计中,流型直接支配传热系数。例如蒸发器从“核态沸腾”过渡到“膜态沸腾”时传热会急剧恶化(烧干现象),而这种转变与流动特性紧密相关。

其次是与“多相流体力学”及“粒子法模拟”的关联。本工具采用相对简单的关联式(属于“系统代码”级方法),若需微观追踪单个气泡的运动,则需使用VOF法或水平集法等界面追踪技术。反之,在针对大型系统(如全厂装置)的“工厂动态特性模拟器”中,正是此类实用关联式被嵌入各组件模型。这是典型的不同尺度选用不同工具的案例。

此外还延伸至“测量工程”领域。两相流流量测量难度大,如何测定气液比例(空隙率)及各相流速始终是挑战。若能通过模拟预测流型,便可为测量方案选择提供依据——例如“气泡流宜用超声波法,环状流适合光学测量”。设计与测量实为一体两面。

进阶学习指引

若想深入理解本工具背后的理论,建议首先学习“量纲分析”与“无量纲数”。为何贝克流型图的坐标轴呈现复杂形式?实为巧妙表征各相惯性力、粘性力与表面张力的比值。例如理解韦伯数 $We$ 与雷诺数 $Re$ 在两相流中的拓展形式,会彻底改变你对流型图的解读方式。

接下来可追溯“流型图谱”的发展脉络。贝克图仅是其中一例。针对管道朝向(水平/垂直/倾斜)与流体组合(可凝/不可凝),学界已提出曼达姆图、泰特尔-杜克勒图等多种图谱。比较各图谱的实验背景与侧重维度,能切身感受工程模型构建的现场感。

最终可挑战“动量守恒定律在各相及界面如何表述”这一核心命题。两相流基础方程“双流体模型”正是该思想的数学表达:分别建立气液两相运动方程,并通过“相互作用项”描述界面动量交换(如摩擦力)。洛克哈特-马蒂内利关联式等简洁公式,可视为对此复杂模型进行实用化“约简”的结果。理解至此,阅读新关联式论文时将不再畏惧。