Perhaps the conversation is getting a bit narrow, but I’m hoping it will catch some attention and gather some momentum.

The FDTD method relies on the “staggered” (or nested) grids for the dual quantities. That works well when the grid is Cartesian.

The same outcome can be achieved with FEM grids via a careful choice of basis functions or choice of weak form.