基于流体底层公理的高超声速边界层扰动饱和约束——可直接嵌入CFD的工程修正模型

作者：张苏杭

             洛阳数学学派

摘要

高超声速边界层转捩与湍流扰动模拟长期存在核心工程难题：传统湍流模型、转捩模型依赖大量经验阻尼、人工限幅、经验间歇因子，通用性差、工况鲁棒性弱、易产生虚假梯度发散与伪高热流。为解决工程CFD普遍存在的扰动过度增长、数值爆破、转捩饱和失真问题，本文依托《切向收敛、法向有限》流体底层公理体系，建立一套物理原生、无经验拟合、可直接嵌入求解器的扰动饱和约束工程模型。

本模型以法向有限公理为核心，构造横向扰动梯度全局上界判据，实现湍流扰动自然饱和、自动抑散、物理真实约束，完全替代传统人工阻尼与数值限幅操作；依托切向收敛公理，保证扰动饱和后平均主流稳态不变，符合真实边界层转捩物理机制。

该方法结构极简、兼容性极强、计算开销极低，可直接嵌入现有RANS/LES求解器，对高超声速边界层转捩位置、湍流饱和强度、壁面热流峰值、激波边界层干扰扰动抑制均有显著改善效果。全部公式可直接代码实现，具备即插即用、工程立刻可用的特征。

关键词：高超声速流动；边界层转捩；湍流饱和；CFD数值约束；扰动上界；流体公理；无经验模型

1 引言

高超声速飞行器气动设计的核心瓶颈集中在边界层转捩预测、湍流热流预测、强梯度扰动数值稳定性三大问题。

当前工程CFD普遍存在以下痛点：

1. 传统转捩模型（γ、Reθ、SST-γ等）依赖大量试验拟合参数，跨马赫数、跨工况泛化能力差；

2. 高剪切、强扩张边界层区域极易出现扰动梯度过度增长，导致数值尖峰、伪高热流、局部发散；

3. 工程上只能依靠人工阻尼、经验限幅、局部滤波、非线性压制补丁维持计算稳定，缺乏底层物理依据；

4. 无法从机理上解释“扰动为何自动饱和、湍流为何不会无限增长”，只能现象级拟合。

本质原因：

现有工程模型缺少真实流体自带的梯度有界约束，数学方程允许扰动无限增长，而真实物理存在天然饱和上限。

依托第一篇基础理论工作建立的切向收敛、法向有限双公理体系，真实流体的所有横向扰动天然存在与时间无关的全局上界。

本文将该底层物理公理工程化、公式化、代码化，构建一套可直接植入CFD求解器的扰动饱和约束模型，实现：

- 无经验参数、无人工阻尼、无拟合常数；

- 物理原生饱和、自动抑制虚假梯度；

- 稳定主流、抑制伪脉动；

- 适配所有边界层流动、涡扰动流动、剪切湍流流动。

2 理论依据（极简工程复用版）

为保证工程落地简洁高效，本章仅保留必要可用结论，完整公理证明与体系论域参见上篇理论论文。

2.1 法向有限公理（工程核心）

全域横向扰动梯度存在全局上界：

\exists C>0,\quad \sup_{t}\|\nabla \boldsymbol{v}_\perp(t)\|_{L^\infty(\Omega)}\le C

横向扰动包含：

1. 边界层壁面法向扰动

2. 涡结构径向扩散扰动

两类扰动工程物理机制统一、约束统一、饱和机制统一。

2.2 切向收敛公理（稳定性保障）

切向平均流长期趋于稳态：

\lim_{t \to \infty} \overline{\boldsymbol{v}_\parallel}(t) \rightarrow \text{稳态主流}

物理效果：

湍流扰动饱和只压制局部无序脉动，不破坏全局主流，完全符合真实转捩物理。

3 工程可执行扰动饱和模型（本文创新、可直接编码）

3.1 统一扰动判定定义

定义工程横向扰动梯度变量：

G_\perp = |\nabla \boldsymbol{v}_\perp|

所有垂直流线方向的剪切、抬升、涡扩散、边界层脉动全部纳入该变量。

3.2 物理上界饱和判据（核心公式）

本文建立工程可用饱和机制：

1. 当 G_\perp < C：扰动自由发展，模型不干预；

2. 当 G_\perp \ge C：触发物理饱和约束，扰动梯度停止增长、自动封顶；

3. 饱和上限 C 由流动拓扑与雷诺数、马赫数决定，无人工调参。

工程简化可直接迭代形式：

