瞬子是什么
三维不可压缩湍流,由大尺度随机力驱动。力在空间光滑,时间上是白噪声。
问:原点、最终时刻轴向涡量ω_z等于某个大值a,概率多大?概率随a指数衰减:P ~ exp-S(a)。S(a)是作用量。
最可能产生涡量的方式,对应最小作用量。对应的最优流场,就叫“瞬子”。它不是人造的,是实际湍流中主导涡结构自然出现的物理配置。
老办法:旋转加挤压
先加方位角力F_θ让流体绕轴转;再加径向和轴向力,让中心平面径向内流、轴向外流。
涡管被径向压缩、轴向拉伸。中心径向应变为负,轴向应变为正。涡量单调增长。
单调增长有阈值
作用量随a二次增长:log P = -S ∝ -a²。概率是高斯分布。
中心涡量满足 ∂tω_z = σω_z + νΔω_z + ∂F_θ/∂r。σ = ∂z u_z,是轴向应变。
低a时,σ在z=0整段时间为正。涡量和应变还没太大,科里奥利力来不及阻止径向挤压。
径向流 u_r = -σr,会在背景涡量ω̄区域产生方位科里奥利力 ω̄σr。外力旋度 ∂F_θ/∂r 若压不过ω̄σ,作用量就随a二次增长。
科里奥利力变盟友
a变大后,系统找到更优策略,F随a增长更慢。
数据:a=34时,∂Fθ/∂r≈260,ω̄σ≈105。a=73时,∂Fθ/∂r≈420,ω̄σ≈525。两者已接近。
a>80后,S(a)增长慢于二次。
科里奥利力是保守恢复力。它让涡量最大值分裂,并让涡量波沿旋转轴传播。
最优策略是在这些传播脉冲内部增长涡量。大尺度F_θ模式跟着脉冲跑,F_z产生相反的轴向流,减慢脉冲。
每次分裂前,σ在最大值处为负;分裂后转正。最终阶段,z=0处σ为正。中心最终的大涡量,靠两个大脉冲碰撞,撞进轴向拉伸、径向挤压的区域。
和奇点构造像在哪
OpenAI最近提出,用“空间振荡脉冲”在有限时间Navier-Stokes爆破中制造无界涡量增长。
关键共同点:涡丝必须通过非线性运动方程自己长出来。强迫只能间接作用。
中国为什么要在意
航空发动机、燃气轮机、火箭发动机,里面全是湍流。涡量怎么生成、集中、耗散,直接关系燃烧稳定、气动加热、噪声。
高超声速边界层转捩、核聚变托卡马克边缘湍流、气象海洋预报,本质都是旋转流体里的涡量演化。
中国要走向前沿,不能只会买软件、跑算例。天河、神威、曙光算得动,绵阳风洞群、JF-12测得到,北大、中科院力学所、天津大学有人做理论。Navier-Stokes是千禧年难题,中国数学界也在啃。
这条路值得算下去
瞬子方法把极端事件变成可计算的优化问题。这是基础研究,也是工程湍流的根。不喊口号,算下去。
— 感谢阅读 · 欢迎点赞关注 · 小编 —
热门跟贴