This question is of pure interest. I would like to know, how a mixed boundary value problem like the following can be solved numerically:
Lets say I have two conducting plates (not necessarily parallel) in two dimensions (infinite in one, so 3D reduces to 2D). On the plates are total charges $Q_1$ and $Q_2$. Now I place a point charge $q$ somewhere between. How is such problem solved? I mean numerically, not analytically.
If potentials $\Phi_1$ and $\Phi_2$ on the plates would be given, and there is no additional charge, I know, that this can be solved by relaxation methods on a grid.
But in this situation the potentials on the plates are not known in advance, the only thing is to know, that they are constant on each plate and that the surface integral of the fields directly at the plates must yield charges $Q_1$ and $Q_2$. So I'm not able to start some kind of iteration with fixed values at the plates. Additionally, how are point charges considered in such case? The surface integral around the point charge must yield $q$, but how can this be transformed into a numerical method?
What literature can be recommended to learn such things? I have good books on classical electrodynamics (like Jackson, ...) but here the scope is to solve analytically.
Of course, there is FEM software available. But I want to know HOW available software solves such problems in principle. I don't need it for an urgent solution.