G_\perp^{\text{new}} = \min\big(G_\perp^{\text{calc}},\, C(Re,Ma,\mathcal{K})\big)

该式即为全文最核心工程公式，可直接写入CFD迭代步。

3.3 与传统模型的本质区别

1. 传统CFD

增长无物理上限 → 必须人工阻尼压制 → 参数凭经验调参 → 跨工况失准。

2. 本文公理模型

物理天然存在上限 → 自动饱和、自动封顶 → 零经验参数、零人工干预。

3.4 饱和后流场稳定性保障

饱和约束仅压制横向过度扰动，

切向速度场始终服从稳态收敛规律：

\overline{\boldsymbol{v}_\parallel} \quad \text{保持稳定、不溃散、不衰减}

实现工程关键效果：

抑制虚假热流尖峰，但不破坏真实主流气动特性。

4 CFD植入流程（标准工程步骤、可直接交付开发）

步骤1：每迭代步提取横向扰动梯度

计算节点法向/径向横向速度梯度模 G_\perp。

步骤2：读取当地流动参数

雷诺数、马赫数、壁面曲率 \mathcal{K}。

步骤3：计算物理上界 C

（可先采用标准工程标定区间，后续精修标度律）

步骤4：梯度物理封顶

G_\perp = \min(G_\perp,C)

步骤5：回代粘性项、湍动能源项、压力梯度项

完成本轮物理约束修正。

植入层级：求解器外层物理约束模块，不改动N-S内核、不改动差分格式、不破坏原有湍流模型。

5 工程效果与解决的具体问题

5.1 彻底消除数值爆破与伪奇点

高超声速强压缩、强剪切区域常出现梯度飞涨、局部发散，属于数学解越界，非物理现象。

本模型依据真实流体公理强行锁定物理论域，完全杜绝非物理解。

5.2 转捩饱和更贴合真实物理

传统转捩模型靠经验间歇因子压制扰动；

本文模型：饱和是自然物理结果，不是人工拟合结果。

转捩起始、转捩充分发展位置更自然、鲁棒性更强。

5.3 高热流尖峰抑制（工程重大价值）

飞行器热防护设计最怕数值伪热流过高估计导致结构超重，或低估湍流峰值导致烧蚀风险。

本模型给出物理严格的热流扰动上限，为热防护裕度设计提供理论硬边界。

5.4 统一处理边界层扰动与涡扰动

传统需要边界层模型、涡模型、分离模型多套补丁；

本文一套约束通管全部横向扰动，模型结构大幅简化。

6 模型优势对比（工程对标）

1. 无经验参数：摆脱海量拟合系数、阻尼系数、修正函数；

2. 全域稳定：高Re、高Ma、强剪切、激波干扰区域稳定性显著提升；

3. 物理保真：只剔除虚假数值解，保留真实湍流脉动；

4. 轻量化极低开销：单轮迭代仅一次最值判定，几乎不增加计算成本；

5. 全工况兼容：可嵌入RANS、DES、LES、混合求解器。

7 讨论与后续工程迭代方向

本文模型为公理首次工程落地版本，结构极简、稳定性极强、立即可用。

后续可继续量化升级：

1. 完善 C(Re,Ma,\mathcal{K}) 显式标度律，实现全工况精准预测；

2. 嵌入激波–边界层干扰专用约束；

3. 拓展至多相流、气泡动力学流动，对接相关方程体系；

4. 构建基于公理的全新无参数转捩预测模型。

8 结论

1. 依托真实流体底层公理体系，本文构建了物理原生、无拟合参数、可直接工程落地的湍流扰动饱和约束模型。

2. 通过横向扰动梯度物理上界约束，替代传统人工阻尼与数值限幅，从机理上解决高超声速边界层扰动过度增长、虚假梯度发散、伪热流尖峰问题。

3. 模型保证局部扰动有限饱和的同时维持全局主流稳态，贴合真实转捩与湍流物理机制。

4. 方法极简、开销极低、兼容性强，可直接嵌入现有CFD求解器，具备极强的工程实用价值与推广价值。

参考文献

[1] 张苏杭. 流体力学底层公理：切向收敛、法向有限——湍流有限性、全局流向稳态与N-S爆破解论域辨析, 2026.

[2] 张苏杭. 流体运动形态的几何统一, 2026.

[3] 张苏杭. 纳维-斯托克斯方程全局光滑解在DOG/MOC/ECS/MIE框架下的证明, 2026.

[4] 傅德薰, 马延文. 计算流体力学原理与应用. 科学出版社.

[5] Spalart P R. Modern progress in CFD, 2020.