后向台阶流
后向台阶流的理论基础
概述
老师，后向台阶流是什么问题？
流道中间有台阶，流在那里分离并重新附着。这是分离与重附着流的最基本基准问题，长期以来被用于CFD代码验证。分离Armaly et al. (1983)的实验数据很著名。
台阶后面会产生涡吗？
是的。台阶下游会产生再循环区域（分离泡）。其长度 $x_r$ 除以台阶高 $h$ 得到的重附着长度 $x_r/h$ 是Reynolds数的函数。这是最重要的验证指标。
支配方程
控制方程是Navier-Stokes吧？
不可压缩Navier-Stokes方程和连续方程。
Reynolds数用台阶高 $h$ 和入口平均流速 $U$ 定义。
重附着长度与Re的关系怎样？
层流区域（$Re < 400$左右）重附着长度基本与Re成正比。有一个近似公式 $x_r/h \approx 0.06 \times Re$。Armaly等人的实验中展开比 $ER = (H+h)/H = 1.94$，Re=100时 $x_r/h \approx 5$，Re=400时 $x_r/h \approx 14$。
Re越大重附着点越远啊。
但是当 $Re > 400$ 前后三维效应变得明显，仅用2D计算就和实验对不上了。这是著名的"Armaly问题的2D-3D转变问题"。
流动结构
再循环区域外还有其他结构吗？
Re增大时，台阶对面的上壁一侧也会出现二次分离泡。再往后台阶角部还会形成小涡。整个流场结构对Re强烈依赖。
| Re范围 | 流动特性 |
|---|---|
| Re < 200 | 仅有下壁主再循环 |
| 200 < Re < 400 | 上壁出现二次分离泡 |
| Re > 400 | 三维不稳定性、跨向变动 |
| Re > 1000 | 湍流转变、非定常涡脱落 |
2D计算足够到Re=400左右就好了。受教了。
重附着点位置为何成为"湍流模型选择的试金石"
后向台阶流中"分离后流体重新附着在台阶下游壁面的距离（重附着长度）"是对比湍流模型性能的标准指标。实验值约为台阶高的6-8倍。但用标准k-ε模型求解常常得到9-11倍的重附着位置，换成SST模型精度就上升。"同一问题却有不同答案"这个差异是选择湍流模型的实际根据。后向台阶看似简单，却包含了"分离、再循环、重附着"三个关键现象，所以成为新湍流模型开发时的标准测试台。
后向台阶流的数值计算方法
数值方法
后向台阶流的数值求解要注意什么？
这个问题中压力和速度的耦合很重要。要用SIMPLE系或投影法这样的标准不可压缩流算法。
压力-速度耦合
SIMPLE法是怎样的机制？
Semi-Implicit Method for Pressure-Linked Equations的缩写。Patankar & Spalding (1972)提出的。步骤是这样的。
1. 从动量方程求得暂时速度场
2. 求解压力修正方程（Poisson型）
3. 修正速度和压力
4. 重复直至收敛
派生方法有SIMPLEC和PISO（针对非定常）。定常计算用SIMPLE/SIMPLEC，非定常用PISO是标准做法。
空间离散化
对流项用什么格式比较好？
后向台阶有再循环，一阶迎风差分数值扩散太大会把重附着长度算长。至少需要2阶精度。
| 格式 | 精度 | 稳定性 | 对重附着长度的影响 |
|---|---|---|---|
| 1阶迎风 | 1阶 | 高 | 高估（数值扩散大） |
| 2阶迎风 | 2阶 | 中 | 合适 |
| QUICK | 3阶 | 略低 | 合适 |
| 中心差分 | 2阶 | 低 | 有振荡风险 |
网格设计
网格要在哪些地方加密？
壁面第一层的 $y^+$ 在层流可以不管，但湍流计算建议 $y^+ < 1$（不用壁面函数时）。结构网格从壁面离开时膨胀比控制在1.1-1.2为好。
出口要这么远啊，30h，相当长了。
否则出口边界条件的影响会改变重附着长度。验证时一定要确认出口位置的敏感性。
台阶后网格的"密度梯度"问题
后向台阶数值计算中易被忽视的是"台阶角部后网格转变"问题。分离点（台阶边缘）是速度梯度最大的地方，需要细网格，但如果从那里到重附着点急速粗化，数值扩散增加会导致再循环区高估或低估。经验法则是"在台阶端周围用不超过台阶高1/10的单元尺寸，向下流方向用膨胀率≤1.1逐步放大"比较安全。另外，即使2D模型中重附着长度与实验吻合，由于忽视了3D效应（角流），这个设置用到3D模型时常常不适用。"2D验证→3D实际"时设置不要直接沿用。
后向台阶流的工程应用
实战指南
实际开展后向台阶分析时，怎样推进比较好？
按以下步骤进行。
分析步骤
1. 几何定义
- 确定展开比 $ER = (H+h)/H$。Armaly问题取 $ER = 1.94$（$h = 0.0049$ m，$H = 0.0052$ m）
- 入口助走段：$10h$ 以上（若用完全发达流入口条件则不需要）
- 出口段：$30h$ 以上
2. 边界条件
- 入口：完全发达抛物线速度分布 $u(y) = U_{max}[1-(2y/H)^2]$ 或均匀流速
- 出口：自然流出条件（$\partial u/\partial x = 0$, $p = 0$）
- 壁面：无滑条件
3. 网格收敛性确认
| 网格水平 | 单元数（2D） | $x_r/h$（Re=200） |
|---|---|---|
| 粗 | 5,000 | 10.8 |
| 中 | 20,000 | 10.2 |
| 细 | 80,000 | 10.05 |
| 参考值 | - | 10.0 |
粗网格就偏差8%啊。
重附着长度取决于壁面剪切应力为零的位置，对网格敏感度高。至少用3种网格规模做Richardson外推比较好。
验证要点
计算结果的合理性怎样确认？
看壁面摩擦系数的符号变化位置就对了。
没错。$C_f > 0$ 是顺流，$C_f < 0$ 是逆流。符号转变点就是重附着点。
燃烧器中的"台阶"——有意制造再循环来稳定火焰
后向台阶流经常被说成"应该避免的流动现象"，但在燃气轮机和工业炉的燃烧器设计中反而"刻意用台阶产生再循环区"。再循环域内滞留高温气体，成为新进入的燃料、空气的点火源——即"火焰的种子"。没有这个再循环区，冷的主流会吹灭火焰（熄火）。现代航空发动机燃烧器设计中CFD之所以不可或缺，部分原因就是要优化台阶后的再循环形状。这是"问题现象"变成"有用现象"的拓扑学好例子。
后向台阶流的软件对比
商用工具实现
主流CFD软件中的设置方法能教我吗？
后向台阶几乎所有CFD软件都有教程。按软件总结要点。
Ansys Fluent
OpenFOAM
simpleFoam（定常）、icoFoam/pimpleFoam（非定常）blockMesh生成结构网格。grading参数细分壁面附近divSchemes中设Gauss linearUpwind grad(U)fixedValue给抛物线分布。codedFixedValue很方便OpenFOAM用blockMesh造台阶形状啊。
两个块组合就成。上游的高度是$H$，下游是$H+h$，按这个定义块的高度就对了。
STAR-CCM+
工具对比
| 项目 | Fluent | OpenFOAM | STAR-CCM+ |
|---|---|---|---|
| 设置便利性 | GUI操作 | 文本字典 | GUI+Java |
| 网格灵活性 | 非结构/结构 | blockMesh/snappyHex | 多面体/Directed |
| 重附着长度精度 | 相同 | 相同 | 相同 |
| 并行性能 | 良好 | 良好 | 良好 |
哪个软件精度都一样啊。网格和设置才是关键吧。
对。求解器本身的差别不如网格品质和对流格式选择对结果的影响大。
厂商对比中"后向台阶验证"成为信任证明的原因
CFD求解器厂商经常在宣传资料和网站上发布后向台阶的验证结果。原因是"这个题目在业界被广泛认可为湍流计算能力的综合评价指标"。因为分离、再循环、重附着都包含，实验数据也充分（Driver & Seegmiller 1985、Armaly et al. 1983等），容易对标。选择求解器时不妨问"能展示一下后向台阶的基准结果吗"，用哪个模型、误差多少、公开讨论的厂商往往技术实力和诚信度都靠得住。
后向台阶流的前沿研究
前沿课题
后向台阶流的研究还在进展吗？
基础问题所以一直被当作新方法的验证工具在用。
湍流模型验证
高Re（Re > 5000）时RANS模型的性能差异明显。Driver & Seegmiller (1985)的实验（$Re_h = 37,400$）是湍流验证的标准。
| 湍流模型 | 重附着长度预测 | 与实验值偏差 |
|---|---|---|
| $k$-$\varepsilon$ 标准 | 偏小（$x_r/h \approx 5.5$） | -15% |
| $k$-$\omega$ SST | 良好（$x_r/h \approx 6.2$） | -4% |
| RSM（Reynolds应力） | 良好（$x_r/h \approx 6.4$） | -1% |
| LES（Smagorinsky） | 良好（$x_r/h \approx 6.5$） | 0% |
| 实验值 | $x_r/h \approx 6.5$ | - |
标准$k$-$\varepsilon$误差这么大啊。
逆压力梯度强的分离流中，$k$-$\varepsilon$标准模型往往高估了涡粘性。该用SST及以上级别的模型。
LES/DNS研究
Le, Moin & Kim (1997)用DNS算了$Re_h = 5,100$的情况。再循环区内湍流统计量（Reynolds应力张量各分量）都求出来了，成了RANS模型改进的基础数据。
近来Wall-Modeled LES（WMLES）能算工程Re数（$Re_h > 10^5$）的问题了。
三维效应与跨向不稳定性
刚才说Re=400左右出现3D效应，能详细说下吗？
Barkley等人(2002)的线性稳定性分析表明，$Re \approx 748$时一个3D模式（跨向波长$\lambda \approx 6.9h$）首次失稳。这是定常分岔，不是Hopf分岔。也就是说出现了3D但仍是定常的图案。
不是非定常，而是定常的3D图案啊。有意思。
这个分岔结构的理解，对判断2D RANS计算的适用范围很关键。
LES中看到的"再循环点非定常摇摆"改变设计
用RANS解后向台阶时再附着长度收敛到一个定值，但LES和实验中重附着点在上下游方向随机摇晃个不停。这不是噪声，而是产生对物体的动态压力载荷的真实现象。在城市风工学中，建筑壁面上流的重附着产生的这种变动载荷会导致疲劳破坏，所以RANS定常解析会低估设计荷载。"台阶后要安装东西"这种设计题时，强烈推荐用LES取时间序列数据，而不仅仅依赖RANS定常解。
后向台阶流的故障排除
故障排除
后向台阶分析时常见什么问题？
典型问题和对策总结一下。
重附着长度与参考值不符
原因和对策：
- 网格不足：壁面附近分辨率不够。增加壁面法向单元数
- 1阶精度格式：数值扩散大，再循环被压平。改2阶或以上格式
- 出口太近：下游至少30h。和40h对比检查敏感性
- 入口条件不符：是否用了完全发达流，或者均匀流？要和文献条件一致
入口条件对结果影响这么大？
很大。入口完全发达抛物线分布和均匀流速区别很大。如果入口不够长，速度分布还没发展完全就到了台阶，会导致重附着点比参考值远。文献的入口条件一定要复现。
计算不收敛
定常计算发散的情况：
- Re太高，实际是非定常→改非定常计算
- 松弛系数太大→压力降到0.2-0.3、速度到0.5-0.7
- 网格纵横比过极端→调整壁面附近分级
非定常计算的CFL条件：
2D计算与实验不一致
Re=600左右2D和实验完全对不上...
那就是3D效应。Re > 400就得用3D。2D会高估重附着长度的倾向。改3D模型，加周期边界条件。跨向宽度取8h以上较好。
检查清单
| 检查项 | 对策 |
|---|---|
| 重附着长度的网格收敛性 | 3种以上网格规模做Richardson外推 |
| 出口距离是否足够 | 至少30h，感度检查 |
| 入口速度分布的妥当性 | 抛物线或足够长的助走段 |
| 对流格式的精度 | 2阶以上 |
| 适当的Re范围 | 2D: Re < 400、3D: 更大 |
先检查网格收敛和出口距离，大多问题都能解决吧。
没错。基础扎实就能成。后向台阶是非常教育性且可靠的基准问题。
"重附着点比实验远"不一定是错
后向台阶中"重附着长度超过实验值"时易误认为计算有误，但先要确认：实验的Re数、入口边界层状态（层流还是湍流、厚度多少）是否精确再现了。入口不是完全发达湍流，而是均匀流，短助走段内速度分布没发展完全就到台阶，剪切层发育晚了，重附着点会比实验远。这种情况是"计算正确，入口边界条件再现有偏"。入口处设"驾驶段"让湍流充分发展后再进台阶，是后向台阶分析的讲究地方。
相关课题
更多
错误