“…Other articles have offered methods for simulating isotropic dielectrics with sharp interfaces: for example [10][11][12] describe methods that use effective diagonal dielectrics, while [13] uses an effective tensor dielectric in 2D simulations. There are also methods that, rather than using an effective dielectric in each grid cell, alter the shapes of the cells to coincide with the dielectric shape, such as [14], which has interesting parallels to this paper, running into an accuracy versus stability trade-off, thwarted by the need to collocate different field components.…”