Solution (source code)

= Solution

Taking the curl of Euler and using incompressibility gives $\partial_t\omega+(u\cdot\nabla)\omega=(\omega\cdot\nabla)u$, equivalently $D\omega/Dt=(\omega\cdot\nabla)u$.