Second order approximations for kinetic and potential energies in Maxwell's wave equations

In this paper we propose a numerical scheme for wave type equations with damping and space variable coefficients. Relevant equations of this kind arise for instance in the context of Maxwell's equations, namely, the electric potential equation and the electric field equation. The main motivation to study such class of equations is the crucial role played by the electric potential or the electric field in enhanced drug delivery applications. Our numerical method is based on piecewise linear finite element approximation and it can be regarded as a finite difference method based on non-uniform partitions of the spatial domain. We show that the proposed method leads to second order convergence, in time and space, for the kinetic and potential energies with respect to a discrete L2-norm.