TGV Dissipation/Enstrophy 归一修正后的数据迁移(旧数据 ×2 / ×4) #5

开放中
TiMidlY2026-06-14 17:32:37 +08:00创建 · 0评论
所有者

背景:代码侧已修

CFlowOutput 计算 TGV Dissipation/Enstrophy 归一所用的涡尺度 tgvLength,原先从 REYNOLDS_LENGTH(=1) 派生:

const su2double tgvLength = config->GetLength_Reynolds() / (2.0 * acos(-1.0));  // = 1/(2π)

而初始条件 CUserDefinedSolution 与所有 Python 后处理(plot.py / extract.py / runner.py)用的是 MESH_BOX_LENGTH/(2π) = 2/(2π) = 1/π。同一个涡尺度取了不同的源,差 2×(涡尺度本就必须 = 盒长/(2π) 才能让 u=sin(x/tgvLength)[-1,1] 周期)。

已改 SU2_CFD/src/output/CFlowOutput.cpp:870 为与 IC 一致的盒长派生:

const su2double tgvLength = config->GetMeshBoxLength(0) / (2.0 * acos(-1.0));  // = 1/π

由此 t_c = tgvLength/U = 1/π,用于 Dissipation = -t_c·dEk/dt:937)与 Enstrophy *= t_c²/(ρ·Vol):949)。(代码改动在工作区,待提交。)

本 issue 留作数据迁移

修复改变了新产出的归一化幅值。旧数据(用 1/(2π) 产出)换算到新约定:

换算
Dissipation new = 2 × old
Enstrophy new = 4 × old
Kinetic_Energy 不变(无 t_c 归一)

待办(按既有 job 的重要性逐批处理):

  • 盘点哪些已发布 / 在用的 Diss/Enstrophy 结果是旧归一产出
  • 每批二选一:后处理重归一(×2 / ×4)或重跑
  • 同步更新对应图 / CSV / 文字里的绝对幅值

影响范围

仅常数因子,不改变标度 / 相对结论(耗散律幂次、KE 相对关系、谱熵标度等);只影响这两条曲线的绝对幅值与 KE–耗散预算闭合。

## 背景:代码侧已修 `CFlowOutput` 计算 TGV `Dissipation`/`Enstrophy` 归一所用的涡尺度 `tgvLength`,原先从 `REYNOLDS_LENGTH`(=1) 派生: ```cpp const su2double tgvLength = config->GetLength_Reynolds() / (2.0 * acos(-1.0)); // = 1/(2π) ``` 而初始条件 `CUserDefinedSolution` 与所有 Python 后处理(`plot.py` / `extract.py` / `runner.py`)用的是 `MESH_BOX_LENGTH/(2π) = 2/(2π) = 1/π`。同一个涡尺度取了不同的源,差 2×(涡尺度本就必须 = 盒长/(2π) 才能让 `u=sin(x/tgvLength)` 在 `[-1,1]` 周期)。 **已改** `SU2_CFD/src/output/CFlowOutput.cpp:870` 为与 IC 一致的盒长派生: ```cpp const su2double tgvLength = config->GetMeshBoxLength(0) / (2.0 * acos(-1.0)); // = 1/π ``` 由此 `t_c = tgvLength/U = 1/π`,用于 `Dissipation = -t_c·dEk/dt`(`:937`)与 `Enstrophy *= t_c²/(ρ·Vol)`(`:949`)。(代码改动在工作区,待提交。) ## 本 issue 留作数据迁移 修复改变了**新产出**的归一化幅值。**旧数据**(用 `1/(2π)` 产出)换算到新约定: | 量 | 换算 | |---|---| | `Dissipation` | `new = 2 × old` | | `Enstrophy` | `new = 4 × old` | | `Kinetic_Energy` | 不变(无 t_c 归一) | 待办(按既有 job 的重要性逐批处理): - [ ] 盘点哪些已发布 / 在用的 Diss/Enstrophy 结果是旧归一产出 - [ ] 每批二选一:后处理重归一(×2 / ×4)或重跑 - [ ] 同步更新对应图 / CSV / 文字里的绝对幅值 ## 影响范围 仅常数因子,**不改变标度 / 相对结论**(耗散律幂次、KE 相对关系、谱熵标度等);只影响这两条曲线的绝对幅值与 KE–耗散预算闭合。
TiMidlY 修改标题CFlowOutput: TGV Dissipation/Enstrophy 归一用 REYNOLDS_LENGTH(t_c=1/(2π)),与 IC/论文/Python 的 1/π 差 2×/4×TGV Dissipation/Enstrophy 归一修正后的数据迁移(旧数据 ×2 / ×4) 2026-06-14 18:13:21 +08:00
登录并参与到对话中。
1位参与者
通知
到期时间
到期日期无效或超出范围。请使用“yyyy-mm-dd”格式。

未设置到期时间。

依赖议题

未设置依赖项。

引用
TiMidlY/SU2-Quantum#5
没有提供说明。