可压缩层流剪切流中球形颗粒的升力模型研究

《Journal of Fluid Mechanics》：A lift model for spherical particles in unbounded compressible laminar shear flows

本文推荐研究针对可压缩剪切流中球形颗粒的升力预测难题，通过直接数值模拟和尺度分析，揭示了动态压力梯度是升力系数的关键影响因素，建立了考虑马赫数、雷诺数和尾流结构的升力模型，为火箭羽流和边界层等工程应用中的颗粒运动模拟提供了重要理论工具。

　　

在火箭发动机羽流、高超声速边界层等强剪切流动中，固体颗粒会受到显著的升力作用。这种升力可能来源于颗粒表面的不对称压力和剪切应力分布、粘性力不对称引起的旋转（马格努斯效应）以及靠近壁面时的不对称压缩（壁面效应）。然而，目前对于可压缩剪切流中球形颗粒的升力预测，特别是当流动同时存在密度梯度和速度梯度时，仍缺乏普适的模型。这限制了我们准确模拟颗粒在诸如燃烧、颗粒输运等复杂场景中的运动轨迹。

为了攻克这一难题，明尼苏达大学的Edward Feist和Graham Candler在《Journal of Fluid Mechanics》上发表了他们的最新研究。他们通过精细的直接数值模拟（DNS），系统研究了无界可压缩层流线性剪切流中固定球形颗粒的升力特性，旨在开发一个能够预测升力方向和量级的实用模型。

研究人员采用了一套先进的计算流体力学方法。他们使用US3D代码求解完美气体的可压缩Navier-Stokes方程。对于稳态流动，对流通量采用Modified-Steger-Warming (MSW) 迎风通量格式计算，粘性通量则采用带延迟修正的加权最小二乘误差梯度法。时间推进在近壁区域使用数据并行线松弛（DPLR）方法，在流场其他区域使用全矩阵数据并行（FMDP）方法。对于非稳态流动，他们采用了低耗散数值方法，结合中心动能一致（KEC）通量和耗散性MSW通量，并通过Larsson开关和数值海绵层来控制数值耗散，以准确捕捉激波和涡结构。计算网格经过精心设计，对于高超声速流动，网格还进行了激波对齐以减小数值噪声。研究系统考察了无量纲动态压力梯度（?Q）、马赫数（M_s）、雷诺数（Re_p）、壁面温度比（T_w/T_aw）以及尾流结构对升力系数的影响。

研究结果揭示了以下几个关键发现：

升力系数与动态压力梯度的尺度关系

研究发现，在稳态流动中，升力系数与无量纲动态压力梯度（?Q = D‖?Q_∞‖/Q_s）呈线性比例关系。升力的方向总是指向动态压力梯度相反的方向，即颗粒倾向于从动态压力较高的区域被“抬升”到动态压力较低的区域。这种关系是通过直接数值模拟和尺度分析共同验证的。

尾流结构对升力系数的影响

研究观察到，尾流结构（稳态轴对称SA、稳态平面对称SP或非稳态）对升力系数有显著影响。稳态平面对称尾流会产生比稳态轴对称尾流更大的升力系数，尤其是在低马赫数下。而非稳态尾流（如带有方位角振荡的hairpin尾流HaWAO）会产生难以预测方向的升力，甚至可能出现与梯度方向相反的升力。这表明预测升力必须先能预测尾流结构。

不同马赫数和雷诺数下的升力特性

研究还系统评估了马赫数、雷诺数和壁面温度比对升力的影响。结果表明，雷诺数对升力系数的影响最为显著，特别是在低马赫数下，因为它直接改变了尾流结构。马赫数的影响次之，而壁面温度比的影响相对较小。在较高雷诺数（Re_p≥ 200）和较高马赫数（M_s≥ 0.8）下，升力系数主要取决于动态压力梯度本身，而不依赖于该梯度是由密度梯度还是速度梯度单独产生。

稳态与非稳态流动的界定

基于大量模拟数据，研究提出了一个近似界定方法来区分稳态轴对称、稳态平面对称和非稳态尾流。该界定使用临界马赫数（M_ts和 M_tp）作为雷诺数的函数。当流动马赫数高于临界值时，尾流更可能呈现稳态行为。

可压缩剪切诱导升力模型的建立

综合以上发现，研究者最终建立了一个可用于颗粒跟踪模拟的升力模型。该模型首先判断流动是否为稳态（使用M_ts界定）。若非稳态，则假定升力为零；若为稳态，则计算升力方向（与动态压力梯度方向相反）和升力系数大小。升力系数被分解为压力升力系数（C_L,P）和粘性升力系数（C_L,τ）两部分，分别用包含马赫数、雷诺数、无量纲动态压力梯度以及尾流结构修正项的公式进行拟合。总升力系数为两者之和。该模型在多个盲验证案例中表现出良好的预测精度。

研究的结论强调，在无界可压缩层流剪切流中，作用于球形颗粒的升力主要由动态压力梯度驱动，方向与梯度相反。升力系数的大小不仅与此梯度线性相关，还强烈依赖于由马赫数和雷诺数共同决定的尾流结构。所开发的升力模型结合了这些关键物理因素，为可压缩剪切流中的颗粒运动模拟提供了重要的理论工具。然而，作者也指出了模型的局限性，例如它目前仅适用于连续流和线性剪切流，未考虑稀薄气体效应和强非线性剪切的影响。未来的研究可以在此基础上，进一步扩展模型以涵盖更广泛的流动条件。这项研究深化了对可压缩流中颗粒动力学的理解，对航空航天、能源环境等领域的颗粒输运过程具有重要的指导意义。