Perturbation of a uniform electrostatic field by a dielectric cube

[Apteronotus]

Hi,

Is there any way to analytically calculate the perturbation of a uniform electrostatic field by a dielectric cube.
I know a solution exists for dielectric spheres but I haven't been able to come across the solution, when dealing with a cube.

Ohh.. and I'm assuming the simplest case, where the field is parallel to the normal of one of the cube's sides.

Thanks in advance.

 

[Vanadium 50]

It's going to be a mess, because you won't have Legendre polynomials in your solution. It's equivalent to solving the electric field from a cube-shaped dipole, and I would attack that by starting with the electric field from a small dipole of length a, and making a cube of length L out of many aligned dipoles, and taking the limit as a/L -> 0.

 

[vanhees71]

Instead of spherical harmonics, I'd use the eigenfunctions of the Laplace operator from the separation ansatz in Cartesian coordinates with the appropriate boundary conditions for the dielectric.

 

[Vanadium 50]

But the variables aren't going to separate nicely in Cartesian space.

 

[Gordianus]

Check Smythe's "Static and Dynamic Electricity". I've found the answer to many complex situations there.

 

[vanhees71]

Of course, the Laplacian separates in Cartesian coordinates. For the Laplace Equation,

\Delta \Phi=0.

you make the ansatz \Phi(\vec{x})=X(x) Y(y) Z(z). Then in the usual way, you can show that

\frac{X''(x)}{X(x)}=-k_1^2, \quad \frac{Y''(y)}{Y(y)}=-k_2^2, \quad \frac{Z''(z)}{Z(z)}=k_1^2+k_2^2

Then you get the set of solutions

\Phi_{k_1,k_2,\pm}(\vec{x})=\exp \left (\mathrm{i} k_1 x+\mathrm{i} k_2 y \pm \sqrt{k_1^2+k_2^2} z \right ).

A general harmonic function can be built out of these also in the usual way. Of course, for the present problem, you have to find the restriction due to the boundary conditions imposed by the presence of the dielectric cube first and then you can linearly superpose the remaining set of functions to find the solution for your situation.

 

[Vanadium 50]

Of course, for the present problem, you have to find the restriction due to the boundary conditions imposed by the presence of the dielectric cube

But that's the point - that's the part that doesn't end up looking so nice. The boundary conditions are ugly, because they are all finite. With spherical or cylindrical symmetry things are much cleaner.

 

[Apteronotus]

Following the description above and solving \frac{X''}{X}=-k_1^2, I find
X(x)=A_1e^{k_1 ix}+A_2e^{-k_1 ix}
similarly
Y(y)=B_1e^{k_2 iy}+B_2e^{-k_2 iy} and
Z(z)=C_1e^{k_3 z}+C_2e^{-k_3 z}, as general solutions.

So,

\Phi(\vec{x})=<br /> \left[(A_1+A_2)\cos k_1x+i(A_1-A_2)\sin k_1x\right]\cdot<br /> \left[(B_1+B_2)\cos k_2x+i(B_1-B_2)\sin k_2x\right]\cdot<br /> \left[C_1e^{k_3 z}+C_2e^{-k_3 z}\right]
rather than what is stated above by vanhees71.

Also for the first boundary condition I have found

\Phi(\vec{x})\rightarrow -E_0z as any of x,y \mbox{ or } z \rightarrow \pm \infty
which simply states that as we move far enough away from the cube, the potential approaches that which would exist without the cube present.

But I'm not sure how this BC would simplify the electric potential equation.

 

[Apteronotus]

Gordianus, thank you for the reference. I just went through Smyth and unfortunately there was nothing on the subject.

Can you think of any other references which may be useful?

 

[Apteronotus]

For any unfortunate sole that may come across this post hoping for a hint as to how to approach this problem, I would ignore the above posts.

The problem can not be solved using the Cartesian form of Laplace. The solution is attained using Green's theorem.