Hydrodynamic Behavior of Non-spherical Particles in Confined Vertical Flows: A Resolved CFD-DEM Study

本研究采用解析的CFD-DEM模拟来证明，由于形状诱导的尾迹不对称性，非球形多金属结核与等体积球体相比，其阻力显著增强且终点速度降低，同时揭示了颗粒尺寸和约束条件如何支配垂直水力输送过程中不同的阻力变化行为。

要点

想象一下，你正试图通过一根巨大的垂直吸管，将一堆奇怪、凹凸不平的岩石（被称为多金属结核）从海底运送到船上。这就是深海采矿的工作原理。工程师们面临的一个大问题是：当水流冲过吸管携带这些岩石向上移动时，这些形状怪异的岩石会如何表现？

大多数用于设计此类系统的计算机模型都将每块岩石视为一个完美的、光滑的弹珠。但在现实中，这些岩石是凹凸不平、形状不规则的，一点也不像弹珠。这篇论文探讨的是：将这种凹凸不平的岩石当作光滑的弹珠来处理，真的有效吗？还是说会给我们错误的答案？

为了找出答案，研究人员构建了一个超级精细的计算机模拟（就像一个高科技的视频游戏物理引擎），它没有使用任何捷径。他们没有去猜测水对岩石的推力，而是计算了水对岩石每一个凸起和缝隙的推力。

以下是他们发现的研究结果，用简单的语言进行了解释：

1. “凹凸不平的岩石” vs. “光滑的弹珠”

当研究人员将这些凹凸不平的岩石放入静水中观察其下沉速度时，发现这些凹凸不平的岩石比重量和大小完全相同的光滑弹珠下沉速度慢了约 28%。

• 类比： 想象你在水中游泳。如果你是一个线条流畅、流线型的海豚，你可以轻松滑行；如果你是一块凹凸不平、锯齿状的漂浮木头，你在向下移动的过程中会阻挡更多的水。
• 发生原因： 凹凸不平的岩石拥有更大的“迎风面积”（它们捕捉了更多的水），并且在身后留下了一个混乱、不对称的尾迹（就像一串混乱的气泡轨迹）。这种额外的阻力显著减慢了它们的移动速度。
• 关键点： 尽管它们下沉得更慢，但作用在它们身上的总向上推力（浮力）与弹珠是相同的。它们只是需要以更慢的速度运动来平衡那个力。

2. 管中的“交通拥堵”

接下来，他们模拟了水流冲过管道以携带这些岩石的过程。他们研究了两种尺寸：“小”岩石和“大”岩石。

• 光滑的弹珠： 当水速增加时，光滑的弹珠表现得很可预测。在低速时，它们会摇晃并沉降；在高速时，它们会像车辆汇入高速公路一样，笔直地向上飞驰。
• 凹凸不平的岩石： 它们的表现要混乱得多。
  • 在低速时： 小型的凹凸不平岩石甚至无法进入管道！它们在底部附近徘徊，原地旋转和摇晃，无法克服重力。然而，光滑的弹珠却能成功向上移动。
  • 在高速时： 即使水流速度足以携带它们，凹凸不平的岩石到达顶部所需的时间也更长，且移动路径更加飘忽不定。它们就像一群试图在旋转的扶梯上向上奔跑的人，而光滑的弹珠只是在直直地向上跑。

3. “陀螺效应”

研究人员观察了岩石旋转方式的最显著差异。

• 光滑的弹珠： 它们大多只是向上移动，几乎不怎么旋转。
• 凹凸不平的岩石： 由于它们表面凹凸不平，水流撞击它们会使它们剧烈旋转。这种旋转（转动）与其上下移动的运动紧密相关。
• 类比： 把光滑的弹珠想象成从枪里射出的子弹——它直行前进。把凹凸不平的岩石想象成被扔进风洞的飞镖或旋转陀螺。由于其形状的原因，它们会不断扭转、转向并改变方向。这种旋转产生了额外的“摩擦力”，使得运输变得更加困难。

4. “力波动”（颠簸的旅程）

研究人员测量了水给岩石带来的“推力”（阻力）。

• 小岩石： 无论是光滑的还是凹凸不平的，推力都相对稳定。
• 大岩石： 这里的情况变得非常疯狂。
  • 大型光滑弹珠： 当水流经过时，推力会有一些变化，产生一种可预测的“波动”模式。
  • 大型凹凸不平岩石： 推力是极度不可预测的。因为岩石在旋转，且相对于水流不断改变形状，推力会突然激增。这就像是在平坦的道路上开车（光滑弹珠）与在一条路面颠簸程度随车身倾斜角度每秒都在变化的道路上开车（凹凸不平的岩石）之间的区别。

核心结论

研究得出结论：虽然你可以使用光滑弹珠模型来获得一个大致的运行情况，但它忽略了细节。

• 如果你使用光滑弹珠模型，你可能会认为岩石向上移动的速度会比实际更快、更容易。
• 凹凸不平的岩石需要更高的水流速度才能开始移动，而且一旦它们开始移动，由于它们会旋转和摇晃，它们会更不稳定且更难控制。

简而言之：自然界是杂乱无章的。 如果你想设计出一台真正有效的机器，你就不能假装一个锯齿状的岩石是一个完美的球体。这种“凹凸不平”增加了许多额外的阻力和混乱，而简单的模型忽略了这一点。

技术摘要

技术摘要：受限垂直流中非球形颗粒的水动力行为

问题陈述
深海平原多金属结核（PMNs）的回收依赖于通过垂直上升管输送粗大、不规则且非球形颗粒的水力系统。这些系统的效率取决于预测在强重力作用和几何限制下的颗粒夹带、悬浮稳定性及输运行为。目前的工程设计实践通常依赖于未解析的计算流体力学-离散元法（CFD-DEM）模型，这些模型使用针对球形颗粒校准的经验阻力关联式，或通过球形度等标量描述符进行参数化。由于这些方法无法捕捉非球形颗粒固有的取向相关阻力、不对称壁面偏转尾迹以及旋转-平移耦合运动，因此在受限流（其中颗粒与管道直径比 $d/D \ge 0.2$）中表现不佳。此外，目前尚不清楚使用等体积球体替代真实的 PMN 形貌是否能保持系综层面的输运指标，或者是否会使对夹带阈值和压降的预测产生偏差。

研究方法
本研究采用全解析 CFD-DEM 框架，研究了代表深海采矿上升管的垂直管道中非球形 PMN 的水动力学特性。

• 数值框架： 该方法将不可压缩 Navier-Stokes 方程（通过 OpenFOAM 求解）与牛顿运动方程（通过 LIGGGHTS 求解）相结合，用于处理离散颗粒。
• 流固耦合： 使用浸没边界法（IB）直接解析界面处的流固相互作用，而不使用经验阻力定律。通过直接强迫法实现无滑移条件，并通过对颗粒表面进行流体应力张量积分来计算水动力。
• 颗粒表示： 基于 X 射线计算机断层扫描（CT）重建的复杂 PMN 几何形状，采用多球近似法（由重叠球形子颗粒组成的刚性集合体）进行建模。这使得在保持计算可行性的同时，能够高保真地表示不规则表面。
• 验证： 该框架通过了单球沉降（Ten Cate et al., 2002）和多球曳动-接触-翻滚（DKT）基准测试（Glowinski et al., 2001; Sharma and Patankar, 2005）的验证。此外，通过将球体的多球近似模型与光滑球体进行对比，验证了几何表示的准确性。
• 模拟设置： 模拟针对有效直径为 20 mm 和 54 mm 的 PMN 及其等体积球体（等体积直径分别为 16.4 mm 和 44 mm）在直径 $D=200$ mm 的管道中进行。研究覆盖了粒子雷诺数（$Re_p$）从 98 到 2904 以及限制比（$d/D$）为 0.082 和 0.22 的情况。输运行为在流速范围 $1.0u_t$ 至 $3.0u_t$（其中 $u_t$ 为终端沉降速度）内进行了分析。

关键结果

1. 沉降与阻力增强：

   • 尽管具有相同的质量和浮力，非球形 PMN 的沉降速度比等体积球体慢 27–29%。
   • 这种终端速度的降低对应于阻力系数较球体增强了 1.8 到 2.0 倍。
   • 阻力的增加源于两个因素：(1) 更大的瞬时迎风面积（由于外接球比率，比等体积面积大至 1.5 倍）以及 (2) 形貌诱导效应，包括不对称尾迹结构和旋转-平移耦合。
   • 至关重要的是，对于两种形状，终端阻力始终等于其淹没重量；PMN 较低的速度是由于在较低的滑移速度下即可产生所需的流体动力阻力。

2. 停留时间与输运机制：

   • 输运行为表现出随速度驱动的从间歇性、沉降主导运动到稳定、对流主导夹带运动的转变。
   • 在低流速（$u_f = 1.0u_t$）下，小型 PMN 处于入口附近的边缘悬浮状态，表现为无净向上输运的振荡运动，而球体则表现为垂直扩散。
   • 随着速度增加，所有颗粒的停留时间分布均变窄。然而，由于较高的阻力和更强的旋转-平移耦合，PMN 的平均停留时间和方差始终高于球体。
   • 在高流速（$u_f = 3.0u_t$）下，所有颗粒均达到空间均匀分布，但 PMN 保持着增强的旋转运动（角速度高达 80 rad/s）和螺旋轨迹。

3. 阻力统计与波动：

   • 平均归一化阻力（$\hat{f}$）在所有条件下均接近于 1，证实了阻力与淹没重量之间的稳态平衡。
   • 小型颗粒（$d/D = 0.082$）： 球体和 PMN 均表现出窄阻力分布。对于小型 PMN，快速的重新定向有效地平均了取向相关的投影面积波动，使其行为在定量上与球体相似。
   • 大型颗粒（$d/D = 0.22$）：
     • 大型球体表现出随速度增加而单调加宽的阻力分布，这是由尾迹不稳定性和限制效应驱动的。
     • 大型 PMN 则表现出显著加宽且偏斜的分布，并具有延长的强阻力尾部。这是由于尾迹历史效应与取向相关强迫作用的叠加所致。与小型颗粒不同，大型 PMN 的旋转时间尺度与尾迹演化时间尺度相当，这阻止了投影面积变化的平均化过程。

意义与结论
本研究阐明了颗粒形状如何影响受限垂直流中的阻力、尾迹动力学以及旋转-平移耦合。其主要贡献在于确立了应用等体积球形模型于非球形 PMN 输运时的定量界限。

• 定性相似性： 底层输运物理机制（从沉降到对流主导机制的演进）在 PMN 和球体之间具有定性相似性。这表明，使用校准后的阻力定律，球形模型可以捕捉到一阶输运行为。
• 定量差异： 存在显著的定量差异。PMN 的停留时间比球体长 40–90%，且阻力系数高出 1.8–2.0 倍。这些差异对于预测深海采矿上升管的最小悬浮速度和压降至关重要。
• 建模启示： 虽然球形近似可能足以用于系综行为的降阶建模，但它们无法捕捉与高限制比下的非球形大颗粒相关的特定间歇性和力波动。本文提出的解析 CFD-DEM 方法为开发此类模型的形状特定修正提供了必要的基准。

作者指出，这些发现适用于稀释系综，并作为校准输入；涉及阻碍沉降或塞流的稠密浆料的推断仍超出本研究范围。