Can Laplace's Equation Be Solved in Cylindrical Coordinates on a 3-Sphere?

[Old Smuggler]

I wish to find exact solutions of Laplace's equation in cylindrical coordinates on (a subset of) the 3-sphere.
This pde is linear but not separable. The potential {\Phi}(x,z) must fulfil the following pde:

<br /> (1-{\frac{x^2}{a^2}}){\frac{{\partial}^2}{{\partial}x^2}}{\Phi}(x,z)+<br /> (1-{\frac{z^2}{a^2}}){\frac{{\partial}^2}{{\partial}z^2}}{\Phi}(x,z)+<br /> {\frac{1}{x}}(1-{\frac{3x^2}{a^2}}){\frac{{\partial}}{{\partial}x}}{\Phi}(x,z)-<br /> {\frac{2xz}{a^2}}{\frac{{\partial}^2}{{\partial}x{\partial}z}}{\Phi}(x,z)-<br /> {\frac{3z}{a^2}}{\frac{{\partial}}{{\partial}z}}{\Phi}(x,z)=0<br />

Here a is a constant (and x,z&lt;a, z{\neq}0). Does anyone know how to solve this equation?
(I'm aware that a transformation of this equation to spherical coordinates yields a separable pde, but
this gives a bunch of useless solutions blowing up near the origin.)

 

[HallsofIvy]

No, it doesn't. One of the "boundary conditions" for a problem like that is that it NOT go to infinity as r goes to 0.

The general solution to the Laplace's equation in spherical coordinates involves a double sum, over l and m, of things like (a_mlr^l+ b_ml/r^l+1)Y_l^m(theta, phi) where Y_l^m(theta, phi) are the spherical harmonics. If your region includes the origin, then one of the conditions must be that all b_ml are 0.

 

[Old Smuggler]

No, it doesn't. One of the "boundary conditions" for a problem like that is that it NOT go to infinity as r goes to 0.

The general solution to the Laplace's equation in spherical coordinates involves a double sum, over l and m, of things like (a_mlr^l+ b_ml/r^l+1)Y_l^m(theta, phi) where Y_l^m(theta, phi) are the spherical harmonics. If your region includes the origin, then one of the conditions must be that all b_ml are 0.

Let me state the problem in more detail. I wish to find the Newtonian potential outside a flat, very thin disk lying in
the surface z=0 in cylindrical coordinates. If one first solves Laplace's equation in cylindrical coordinates in ordinary
Euclidean space, the solution is given by (from MAPLE)
<br /> {\exp}(-kz)[C_1J_0(kx)+C_2Y_0(kx)] , \qquad (z&gt;0), <br />
where k, C_1, C_2 are constants and J_0, Y_0 are Bessel functions of the first and second kind, respectively (C_2 is set to zero for physical reasons). Since k is a continuous parameter, solutions can be summed by integration over k, and using Gauss theorem over the disk, and moreover Hankel transforms, we can find the potential in the disk for any given density distribution. All this is standard stuff. On the other hand, if one changes to spherical coordinates, the solution is
<br /> {\Big [}C_3r^{-1/2+{\sqrt{1+4C}}/2}+C_4r^{-1/2-{\sqrt{1+4C}}/2}]{\Big ]}<br /> {\Big [}C_5P^0_{-1/2+{\sqrt{1+4C}}/2} ({\cos}{\theta})+C_6Q^0_{-1/2+{\sqrt{1+4C}/2}}({\cos}{\theta}){\Big ]} <br />,
where C, C_3,C_4, C_5, C_6 are constants, and P^{\mu}_{\nu}, Q^{\mu}_{\nu} are Legendre functions of the first and second kind, respectively. However, if one reinserts cylindrical coordinates into this solution, it is not clear to me how to get back the above solution involving Bessel functions in a simple manner, and how the constants k and C are related.

Now for the similar problem on a subset of the 3-sphere. Laplace's equation on a subset of the 3-sphere is separable in spherical coordinates, and again the solution involves Legendre functions (but this time with non-zero order).
However, in cylindrical coordinates the equation is non-separable. If there is a simple way of taking the solution for the equation in spherical coordinates, reinsert cylindrical coordinates into this solution and then get something that resembles the flat space solution involving Bessel functions, that would be great. However I don't know how to do that even in the Euclidean case. So it seems that I have to find solutions of the non-separable pde after all.

 

[Physicslad78]

Guys I need help to solve this differential equation:

(1-x^2)~\frac{\partial ^2 F(x,y)}{\partial x^2}+\left(\frac{1-2x^2}{x}\right)~\frac{\partial F(x,y)}{\partial x}-\left[ax^2(1-2y^2)-c\right]~F(x,y)=0 where a and c are constants and x and y are the variables

Does this have a solution in terms of well known functions? Mathematica could not find any so I proposed a non separable series solution of the form :

\sum_{n}\sum_{p}C_{n,p} x^n y^p

But the problem is when I tried to find the eigenvalues of the formed matrix it got too complicated and c (which is a constant in terms of the energy of the system) turned out to be complex in some cases! I want real values as the system would not be physical..can anyone help me in this matter please.

Thanks