SCHEMES AND ESTIMATES FOR THE LONG-TIME NUMERICAL SOLUTION OF MAXWELL'S EQUATIONS FOR LORENTZ METAMATERIALS

We consider time domain formulations of Maxwell's equations for the Lorentz model for metamaterials. The field equations are considered in two different forms which have either six or four unknown vector fields. In each case we use arguments tuned to the physical laws to derive data-stability estimates which do not require Gronwall's inequality. The resulting estimates are, in this sense, sharp. We also give fully discrete formulations for each case and extend the sharp data-stability to these. Since the physical problem is linear it follows (and we show this with examples) that this stability property is also reflected in the constants appearing in the a priori error bounds. By removing the exponential growth in time from these estimates we conclude that these schemes can be used with confidence for the long-time numerical simulation of Lorentz metamaterials.