数值算法与湍流仿真规范曲率约束的NS数值格式、内蕴时间离散方案作者方见华单位世毫九实验室世毫九·内蕴时空正则化纲领数值实现篇面向三维不可压Navier‑Stokes高雷诺数湍流DNS解决传统格式有限时间数值爆破、BKM准则触发的时间刚性、小尺度间歇性奇异性、人工耗散污染物理级串不修改空间微分算子、保留Galilean协变、不引入全局黏性通过几何约束时间测度重参数化实现正则化。1 物理背景与核心矛盾三维不可压NS\begin{cases}\partial_t \boldsymbol u (\boldsymbol u\cdot\nabla)\boldsymbol u -\nabla p\nu\Delta \boldsymbol u\\\nabla\cdot\boldsymbol u0\end{cases}传统数值困境1. 涡量拉伸 \boldsymbol\omega\cdot\nabla \boldsymbol u 驱动局部涡量爆炸BKM判据 \int_0^T\|\boldsymbol\omega\|_{L^\infty}\mathrm dt\infty 决定正则性临近疑似奇点显式格式时间步长被压到极小计算刚性爆炸。2. 常规正则化LES滤波、超黏性、人工耗散会扭曲能量级串、破坏不变性引入非物理耗散。3. 几何视角湍流小尺度结构对应流场诱导的内蕴联络‑曲率传统直角坐标离散忽略流场自身生成的几何结构把几何奇异性直接转化为数值发散。本方案两条核心支柱1. 规范曲率约束NS弱格式借鉴Yang‑Mills规范丛由Clebsch分解建立流场‑规范联络对偶在变分弱形式引入曲率最小二乘惩罚约束局部自适应压制曲率爆炸维持规范不变性。2. 内蕴分形时间离散流场依赖的Sundman‑型测度重参数化 d\tau\Phi(\boldsymbol\omega)\mathrm dt将物理时间潜在爆破区间映射为内蕴时间的有限光滑演化区间构造适配该测度的时间离散算子满足离散几何守恒律。2 规范曲率场构造流体‑规范对偶2.1 Clebsch表示与SU(2)伴随丛曲率速度Clebsch参数化\boldsymbol u\nabla\phi\alpha\nabla\beta把 (\alpha,\beta) 视作SU(2)丛上标量场构造流诱导联络 \mathcal A_i对应的场强规范曲率二形式\mathcal F_{ij}\partial_i\mathcal A_j-\partial_j\mathcal A_i[\mathcal A_i,\mathcal A_j]\mathcal F_{ij} 直接编码涡量、应变率、涡拉伸拓扑曲率模 |\mathcal F|^2 作为局部奇异强度指示器规范变换下不变不依赖坐标系选择。物理含义流场涡旋拓扑畸变等价于底流形上规范丛的曲率增长奇点对应 |\mathcal F|\to\infty。2.2 带规范曲率约束的NS弱形式求解速度‑压强混合变分增加曲率惩罚项最小二乘弱约束不是方程修改是变分空间约束\begin{aligned}\int_\Omega\big(\partial_t^\tau \boldsymbol u\cdot\boldsymbol v (\boldsymbol u\cdot\nabla)\boldsymbol u\cdot\boldsymbol v\big)\mathrm d\Omega\nu\int_\Omega\nabla\boldsymbol u:\nabla\boldsymbol v\,\mathrm d\Omega-\int_\Omega p\,\nabla\cdot\boldsymbol v\,\mathrm d\Omega\\\boldsymbol\lambda_c \int_\Omega \Big(|\mathcal F[\boldsymbol u]|^2-\mathcal F^2_\mathrm{th}\Big)_\,\mathcal D_{\mathcal F}\mathcal F[\boldsymbol u;\delta \boldsymbol u]\,\mathrm d\Omega0\\\int_\Omega q\,\nabla\cdot\boldsymbol u\,\mathrm d\Omega0\end{aligned}• \partial_t^\tau对内蕴时间 \tau 的时间导数• \big(\cdot\big)_正部算子仅当局部规范曲率超过阈值 \mathcal F_\mathrm{th} 才激活约束平滑区惩罚自动关闭• \mathcal D_{\mathcal F}\mathcal F[\boldsymbol u;\delta\boldsymbol u]曲率场对速度的Fréchet导数• \boldsymbol\lambda_c自适应惩罚系数由RG不动点条件、局部雷诺数动态调节不是固定常数。关键性质1. 约束是局部激活不会全局抹掉湍流小尺度2. 满足规范不变、Galilean不变3. 不改变NS原始微分算子仅在变分层面限制解所在函数空间过滤曲率爆炸的伪解。3 内蕴时间重参数化与离散方案3.1 连续内蕴时间变换物理时间 t内蕴时间 \tau(x,t)测度变换\mathrm d\tau\Phi\left(\|\boldsymbol\omega\|_{L^p(\mathrm{local})}\right)\mathrm dt,\quad\Phi(s)\theta_0\frac{s^2}{1s^2/s_*^2}\theta_00 保证变换因子严格有正下界规避Sundman变换可微性缺口s_* 涡量特征尺度。逆关系\frac{\partial}{\partial t}\Phi(\boldsymbol\omega)\,\frac{\partial}{\partial\tau}物理时间潜在爆破时刻 T_*映射到 \tau 上有限区间在内蕴时间坐标系演化算子不再出现局部时间刚性。3.2 全离散内蕴时间离散格式记第 n 物理时刻流场 \boldsymbol u^n,p^n,\boldsymbol\omega^n计算局部权重 \Phi^n\Phi(\boldsymbol\omega^n)。内蕴时间步\Delta\tau^n\Phi^n\,\Delta t^n采用内蕴时间上的二阶BDF向后差分保证几何守恒\partial_t^\tau \boldsymbol u\approx\frac{3\boldsymbol u^{n1}-4\boldsymbol u^n\boldsymbol u^{n-1}}{2\Delta\tau^n}代回物理时间\partial_t^\tau \boldsymbol u\frac{1}{\Phi^n}\partial_t \boldsymbol u\Rightarrow\partial_t \boldsymbol u \approx \Phi^n\cdot \frac{3\boldsymbol u^{n1}-4\boldsymbol u^n\boldsymbol u^{n-1}}{2\Delta\tau^n}迭代求解流程每一步1. 由 \boldsymbol u^n 计算涡量 \boldsymbol\omega^n、规范曲率场 \mathcal F^n2. 计算权重场 \Phi^n(x)评估局部曲率是否超过阈值确定惩罚强度 \lambda_c^n(x)3. 组装带曲率约束的内蕴‑BDF全离散弱形式4. 求解速度‑压强鞍点系统Stokes‑type隐式迭代5. 后验校验离散几何守恒律若局部曲率残余过高局部回退并自适应缩小 \Delta\tau6. 更新物理时间\Delta t^{n1}\Delta\tau/\Phi^{n1}进入下一轮。几何守恒律要求坐标‑测度变换的雅可比必须在离散层面严格匹配质量、动量积分避免变换引入虚假源项。3.3 与传统自适应时间步长本质区别传统CFL\Delta t\propto 1/|\boldsymbol u|只考虑对流速度本方案\Delta t(x)\propto 1/\Phi(\boldsymbol\omega)时间步长由涡量‑规范曲率拓扑畸变主导涡拉伸剧烈区域自动压缩物理时间步平滑大涡区可以放大步长整体仿真效率提升。4 格式关键理论性质1. 规范不变性规范变换下曲率约束项不变不引入规范伪物理解2. 离散H¹稳定性在内蕴时间框架能量耗散具备一致下界规避BKM准则带来的刚性3. 无全局人工黏性约束只作用在曲率超阈值局部区域湍流级串在光滑区域完整保留4. 几何守恒律兼容测度变换的雅可比离散满足质量、动量守恒5. RG尺度匹配惩罚阈值 \mathcal F_\mathrm{th} 可以耦合重正化群不动点条件匹配湍流惯性区标度实现与物理多尺度自洽。5 数值验证算例体系1. 三维Taylor‑Green涡对比DNS基准检验能量、拟能演化观察格式抑制伪爆破的能力2. 各向同性高雷诺数湍流检验能谱 E(k)\propto k^{-5/3} 是否保留间歇性高阶统计3. 涡管拉伸奇点测试构造BKM临界初值对比传统Runge‑Kutta爆破行为本格式维持稳定演化4. 收敛阶测试空间采用DG/DEC离散检验速度、压强收敛阶不受曲率约束破坏。6 现存开放问题与后续拓展1. 曲率场Fréchet导数 \mathcal D_{\mathcal F}\mathcal F 的高效稀疏组装降低计算开销2. 阈值 \mathcal F_\mathrm{th} 的自洽RG演化实现完全无参数3. 拓展至可压缩NS建立可压缩流的规范‑曲率对偶4. 耦合几何机器学习每步由流场几何特征预判奇点风险预调节惩罚与步长。7 简短对比与主流正则化方法方案 核心机制 是否改变NS算子 不变性保留超黏性/滤波LES 高阶耗散/空间滤波 修改微分算子 部分破坏级串SUPG/流线扩散 稳定化残差 添加稳定项 Galilean近似成立本规范曲率内蕴时间 变分曲率约束时间测度重参数化 不改空间微分算子 Galilean、规范不变严格保留