# Energy conservation in the solution of the Helmholtz equation

This might be a silly question, but I know very little about the theoretical properties finite elements, so here goes. Suppose you were to solve the Helmholtz equation (let's say in 2D) with a spatially varying wave speed using finite elements. For reference the equation is $$[\nabla^2 + k^2(\vec{r})]\psi(\vec{r}) = f(\vec{r})$$ This being a scattering problem, the source term is limited to, say, Dirichlet values on the boundary of the domain. My question is, when the domain is large (say ten or more of wavelength (I consider this large, it might not be)), how well does energy conservation hold? More precisely, what kind of guarantees can you place on conservation of energy as a function of domain size?

Assume for now linear elements (in practice, we use finite differences, is there a huge difference?) since $k(\vec{r})$ is piecewise constant, so higher orders aren't terribly helpful. I hope I have not revealed too much of my ignorance on this subject!

Edit: I am referring to energy in the "intuitive" sense. I work mostly with Maxwell's equations, which in 2D at constant frequency becomes the Helmholtz equation. I believe the mathematical definition of energy flow at a point is $\Psi\nabla \Psi$, or something to that effect (it should be the Poynting vector). Energy conservation would say that in a source-free region, the integral of the energy in the surface normal direction along the boundary of the region should be identically zero. Intuitively, I should see waves obeying the inverse square law (or whatever its analog is in 2D) from a point source exactly.

Also, I realize there is a difference between finite differences and finite elements. If you could comment on both, that would be even better.

• On non-rectangular grids, in general, FEM and FDM should differ. Mar 6, 2013 at 12:17
• How is your energy defined? Mar 6, 2013 at 12:18
• @vanCompute: I have updated the question with some more details about energy. Mar 6, 2013 at 23:09
• I'm not a physicist. But should not the space-time density of energy be something like $(\Psi_t)^2+|\nabla_x \Psi|^2$ in time domain? Mar 7, 2013 at 8:15
• @Hui: perhaps it's something like that. I only mentioned energy flow in the question. Also note that we are in the frequency domain, so it would be like $|\Psi|^2 + |\nabla\Psi|^2$. Mar 7, 2013 at 8:49