本系统是一个基于 OpenGL 的 3D 粒子可视化教学平台，专门用于解释和展示三维不可压纳维-斯托克斯方程光滑解存在性问题（克雷数学研究所千禧年大奖难题之一）。

纳维-斯托克斯千年难题教学系统（3D 可视化）介绍

一、本系统主要做什么

本系统是一个基于 OpenGL 的 3D 粒子可视化教学平台，专门用于解释和展示 三维不可压纳维-斯托克斯方程光滑解存在性问题（克雷数学研究所千禧年大奖难题之一）。它将抽象的数学概念（平滑初值、有限时间爆破、能量级联、壳模型、数值证据链）通过动态粒子云、柱状图、曲线和标注箭头等直观方式呈现出来。系统包含自动讲解流程（从“什么是 N-S 难题”到“当前研究进展”），并提供多种交互模式（总览、平滑初值、两条结局、二维 vs 三维、能量级联、Sabra 证据链、统计标度、Clay 标准、陈氏同步证据链、陈氏实时回放），方便教学和研究表达。

二、研究是什么——三维 N-S 光滑解存在性问题

三维不可压纳维-斯托克斯方程是流体力学的基本方程，用于描述粘性流体的运动。数学上，给定光滑、无散度、有限能量的初始速度场，是否存在一个全局光滑（即所有导数都有限）的速度场和压力场满足方程？或者，是否会在有限时间内出现奇点（导数发散）？

核心问题：
设 u₀(x) 为 ℝ³ 上光滑、无散度、无穷远衰减的初始速度场。问：是否存在定义在 (0, ∞) 上的光滑速度场 u(x,t) 和压力 p(x,t) 满足：
    ∂u/∂t + (u·∇)u = -∇p + ν Δu,   ∇·u = 0,
    u(x,0) = u₀(x),
    
且能量有限？
如果存在，则称为“全局光滑解”；如果不存在，则称“有限时间爆破”（即某个量在有限时间内趋于无穷大）。

目前已知：

二维情形：已经严格证明全局光滑解存在（Ladyzhenskaya 等）。
三维情形：至今未解决。人类已知弱解（Leray 弱解）存在，但不知道弱解是否光滑。既没有证明全局光滑，也没有构造出公认的严格爆破解。

三、陈恩华的研究是什么

在本系统中，“陈氏同步证据链”（模式8）和 “陈氏实时回放”（模式9）是开发者（陈恩华）个人对 N-S 奇点逼近过程的数值直觉表达和教学演示。

陈恩华利用 Sabra 壳模型（一种简化的湍流级联模型）进行数值模拟，提取关键指标：

E(t)：总能量，随时间变化。
max(k|u|)：高频壳的梯度量，反映小尺度梯度危险。
maxPi：壳间能量通量，反映非线性级联强度。
dt：时间步长，受 CFL 限制，当危险增大时时间步长被压小。

这些指标在 Sabra 模拟中表现出“同步失稳”现象：当模拟逼近可能的奇点时，四个指标同时快速增长或剧烈变化。陈恩华将其可视化为“同步锁定”结构，认为奇点风险不是单一变量异常，而是多指标协同的“结构事件”。

重要说明： 该演示是 数值直觉实验台，旨在帮助理解奇点逼近的形态，不等于三维 N-S 方程的严格数学证明。它用于教学和研究表达，不能替代 Clay 研究所要求的严格解析证明。

四、相关公式排序

以下按从一般到具体、从基本方程到壳模型的顺序列出主要公式：

三维不可压纳维-斯托克斯方程（粘性系数 ν > 0）：

∂u/∂t + (u·∇)u = -∇p + ν Δu,   ∇·u = 0,   u(x,0) = u₀(x).

能量方程（由上式点乘 u 积分得到）：
```
dE/dt = -ν ∫|∇u|² dx,   E = (1/2)∫|u|² dx.
```
能量耗散率与粘性有关，不增。
涡量方程（对 N-S 取旋度）：
```
∂ω/∂t + (u·∇)ω = (ω·∇)u + ν Δω,   ω = ∇×u.
```
涡量拉伸项 (ω·∇)u 是三维特有的危险来源，在二维中此项为零，故二维光滑。
能量级联与壳模型（Sabra 壳模型）：
将波数空间分成对数均匀的壳层 k_n = k₀ 2ⁿ，定义速度壳变量 u_n（复数），满足：
```
        (d/dt + ν kₙ²) uₙ = i kₙ [ a₁ uₙ₊₁ uₙ₊₂ + a₂ uₙ₋₁ uₙ₊₁ + a₃ uₙ₋₁ uₙ₋₂ ]* + fₙ,
        
```
其中 a₁, a₂, a₃ 是耦合系数（通常取 a₁ = 1, a₂ = -β, a₃ = -(1-β) 保证守恒），* 表示共轭，fₙ 为外力。能量 E(t) = Σ |uₙ|²，梯度量可用 kₙ|uₙ| 表征。
壳模型中的几个关键观测指标：
- 总能量：
```
E(t) = Σ_n |uₙ|²
```
- 最大加权梯度：
```
max(k|u|) = max_n (kₙ |uₙ|)
```
- 能量通量（壳间非线性转移）：
```
Π_n = Im[ uₙ* · (非线性项) ]
```
  ，常用最大值 maxΠ 表示通量强度。
- CFL 条件（数值稳定性要求时间步长 dt 足够小）：
```
dt ≤ C / (max(k|u|))
```
  ，因此当 max(k|u|) 暴涨时，dt 被急剧压小。
陈氏同步证据链的判别准则：
当以下四项同时出现显著变化时，视为“同步失稳”：
```
        E(t) ↑,   max(k|u|) ↑↑,   maxΠ ↑↑,   dt ↓↓.
        
```
该准则用于教学演示，不作为严格数学证明。
Clay 研究所要求的严格数学命题：
要么证明：
```
∀ 光滑无散度有限能量 u₀,  ∃ 全局光滑解 u ∈ Cᵀ(ℝ³×(0,∞))，
```
要么构造一个具体的 u₀ 使得存在 T<∞ 满足：
```
lim_t→T max_x |∇u(x,t)| = ∞
```
（或其他等价奇点判定）。

总结： 本系统通过 3D 粒子可视化，将抽象的 N-S 难题分解为直观的教学模块，帮助用户理解问题的本质、可能的研究路径以及数值模拟提供的直觉线索。陈恩华的贡献在于将 Sabra 模型中的同步失稳现象提炼为可视化的“证据链”，用于教学和启发思考，但始终强调这不等于严格证明。

总结： 本系统通过交互式 3D 粒子可视化，将上述抽象公式和概念转化为直观的图形，帮助学习者和研究者更深入地理解 BSD 猜想的结构与挑战。开发者陈恩华在此框架内加入了自己对证明路径的视觉化表达，旨在促进数学教育与研究交流。

代码下载：main.cpp