I am trying to use an energy argument for to show that the global Cauchy problem for the three-dimensional wave equation has a unique solution.
The wave equation is $$\partial^2u/\partial t^2=\nabla^2u$$
I looked up the energy functional and I want to use this $$E[u]=\frac1{2}\int_{R^3}((\partial u/\partial t)^2+\nabla u \nabla u))d^3 r$$
I have been reading what may be a similar question http://www.sgo.fi/~j/baylie/Partial%20Differential%20Equations%20in%20Action%20-%20From%20Modelling%20to%20Theory%20-%20S.%20Salsa%20(Springer,%202008)%20WW.pdf on page 265-266 but I can't figure out how to go to the 3D case.

I know the idea is the following:
1) define the energy integral as non negative non increasing function
2) for t =0 the energy is zero and so the energy is 0 for all t greater than or equal to 0. Then due to positivity of energy and the zero ic and bc we get solution as zero.
I just am struggling with how to find the derivative of the Energy integral above.
Thanks for your time.
The author uses the Leibniz integral rule, and this is possibly why he chose to expand the energy integral. For \begin{align*} e(t) & = \frac{1}{2}\int_0^{c(t_0 - t)}dr \int_{\partial B_r(x_0)}\left(u_t^2 + c^2|\nabla u|^2\right)d\sigma \\ & = \frac{1}{2}\int_0^{c(t_0 - t)}\left(\int_{\partial B_r(x_0)}\left(u_t^2 + c^2|\nabla u|^2\right)d\sigma\right)dr \end{align*} applying Leibniz integral rule and product rule gives \begin{align*} \dot{e}(t) & = \frac{1}{2}\left(\int_{\partial B_r(x_0)}\left(u_t^2 + c^2|\nabla u|^2\right)d\sigma\right)\bigg|_{r = c(t_0 - t)} \left(\frac{d}{dt}\left[c(t_0 - t)\right]\right) \\ & \qquad + \frac{1}{2}\int_0^{c(t_0 - t)}\frac{\partial}{\partial t}\left(\int_{\partial B_r(x_0)}\left(u_t^2 + c^2|\nabla u|^2\right)d\sigma\right)dr \\ & = -\frac{c}{2}\int_{\partial B_{c(t_0 - t)}(x_0)}\left(u_t^2 + c^2|\nabla u|^2\right)d\sigma + \frac{1}{2}\int_0^{c(t_0 - t)}\int_{\partial B_r(x_0)}\left(2u_tu_{tt} + 2c^2\nabla u\cdot\left(\nabla u\right)_t\right)d\sigma dr \\ & = -\frac{c}{2}\int_{\partial B_{c(t_0 - t)}(x_0)}\left(u_t^2 + c^2|\nabla u|^2\right)d\sigma + \int_{B_{c(t_0 - t)}(x_0)}\left(u_tu_{tt} + c^2\nabla u\cdot\nabla (u_t)\right)d\mathbf{x} \end{align*} The rest should be pretty straightforward I hope. The crucial difference compared to the case of bounded domain is that divergence theorem doesn't apply to infinite domain, so we need to consider the energy integral $e(t)$ over some bounded domain and do a limiting argument. In the case of wave equation, instead of considering $e(t)$ over the sphere we exploit the finite propagation property and consider the backward cone.