[论文解读] On a separation criterion for symmetric elliptic bluff body flows

Kaushik Srinivasan|ArXiv.org|Nov 30, 2005

Fluid Dynamics and Vibration Analysis被引用 5

一句话总结

该论文提出了一种新的分析准则，通过分析圆柱体驻点附近流速（$u$）的曲率，以确定对称椭圆钝体的精确分离雷诺数（$Re_s$）。基于高分辨率的N-S方程模拟，发现$∂^2 u / ∂ ξ^2$在$Re_s$附近与$Re$呈线性关系，从而可精确外推得到圆柱体的$Re_s = 4.8$，解决了文献中长期存在的模糊性问题。

ABSTRACT

A new analytical criterion that captures the onset of separation of flow past elliptic cylinders is derived by considering the variation of the wall normal velocity in Reynolds number parameter space. It is shown that this criterion can be used to calculate the separation Reynolds number Re_s for the classical problem of flow past a circular cylinder, a contentious and unresolved issue till date. The two dimensional Navier-Stokes equations are solved computationally and an exact value of Re_s is obtained by applying the aforementioned criterion.

研究动机与目标

• 解决长期以来关于圆柱体绕流中分离雷诺数（$Re_s$）精确值的文献争议。
• 基于驻点处的局部流动行为，为一类对称椭圆圆柱体建立通用的分析准则。
• 识别出一个计算上易于处理的物理量——$∂^2 u / ∂ ξ^2$——其在$Re_s$附近与$Re$呈线性关系，从而可无需依赖视觉气泡检测，精确确定$Re_s$。
• 证明普朗特剪切应力准则不足以预测$Re_s$，因其仅在分离发生后适用。

提出的方法

• 基于雷诺数空间中法向速度的变化，提出以圆柱体驻点$P$为关注点的判据。
• 在椭圆柱坐标系（$ξ$-$\eta$）中，采用贴体、边界层加密的网格，求解二维不可压N-S方程。
• 对非线性项采用三阶迎风格式，对导数项采用中心差分格式，通过稳定化的双共轭梯度法求解流函数-涡度泊松方程。
• 在一系列雷诺数下计算驻点$P$处的$∂^2 u / ∂ ξ^2$，并观察其在$Re_s$附近的线性依赖关系。
• 通过将$∂^2 u / ∂ ξ^2$的线性趋势外推至零，确定$Re_s$，并利用在$P$点成立的导出关系$∂^2 u / ∂ ξ^2 \propto ∂ p / ∂ ξ$。
• 通过对比多种网格分辨率和外域尺寸的结果，验证方法的鲁棒性，确认线性关系与$Re_s$估计的稳定性。

实验结果

研究问题

• RQ1 对于圆柱体绕流，其分离雷诺数（$Re_s$）的精确值是多少？尽管文献中存在长期争议。
• RQ2 能否基于驻点处的局部流动行为，为对称椭圆钝体家族推导出一个通用的分析准则，以预测$Re_s$？
• RQ3 在$Re_s$附近，流速的二阶导数（$∂^2 u / ∂ ξ^2$）是否与$Re$呈线性关系，从而实现精确外推？
• RQ4 为何传统数值与实验方法在$Re \approx 6$以下难以一致地分辨分离的起始？
• RQ5 所提出的准则与普朗特剪切应力准则有何本质区别？为何后者无法用于预测$Re_s$？

主要发现

• 在驻点$P$处，流速的二阶导数$∂^2 u / ∂ ξ^2$在$Re_s$附近与雷诺数呈线性关系，从而可实现精确外推。
• 确定了圆柱体的$Re_s$精确值为4.8，解决了文献中长期存在的模糊性问题。
• $∂^2 u / ∂ ξ^2$在$Re_s$附近的线性行为在不同网格分辨率和外域尺寸下均保持稳健，$Re_s$值仅略有变化（在30倍与40倍圆柱直径的外域尺寸间差异约7%）。
• 尽管基本方程成立，驻点处的压力梯度$∂ p / ∂ ξ$在$Re_s$附近与$Re$并不呈线性关系，因此不适用于外推。
• 所提出的准则与普朗特剪切应力准则有本质不同，后者仅在分离发生后适用，无法预测$Re_s$。
• 该方法为避免在密集$Re$网格上进行N-S方程的暴力数值求解以检测分离起始，提供了一种计算高效的替代方案。