BEHIND THE ART
先构造流函数,再求速度。
流函数由两组不同尺度的正弦模式相加:ψ=sin(kx+0.3t)sin(ky−0.2t)/k + 0.35b sin(2kx−0.4t)cos(2ky+0.1t)/(2k)。解析求导得到速度,两项混合偏导相消,因而理论散度为零。
粒子轨迹由离散时间积分近似;粒子重生与有限采样会改变屏幕上的点密度,不能据此直接判断连续速度场是否有散度。本作品不是求解带压力和黏性的完整流体方程。
离开取样窗口或达到寿命的粒子会重新播撒。显示轨迹会在重生时清空,避免把两次出生位置错误连接;亮度只用于呈现轨迹新旧,不是物理能量。