Local Biot-Savart law in $B(x_o,r) \subset \mathbb R^2$ Let $u: \mathbb R^2 \to \mathbb R^2$ and let $\omega = \text{curl } u$ be the 2D vorticity of $u$, where $u, \omega \in L^2(\mathbb R^2)$ and $\nabla \cdot u = 0$. The classical Biot-Savart law states that
$$
u_1(x) = -\frac1{2\pi}\int_{\mathbb R^2}\frac{x_2-y_2}{|x-y|^2}\omega(y) \, dy
$$
and
$$
u_2(x) = \frac1{2\pi}\int_{\mathbb R^2}\frac{x_1-y_1}{|x-y|^2}\omega(y) \, dy.
$$
I am interested in a local version of the law in $B(x_o, r) \subset \mathbb R^2$. More precisely, suppose $u: \Omega \to \mathbb R^2$, where $\Omega \subset \mathbb R^2$ is bounded, and $u, \omega \in L_{loc}^2(\Omega)$ with $\nabla \cdot u = 0$. Is it true that
$$
u_1(x) = -\frac1{2\pi}\int_{B(x_o,r)}\frac{x_2-y_2}{|x-y|^2}\omega(y) \, dy+A_1^{x_o}(x)
$$
and
$$
u_2(x) = \frac1{2\pi}\int_{B(x_o,r)}\frac{x_1-y_1}{|x-y|^2}\omega(y) \, dy+A_2^{x_o}(x).
$$
for all $x \in B(x_o,r) \subset \Omega$ and for some smooth-ish function $A^{x_o}$? If not this, is there something similar which is true?
Finally, I am also interested in the general case of $\mathbb R^d$ for $d >2$, but for simplicity we may first restrict to the case $d=2$.
 A: Remember how the classical Biot-Savart law is obtained : first, you take the curl of the defining equation $\omega = \nabla \wedge u$, and then, you use the fact that you know the solution of the Poisson equation in terms of the source. So, your problem boils down to computing explicitly some kernel.
Let's go through the proof in the 2D case, the necessary modifications for 3D and higher being easy.
First, extend $\omega$ by $0$ outside $\Omega$ and denote this extension by $\tilde{\omega}$. Now, $\tilde{\omega}$ is defined on the plane.
Let $\varphi$ be the unique solution in $L^2(\mathbb{R}^2)$ (actually also in $H^2$) of the Poisson equation $\Delta \varphi = \tilde{\omega}$. Define a new vector field $v := \nabla^{\perp} \varphi$, i.e. 
$$v_1 = -\partial_2 \varphi \quad , \quad v_2 = \partial_1 \varphi .$$
Obviously, $\nabla \cdot v = 0$ and $\nabla \wedge v = \tilde{\omega}$. Thus, on $\Omega$, $u$ and $v$ agree up to a harmonic polynomial, which is smooth anyway.
Now, take a compact subset $K$ of $\Omega$ and choose a fattening $K'$ of $K$, that is, another compact subset of $\Omega$ whose interior contains $K$. $K'$ may very well be a set of the type $K_{\varepsilon} := \{x \in \Omega ; \text{dist}(x,K) \leq \varepsilon \}$ for some small $\varepsilon$.
Write $v = v_1 + v_2$ with $v_1 := \int_{K'} \nabla_x \left( \ln |x-y| \right) \omega (y) \: dy$ and $v_2$ is the same integral, but over the complement of $K'$.
As $v_2$ is harmonic over $K$, with the boundary of $K$ uniformly far from that of $K'$, $v_2$ enjoys both smoothness and integrability over $K$.
Up to some possible harmonic polynomial, this may satisfy you, I hope.
(On a side note : are you working with or around Navier-Stokes and fluid mechanics ?)
