审核日期:2026-10-23。审核对象为 live_core.py 的增量求解、live_server.py 的状态发布和 live-view.js 的实际节点显示。此目录保存审核证据,网页不会读取它们作为回放数据。
独立启动快速网格(459 个顶点、844 个三角形)、Δt=.05、f_n=.16、m*=10、ζ=.01。从 U=.6 持速 50 个参考时间单位,接着用 160 个时间单位连续升至 U=1.5,再持速 40;完成 5000 个数值时间步,共计算时间 250,耗墙钟约 351 秒。时间连续递增、风速连续非减,每步最大入口变化 .00028125。运行时设置了禁止读取 ale/runs 与 fsi-runs 的审计钩子。
流场直接来自每次新装配的 ALE 速度—压力系统。结构由上一时刻已经算出的压力和黏性牵引驱动,运动更新局部网格与无滑移边界。程序没有按风速挑选存储工况,也没有处方正弦力、指定脱涡频率或规定结构振幅。
用不依赖 U、St 或 f_n 的独立 Hann 窗频谱,在 40 时间单位窗口中观察到:t=80 时尾流/升力主频约 .097;t=140 时约 .159,接近设置的 f_n=.16;t=240 时约 .243,已高于结构频率。40 窗半峰峰位移依次由约 .0011D 增至过程中的 .238D,再降为 .096D。这是有限速率扫速过程,不能当成稳定工况的幅值曲线或已经确认的锁频范围。
与网页相同的未零填充有限窗估计在 t=80、140、250 分别给出约 .100、.150、.250,真实分辨率约 .025。独立审核中零填充仅便于定位谱峰,未提升真实分辨率。初始小扰动和流场发展阶段不应按成熟周期流解释。
全部 5000 步的最小有向面积比为 .818917,最大绝对位移为 .238785D。50 组直接保存的全 P2 速度、P1 压力与实际节点坐标接受了不复用求解器形式的独立梯度和边积分核对:
- 外边界位移严格为 0,圆柱几何跟随误差约 1.7×10⁻¹⁶D,所有单元保持原朝向。
- 壁面无滑移速度误差为 0,网格速度与实际位移更新一致。
- 独立压力加黏性牵引积分与记录力的最大差约 6.7×10⁻¹⁶。
- 弱散度残差约 1.4×10⁻¹⁴,边界总流量差约 2.6×10⁻¹⁴;强散度空间 RMS 最大约 .137,不能将弱残差等同于连续解误差。
- 结构的上一时刻流体力平衡相对 RMS 约 1.9×10⁻¹⁴;同时刻力并非本显式分区算法所用力,其相对残差约 7.0%。
逐步离散能量包括物理耗散和隐式 Euler 数值耗散。此次累计数值耗散约为物理耗散的 2.68 倍,纳入两者后的离散能量残差约 3.9×10⁻¹⁴。因而输入的 ζ=.01 只表示物理阻尼,不能说整个离散模拟只有 1% 阻尼,也不能据代数残差宣称时间或网格收敛。
参数关系与文献定义见 参数核对。直接调 Sc 保持 m* 并回算 ζ=Sc/(π²m*)。改变 f_n 不改变此定义下的 Sc;k、c 必须根据新参数同步更新。支承切换锁住当前位移并将速度设零,再释放保持当前网格位置。约束切换的瞬时快照先表达参数/支承变化,完整流场在下一数值步重新求解。
后台把命令排队,在两次数值求解之间统一应用,普通调参和扫速保持同一会话与历史;只有显式重新开始才创建新初始流场。前端使用求解器快照中的真实节点,不在浏览器中另行制造网格变形。
独立合成信号测试确认:已知 .125 频率不随 f_n 的更改而被强制改变;零信号、纯趋势、趋势叠加低于阈值的振动、不足 2.5 周期、超出支持频带、参数修改后过短的新窗口均不报告有效频率。该检查验证了“不以结构频率代替测量”的行为,不代表扫速频谱具有稳态精度。
机器运行速度约为每墙钟秒推进 .71 个参考时间单位;这说明“在线按当前输入求解”,并不保证计算时间和墙钟时间 1∶1。
审核专门从上述在线计算的 t=140 状态复制出两个完全相同的起点,保持 U=1.10625、m*=10、f_n=.16,分别继续求解 40 个时间单位。一个保持 ζ=.01(Sc=.986960),另一个只改为 ζ=.05(Sc=4.934802)。压力、速度、网格坐标、结构位移和速度的起点完全相同;此复现实验只在审核脚本中恢复快照,网页在线求解没有读取该快照的功能。
在继续计算的后 20 时间窗,低 Sc 分支位移半峰峰幅值为 .263055D,高 Sc 分支为 .152127D,比值约 .5783;此有限窗观察中提高 Sc 后振幅降低约 42.2%。这个数值只描述此次相同起点、相同风速和当前离散精度的参数变化实验,不作为稳态减振率或普遍阈值。高 Sc 分支的升力 RMS 反而由 .23562 增为 .26673,说明响应减小不能被解释成流体力必然减小。
两分支各完成 800 步,最小有向面积比分别为 .80008 和 .83499,所有数值有限,外边界仍固定。检查并未预设高 Sc 分支必须更小,只记录实际求解结果。配置和有限窗结果见 Scruton 对照原始报告。