后流（尾流）
后流（尾流）的理论基础
概要
老师，后流是指物体后面形成的流动对吧？
完全正确。当物体在流体中运动（或流体流过物体）时，在物体下游侧形成的速度缺失区域被称为后流（wake）。飞机尾流乱气流、风力发电的尾流干扰、汽车空气阻力、运动员的滑流等，应用范围极其广泛。
后流的构造
让我们考虑物体充分远下游的远场后流。后流的速度分布是均流速度 $U_\infty$ 减去速度缺失 $u_{def}(x, y)$。
在远场后流中，可以假设 $u_{def} \ll U_\infty$，因此可以应用线性化的边界层方程。
自相似解
后流也存在自相似解吗？
远场后流的自相似性已由Townsend理论确立。
二维后流（圆柱等二维物体）：
轴对称后流（球等三维物体）：
喷流是 $u_c \propto x^{-1}$，后流的衰减比喷流慢呢。
正是如此。后流是均匀来流中的一个"空洞"，不像喷流那样具有强的自诱导作用，因此扩散较慢。
动量积分与阻力
我听说可以从后流推导出物体的阻力。
这是一个极其重要的关系。对后流的动量缺失进行积分可以得到物体的阻力。
二维物体情况（单位展长）：
轴对称物体情况：
这个关系也被称为Jones公式，是风洞实验中非接触测量物体阻力的基础。通过皮托管遍历获得速度分布，然后用动量积分计算阻力。
不需要直接测量力就能知道阻力。
正确。这是动量守恒定律的必然结果，也在CFD中作为一种计算技巧被使用——从控制体积的动量收支计算阻力。壁面压力与摩擦力的积分计算结果应当相一致。
后流的稳定性
后流的稳定性分析也很重要吗？
通过后流速度分布的稳定性分析，可以解释卡门涡街的特性。对后流速度分布进行时间稳定性分析，
- 反对称模式（sinuous mode）：涡街的蛇形运动。对应卡门涡
- 对称模式（varicose mode）：后流宽度的脉动。通常比反对称模式的不稳定性弱
Monkewitz (1988) 表明，当后流处于"绝对不稳定"状态（速度缺失足够大）时，后流会自励振荡，卡门涡街会自发形成，不需要来自上游的扰动。
卡车队列行驶与后流的节能——滑流效应的计算
大型卡车在高速公路上编队行驶时，后续车辆进入前车的后流（尾流）区域，空气阻力会降低20～30%。这就是"滑流"效应，也是自动驾驶车队行驶的节能基础，实际应用正在推进。后流速度缺失能恢复到多远（后流回复长度）取决于雷诺数和物体形状，直接关系到卡车间最优间距的设计。CFD采用"将前车后流的速度分布作为后续车的入口边界条件"进行耦合解析，计算表明随着连接车数增加，后续车的节能效果也会增加。
后流（尾流）的数值计算手法
数值手法的选择
后流的CFD用什么手法？
后流的解析有两个方面：物体自身的分析和后流区域的分析。
| 目的 | 手法 | 备注 |
|---|---|---|
| 物体近旁的分离与近场后流 | RANS / DES / LES | 需要壁面解析 |
| 远场后流的扩散与恢复 | RANS / LES | 需要广阔的计算域 |
| 后流的稳定性分析 | DNS + Floquet / BiGlobal | 以精密计算的基态为前提 |
| 风力发电尾流干扰 | Actuator Line/Disk + LES | 风机用模型表示，集中计算后流 |
后流域的网格设计
后流区域的网格应注意什么？
后流向下游扩展，网格也必须追踪这种扩展。
- 物体直后：最细的网格。直到再循环区长度范围内的解析度与物体表面网格相当
- 中间后流（$5D\text{--}20D$）：涡结构的崩壊过程。流向上网格逐步粗化（增长率 $< 1.1$）
- 远场后流（$> 20D$）：自相似区域。后流宽度处至少配置10个网格单元
- 横向：后流宽度的3倍以上的区域
数值扩散导致后流过早消散怎么解决？
后流的速度缺失向下游逐步减小，容易受数值扩散影响。对策是，
1. 高阶精度格式：至少二阶精度。LES用中心差分系
2. 网格等向性：流向不要过度拉伸网格。长宽比 $< 5$
3. 充分解析：即使速度缺失低于 $1\%$ 的区域，也应保证能分解缺失分布
4. AMR（自适应网格加密）：基于涡度或速度梯度动态细化网格
动量积分法计算阻力
请教我CFD中从后流计算阻力的方法。
在物体下游充分远处的截面（如下游 $10D$ 处）提取速度分布，执行动量积分。
第二项压力项在远场较小，但在物体近场不可忽视。用这种方法获得的阻力应与壁面压力和摩擦力直接积分的结果一致，这种一致性是CFD验证的良好实践。
OpenFOAM 下的后流解析
用 OpenFOAM 获取后流统计量怎么做？
使用 fieldAverage function object 计算时间平均场。
```
functions
{
fieldAverage1
{
type fieldAverage;
libs ("libfieldFunctionObjects.so");
writeControl writeTime;
fields
(
U { mean on; prime2Mean on; base time; }
p { mean on; prime2Mean on; base time; }
);
}
}
```
这样会输出 UMean（时间平均速度）和 UPrime2Mean（雷诺应力张量）。后流的速度缺失分布从 UMean 减去 $U_\infty$ 得到。
后流的"自相似性"——走再远形状还是一样的流动
后流（尾流）最有趣的性质之一是"自相似性"。在物体充分下游的远场后流区域，速度缺失分布收敛为高斯分布形状，改变尺度后与前面的截面形状相同。这种自相似解称为Gaussian far-wake solution，可以从理论推导得出。测量位置改变时，用相似参数重新标度，所有曲线都会重合在一条线上。要精确追踪后流的数值计算需要"足够长的计算域（物体后方直径的50～100倍）"。如果只在少数几个截面评估后流，可能在解析完成前都无法确认自相似解的收敛。
后流（尾流）的实务应用
实务下的后流解析
什么时候会在实务中需要后流解析？
主要应用场景列举如下。
| 应用 | 关注的量 | Re范围 | 典型手法 |
|---|---|---|---|
| 飞机尾流乱气流 | 翼尖涡强度与持续距离 | $10^7$ | RANS + 涡追踪 |
| 风力发电尾流干扰 | 速度缺失率、湍流附加 | $10^7$ | Actuator Line + LES |
| 汽车空力 | 后流区压力回复、$C_D$ | $10^6$ | RANS (SST) / DES |
| 建筑物周风环境 | 行人高度风速 | $10^6\text{--}10^7$ | RANS (k-ε) / LES |
| 海洋结构物（VIV） | 涡脱频率、横向力 | $10^4\text{--}10^6$ | URANS / DES |
风力发电尾流的解析
风力发电的尾流干扰是大问题吗？
这是一个巨大的问题。上游风机的尾流导致下游风机接收的风速降低 $20\text{--}40\%$。输出几乎与 $U^3$ 成正比，所以速度降低以3次方影响发电量。
尾流解析方法有3个阶段。
1. 工程模型：Jensen/Frandsen模型。假设尾流线性扩展。用于风电场布置优化
2. Actuator Disk/Line + RANS：用源项表示风机。单个风机尾流预测
3. Actuator Line + LES：用线源表示叶片，解析后流的湍流结构。最精确但成本大
Jensen模型是什么公式？
最简单的尾流模型。
其中 $C_T$ 是推力系数，$k_{wake}$ 是尾流扩展率（陆地约 $0.075$，海上约 $0.04$），$D$ 是转子直径。
后流的回复距离
后流恢复到多远？
取决于物体类型和流动条件。参考值如下。
| 物体 | 后流回复95%的距离 | 备注 |
|---|---|---|
| 圆柱 | $50\text{--}100D$ | 亚临界区 |
| 球 | $30\text{--}50D$ | 三维因此恢复快 |
| 风力涡轮机 | $10\text{--}20D$ | 大气湍流加速混合 |
| 飞机（翼尖涡） | 数km～十几km | 持续性极强 |
飞机尾流这么长啊。
翼尖涡的涡核直径小，涡度集中，粘性散逸需要很长时间。大型飞机的尾流对后续飞机构成严重威胁，因此管制间隔与翼尖涡崩溃预测（Crow不稳定性）的研究很活跃。
雷诺的实验（1883年）——发现湍流的时刻
奥斯博恩·雷诺做过一个经典实验：在管中的水里加入墨水，随着流速增加，墨水线在某一刻突然变得混乱。这就是层流向湍流的过渡。雷诺用数学方式表达这一时刻为无量纲数 $Re = \rho uD/\mu$。过了一百多年，CFD工程师仍然最先检查的就是这个雷诺数。
后流（尾流）的软件比较
工具别的特点
哪些CFD工具适合后流解析？
根据应用领域的工具选择。
| 应用 | 推荐工具 | 原因 |
|---|---|---|
| 汽车空力 | Fluent, STAR-CCM+, PowerFLOW | DES/LES对应、自动网格 |
| 风力发电尾流 | SOWFA (OpenFOAM), EllipSys3D | 已实装Actuator Line |
| 飞机尾流 | OVERFLOW, FUN3D, elsA | NASA/航宇高精度代码 |
| 通用后流解析 | OpenFOAM, Fluent | 适用范围广 |
| 建筑风工学 | WindSim, Fluent, OpenFOAM | 大气边界层建模 |
STAR-CCM+ 下的汽车后流解析
听说汽车后流解析多用STAR-CCM+。
STAR-CCM+ 在汽车CFD中实绩丰富。特点是，
- 自动网格：多面体网格 + 棱柱层自动生成。适应复杂车身形状
- IDDES：Improved Delayed DES。直接解析后流大涡结构
- Surface Wrapper：脏CAD数据清理。实务中很重要
- Design Manager：形状参数自动优化
汽车的 $C_D$ 预测关键在于后流区压力恢复。用DES可以捕捉RANS无法再现的后流非定常结构（后柱涡、C柱涡等），提高 $C_D$ 预测精度。
OpenFOAM (SOWFA) 下的风力发电尾流
SOWFA是什么？
Simulator fOr Wind Farm Applications 的缩写，由美国国家可再生能源实验室（NREL）开发，基于 OpenFOAM 的工具。在大气边界层LES中内置Actuator Line模型，可以解析整个风电场的尾流干扰。
主要特点是，
- 大气稳定度（稳定/中立/不稳定）的模拟
- 科氏力和地表粗糙度效应
- Actuator Line/Disk模型的转子表示
- 上游风机尾流对下游风机影响的定量评估
Fluent 下的后流解析设置
用Fluent做物体后流解析的推荐设置是什么？
以汽车Ahmed车体（Re $\sim 10^6$）为例，
- 网格：Poly-hexcore。后流区用Body of Influence追加细化
- 湍流模型：SST $k$-$\omega$ DDES（解析后流非定常结构）
- 空间离散：Bounded Central Differencing（LES区）、Second Order Upwind（RANS区）
- 时间步长：$\Delta t \approx 5 \times 10^{-4} L/U_\infty$
- 统计收集：$100 L/U_\infty$ 以上的统计平均时间
Ahmed车体的后斜角度会大幅改变后流结构对吧。
正确。后斜角 $\phi = 25°$ 时后流结构从C柱涡主导急速转变为全压缩型，$C_D$ 不连续变化。这是非常著名的CFD基准。
后流CFD工具比较——SimScale的云CFD对精密后流解析够用吗
后流（Wake）解析工具的选择近年大幅增加。SimScale是浏览器基础的云CFD，利用OpenFOAM后端，无需专业HPC环保就能执行LES解析。对汽车后流简易解析或草稿走行概念检讨级别实用，但高精度后流回复距离预测可能精度不足。反观NUMECA FINE/Open与开源联动容易，内置半解析LES (SALES) 模型用于风力后流。作为工业标准参考，OpenFOAM + windFarmer的组合是研究机关的标准选择。工具选定通常由计算目的精度要求与许可证成本平衡决定。
后流（尾流）的先端研究
后流的绝对/对流不稳定性
请详细讲讲后流的稳定性理论。
Huerre & Monkewitz (1985, 1990) 的框架是标准。
- 对流不稳定性：扰动随下游流动增长。上游无信息传播。外部扰动的放大器
- 绝对不稳定性：扰动在发生地点留存并增长。引发自励振荡。卡门涡列的发生源
后流速度分布的参数 $\Lambda = (U_{min} - U_\infty)/(U_{min} + U_\infty)$ 小于临界值 $\Lambda_{cr} \approx -0.9$ 时（速度缺失足够大），后流转为绝对不稳定。
圆柱直后绝对不稳定，充分下游变成对流不稳定对吧。
完全正确。绝对不稳定区充当"全局振荡子"，决定卡门涡列的频率。这就是Strouhal数的物理来源。
后流的涡追踪手法
飞机翼尖涡这样的长距离后流怎么计算？
全Navier-Stokes计算会导致计算域庞大，因此采用混合手法。
- 涡粒子法：用离散涡粒子表示涡度场。无网格，追踪涡的运输与扩散
- 时间方法：飞机直后的截面作为初始条件，时间演化追踪涡衰减（时空坐标变换）
- CFD → 涡粒子法承接：飞机近旁用CFD，远场用涡粒子法高效计算
Proper Orthogonal Decomposition (POD)
听说后流解析中用POD。
POD从快照（瞬时流场集合）进行能量最优的模式分解。在后流解析中，
第1模式对应涡脱的对称/反对称模式，能量占比通常超过50%。被用于低维模型（ROM）构建和流动控制设计。
风力发电的尾流最新研究
风力尾流研究的最前线是什么？
以下话题很活跃。
- 尾流转向：故意偏转上游风机的偏角，将尾流横向偏转以减轻下游风机影响。Howland et al. (2019, Nature Energy)
- 动态诱导控制：时变风机推力促进后流混合
- LiDAR计测集成：前方LiDAR实时检测尾流，反馈到控制系统
- 机器学习尾流模型：GAN和自编码器实现高速尾流预测
尾流转向在实际风电场应用了吗？
西门子歌美飒和维斯塔斯进行了实机测试，报告称风电场总发电量提高 $1\text{--}3\%$。虽然看起来小，但对大规模海上风力能产生数十亿元的年度差异。
风力发电机后流的叠加效应——海上风电场CFD优化
海上风力发电所（Offshore Wind Farm）中，上游涡轮机的后流（Wake）冲向下游涡轮机，导致发电量下降15～40%的"尾流效应"是大问题。后流回复距离取决于大气湍流强度和风机间距（D：转子直径），通常需要8D～10D。RANS-CFD风电场级解析能同时处理数十台涡轮机，但各后流结构只能近似捕捉。最前线采用LES解析涡轮机间的复杂涡干扰，与机器学习（GP回帰）结合进行风电场优化（偏角控制、倾斜），已展现出仅布置优化就提高年发电量5～10%的实绩。
后流（尾流）的故障应对
常见故障
后流CFD常见问题有哪些？
1. 后流过早消散（速度缺失过小估计）
原因与对策：
| 原因 | 对策 |
|---|---|
| 数值扩散（一阶风上格式） | 二阶及以上精度。LES部分用中心差分 |
| 后流区网格粗糙 | 后流宽度至少10个单元。使用Body of Influence |
| 网格长宽比大 | 后流区长宽比控制在5以下。追求等向网格 |
| RANS涡粘性过大 | 切换到DES/DDES。后流LES区降低涡粘性 |
2. 后流非对称（物理上应该对称）
对称物体的后流却偏斜了。
检查点：
- 网格对称性：非结构网格的对称面网格可能不对称。应沿对称面对齐网格生成
- 非定常性：也许实际上是物理非对称。卡门涡列瞬间非对称，时间平均后应对称
- 迟滞：高Re数下可能存在多个稳定状态。Ahmed车体 $\phi = 25°$ 附近等
3. 壁面分离位置偏差
后流结构由分离位置主导，准确预测壁面分离位置最为关键。
对策：
- $y^+$ 确认：壁函数使用时 $30 < y^+ < 300$，壁面解析时 $y^+ < 1$
- 过渡模型：层流边界层过渡位置重要时（球的阻力危机等），使用γ-Re_θ过渡模型
- 网格诱导分离 (GIS)：DDES中将RANS边界层切换为LES时的非物理分离防止
4. DES/DDES的灰度区问题
DES 不是万能的？
DES的最大问题是"灰度区"。RANS和LES过渡区中两种模型都发挥不了作用。
症状：分离直后的涡无法解析，应有的湍流结构发育不全。也叫Modelled Stress Depletion (MSD)。
对策：
5. 动量积分与壁面力积分不一致
后流动量积分得出的阻力与壁面压力摩擦力积分结果不一致时，
- 动量积分截面位置不当：过远会超出计算域，过近则压力项大而不确定
- 出口边界条件影响：出口太近导致动量不守恒
- 网格未收敛：两个值随网格加细都无法一致，说明有更根本问题
两者一致是解的信任度好的检验。
完全正确。壁面积分与后流动量积分一致是"守恒检验"，非常有用。差异在 $5\%$ 以内为可接受，$1\%$ 以内为极好结果。
后流CFD收敛不良——RANS定常解析非定常卡门涡的矛盾
圆柱后流CFD解析中"残差无法收敛"的问题，实际上往往不是CFD的问题，而是物理问题。Re>46的圆柱后流会产生卡门涡街，流场本质上是非定常的。试图用定常RANS求解时，求解器找不到稳定解，残差不断振荡。解决方案有两个：①切换为非定常RANS（URANS），将结果作为时间平均量评估；②采用SAS（尺度自适应模拟）或DES直接计算大涡。工程实践中"残差振荡 = 物理非定常 = 定常解析假设成立"是信号，需要从根本上重新检视解析手法。
价值
更详细
错误