Meep supports the simulation of Maxwell's equations in cylindrical coordinates, albeit only for structures that have continuous rotational symmetry around the z axis. This reduces problems in 3d to 2d, and 2d to 1d, if you have sufficient symmetry.
See: