流向旋转角形颗粒在惯性流中的尾迹转捩与标度律研究

《Journal of Fluid Mechanics》：Wake transitions and scaling laws of a streamwise-rotating angular particle in an inertial flow

　　

： 通过颗粒分辨直接数值模拟（PR-DNS），研究了中等雷诺数惯性流中流向旋转角形颗粒的尾迹动力学与流态转捩，重点关注旋转立方体。研究人员首先识别了前缘与后表面上涡生成的机制，表明旋转诱导的涡合并主导了若干关键流态转捩。研究揭示了一种此前未见报道的矛头状流态（spearhead regime），其特征为螺旋近尾迹与稳态丝状远尾迹结构。研究人员证明该流态源于边缘涡与中心旋转诱导涡的部分合并。流向旋转被发现能适度稳定尾迹，并诱导阻力系数\overline {C_d}与流向力矩\overline {T_x}随旋转速率\varOmega呈二次标度关系。升力系数的演化与回流区长度L_r的变化密切相关。对于旋转立方体，流向力矩数据以旋转雷诺数{\textit{Re}}{\varOmega }表达时坍缩至单一曲线。频率分析揭示，横向力的基频由旋转对称折数n{\varOmega }决定，而高次谐波则反映流-颗粒相互作用的非线性本质。研究人员进一步将分析推广至双角位形下的一般角形颗粒。处于对偶构型（sRCV，即一角迎流的旋转立方体；sROF，即一面迎流的旋转八面体）的颗粒可呈现几乎相同的流态转捩，因为它们具有相同的旋转对称折数与相当的涡对强度。水动力响应仍强烈依赖于颗粒形状与取向。其频率特征按n_{\varOmega }归类时同样坍缩。

论文解读文章：

一、研究背景与问题

颗粒负载流广泛存在于自然与工程系统中，从沉积物与污染物输运、流化床反应环境到生物医学应用均涉及颗粒与载流体的动量交换。对于固定钝体，当雷诺数超过临界值时，尾迹失去对称性并发生涡脱落，产生非定常载荷，进而诱发振动、疲劳与结构失效。旋转从根本上改变了这些动力学过程，可减弱或抑制涡脱落，而自由沉降颗粒的旋转响应则主导路径失稳并调节惯性悬浮体的集体行为。然而，尽管球体研究已相当充分，流场与旋转角形颗粒之间的相互作用仍相对未被探索。尖边与尖角引入了独特的涡量生成机制与尾迹失稳模式，增强了力与力矩，并产生了光滑形状中不存在的不同流态。

对于流向旋转的球体，已有研究识别了从轴对称流到低螺旋态、高螺旋态乃至近混沌态的转捩序列，并揭示了旋转对尾迹的稳定化作用。对于非轴对称体，椭球、扁盘及子弹形体的流向旋转研究均表明，适度旋转促进尾迹轴对称化与稳态化，而强旋转则触发次级失稳与更复杂的非定常动力学。相比之下，具有尖边的角形颗粒研究仍显不足。立方体作为典型的角形几何体，其静止绕流已建立了较完整的流态图谱，但流向旋转下其涡结构如何重组、力与力矩如何标度、频率内容由何决定，以及旋转对称性与多面体对偶性如何影响这些结论的普适性，均缺乏系统性认识。

本研究旨在回答一个核心问题：流向旋转如何重塑角形颗粒的尾迹拓扑，其机制是什么？并辅以三个支撑问题：尾迹重组如何体现在水动力响应及其标度关系中；非定常力的频率内容由什么决定、能否先验预测；旋转对称性、角形度与多面体对偶性如何决定这些发现的跨形状普适性。

二、研究内容与主要结论

研究人员采用颗粒分辨直接数值模拟（PR-DNS），系统研究了中等雷诺数下流向旋转立方体的尾迹动力学，并将分析推广至一般等距角形多面体及旋转球体。研究发现，流向旋转能适度稳定尾迹，并诱导阻力系数与流向力矩随旋转速率呈二次标度关系。研究揭示了一种此前未见报道的矛头状流态（spearhead regime），其近尾迹呈螺旋结构、远尾迹呈稳态丝状结构，源于边缘涡与中心旋转诱导涡的部分合并。升力系数的演化与回流区长度密切相关。旋转立方体的流向力矩数据以旋转雷诺数（{\textit{Re}}{\varOmega }）表达时坍缩至单一曲线。频率分析表明，横向力基频由旋转对称折数（n{\varOmega }）决定，高次谐波反映流-颗粒相互作用的非线性。对偶构型颗粒（sRCV与sROF）可呈现几乎相同的流态转捩，因其共享相同的旋转对称折数与相当的涡对强度。

