Consider nice domain $D\subset \mathbb R^d$ and $\Delta u =0$ with $u\big|_{\partial D}=g$. It is well known that $u(x)=E^x[g(B(\tau))]$ where $\tau$ is exit time of $B$ from the domain $D$.

What if we consider Neumann boundary conditions? Is there a stochastic representation?