After a couple of months without a new post I've started writing about simulating charge distributions, rather than just boundary value problems. I'll explain how #FEniCS can calculate the electrostatic fields of arbitrarily shaped charge densities, which are specified with mesh subdomains. These will have hard edges for now but I want to look into having more general charge distributions.
I've also run 3D simulations for calculating capacitances between arbitrarily shaped conductors and 2D linear dielectrics. I need to do some more testing & comparison with other FEA tools/ analytical methods before I publish any posts. Seems to be going well though!