三、关键技术方法

研究采用分布式拉格朗日乘子/虚拟域方法（DLM/FD）模拟刚性颗粒在牛顿流体中的运动。计算域为边长L_0=40的立方体，颗粒质心位于x_p=(10,20,20)，绕流向x轴旋转。所有角形颗粒具有相同体积，雷诺数（{\textit{Re}}）基于体积等效球直径（D{\textit{sph}}^*）定义。涡结构采用\lambda 2=-1准则识别。旋转速率以无量纲旋转率\varOmega及旋转雷诺数（{\textit{Re}}{\varOmega }=\varOmega {\textit{Re}}）表征。力与力矩通过拉格朗日乘子在颗粒域上的积分计算，频率特征以斯特劳哈尔数（\mathcal{S}tr）量化。研究考察了旋转对称折数n{\varOmega }、球形度\phi等参数的影响。

四、研究结果

1. 1.
   流向旋转立方体的流态图谱
   研究人员建立了流向旋转立方体的流态图，识别了若干关键流态转捩。旋转诱导的涡合并被证明是主导这些转捩的核心机制。随着旋转速率增加，尾迹经历从稳态轴对称到螺旋态、再到更复杂非定常态的转捩序列。流向旋转总体上有助于稳定尾迹，抑制或延迟涡脱落。
2. 2.
   矛头状流态的发现与机制
   研究揭示了一种此前未见报道的矛头状流态，其近尾迹呈螺旋结构、远尾迹呈稳态丝状结构。该流态源于边缘涡与中心旋转诱导涡的部分合并，这一部分合并机制是理解该流态形成的关键。该发现拓展了流向旋转角形颗粒尾迹拓扑的已知范围。
3. 3.
   水动力响应的标度行为
   阻力系数\overline {C_d}与流向力矩\overline {T_x}随旋转速率\varOmega呈二次标度关系。升力系数的演化与回流区长度L_r的变化紧密对应。旋转立方体的流向力矩数据以旋转雷诺数{\textit{Re}}_{\varOmega }表达时坍缩至单一曲线，表明该参数可有效统一不同旋转条件下的力矩响应。
4. 4.
   频率特征与旋转对称性
   频率分析表明，横向力的基频由旋转对称折数n_{\varOmega }决定，高次谐波则反映流-颗粒相互作用的非线性本质。这一发现为预测非定常力的频率内容提供了先验依据，即根据颗粒的旋转对称性即可预知其主导频率。
5. 5.
   对偶构型颗粒的普适性
   研究人员将分析推广至一般角形颗粒的双角位形。对偶构型（sRCV与sROF）颗粒因共享相同的旋转对称折数与相当的涡对强度，可呈现几乎相同的流态转捩。然而，水动力响应仍强烈依赖于颗粒形状与取向。其频率特征按n_{\varOmega }归类时同样坍缩，验证了旋转对称性作为统一组织原则的有效性。

五、讨论与结论

讨论部分强调，旋转诱导的涡合并是理解流向旋转角形颗粒尾迹转捩的统一机制框架。矛头状流态的发现表明，部分涡合并可产生独特的尾迹拓扑，丰富了已知流态谱系。二次标度关系为工程预测提供了简洁的定量工具。频率特征与旋转对称折数的关联则为非定常载荷预测提供了先验途径。对偶构型颗粒的相似行为揭示了多面体对偶性在流体动力学中的物理后果，即通过互换面与顶点改变迎流几何元素，可产生相似的尾迹组织。

研究结论指出，流向旋转适度稳定尾迹，诱导阻力与力矩的二次标度，矛头状流态源于边缘涡与旋转诱导涡的部分合并，横向力基频由旋转对称折数决定，对偶构型颗粒可呈现相似的流态转捩。这些发现为理解角形颗粒在一般旋转运动下的行为提供了综合框架，对颗粒负载流建模、沉降动力学预测及工程流体系统设计具有重要参考价值。未来工作可进一步探索更高雷诺数、更复杂旋转模式及多颗粒相互作用下的尾迹动力学。