[论文解读] Polymer Maximum Drag Reduction: A Unique Transitional State

Yves Dubief, Christopher White+2 more|arXiv (Cornell University)|Jun 22, 2011

Fluid Dynamics and Turbulent Flows参考文献 20被引用 10

一句话总结

本文提出，聚合物最大阻力降低（MDR）是湍流中非线性失稳临界点处的一种独特过渡态，此时由于聚合物-涡旋相互作用破坏了壁湍流的自持循环，阻力降低达到峰值。通过直接数值模拟（DNS），研究发现MDR对应于失稳前与失稳后阶段之间的振荡，解释了实验中雷诺应力的差异，并为阻力降低提供了预测性框架。

ABSTRACT

The upper bound of polymer drag reduction is identified as a unique transitional state between laminar and turbulent flow corresponding to the onset of the nonlinear breakdown of flow instabilities.

研究动机与目标

• 解决关于最大阻力降低（MDR）是否为唯一确定状态，还是在不同实验中存在差异的争议。
• 研究MDR是否对应于特定的过渡流态，而非稳态湍流。
• 通过过渡边界层的DNS数据，检验MDR源于流动不稳定性非线性失稳临界点的假设。
• 通过将MDR处实验测量的雷诺应力差异与流动中的瞬态振荡行为关联，解释其原因。
• 通过将MDR识别为受聚合物-涡旋相互作用控制的过渡态，建立聚合物阻力降低的预测基础。

提出的方法

• 利用Wu & Moin（2009）关于边界层旁路过渡的直接数值模拟（DNS）数据，分析不同雷诺数下的速度与应力统计特性。
• 比较$Re_{\theta} = 200$、300和900下的平均速度剖面与雷诺应力，识别与Virk剖面（公式1）匹配的状态。
• 计算时间与空间平均统计量（$\overline{\bullet}$）以及同质方向上的时间依赖平均量（$\left\langle{\bullet}\right\rangle_{(t)}$），以评估MDR附近的振荡行为。
• 使用$Q$-准则（$Q > 0.1$）识别涡旋，并将涡旋体积分数$N_{Q>0.1}$作为空间平均壁面剪切$S(t)$的函数进行追踪。
• 分析阻力（通过$-dP/dx$）、速度与雷诺应力的时间演化，关联波动与涡旋活动及壁面剪切变化。
• 通过在充分发展流与入口流中进行思想实验，模拟MDR振荡并验证过渡态假设。

实验结果

研究问题

• RQ1 最大阻力降低（MDR）是否为唯一确定的状态，还是由于瞬态流动行为而存在差异？
• RQ2 MDR是否对应于向湍流过渡过程中的特定阶段，特别是非线性失稳的临界点？
• RQ3 为何在MDR处，尽管平均速度剖面一致，实验测量的雷诺应力仍存在显著差异？
• RQ4 MDR湍流的振荡特性是否可与流动不稳定性失稳前/后阶段相关联？
• RQ5 在过渡过程中，聚合物-涡旋相互作用如何导致阻力降低峰值，并解释观测到的实验变异性？

主要发现

• MDR状态下平均速度剖面与Virk剖面（$\overline{u}^{+}=11.7\log(y^{+})-17$）高度一致，证实其普适性。
• MDR处的雷诺剪切应力几乎可忽略（$\lesssim 10\%$的湍流水平），与Warholic等（1999）一致，但时间依赖统计量显示其偏差可达牛顿型过渡状态水平的$\sim 50\%$。
• 摩擦系数$C_f$在$Re_{\theta} \approx 200-300$时达到最小值，标志着从线性失稳到非线性失稳阶段的转变。
• MDR对应于非线性失稳临界点附近的狭窄过渡区域，此时涡旋开始形成并开始与聚合物相互作用。
• 时间序列分析显示，在MDR附近存在低阻力与高阻力状态之间的振荡，$-\left\langle{uv}\right\rangle_{(t)}$的最大负向与正向偏差与阻力波动密切相关。
• 涡旋体积分数$N_{Q>0.1}$与壁面剪切$S(t)$呈现顺时针轨道，证实了循环行为：涡旋形成导致剪切增加，随后被阻尼，剪切随之降低。