细长旋成体绕流全局稳定性与流动不对称性机理研究

《Journal of Fluid Mechanics》：Investigation of the asymmetric flow field around a spheroid using global stability analysis

　　

研究人员针对细长旋成体绕流在低雷诺数条件下的流动失稳与非对称现象，采用全局线性稳定性分析方法，系统研究了长径比为6:1的椭球体在不同攻角下的流动演化规律。通过数值求解线性化Navier-Stokes方程的特征值问题，获得主导不稳定模态及其伴随模态，并结合结构灵敏度分析与基流涡量场特征，揭示了静止非对称模态与振荡脱落模态的竞争机制。研究发现，低攻角下流动失稳由静止非对称模态主导，其临界雷诺数随攻角增大而降低；中攻角区间存在两种对称基流解，导致流动对不对称的敏感性发生突变；高攻角下振荡脱落模态成为主导，非对称模态受到抑制。伴随模态分析表明，横向扰动可显著激发非对称响应，解释了实验中观察到的稳健流动不对称现象。该研究建立了钝体尾迹与攻角涡的不稳定性关联，为流动控制提供了理论依据。

研究背景与意义

细长旋成体（如潜艇、飞行器后体）的绕流非对称现象广泛存在于航空航天与海洋工程领域，表现为侧力突变与涡结构不对称，严重影响航行器的操纵性与稳定性。传统实验与模拟观察到此类现象，但其物理机制长期存在争议：一方认为源于分离线位移的几何效应，另一方主张与涡动力学不稳定性相关。由于高雷诺数下非线性效应复杂，线性理论能否解释此类现象尚不明确。本研究旨在通过全局线性稳定性分析，阐明不同攻角下流动失稳的物理机制，解决机制争议，并为流动控制提供理论支撑。研究成果发表于流体力学顶刊《Journal of Fluid Mechanics》。

关键技术方法

研究人员采用有限元法离散线性化Navier-Stokes方程，结合Newton-Raphson法求解基流，利用Krylov-Schur算法与shift-invert技术计算特征值谱。通过FreeFEM开源平台实现数值求解，采用MINI单元空间离散，MUMPS库进行线性系统求解。为处理下游边界影响，引入空间变化粘度模型抑制数值振荡。伴随模态通过离散伴随方法获得，结构灵敏度场定义为直接模态与伴随模态的内积模。计算网格针对各攻角与长径比独立生成，以精确捕捉尾迹与涡结构。

研究结果

长径比影响

零攻角下，长径比从1:1增至6:1，主导模态均为静止非对称模态，其临界雷诺数单调上升。流动可视化显示，随着长径比增大，回流区缩短，基流拖曳力修正因子减小，表明流动趋于稳定。结构灵敏度分析表明，失稳区域始终位于尾迹再循环区，与涡量梯度强相关。

攻角效应

建立6:1椭球体在攻角α∈[0°,90°]与雷诺数Re_D∈[90,3000]的稳定性图谱。低攻角（α≤45°）下，静止非对称模态主导失稳，临界雷诺数随攻角增大先降后稳。中攻角（45°<α<70°）出现两种对称基流解，分别对应有无静止对称模态，导致流动对不对称的敏感性发生阶跃式变化。高攻角（α≥70°）下振荡脱落模态成为主导，静止非对称模态被抑制。

低攻角流动结构

皮肤摩擦线显示，α=0°时为环形分离，α=5°时出现双侧分离线，α=10°时出现稳定焦点分离。轴向速度场与伴随模态分析表明，横向扰动可通过伴随模态高效激发非对称响应，解释了DARPA SUBOFF实验中观测到的稳健不对称现象。

中攻角双稳态机制

α=49.25°与49.5°的基流差异恰好对应静止对称模态形态，证实存在两种对称基流解。其中一种含对称模态的基流对不对称扰动更具鲁棒性，导致临界雷诺数突变。伴随对称模态的空间分布表明，沿流向加速流动可增强流动对对称扰动的稳定性。

高攻角涡结构演变

α=65°时背风面存在强下洗流，形成一对反向旋转的相干涡，结构灵敏度场与涡核重合。α=70°时下洗流减弱，无法维持相干涡对，静止非对称模态失去支撑，振荡脱落模态主导流动。

讨论与结论

研究人员指出，静止非对称模态本质上是涡动力学不稳定性，其驱动区域与涡量梯度及集中涡位置一致，支持了该模态源于涡相互作用而非分离线位移的观点。低攻角下，伴随模态的大幅值反映了流动的强对流型非正规性，解释了微小扰动即可触发显著不对称的实验现象。中攻角区的双稳态特性预示了雷诺数与攻角变化中可能出现的迟滞行为，为流动控制提供了新思路：通过注入对称模态可增强流动抗不对称能力。高攻角下相干涡对的消失导致非对称模态受抑，进一步验证了涡结构对非对称失稳的关键作用。本研究首次将钝体尾迹的非对称性与攻角涡不稳定性统一于同一理论框架，解决了长期存在的机制争议，并为水下航行器与飞行器的流动控制设计提供了定量依据。