Take $\Omega$ to be a bounded smooth domain with boundary $\partial\Omega = \Gamma_1 \cup \Gamma_2$, where $\Gamma_1$ and $\Gamma_2$ are disjoint.
Consider the problem
$$\Delta u = f \quad\text{in $\Omega$}$$ $$u = g \quad\text{on $\Gamma_1$}$$ $$\partial_\nu u = h \quad\text{on $\Gamma_2$}$$ where $f, g, h$ are given data such that $f,g,h \in L^\infty\cap H^1$ on the respective domains and $\nu$ is the outward-pointing unit normal vector. The data can be taken to be more regular Sobolev functions if needed.
Does anyone know a reference or method of proof to obtain an $L^\infty(\Omega)$-estimate on the weak solution $u$, in terms of the data?