3. 欧拉-欧拉双流体模型基础

欧拉-欧拉双流体模型（Eulerian-Eulerian Two-Fluid Model）是鼓泡塔、流化床、气升式反应器等气液或气固两相工业装置仿真的主力方法。它的核心思想是：不追踪单个气泡或颗粒，而是把每一相都视为充满整个计算域的连续介质，两相在同一空间网格上共存，各自拥有独立的体积分数、速度、压力和温度场，通过相间作用力耦合。理解这套模型，关键要抓住四个概念：相体积分数、相间动量传递、曳力与升力/壁面润滑力、以及面向鼓泡塔和流化床的建模要点。

相体积分数：两相共存的“占位”约束

相体积分数 α_q 表示第 q 相在控制体内占据的体积比例，满足 Σα_q = 1。气相体积分数 α_g 在鼓泡塔中就是气含率，在流化床中则直接决定床层处于固定床、鼓泡流化还是快速流化状态。Fluent 中每相都求解一个连续性方程：

∂(α_q ρ_q)/∂t + ∇·(α_q ρ_q u_q) = S_q

其中 Σα_q = 1，S_q 为相变或反应引起的质量源项

体积分数不是被动结果，它反过来影响物性：混合密度 ρ = Σα_qρ_q，混合黏度同理。因此 α_g 的空间分布（塔中心高、近壁低，或流化床的床层膨胀高度）本身就是最重要的仿真输出之一。工程上常把 α_g 的径向分布与实验测得的“壁面峰值”或“中心峰值”对比，作为判断流型是否合理的依据。

相间动量传递：两相如何“互相拖拽”

两相各自求解动量方程，相间通过界面作用力交换动量。以气相为例，其动量方程可写为：

∂(α_g ρ_g u_g)/∂t + ∇·(α_g ρ_g u_g u_g)
  = -α_g ∇p + ∇·(α_g τ_g) + α_g ρ_g g + F_inter

F_inter = F_drag + F_lift + F_wall + F_vm + F_td

其中 F_drag 是曳力，F_lift 是升力，F_wall 是壁面润滑力，F_vm 是虚拟质量力，F_td 是湍流扩散力。对鼓泡塔和流化床，曳力通常占主导，升力与壁面润滑力决定径向分布，虚拟质量力在气液体系中不可忽略，而湍流扩散力影响相分布弥散程度。这些力不是全部都要开，而是按物理重要性取舍——这是欧拉模型调参的核心。

曳力模型：决定相间滑移的关键

曳力是两相速度差的函数，F_drag = K_gl(u_g − u_l)，其中 K_gl 是相间动量交换系数。K_gl 由曳力系数 C_D 和相界面面积浓度决定：

K_gl = (3/4) · C_D · (α_g α_l ρ_l |u_g - u_l|) / d_b

C_D 由曳力模型给出，d_b 为气泡/颗粒直径

Fluent 提供多种曳力模型，选择依据是流型与相含率：

以鼓泡塔为例，若表观气速 0.05 m/s、气泡直径 5 mm、α_g 约 0.1，用 Grace 或 Ishii-Zuber 通常比 Schiller-Naumann 更接近实验气含率。流化床则常用 Gidaspow，因为床层密相区颗粒浓度高，Ergun 项主导压降。

升力与壁面润滑力：径向分布的控制者

升力 F_lift = C_L α_g ρ_l (u_g − u_l) × (∇ × u_l)，方向垂直于相对速度与液相涡量。在鼓泡塔中，升力把气泡从高速区推向低速区，是形成“壁面峰值”气含率分布的主要原因。升力系数 C_L 默认 0.5，但对变形气泡可降至 0.1 甚至变号，需按流型调整。

壁面润滑力 F_wall 是近壁排斥力，防止气泡穿透壁面，形式为 F_wall = −C_w α_g ρ_l |u_r|² n_w，其中 n_w 是壁面法向。它把气泡推离壁面，与升力共同塑造近壁 α_g 剖面。若忽略壁面润滑力，近壁气含率会异常偏高，塔内循环结构失真。

鼓泡塔与流化床建模要点

鼓泡塔：气相为离散相、液相为连续相，α_g 通常 0.05～0.2。建议开启曳力（Grace/Ishii-Zuber）、升力、壁面润滑力；虚拟质量力在瞬态启动阶段开启；湍流用 k-ε 或 RNG k-ε 加分散相湍流模型。气泡直径可用常数或群体平衡模型（PBM）描述聚并破碎。入口用速度入口给定表观气速，出口用压力出口，壁面用无滑移。

流化床：气固两相，颗粒相需设置颗粒压力与颗粒黏度（如 Gidaspow 或 Syamlal-O'Brien 本构）。曳力用 Gidaspow 或 Syamlal-O'Brien，升力与壁面润滑力通常影响较小可关闭。重点监控床层压降、最小流化速度 U_mf、床层膨胀高度与气泡行为。时间步长建议 1e-3～1e-4 s，先稳态再瞬态。

掌握以上四点，就能为鼓泡塔和流化床搭建物理上合理的欧拉双流体模型，并为后续的曳力模型对比、网格无关性验证和实验标定打下基础。