低正则弱解与实验流场的规范判定：无导数、抗噪声的工程奇点检测

低正则弱解与实验流场的规范判定：无导数、抗噪声的工程奇点检测
作者：方见华
单位：世毫九实验室
摘要
流体实验PIV/PTV测得流场天然带有测量噪声、截断误差，直接速度梯度求导（涡量、应变率、不变量、λ₂准则等）属于导数型算子，对噪声高度敏感；在低正则弱解框架下，欧拉/NS方程弱解允许速度场本身连续、导数局部无界，即工程意义流场奇点（剪切层撕裂、涡核破裂、间断雏形、集中耗散区）。传统奇点识别依赖一阶、二阶导数，在弱解场景直接失效。本文提出一套无导数、仅使用速度场原始离散数据的规范判定体系，不需要计算梯度，实现抗噪声的流场奇点定位，适配实验流场、大涡模拟粗分辨率弱解场。
关键词：低正则弱解；Navier‑Stokes；流场奇点；PIV噪声；无导数检测；集中耗散；规范判据
1 背景问题
1. NS方程数学理论：弱解u\in L^2\cap L^4，速度本身可积连续，空间导数可以局部发散，即解正则性局部崩塌，对应数学奇点候选。
2. 工程现实：实验流场只有离散采样速度\boldsymbol u(\boldsymbol x_i)，叠加高斯/椒盐测量噪声。
◦ λ₂、Q判据、Δ判据、Okubo‑Weiss：全部建立于速度梯度\nabla \boldsymbol u，需要做数值微分，噪声会虚假生成大量“伪涡、伪奇点”。
◦ 对低正则弱解：奇点恰恰发生在\nabla u趋向无穷位置，数值微分直接爆炸，计算梯度本身在奇点邻域数学上无良好定义。
矛盾：要检测导数发散的区域，却依赖求导数的工具。
现有降噪方案（高斯光滑、POD滤波）代价：抹平真实弱解奇点信号，把集中耗散做弥散，丢失物理。
因此需要规避微分运算，仅使用速度场点取值，构造尺度依赖的积分型判定量。
2 低正则弱解的局部特征（无梯度语言重述）
设空间邻域B_r(\boldsymbol x_0)，半径r，中心点\boldsymbol x_0为待检测点。
对经典光滑解：邻域内速度差分随尺度衰减
\iint_{B_r}\big|\boldsymbol u(\boldsymbol x)-\boldsymbol u(\boldsymbol x_0)\big|^2\mathrm d\boldsymbol x \sim r^{2\alpha+3}
\alpha为Hölder正则指数。
• 光滑区域：\alpha\approx1；
• 剪切层、弱解奇点邻域：正则性下降 \alpha<1；临界奇点候选：\alpha\to0，差分积分随r增长极慢。
核心观测：不需要求导，速度自身的局部差分统计，编码了导数爆破信息。导数发散等价于速度场局部Hölder指数降低。
3 无导数奇点检测规范判据：离散工程版本
3.1 局部速度振荡积分算子S_r(\boldsymbol x)
离散形式（PIV网格，采样点\{\boldsymbol x_k\}\in B_r(\boldsymbol x_0)）
S_r(\boldsymbol x_0)=\frac{1}{N_r}\sum_{\boldsymbol x_k\in B_r(\boldsymbol x_0)}\big\|\boldsymbol u(\boldsymbol x_k)-\boldsymbol u(\boldsymbol x_0)\big\|^2
N_r为邻域内采样点数量。
S_r仅做速度相减、模平方、平均，全程无任何空间导数。
尺度扫描：取一组嵌套尺度r_1<r_2<r_3，计算S_{r_1},S_{r_2},S_{r_3}，做对数‑对数回归：
\log S_r \sim 2\alpha \log r + C
拟合得到局部Hölder指数\hat \alpha(\boldsymbol x_0)。
物理阈值规范
1. 光滑流场区：\hat\alpha \ge 0.7；
2. 强剪切/涡边界：0.3<\hat\alpha<0.7；
3. 奇点候选区（低正则崩塌）：\boldsymbol{\hat\alpha \le 0.3}。
噪声区分机制：
纯测量噪声，和空间尺度r无关，S_r几乎不随r变化，但噪声贡献一个常数基底S_\mathrm{noise}。预处理：用远离剪切的静止区标定噪声基底S_0，做补偿S_r'=S_r-S_0。
• 伪奇点（纯噪声）：扣除噪声基底后S_r'\approx0；
• 真实低正则奇点：扣除基底后S_r'随r显著上升。
3.2 补充不变量：速度环流无导数替代C_r(\boldsymbol x_0)
传统环流\oint \boldsymbol u\cdot \mathrm d\boldsymbol l，离散可以直接沿邻域多边形采样点求和，依然不需要微分：
C_r=\sum_{k}\boldsymbol u(\boldsymbol x_k)\cdot (\boldsymbol x_{k+1}-\boldsymbol x_k)
|C_r|局部异常抬升，代表邻域有旋转结构。
组合判定规则（工程规范）
\textbf{奇点标记：}\quad
\hat\alpha\le0.3 \quad \textbf{且}\quad S'_r>0 \quad \textbf{且}\quad |C_r|>\epsilon_\mathrm{circ}
三条同时满足判定为流场奇点候选。只满足\hat\alpha\le0.3可能只是强剪切层；C_r用来区分旋转主导奇点（涡核破裂）和纯剪切奇点（间断层）。
4 和传统导数判据数学对应关系
对充分光滑场，泰勒展开
\boldsymbol u(\boldsymbol x)=\boldsymbol u(\boldsymbol x_0)+\nabla\boldsymbol u \cdot(\boldsymbol x-\boldsymbol x_0)+o(|x-x_0|)
代入S_r：
S_r \propto \|\nabla \boldsymbol u\|^2 \, r^2
对数拟合直接给出\alpha\approx1，和光滑理论自洽。
当\|\nabla u\|\to\infty（数学奇点）：有限采样下，邻域内速度剧烈跳变，S_r增长变缓，拟合\alpha压低，对应低正则。
关键：传统判据\nabla u是正向直接观测；本套算子S_r,\,C_r是间接积分观测，当\nabla u本身数学上不存在时依然有定义。
5 抗噪声工作流程（可直接用于PIV后处理）
1. 噪声标定：选取流场静止区域，统计S_0，即纯噪声的S_r基底；
2. 多尺度环形邻域滑动扫描全场每个网格点；
3. 每个点计算多尺度S_r，扣除噪声基底S_r'=S_r-S_0；
4. 对数线性回归得到局部Hölder指数\hat\alpha；
5. 同时计算邻域离散环流C_r；
6. 应用组合阈值输出奇点掩码；
7. 后处理：掩码连通域聚类，得到奇点几何位置、空间尺度。
注意事项：
• r不能小于PIV空间分辨率；最小尺度取2‑3倍PIV网格步长，规避亚网格混叠；
• 不能单靠单一尺度r判定，必须多尺度回归，单点单尺度极易被噪声欺骗。
6 弱解理论层面规范边界
1. 该方法检测的是正则性局部下降区域，不等同严格数学NS爆破奇点。工程上称之为奇点候选。有限分辨率实验永远无法严格证明数学爆破，只能给出证据。
2. 物理区分两类低正则区域：
◦ A类：粘性正则化：极小尺度粘性抹平，大尺度看表现低正则，本质光滑（强涡、薄剪切层）；
◦ B类：NS弱解本身集中耗散，即真正数学意义奇异行为。
本套无导数判定只能定位A+B，进一步区分A/B需要多分辨率收敛检验：不断增大采样分辨率，观察\hat\alpha是否持续压低。
收敛检验：加密采样，如果\hat\alpha持续下降→指向B类弱解奇异行为；分辨率提高后\hat\alpha回升→属于粘性薄层A类。
7 对比总结表
判据 是否需要导数 噪声鲁棒性 适用对象 局限
是（） 差 光滑涡结构 低正则弱解邻域失效，噪声生成伪结构
Okubo‑Weiss 是 差 二维流场 三维不可用，微分噪声敏感
本文组合 无导数，积分差分 高 实验PIV、LES粗场、NS低正则弱解 不能直接给出涡量，需要多尺度分析，计算开销更大
8 开放工程问题与延伸
1. 三维离散流场下球面邻域采样优化、各向异性修正：拉伸涡会使得球形邻域估计\alpha产生偏差，需要椭球加权邻域。
2. 阈值\hat\alpha=0.3是经验工程阈值，可以通过DNS基准数据集标定概率分布，给出奇点置信度而不是硬0‑1掩码。
3. 拓展：把这套无导数正则性检测思路推广到其它场：标量浓度场、压力场，检测标量间断、混合层奇点。
4. 和集中耗散测度联系：\int_{B_r}|\nabla u|^2\mathrm d x耗散积分在奇点发散，在本框架，不需要求梯度，S_r就是耗散积分无导数代理量。