Let $\Omega$ be a smooth bounded domain of $\mathbb{R}^N$ where $N\geq 2$.
Consider the Laplace equation with a Neumann boundary condition
$$
-\Delta u = 0 \quad\mbox{in } \Omega, \qquad
\frac{\partial u}{\partial \mathbf{n}} = \varphi(x) \quad\mbox{on } \partial \Omega,
$$
where $\varphi$ is a sufficiently smooth function defined on $\partial\Omega$
for which the solutions are in $C^{2}(\overline{\Omega})$.
Assume the existence of a positive solution $u_0 > 0$ in $\overline{\Omega}$,
meaning that the solution $u_0$ is strictly positive in $\overline{\Omega}$.
My question is that: can $u_0$ be extended to some domain $\Omega^\prime \supset \overline{\Omega}$ in the manner that $u_0 > 0$ in $\Omega^\prime$, and $-\Delta u_0 = 0$ in $\Omega^\prime$ ?
If yes, then I think the Harnack inequality provides us with the assertion that $\sup_{\Omega}u \leq c \inf_{\Omega}u$ for any positive solution $u$ of this problem, where $c>0$ is a constant which does not depend on $\varphi$. Is it true ?