gt.geometric topology – Harmonic functions on knot complements

In Axler’s Harmonic Function Theory, he and his coauthors develop the theory of harmonic polynomials on $mathbb{R}^n$ by considering the restrictions of arbitrary polynomials on the sphere $S^{n-1} = {x in mathbb{R}^n : ||x||^2 = 1 }$ and taking the Poisson integral to get a harmonic polynomial in the interior ball. One can then take the Kelvin transform to get a harmonic polynomial on the exterior of the sphere. This process yields a canonical projection $mathscr{P}(mathbb{R}^n) to mathscr{H}(mathbb{R}^n)$, from the space of polynomials to the space of harmonic polynomials, factoring through the restriction map to $L^2(S^{n-1})$. By taking polynomial approximations via the Stone-Weierstrass theorem, and verifying some series convergence, I believe one can expand the domain and codomain of the projection to larger, more general function spaces. I have not verified this myself, so I could be wrong.

Does this theory generalize to knot complements? Say we have a knot $K subseteq mathbb{R}^3$, and we take a small tubular neighborhood $V$ around $K$, whose boundary is topologically a torus $T$. Given a function on the knot complement, one could restrict to $T$ and then solve the Dirichlet problem on the knot complement to get a projection like the one above. However, in the sphere case, there are many nice properties of the projection; namely it comes with an efficient algorithm for computation which involves repeatedly differentiating the function $f(x) = |x|^{2-n}$.

Is anyone aware of any theory along this vein? Are there any obstacles to generalizing what happens in the sphere case?