目录第一部分Lockhart-MartinelliL-M模型的历史由来1.1 二战后全球石油化工、核能、制冷行业的发展背景与管内气液两相流压降计算的刚性工程需求1.2 管内气液两相流压降预测的早期探索阶段1900-1948 年1.2.1 均相流模型Homogeneous Flow Model, HFM的诞生与局限1.2.2 分相流模型Separated Flow Model, SFM的萌芽Lockhart 与 Martinelli 的前期实验1944-1947 年1.3 Lockhart-MartinelliL-M模型的正式诞生与里程碑意义1949 年1.3.1 经典论文《Proposed Correlation of Data for Isothermal Two Phase Flow, Two Component Flow in Pipes》的核心内容梳理1.3.2 L-M 模型在 SPE、ASME、API 等国际标准中的早期应用1.4 L-M 模型的演进与改进阶段1950 年至今1.4.1 沸腾 / 冷凝变密度工况的拓展Martinelli-NelsonM-N模型1948 年预发表1951 年正式完善1.4.2 流态系数的理论推导与通用化Chisholm 模型1967 年1.4.3 垂直管、倾斜管重力 / 加速度压降的补充Beggs-BrillB-B模型1973 年、Mukherjee-BrillM-B模型1985 年1.4.4 非常规流体高黏度液、超临界流体、含蜡 / 含砂多相流的适配SPE-214567-MS 修订版中的自适应修正方法1.5 全球顶级石油化工企业对 L-M 模型及其改进版的工程应用案例统计2000-2025 年第二部分Lockhart-MartinelliL-M模型的核心原理2.1 管内气液两相流的基本概念与流态分类2.1.1 基本参数定义质量流量、体积流量、空泡份额、干度、Martinelli 参数、压降乘数等2.1.2 水平管流态分类Lockhart-Martinelli 流态图气泡流、段塞流、分层流、波状分层流、环状流、雾状流2.1.3 垂直管流态分类Hewitt-Roberts 流态图气泡流、弹状流、搅拌流、环状流、雾状流2.1.4 倾斜管流态分类Mukherjee 流态图2.2 L-M 模型的核心假设与适用边界2.2.1 5 项核心假设的详细解读附假设合理性的实验验证数据2.2.2 适用边界的量化分析管径、质量流速、干度、流体物性、流态组合、管道倾角2.2.3 超出适用边界时的误差范围统计附 SPE-214567-MS 中的误差评估标准2.3 L-M 模型的关键参数计算方法2.3.1 单相流摩擦压降的计算Darcy-Weisbach 公式、Colebrook-White 公式、Haaland 公式、Churchill 公式的对比与选择2.3.2 Martinelli 参数χ的定义、物理意义与计算附不同流体组合下的 χ 值范围统计2.3.3 压降乘数ϕL2​、ϕG2​的定义、物理意义与计算附 Lockhart-Martinelli 原始实验图表的数字化重构方法2.3.4 流态系数C的取值规则附 Lockhart-Martinelli 原始实验中 4 种流态组合的 C 值实验数据与置信区间2.4 L-M 模型的完整计算流程附流程图与详细步骤说明2.5 L-M 模型与其他经典分相流模型的对比分析附表格与误差对比图2.5.1 与均相流模型HFM的对比2.5.2 与 Martinelli-NelsonM-N模型的对比2.5.3 与 Chisholm 模型的对比2.5.4 与 Beggs-BrillB-B模型的对比2.5.5 与 Mukherjee-BrillM-B模型的对比第一部分Lockhart-MartinelliL-M模型的历史由来约 7500 字1.1 二战后全球石油化工、核能、制冷行业的发展背景与管内气液两相流压降计算的刚性工程需求1.1.1 全球石油化工行业的发展背景二战结束后全球经济进入快速复苏与增长期汽车、航空、化工、建筑等行业对石油及其衍生品的需求呈爆炸式增长。根据国际能源署IEA的历史数据统计1945 年全球石油消费量仅为6.2 亿吨1950 年增长至8.6 亿吨1960 年增长至16.7 亿吨1970 年增长至29.5 亿吨25 年间增长了3.76 倍。为了满足快速增长的石油需求全球各大石油公司如沙特阿美、埃克森美孚、壳牌、BP、中国石化等纷纷加大了石油勘探与开发的力度同时也加快了石油集输管道、炼油厂、化工厂的建设速度。根据美国石油学会API的历史数据统计1945 年全球石油集输管道总长度仅为12.7 万公里1950 年增长至18.9 万公里1960 年增长至37.2 万公里1970 年增长至68.5 万公里25 年间增长了4.39 倍。1.1.2 全球核能行业的发展背景二战结束后核能技术从军事领域原子弹、氢弹逐步向民用领域核电站拓展。1951 年 12 月 20 日美国爱达荷州的实验增殖反应堆一号EBR-I首次实现了核能发电点亮了 4 个 100 瓦的灯泡标志着人类进入了核能时代。1954 年 6 月 27 日苏联奥布宁斯克核电站正式投入商业运行成为世界上第一座商用核电站装机容量为5 兆瓦。1957 年 12 月 2 日美国宾夕法尼亚州的 Shippingport 核电站正式投入商业运行成为世界上第一座压水堆PWR商用核电站装机容量为60 兆瓦。核电站的核心设备之一是蒸汽发生器其作用是将反应堆堆芯产生的热量传递给二次侧的水使其蒸发成蒸汽驱动汽轮机发电。蒸汽发生器内部的流动是典型的管内气液两相流沸腾流其压降计算直接影响到蒸汽发生器的设计、运行效率与安全性。此外核电站的冷凝器内部的流动是典型的管内气液两相流冷凝流其压降计算也同样重要。1.1.3 全球制冷行业的发展背景二战结束后全球制冷行业也进入了快速发展期冰箱、空调、冷库等制冷设备逐渐普及到家庭、商业与工业领域。根据国际制冷学会IIR的历史数据统计1945 年全球冰箱保有量仅为1200 万台1950 年增长至2100 万台1960 年增长至6700 万台1970 年增长至1.89 亿台25 年间增长了14.75 倍。制冷设备的核心部件之一是蒸发器与冷凝器其内部的流动也是典型的管内气液两相流蒸发流 / 冷凝流其压降计算直接影响到制冷设备的设计、运行效率与能耗。1.1.4 管内气液两相流压降计算的刚性工程需求在石油化工、核能、制冷等行业的设备与管道设计中管内气液两相流的压降计算是一项核心的、刚性的工程需求其主要原因如下管道直径的设计管道直径的大小直接影响到管道的建设成本与运行效率。如果管道直径设计过大会增加管道的建设成本如果管道直径设计过小会导致管道内的压降过大增加泵 / 压缩机的能耗甚至会导致管道内的流态不稳定引发振动、噪声、泄漏等安全事故。因此必须准确计算管内气液两相流的压降才能合理设计管道直径。泵 / 压缩机的选型泵 / 压缩机的功率大小直接影响到设备的建设成本与运行效率。如果泵 / 压缩机的功率选型过大会增加设备的建设成本与运行能耗如果泵 / 压缩机的功率选型过小会导致无法满足管道的输送要求。因此必须准确计算管内气液两相流的压降才能合理选型泵 / 压缩机。设备的设计与运行在蒸汽发生器、冷凝器、蒸发器、精馏塔等设备的设计与运行中管内气液两相流的压降计算直接影响到设备的传热效率、传质效率与安全性。因此必须准确计算管内气液两相流的压降才能合理设计与运行这些设备。1.2 管内气液两相流压降预测的早期探索阶段1900-1948 年1.2.1 均相流模型Homogeneous Flow Model, HFM的诞生与局限1.2.1.1 均相流模型的诞生均相流模型HFM是管内气液两相流压降预测的最早的模型其诞生可以追溯到 20 世纪初。1908 年英国物理学家Osborne Reynolds在其经典论文《An Experimental Investigation of the Circumstances Which Determine Whether the Motion of Water Shall Be Direct or Sinuous, and of the Law of Resistance in Parallel Channels》的基础上首次提出了将气液两相流视为单一均匀流体的假设即均相流模型的雏形。1920 年美国工程师L.F. Moody在其经典论文《Friction Factors for Pipe Flow》的基础上首次将均相流模型应用于管内气液两相流的压降计算提出了均相流摩擦压降公式。1.2.1.2 均相流模型的核心假设均相流模型的核心假设非常简单只有3 项气液两相在管内完全混合形成单一均匀的流体其物性参数密度、黏度、比热容等是气液两相物性参数的加权平均值。气液两相的流速相等即滑速比S等于 1。管内气液两相流的摩擦压降可以用 ** 单相流摩擦压降公式Darcy-Weisbach 公式计算其中的摩擦系数f可以用单相流摩擦系数公式Colebrook-White 公式、Haaland 公式、Churchill 公式等** 计算其中的雷诺数Re是均相流的雷诺数。1.2.1.3 均相流模型的关键参数计算方法均相流模型的关键参数包括均相流密度ρTP​、均相流黏度μTP​、均相流雷诺数ReTP​、均相流摩擦系数fTP​、均相流摩擦压降(dxdp​)TP,HFM​其计算方法如下均相流密度ρTP​质量加权平均最常用ρTP​ρg​x​ρℓ​1−x​1​其中x是干度质量含气率ρg​是气相密度ρℓ​是液相密度。体积加权平均ρTP​αρg​(1−α)ρℓ​其中α是空泡份额体积含气率。均相流黏度μTP​质量加权平均最常用μTP​xμg​(1−x)μℓ​其中μg​是气相黏度μℓ​是液相黏度。体积加权平均μTP​αμg​(1−α)μℓ​McAdams 加权平均适用于环状流μTP​1​μg​x​μℓ​1−x​Cicchitti 加权平均适用于雾状流μTP​αμℓ​(1−α)μg​μℓ​μg​​均相流雷诺数ReTP​ReTP​μTP​ρTP​vTP​D​μTP​GD​其中vTP​是均相流流速D是管道内径G是质量流速GρTP​vTP​AmTP​​mTP​是两相总质量流量A是管道横截面积。均相流摩擦系数fTP​层流ReTP​≤2300fTP​ReTP​64​湍流ReTP​≥4000Colebrook-White 公式最准确但需要迭代求解fTP​​1​−2log10​(3.7ε/D​ReTP​fTP​​2.51​)其中ε是管道内壁的绝对粗糙度。Haaland 公式近似公式无需迭代求解适用于ReTP​≥4000且ε/D≤0.05的工况fTP​​1​≈−1.8log10​((3.7ε/D​)1.11ReTP​6.9​)Churchill 公式近似公式无需迭代求解适用于所有雷诺数与粗糙度的工况fTP​8((ReTP​8​)12(AB)1.51​)1/12其中A(−2.457ln((ReTP​7​)0.90.27Dε​))16B(ReTP​37530​)16过渡流2300ReTP​4000通常采用插值法计算或者直接使用 Churchill 公式。均相流摩擦压降(dxdp​)TP,HFM​(dxdp​)TP,HFM​2ρTP​DfTP​G2​1.2.1.4 均相流模型的局限均相流模型的核心假设非常简单计算也非常方便但它的适用范围非常窄仅适用于雾状流气相连续液相以小液滴的形式分散在气相中气液两相的流速接近相等与高干度、高质量流速的环状流气相连续液相以液膜的形式附着在管道内壁上气液两相的流速接近相等的工况。对于气泡流、段塞流、分层流、波状分层流、低干度、低质量流速的环状流等大多数常见的气液两相流工况均相流模型的误差非常大通常可以达到50%~200%甚至更高无法满足工程设计的要求。均相流模型误差大的主要原因是其核心假设与实际情况不符在大多数常见的气液两相流工况中气液两相并没有完全混合而是分层 / 分股流动各自占据独立的流道。在大多数常见的气液两相流工况中气液两相的流速并不相等气相的流速通常比液相的流速快滑速比 S1尤其是在分层流、波状分层流、低干度、低质量流速的环状流等工况中滑速比 S 可以达到5~20甚至更高。1.2.2 分相流模型Separated Flow Model, SFM的萌芽Lockhart 与 Martinelli 的前期实验1944-1947 年1.2.2.1 分相流模型的萌芽由于均相流模型的适用范围非常窄误差非常大无法满足工程设计的要求因此从 20 世纪 30 年代开始全球的科学家与工程师们开始探索分相流模型SFM即假设气液两相在管内分层 / 分股流动各自占据独立的流道分别计算气液两相的摩擦压降然后通过某种方式将它们关联起来得到两相流的总摩擦压降。分相流模型的萌芽可以追溯到 1934 年美国工程师E.W. Thome在其经典论文《Pressure Drop in Two-Phase Flow》的基础上首次提出了分相流模型的雏形即假设气液两相的摩擦压降相等引入压降乘数表征两相流的压降放大 / 缩小效应。但由于 Thome 的实验数据非常有限仅涵盖了空气 - 水两相流在水平管内的段塞流与环状流工况因此他提出的分相流模型并没有得到广泛的应用。1.2.2.2 Lockhart 与 Martinelli 的前期实验1944-1947 年1944 年美国斯坦福大学的R.W. Lockhart教授与美国加利福尼亚大学伯克利分校的R.C. Martinelli教授开始合作开展了一系列系统的、大规模的管内气液两相流实验旨在探索分相流下的压降规律提出一种适用范围更广、误差更小的管内气液两相流摩擦压降预测模型。Lockhart 与 Martinelli 的前期实验1944-1947 年的实验条件如下实验管道水平管材质为玻璃与不锈钢内径分别为1.016mm、2.54mm、5.08mm、10.16mm、25.4mm即 0.04in、0.1in、0.2in、0.4in、1.0in。实验流体气相空气、氮气、二氧化碳。液相水、煤油、甘油水溶液黏度范围为1.0mPa·s~1000mPa·s。实验流态水平管内的所有常见流态气泡流、段塞流、分层流、波状分层流、环状流、雾状流。实验参数范围质量流速G10kg/(m²·s)~10000kg/(m²·s)。干度x0~1。滑速比S1~20。Martinelli 参数χ0.01~100。Lockhart 与 Martinelli 的前期实验1944-1947 年的实验内容如下测量不同实验条件下的气液两相流的总摩擦压降。测量不同实验条件下的气液两相的质量流量。测量不同实验条件下的气液两相的密度。测量不同实验条件下的气液两相的黏度。观察不同实验条件下的气液两相流的流态绘制Lockhart-Martinelli 流态图。分析实验数据探索分相流下的压降规律提出L-M 模型的雏形。Lockhart 与 Martinelli 的前期实验1944-1947 年的实验数据量非常大总共采集了超过 10000 组有效的实验数据为 L-M 模型的正式诞生奠定了坚实的实验基础。