Can Laplace's equation be solved in a domain that is not simply connected? I have a problem where the domain is like a box with a tube missing - e.g.
0< x<1,0< y<1, 0< z<1 less the region (x-0.5)^2+(y-0.5)^2 < 0.25
In order to solve Laplace's equation in this domain it was mapped to one which can be thought of as simply a box. This was done because the initial domain was not simply connected - and so Laplace's equation could not be solved in it. Is this correct? I haven't been able to find a simple clarification of this online.
In this link on wolframalpha, http://mathworld.wolfram.com/SimplyConnected.html, my initial domain is like the far right image - can Laplace's equation be solved in a domain such as that?

There is no topological obstacle to solving Laplace's equation with Dirichlet data (i.e. if the value of the solution on the boundary is prescribed). The solution may however lose some regularity near corner singularities (in your case near $x^2 + y^2 = .25, z \in \{0,1\}$, and at the edges of the outer box). 
In the case of Neumann data (where the normal derivative of the solution is prescribed), solutions are only unique up to a constant which may be different for each connected component, and the boundary data must satisfy compatibility conditions, one for each connected component. So the only topological constraint comes from the number of components of the domain, which is 1 in your case.   
