衝撃波-境界層干渉
理论与物理
冲击波-边界层干扰是什么
老师，冲击波-边界层干扰（SBLI）在超音速话题中经常出现，具体是什么样的现象呢？
是可压缩流动中冲击波与壁面边界层相交时发生的强烈相互作用。冲击波引起的急剧压力上升传递到边界层时，边界层内的低速流体无法承受逆压梯度而发生分离。这个分离泡与冲击波相互影响，形成非常复杂的流场。
分离和冲击波会同时发生啊。在什么情况下会成为问题呢？
典型的有进气道内部的冲击波/边界层干扰、翼面上的跨音速冲击波引起的分离，以及超燃冲压发动机内的冲击波反射。这些都直接关系到气动性能的急剧恶化、非定常载荷或机体振动，是航空航天领域最重要的设计课题之一。
控制方程与无量纲参数
描述SBLI的方程和指标有哪些呢？
基础是可压缩Navier-Stokes方程，但SBLI中特别重要的是入射冲击波前后的关系式以及边界层的分离判据。首先，斜激波的压力比可以这样表示。
$$ \frac{p_2}{p_1} = 1 + \frac{2\gamma}{\gamma+1}(M_1^2 \sin^2\beta - 1) $$
这里 $\beta$ 是激波角，$M_1$ 是上游马赫数。这个压力跳跃会作用于边界层。
那个压力跳跃越大就越容易分离吗？
没错。根据自由干扰理论（Free Interaction Theory），开始分离所需的压力上升由上游边界层状态决定。层流边界层情况下，分离所需的压力上升非常小，但湍流边界层由于近壁面动量较大，不易分离。干扰强度常用干扰压力参数来表示。
$$ \chi = \frac{p_{plateau}/p_1 - 1}{(C_f/2)^{1/2}} $$
$C_f$ 是上游壁面摩擦系数。$\chi$ 越大，干扰越强，会形成更大的分离泡。
层流和湍流差别好大啊。现实中的飞机还会涉及转捩位置的问题吧...
这正是SBLI的难点所在。有时冲击波会诱发边界层转捩，有时在冲击波上游已经转捩为湍流。冲击波诱发转捩（Shock-Induced Transition）在超燃冲压发动机进气道设计中尤为重要，转捩位置的预测精度直接关系到发动机性能的推算。
干扰模式的分类
SBLI也有好几种类型吗？
Edney的分类很有名，根据冲击波与边界层的相交形态分为Type I到Type VI共6种。我来介绍几个典型的。
分类 干扰模式 特征 实例 Type I 斜激波/边界层 最常见。形成分离泡 进气道壁面 Type II 冲击波/冲击波/边界层 超音速射流撞击壁面 射流偏转 Type III 正激波/边界层 lambda型激波结构 跨音速翼面 Type IV 冲击波/冲击波相交 形成超音速射流，极大加热 前缘形状干扰
Type IV写着极大加热，有多危险呢？
可能会产生非干扰时数倍到十数倍的局部热流密度。航天飞机前缘部位的问题就是个著名案例。如果在设计阶段遗漏这一点，热防护系统可能无法承受，因此必须用CFD进行事先预测。
Coffee Break 杂谈
SWBLI研究的原点——1940年代风洞照片中神秘的分离
冲击波-边界层干扰（SWBLI）的系统性研究始于二战后的1940年代末。当时的研究者们在风洞实验中对超音速翼型进行油膜可视化，发现激波根部出现了意想不到的流动分离图案。“冲击波撞击导致边界层剥离”这一现象的阐明，是超音速飞机设计的当务之急。Lees & Reeves（1956年）等人的“自由相互干扰理论”和Chapman（1958年）的分离长度关联式等，当时的理论至今仍用于入门级的理解。这是在数值计算（CFD）被使用之前，用铅笔和直尺开拓的领域。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水流会不稳定地喷溅，过一会儿才变成稳定的水流对吧？描述这个“变化过程中”的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭导致流动变化，这些都是非定常现象。那么定常分析是什么？就是只看“经过足够时间流动稳定之后”——也就是令此项为零。计算成本会大幅下降，因此先用定常求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶扔进河里会怎样？会被水流带着往下游漂对吧。这就是“对流”——流体运动搬运物质的效果。暖气的热风能到达房间另一端，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——这项包含“速度×速度”，因此是非线性的。也就是说，流速越快，这项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置的经历吗？即使不搅拌，过一会儿也会自然混合对吧。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水对吧。蜂蜜因为粘度（$\mu$）高所以不易流动。粘度越大，扩散项越强，流体的运动就变得“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，Re数大的流动中对流占压倒性优势，扩散则成为配角。
- 压力项 $-\nabla p$：推注射器的活塞，液体就会从针头有力地射出对吧？为什么呢？因为活塞侧压力高，针头侧压力低——这个压差就是推动流体的力。水坝放水也是同样原理。天气图中等压线密集的地方会怎样？没错，会刮强风。“有压差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里的误解点：CFD的“压力”通常指表压而非绝对压力。切换到可压缩分析时结果突然出错，原因可能就是混淆了绝对压力/表压。
- 源项 $S_\phi$：受热的空气会上升——为什么呢？因为比周围更轻（密度更低），被浮力推上去了。这个浮力作为源项添加到方程中。还有，燃气灶的火焰产生化学反应热，工厂的电磁泵对金属熔液施加洛伦兹力…这些都是“从外部向流体注入能量或力”的作用，用源项表示。忘记源项会怎样？自然对流分析中忘记加入浮力，流体就完全不动——变成冬天开了暖气但热空气不上升这种物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3时）：将密度视为常数处理。马赫数0.3以上需考虑压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用的情形：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要激波捕捉）、自由表面流动（需要VOF/Level Set等）
量纲分析与单位制
变量 SI单位 注意事项·换算备忘 速度 $u$ m/s 从入口条件的体积流量换算时，注意截面积单位 压力 $p$ Pa 区分表压和绝对压力。可压缩分析中使用绝对压力 密度 $\rho$ kg/m³ 空气: 约1.225 kg/m³@20°C，水: 约998 kg/m³@20°C 粘性系数 $\mu$ Pa·s 注意与运动粘度系数 $\nu = \mu/\rho$ [m²/s] 混淆 雷诺数 $Re$ 无量纲 $Re = \rho u L / \mu$。层流/湍流转捩的判断指标 CFL数 无量纲 $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性
老师，冲击波-边界层干扰（SBLI）在超音速话题中经常出现，具体是什么样的现象呢？
是可压缩流动中冲击波与壁面边界层相交时发生的强烈相互作用。冲击波引起的急剧压力上升传递到边界层时，边界层内的低速流体无法承受逆压梯度而发生分离。这个分离泡与冲击波相互影响，形成非常复杂的流场。
分离和冲击波会同时发生啊。在什么情况下会成为问题呢？
典型的有进气道内部的冲击波/边界层干扰、翼面上的跨音速冲击波引起的分离，以及超燃冲压发动机内的冲击波反射。这些都直接关系到气动性能的急剧恶化、非定常载荷或机体振动，是航空航天领域最重要的设计课题之一。
描述SBLI的方程和指标有哪些呢？
基础是可压缩Navier-Stokes方程，但SBLI中特别重要的是入射冲击波前后的关系式以及边界层的分离判据。首先，斜激波的压力比可以这样表示。
这里 $\beta$ 是激波角，$M_1$ 是上游马赫数。这个压力跳跃会作用于边界层。
那个压力跳跃越大就越容易分离吗？
没错。根据自由干扰理论（Free Interaction Theory），开始分离所需的压力上升由上游边界层状态决定。层流边界层情况下，分离所需的压力上升非常小，但湍流边界层由于近壁面动量较大，不易分离。干扰强度常用干扰压力参数来表示。
$C_f$ 是上游壁面摩擦系数。$\chi$ 越大，干扰越强，会形成更大的分离泡。
层流和湍流差别好大啊。现实中的飞机还会涉及转捩位置的问题吧...
这正是SBLI的难点所在。有时冲击波会诱发边界层转捩，有时在冲击波上游已经转捩为湍流。冲击波诱发转捩（Shock-Induced Transition）在超燃冲压发动机进气道设计中尤为重要，转捩位置的预测精度直接关系到发动机性能的推算。
SBLI也有好几种类型吗？
Edney的分类很有名，根据冲击波与边界层的相交形态分为Type I到Type VI共6种。我来介绍几个典型的。
| 分类 | 干扰模式 | 特征 | 实例 |
|---|---|---|---|
| Type I | 斜激波/边界层 | 最常见。形成分离泡 | 进气道壁面 |
| Type II | 冲击波/冲击波/边界层 | 超音速射流撞击壁面 | 射流偏转 |
| Type III | 正激波/边界层 | lambda型激波结构 | 跨音速翼面 |
| Type IV | 冲击波/冲击波相交 | 形成超音速射流，极大加热 | 前缘形状干扰 |
Type IV写着极大加热，有多危险呢？
可能会产生非干扰时数倍到十数倍的局部热流密度。航天飞机前缘部位的问题就是个著名案例。如果在设计阶段遗漏这一点，热防护系统可能无法承受，因此必须用CFD进行事先预测。
SWBLI研究的原点——1940年代风洞照片中神秘的分离
冲击波-边界层干扰（SWBLI）的系统性研究始于二战后的1940年代末。当时的研究者们在风洞实验中对超音速翼型进行油膜可视化，发现激波根部出现了意想不到的流动分离图案。“冲击波撞击导致边界层剥离”这一现象的阐明，是超音速飞机设计的当务之急。Lees & Reeves（1956年）等人的“自由相互干扰理论”和Chapman（1958年）的分离长度关联式等，当时的理论至今仍用于入门级的理解。这是在数值计算（CFD）被使用之前，用铅笔和直尺开拓的领域。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水流会不稳定地喷溅，过一会儿才变成稳定的水流对吧？描述这个“变化过程中”的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭导致流动变化，这些都是非定常现象。那么定常分析是什么？就是只看“经过足够时间流动稳定之后”——也就是令此项为零。计算成本会大幅下降，因此先用定常求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶扔进河里会怎样？会被水流带着往下游漂对吧。这就是“对流”——流体运动搬运物质的效果。暖气的热风能到达房间另一端，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——这项包含“速度×速度”，因此是非线性的。也就是说，流速越快，这项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置的经历吗？即使不搅拌，过一会儿也会自然混合对吧。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水对吧。蜂蜜因为粘度（$\mu$）高所以不易流动。粘度越大，扩散项越强，流体的运动就变得“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，Re数大的流动中对流占压倒性优势，扩散则成为配角。
- 压力项 $-\nabla p$：推注射器的活塞，液体就会从针头有力地射出对吧？为什么呢？因为活塞侧压力高，针头侧压力低——这个压差就是推动流体的力。水坝放水也是同样原理。天气图中等压线密集的地方会怎样？没错，会刮强风。“有压差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里的误解点：CFD的“压力”通常指表压而非绝对压力。切换到可压缩分析时结果突然出错，原因可能就是混淆了绝对压力/表压。
- 源项 $S_\phi$：受热的空气会上升——为什么呢？因为比周围更轻（密度更低），被浮力推上去了。这个浮力作为源项添加到方程中。还有，燃气灶的火焰产生化学反应热，工厂的电磁泵对金属熔液施加洛伦兹力…这些都是“从外部向流体注入能量或力”的作用，用源项表示。忘记源项会怎样？自然对流分析中忘记加入浮力，流体就完全不动——变成冬天开了暖气但热空气不上升这种物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3时）：将密度视为常数处理。马赫数0.3以上需考虑压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用的情形：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要激波捕捉）、自由表面流动（需要VOF/Level Set等）
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 从入口条件的体积流量换算时，注意截面积单位 |
| 压力 $p$ | Pa | 区分表压和绝对压力。可压缩分析中使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气: 约1.225 kg/m³@20°C，水: 约998 kg/m³@20°C |
| 粘性系数 $\mu$ | Pa·s | 注意与运动粘度系数 $\nu = \mu/\rho$ [m²/s] 混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转捩的判断指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
RANS分析的局限与湍流模型选择
用CFD求解SBLI时，RANS分析能在多大程度上信赖呢？
这可以说是SBLI在CFD中的最大课题。RANS湍流模型难以准确预测在冲击波引起的急剧逆压梯度下的边界层行为。特别是分离泡的尺寸和再附着位置的精度会成为问题。
不同模型的结果会差很多吗？
差别很大。总结一下典型模型的倾向如下。
| 湍流模型 | 分离预测 | 壁面压力 | 热流密度 | 计算成本 |
|---|---|---|---|---|
| SA（Spalart-Allmaras） | 低估分离 | 稍不准确 | 高估倾向 | 低 |
| SST $k$-$\omega$ | 相对良好 | 良好 | 中等精度 | 低 |
| $k$-$\omega$ GEKO | 参数可调 | 取决于调整 | 取决于调整 | 低 |
| RSM（Reynolds Stress） | 有改善但不稳定 | 良好 | 稍有改善 | 中 |
SST $k$-$\omega$ 在实际工作中平衡性最好，但在强干扰下有时会低估分离长度20-30%。
LES/DES高精度分析
要想进一步提高精度，果然还是得用LES吗？
壁面解析LES精度最高，但近壁面网格要求非常严格，需要 $\Delta x^+ \sim 50$, $\Delta y^+ < 1$, $\Delta z^+ \sim 20$。对于像SBLI问题这样的高Re数壁面流动，网格数量会变得极其庞大。因此在实际工作中，DES（分离涡模拟）或WMLES（壁面模型化LES）是有效的。
DDES和SBES有什么区别呢？
DDES用遮蔽函数区分RANS区域和LES区域，但在切换区域容易产生称为“Grey Area”的模糊区域。SBES是改进了这个Grey Area问题的方法，从RANS应力张量到LES应力张量的过渡更加平滑。Ansys Fluent 2020R1之后可用。
空间离散格式的选择
捕捉冲击波的数值格式该怎么选呢？
相关主题
この記事の評価ご回答ありがとうございます！参考に
なったもっと
詳しく誤りを
報告
要想进一步提高精度，果然还是得用LES吗？
壁面解析LES精度最高，但近壁面网格要求非常严格，需要 $\Delta x^+ \sim 50$, $\Delta y^+ < 1$, $\Delta z^+ \sim 20$。对于像SBLI问题这样的高Re数壁面流动，网格数量会变得极其庞大。因此在实际工作中，DES（分离涡模拟）或WMLES（壁面模型化LES）是有效的。
DDES和SBES有什么区别呢？
DDES用遮蔽函数区分RANS区域和LES区域，但在切换区域容易产生称为“Grey Area”的模糊区域。SBES是改进了这个Grey Area问题的方法，从RANS应力张量到LES应力张量的过渡更加平滑。Ansys Fluent 2020R1之后可用。
捕捉冲击波的数值格式该怎么选呢？
相关主题
なった
詳しく
報告