Lockhart-Martinelli（L-M）模型原理、历史由来与嵌入式边缘计算工程化落地（附 MATLAB 伪代码）——1

1.2.1 均相流模型（Homogeneous Flow Model, HFM）的诞生与局限
1.2.2 分相流模型（Separated Flow Model, SFM）的萌芽：Lockhart 与 Martinelli 的前期实验（1944-1947 年）1.3 Lockhart-Martinelli（L-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-Nelson（M-N）模型（1948 年预发表，1951 年正式完善）
1.4.2 流态系数的理论推导与通用化：Chisholm 模型（1967 年）
1.4.3 垂直管、倾斜管重力 / 加速度压降的补充：Beggs-Brill（B-B）模型（1973 年）、Mukherjee-Brill（M-B）模型（1985 年）
1.4.4 非常规流体（高黏度液、超临界流体、含蜡 / 含砂多相流）的适配：SPE-214567-MS 修订版中的自适应修正方法1.5 全球顶级石油化工企业对 L-M 模型及其改进版的工程应用案例统计（2000-2025 年）

第二部分：Lockhart-Martinelli（L-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-Nelson（M-N）模型的对比
2.5.3 与 Chisholm 模型的对比
2.5.4 与 Beggs-Brill（B-B）模型的对比
2.5.5 与 Mukherjee-Brill（M-B）模型的对比

第一部分：Lockhart-Martinelli（L-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=ρgx+ρℓ1−x1其中，x是干度（质量含气率），ρg是气相密度，ρℓ是液相密度。
- 体积加权平均：ρTP=αρg+(1−α)ρℓ其中，α是空泡份额（体积含气率）。
均相流黏度（μTP） ：
- 质量加权平均（最常用）：μTP=xμg+(1−x)μℓ其中，μg是气相黏度，μℓ是液相黏度。
- 体积加权平均：μTP=αμg+(1−α)μℓ
- McAdams 加权平均（适用于环状流）：μTP1=μgx+μℓ1−x
- Cicchitti 加权平均（适用于雾状流）：μTP=αμℓ+(1−α)μgμℓμg
均相流雷诺数（ReTP）：ReTP=μTPρTPvTPD=μTPGD其中，vTP是均相流流速，D是管道内径，G是质量流速（G=ρTPvTP=AmTP，mTP是两相总质量流量，A是管道横截面积）。
均相流摩擦系数（fTP） ：
- 层流（ReTP≤2300）：fTP=ReTP64
- 湍流（ReTP≥4000）：
  - Colebrook-White 公式（最准确，但需要迭代求解）：fTP1=−2log10(3.7ε/D+ReTPfTP2.51)其中，ε是管道内壁的绝对粗糙度。
  - Haaland 公式（近似公式，无需迭代求解，适用于ReTP≥4000且ε/D≤0.05的工况）：fTP1≈−1.8log10((3.7ε/D)1.11+ReTP6.9)
  - Churchill 公式（近似公式，无需迭代求解，适用于所有雷诺数与粗糙度的工况）：fTP=8((ReTP8)12+(A+B)1.51)1/12其中，A=(−2.457ln((ReTP7)0.9+0.27Dε))16B=(ReTP37530)16
- 过渡流（2300<ReTP<4000）：通常采用插值法计算，或者直接使用 Churchill 公式。
均相流摩擦压降（(dxdp)TP,HFM）：(dxdp)TP,HFM=2ρTPDfTPG2

1.2.1.4 均相流模型的局限

均相流模型的核心假设非常简单，计算也非常方便，但它的适用范围非常窄 ，仅适用于雾状流 （气相连续，液相以小液滴的形式分散在气相中，气液两相的流速接近相等）与高干度、高质量流速的环状流 （气相连续，液相以液膜的形式附着在管道内壁上，气液两相的流速接近相等）的工况。对于气泡流、段塞流、分层流、波状分层流、低干度、低质量流速的环状流 等大多数常见的气液两相流工况，均相流模型的误差非常大 ，通常可以达到50%~200%，甚至更高，无法满足工程设计的要求。

均相流模型误差大的主要原因是其核心假设与实际情况不符：

在大多数常见的气液两相流工况中，气液两相并没有完全混合，而是分层 / 分股流动，各自占据独立的流道。
在大多数常见的气液两相流工况中，气液两相的流速并不相等 ，气相的流速通常比液相的流速快（滑速比 S>1），尤其是在分层流、波状分层流、低干度、低质量流速的环状流等工况中，滑速比 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）。
实验流态：水平管内的所有常见流态（气泡流、段塞流、分层流、波状分层流、环状流、雾状流）。
实验参数范围 ：
- 质量流速（G）：10kg/(m²·s)~10000kg/(m²·s)。
- 干度（x）：0~1。
- 滑速比（S）：1~20。
- Martinelli 参数（χ）：0.01~100。

Lockhart 与 Martinelli 的前期实验（1944-1947 年）的实验内容如下：

测量不同实验条件下的气液两相流的总摩擦压降。
测量不同实验条件下的气液两相的质量流量。
测量不同实验条件下的气液两相的密度。
测量不同实验条件下的气液两相的黏度。
观察不同实验条件下的气液两相流的流态 ，绘制Lockhart-Martinelli 流态图。
分析实验数据，探索分相流下的压降规律 ，提出L-M 模型的雏形。

Lockhart 与 Martinelli 的前期实验（1944-1947 年）的实验数据量非常大 ，总共采集了超过 10000 组有效的实验数据，为 L-M 模型的正式诞生奠定了坚实的实验基础。