The classical point vortex model corresponding to the 2D incompressible Euler equation in vorticity form is the system of $N$ ODES

$$\begin{cases} \dot{x}_i(t) = \sum_{1\leq i\neq j\leq N} a_j K(x_i(t),x_j(t)) \\ x_i(0) = x_i^0\end{cases}\tag{1}$$

where $K(x,y) := -\frac{1}{2\pi}(\frac{(x-y)_2}{|x-y|^2}, -\frac{(x-y)_1}{|x-y|^2})$ is the Biot-Savart kernel.

It is evident from the RHS of equation $\ref{a}$ that the point vortices only have binary (i.e. pairwise) interactions. My question is the following: are there physically motivated models generalizing equation \ref{a} which have binary and ternary interactions? By ternary, I have in mind a term of the form

$$\sum_{1\leq i\neq j\neq k\leq N} a_j a_k \tilde{K}(x_i(t),x_j(t),x_k(t)),$$

where $\tilde{K}$ is an $\mathbb{R}^2$-valued map.