How to solve coupled PDEs for electrostatics?

“Liked” this post because of its similarity to the problem I am working on currently: setting a far-field boundary condition to be zero in a coupled PDE system returns the trivial result for that variable.

Eagerly watching the space for any help/insight.

Thank You
Warm Regards