I am trying to show, with the following PDE: $$\begin{cases}\Delta u=f &\text{on}\>\>\Omega\\ u=g &\text{ on }\partial\Omega \end{cases}$$ that for all $x_0\in\Omega$, we have: $$u(x_0) = \iiint_{\Omega}G(x,x_0)f(x) dx+\iint_{\partial\Omega}g(x)\frac{\partial}{\partial \vec n}G(x,x_0)dS$$ Where $\vec n$ is the normal vector to $\partial\Omega$. Now, we have defined a Green's Function to be: $$G(x,x_0) = H(x)+\Phi(x)$$ Where $\Phi(x)$ is the Fundamental Solution in the respective Dimension. We also have the characteristics of $G$ and $H$, $H$ is a smooth Harmonic function in all of $\Omega$, and $G(x,x_0)=0$ on $\partial\Omega$.
For simplicity let us assume we are in Dimension $3$. Then, I believe we can solve this by Superposition. We had previously solved: $$\begin{cases}\Delta u=0 &\text{on}\>\>\Omega\\ u=g &\text{ on }\partial\Omega \end{cases}$$ The solution was given by: $$u(x_0)=\iint_{\partial\Omega}g(x)\frac{\partial}{\partial\vec n}G(x,x_0)dS$$ Good. This is half of what we need to show. Now, by Superposition I need to find the answer to: $$\begin{cases}\Delta u=f &\text{on}\>\>\Omega\\ u=0 &\text{ on }\partial\Omega \end{cases}$$ I started by applying Green's Second Identity to $H(x)$ and $u$ on $\Omega$: $$\iiint_{\Omega}H\Delta u-\Delta H u dx = \iint_{\partial\Omega}H\frac{\partial u}{\partial\vec n} dS_x-\iint_{\partial\Omega}\frac{\partial H}{\partial \vec n} u dS_x$$ Now, we know that $\Delta H$ is zero on $\Omega$. I also know that $u$ and all its derivatives are $0$ on $\partial\Omega$ by Superposition. So I am left with: $$\iiint_{\Omega}H\Delta u dx=\iiint_{\Omega}G\Delta u dx -\iiint_{\Omega}\Phi(x) \Delta u dx$$ Now I know above that $\Delta u = f$. All I need to show is that the above second term on the right hand side goes to $0$. But how do I show this??
You say that you know that $u$ is zero on the boundary $\partial \Omega$ and that you want to show that:
$$ \int_\Omega \Phi(x) \Delta u \ \mathrm{d}x = 0$$
This is just a classic application of Green's billionth identity. We have:
$$ \int_\Omega \Phi(x) \Delta u \ \mathrm{d} x = \int_\Omega u(x) \Delta \Phi \ \mathrm{d} x+ \int_{\partial\Omega} \Phi (\nabla u \cdot \mathbf{\hat{n}}) - u (\nabla \Phi \cdot \mathbf{\hat{n}}) \ \mathrm{d} S$$
Since $\Delta \Phi = 0$ in $\Omega$ and $u = 0$ on $\partial \Omega$, this leaves us with the following integral:
$$ \int_{\partial\Omega} \Phi (\nabla u \cdot \mathbf{\hat{n}}) \ \mathrm{d} S $$
Note that we have the following identity by the product rule:
$$ \nabla(\Phi u) = \Phi \nabla u + u \nabla \Phi $$
So that $\Phi \nabla u = \nabla(\Phi u ) - u\nabla \Phi $. If we now plug this into our integral we get:
$$ \int_{\partial\Omega} \Phi (\nabla u \cdot \mathbf{\hat{n}}) \ \mathrm{d} S = \int_{\partial\Omega} \nabla(\Phi u ) \cdot \mathbf{\hat{n}} \ \mathrm{d} S - \int_{\partial\Omega} u \nabla\Phi \cdot \mathbf{\hat{n}} \ \mathrm{d} S = 0$$
The first integral vanishes, since we are integrating the conservative vector field $\nabla (\Phi u)$ over a closed curve $\partial\Omega$, and the second integral vanishes since $u = 0$ on $\partial \Omega$. The conclusion is then:
$$ \int_\Omega \Phi(x) \Delta u \ \mathrm{d}x = 0$$
There you go! This is exactly what you wanted to show.