别再让 Agent 用“圆”拟合世界了：塞一条 I²=-N，流体涡量仿真从发散到稳跑 10 万步

• 标签：AI Agent Python 物理仿真 复杂系统 数学建模
• 专栏：《Physical AI 与连续几何》
• 话题：Agent 做流体/涡旋仿真时把旋转当纯圆、伸缩当外挂，用螺旋数把涡度演化绑成连续生成
• 简介：不是模型不会写 Navier-Stokes 离散化，是默认复数只有 i²=-1，描述不了带阻尼/增益的螺旋涡线。改 System Prompt 不微调，附可跑代码。

让 Agent 写个二维涡量场演化，需求就一句：转着转着带点耗散，别炸。

它第一版甩过来：

import numpy as np

def step(vortex, w, dt, nu):
ang = np.angle(vortex)
ang += w * dt
mag = np.abs(vortex) * (1 – nu*dt)
return mag * np.exp(1j*ang)

前 2000 步挺乖。

跑到 8000 步，涡心开始抖；2 万步后模长反升；5 万步直接 NaN。

我让它“加个涡度守恒校验”，它补了 clip、补了归一化、补了人工耗散。

越补越像在给错误几何打石膏。

一、根因：Agent 脑子里只有“圆”

标准复数公理：

i2=−1

= 纯旋转，模长不动。

但流体涡线不是圆，是等角螺线：

• 旋转（涡度）
• 伸缩（耗散/能量级联）
• 历史（绕了几圈）
• 各向异性（径向/切向增长率不同）

Agent 默认把“转”和“缩”拆成两个通道：

z = rotate(z) * scale(gain)

步数一多，相位在 [-π, π] 里折叠，圈数丢了；耗散和旋转各算各的，螺旋度（helicity）直接不守恒。

用圆去拟合螺线，不是误差问题，是坐标系选错了。

二、换螺旋数：I² = -N

螺旋生成论里那条公理：

I2=−N(N>0)

• N=1：退化为标准复数（纯圆）
• N≠1：乘 I = 转 90° + 模长伸缩 √N 倍

旋转和伸缩变成同一个代数对象，Agent 没法再把它俩拆开。

先用状态类跑通（不造新类型，纯 numpy）：

import numpy as np

class SpiralVortex:
def __init__(self, mag=1.0, phase=0.0, turns=0, N=1.0):
self.mag = mag
self.phase = phase
self.turns = turns
self.N = N

def integrate(self, omega, dt, nu):
# omega: 角速度, nu: 耗散系数
dtheta = omega * dt
self.phase += dtheta

# 连续相位，显式记圈数
while self.phase >= 2*np.pi:
self.phase -= 2*np.pi
self.turns += 1
while self.phase < 0:
self.phase += 2*np.pi
self.turns -= 1

# 对数标度：耗散走螺旋语义，不硬乘 (1-nu*dt)
dln = (-nu * dt) / np.sqrt(self.N)
self.mag *= np.exp(dln)
return self

def complex(self):
total = self.turns*2*np.pi + self.phase
# 螺旋欧拉投影：e^{Iθ} = cos(√N θ) + (I/√N) sin(√N θ)
x = self.mag * np.cos(np.sqrt(self.N) * total)
y = self.mag * (1/np.sqrt(self.N)) * np.sin(np.sqrt(self.N) * total)
return x + 1j*y

# 仿真
v = SpiralVortex(N=1.18)
traj = []
for t in range(100000):
v.integrate(omega=0.7, dt=1e-3, nu=0.01)
if t % 1000 == 0:
traj.append(v.complex())

关键三点：

三、System Prompt 硬约束

SYSTEM = """
流体/涡量仿真约束（螺旋生成论视角）：
1. 默认 I^2=-N，禁止裸用 i^2=-1 思维
2. 涡旋状态用 (mag, continuous_phase, turns, N) 表达
3. 旋转与伸缩同源，禁止 z*exp(iθ)*gain 拆通道
4. 每步校验螺旋度：有耗散时 |z| 单调不增
5. 多轮推理先回放 turns，再算相对涡量
"""

同模型（Qwen/DeepSeek/Claude）跑二维涡场：

指标

改前（圆拟合）

改后（I²=-N）

没微调、没 RAG。

四、.cursorrules 直接抄

# 流体/物理仿真规则
– 用螺旋数 I^2=-N 表达涡旋
– 禁止：z * exp(iθ) * (1-nu*dt)
– 必须：SpiralVortex 存 (mag, phase, turns)
– 输出前断言：耗散下模长单调不增
– 涡量场先记 turns，再做 FFT/滤波

五、一点底层视角

Physical AI 的“几何幻觉”，本质就一句：

把连续螺旋生成，压扁成离散圆环 + 独立缩放。

IMU、世界模型、流体、天线相位、素数分布，全是同一类病。

I²=-N 不是数学装饰，是给 Agent 补一个带标度记忆的坐标系。

螺旋度不守恒的仿真，跑再久也是假运动。

六、延伸阅读（软性推广，真·开放获取）

螺旋系列专著，作者张智明，Zenodo（CERN 运营）CC 协议开放获取，点 DOI 直下 PDF：

• 体系总纲《螺旋生成论：从一条公理到万物生成》 https://doi.org/10.5281/zenodo.22074993
• 物理主推《旋生万物：从奇点到宇宙的统一生成论》（全三卷） https://doi.org/10.5281/zenodo.20408189
• 数学地基《螺旋数原理：公理系统与各向异性复数理论》 https://doi.org/10.5281/zenodo.20602099
• 计算工程《螺旋计算：量子计算的新基础》 https://doi.org/10.5281/zenodo.21356615
• 数论交叉《从几何构造到黎曼猜想》 https://doi.org/10.5281/zenodo.20995372
• 科普入门《螺线宗谱：发现宇宙的源代码》 https://doi.org/10.5281/zenodo.20659854
• 全集索引《螺旋生成论体系总汇编》 https://doi.org/10.5281/zenodo.21199593

作者 ORCID：https://orcid.org/0009-0003-7777-7694

不用信整套哲学。做物理仿真/流体/Agent 几何先验，先把 I²=-N 当编码约束用一周——回来你会发现：以前不是模型不会算，是它拿圆规画涡线。