Fenics documentation or demonstration using to solve Poisson-Nernst-Planck equation

Someone else is trying to solve the same problem in FEniCSx so perhaps you can learn from them at How to debug the PETSC ERROR?. I may be solving coupled potential and concentration problem at the turn of the year. Good luck.