Just made a short post about how to simulate electrostatic nonuniform charge density distributions with FEniCS; that is spatially varying charge density distributions. This could be useful for simulating things like particle beams which are commonly assumed to have Gaussian (or similar) distributions. This is a small extension of a previous post about how to simulate uniform charge density distributions.