Solving poisson's equation slowly rotating spherical shell of mass

In summary, the stress energy tensor for a perfect fluid is (I only care about the 0i components here which is where the rotational effects come into play anyways and I have assumed the time component of 4 - velocity is approximately 1 to first order which is fine because this is a slow motion approximation) T_{0i} = -\rho u{^{i}} = -\rho (\boldsymbol{\omega } \times r')^{i} = -\frac{M}{4\pi R^{2}}\delta(r' - R)(\boldsymbol{\omega } \times r')^{i} so in particular T_{0x} = -\
  • #1
WannabeNewton
Science Advisor
5,844
550
Solving "poisson's equation" slowly rotating spherical shell of mass

Homework Statement


We have that [itex]\partial ^{\alpha}\partial _{\alpha}\bar{\gamma _{0\mu}} = -16\pi T_{0\mu}[/itex] which is very similar to Poisson's equation if we treat each component of the metric tensor as a scalar field (the reason we only care about the time - space components is because the to be described mass - energy distribution has negligible pressure and rotates very slowly). The solution, in general, is given by [itex]\bar{\gamma _{0\mu}} = 4\int \frac{T_{0\mu}(x')dS(x')}{|x' - x|}[/itex]. I have a thin hollow slowly rotating spherical shell of uniform mass density [itex]\rho[/itex], radius [itex]R[/itex], and constant angular velocity [itex]\boldsymbol{\omega }[/itex] about the z - axis and I must solve for the components for such a case.

The Attempt at a Solution


The stress energy tensor for a perfect fluid is (I only care about the 0i components here which is where the rotational effects come into play anyways and I have assumed the time component of 4 - velocity is approximately 1 to first order which is fine because this is a slow motion approximation) [itex]T_{0i} = -\rho u{^{i}} = -\rho (\boldsymbol{\omega } \times r')^{i} = -\frac{M}{4\pi R^{2}}\delta(r' - R)(\boldsymbol{\omega } \times r')^{i}[/itex] so in particular [itex]T_{0x} = -\frac{M}{4\pi R^{2}}\delta(r' - R)\omega y' = -\frac{M}{4\pi R^{2}}\delta(r' - R)\omega r'sin\theta 'sin\varphi '[/itex] and [itex]T_{0y} = \frac{M}{4\pi R^{2}}\delta(r' - R)\omega x' = \frac{M}{4\pi R^{2}}\delta(r' - R)\omega r'sin\theta 'cos\varphi '[/itex] (the [itex]\omega[/itex] that shows up here is just the magnitude of [itex]\boldsymbol{\omega }[/itex])

Basically I'm solving for these functions [itex]\bar{\gamma _{0x}}, \bar{\gamma _{0y}}[/itex] that are pretty much like the scalar potential in EM and the solution is given by the usual solution to poisson's equation for a scalar potential as noted above. So, for example, [itex]\bar{\gamma _{0x}} = -4\int\frac{ \frac{M}{4\pi R^{2}}\delta(r' - R)\omega r'sin\theta 'sin\varphi 'dS'}{|r' - r|} = -4\int\frac{ \frac{M}{4\pi R^{2}}\delta(r' - R)\omega r'sin\theta 'sin\varphi 'r'^{2}sin\theta 'dr'd\theta 'd\varphi '}{|r' - r|}[/itex] and I have to evaluate this over the entire sphere. I'm not sure if I even set this up correctly and even if I did set it up correctly I'm also not sure how to go about actually evaluating that thing. [itex]|r' - r| = \sqrt{r'^{2} + r^{2} - 2rr'cos\beta _{r'r}} [/itex] but I don't know how to relate the angle between them, [itex]\beta _{r'r}[/itex] to the spherical coordinate angles so that I could get started off trying to integrate. Again, still not sure if I even set up the whole thing right. Thanks a ton guys!
 
Last edited:
Physics news on Phys.org
  • #2


WannabeNewton said:
[itex]|r' - r| = \sqrt{r'^{2} + r^{2} - 2rr'cos\beta _{r'r}} [/itex] but I don't know how to relate the angle between them, [itex]\beta _{r'r}[/itex] to the spherical coordinate angles so that I could get started off trying to integrate. Again, still not sure if I even set up the whole thing right.

I think you have it set up right although I haven't checked all the numerical factors. To deal with the angle [itex]\beta _{r'r}[/itex] issue, you can set up the problem initially in a coordinate system where the field point (observation point) lies on the z axis and the angular velocity vector lies in the x-z plane. Then [itex]\beta _{r'r}[/itex] is just the angle θ' between the source point r' on the sphere and the z-axis. Later you can convert back to the coordinate system where the angular velocity is along the z axis. Here is a tutorial for a similar calculation of the vector potential of a rotating charged spherical shell.
 
  • #3


I will try that thanks and thank you so much for the tutorial. I will see what I can do and come back with further questions if that is ok. I appreciate the help!
 
  • #4


Hi TSny. So first off can I start off by saying
[tex]-4\int\frac{ \frac{M}{4\pi R^{2}}\delta(r' - R)\omega r'sin\theta 'sin\varphi 'r'^{2}sin\theta 'dr'd\theta 'd\varphi '}{|r' - r|} = -4\int\frac{ \frac{M}{4\pi R^{2}}\omega R^{3}sin\theta 'sin\varphi 'sin\theta 'dr'd\theta 'd\varphi '}{|R - r|} = -\frac{MR}{\pi}\omega\int\frac{ sin^{2}\theta 'sin\varphi 'dr'd\theta 'd\varphi '}{|R - r|}[/tex]?
Secondly, so if I pick my coordinate system so that the field point is on the z - axis I would just have [itex]|R - r| = \sqrt{R^{2} + r^{2} - 2Rrcos\theta '}[/itex] correct? I feel like I'm missing one thing: if I do make this choice of coordinates how exactly will it affect the angular velocity vector? It will no longer be perfectly aligned along the z - axis so I have to go back to the start of my set up and redo the [itex]v^{i} = (\boldsymbol{\omega} \times r')^{i}[/itex] stuff right (so that my above expression for the integral would have to be changed to accommodate this if I'm not mistaken)? Thanks a ton!
 
  • #5


You have two factors of sinθ' in your integrand. But one of these is incorrect. Instead, you should have sinψ' as in equation 11 of the tutorial. To see where this comes from, you have to evaluate the cross product ω x r' where ω makes the angle ψ to the z-axis. ω x r' will have x, y, and z components, but you can show that only the y-component will contribute to the integral.
 
  • #6


I got the first sin factor from replacing the y coordinate from the cross product with it's spherical coordinate representation. The other came from the volume element. Was that not the right procedure? I'll quote it here: [itex]T_{0i} = -\rho u{^{i}} = -\rho (\boldsymbol{\omega } \times r')^{i} = -\frac{M}{4\pi R^{2}}\delta(r' - R)(\boldsymbol{\omega } \times r')^{i}[/itex] so in particular [itex]T_{0x} = -\frac{M}{4\pi R^{2}}\delta(r' - R)\omega y' = -\frac{M}{4\pi R^{2}}\delta(r' - R)\omega r'sin\theta 'sin\varphi '[/itex]. This was for when the angular velocity was aligned with the z - axis. Is what you are talking about when you make it lie on the x - z plane so that [itex]v = \omega \times r' = (-\omega _{z}y',\omega _{z}x' - \omega _{x}z', \omega _{x}y')[/itex]? Thanks.
 
Last edited:
  • #7


WannabeNewton said:
Is what you are talking about when you make it lie on the x - z plane so that [itex]v = \omega \times r' = (-\omega _{z}y',\omega _{z}x' - \omega _{x}z', \omega _{x}y')[/itex]?

Yes, that's right. Note that x' and y' have factors of cos[itex]\phi[/itex]' and sin[itex]\phi[/itex]' respectively. When you integrate over [itex]\phi[/itex]' what happens to these terms?
 
  • #8


Well integrating the [itex]cos\phi[/itex] terms will make them vanish and integrating the [itex]sin\phi[/itex] term will also make them vanish, and the only term without any of these is the [itex]-\omega _{x}z'[/itex] in the y - component, so the only surviving integral should be
[tex]\bar{\gamma _{0y}} = 4\int \frac{T_{0y}dS'}{|r' - r|} = -4\int \frac{\frac{M}{4\pi R^{2}}\delta (r' - R)\omega _{x}z'r'^{2}dr'd\Omega '}{|r' - r|} = -4\int \frac{\frac{M}{4\pi R^{2}}\delta (r' - R)\omega _{x}r'^{3}cos\theta 'sin\theta 'dr'd\theta '}{\sqrt{r'^{2} + r^{2} - 2r'rcos\theta '}}[/tex]
I think.

I am also unsure about where the angle [itex]\Psi[/itex] comes in here. I didn't seem to be able to extract any such angle from the cross product. Thank you!
 
  • #9


WannabeNewton said:
Well integrating the [itex]cos\phi[/itex] terms will make them vanish and integrating the [itex]sin\phi[/itex] term will also make them vanish, and the only term without any of these is the [itex]-\omega _{x}z'[/itex] in the y - component, so the only surviving integral should be
[tex]\bar{\gamma _{0y}} = 4\int \frac{T_{0y}dS'}{|r' - r|} = -4\int \frac{\frac{M}{4\pi R^{2}}\delta (r' - R)\omega _{x}z'r'^{2}dr'd\Omega '}{|r' - r|} = -4\int \frac{\frac{M}{4\pi R^{2}}\delta (r' - R)\omega _{x}r'^{3}cos\theta 'sin\theta 'dr'd\theta '}{\sqrt{r'^{2} + r^{2} - 2r'rcos\theta '}}[/tex]

Yes, that looks good.
I am also unsure about where the angle [itex]\Psi[/itex] comes in here.
The angular velocity ω lies in the x-z plane and makes an angle ψ to the z-axis. You need ωx. See the figure in the tutorial.
 
  • #10


Oh I see what you mean. So yes [itex]\omega _{x} = \omega sin\Psi [/itex] (where [itex]\omega[/itex] is the same angular velocity magnitude as from when it was aligned on the z axis as rotations preserve norm) I think. Also sorry but I think I forgot a factor of [itex]\pi[/itex] in front of the 4 that is outside the integral otherwise I have a leftover pi. I ran the integral through wolfram and then evaluated the bounds and hopefully if I didn't make a mistake, I got [itex]\bar\gamma {_{0y}} = -\frac{2M}{3R}r\omega sin\Psi [/itex]. Am I supposed to rotate back now?

EDIT: I should add that I chose the signs for the square roots resulting from the integration in accordance with r <= R because the problem only asked for the field when within the shell.
 
Last edited:
  • #11


OK. I think that's basically it, although if I'm not mistaken the factor of 2/3 should be 4/3. Anyway, note that you can write ##r\omega sin\psi## (in the negative y-direction) as ##\vec{\omega}\times\vec{r}##. This is now in a form that is independent of coordinate system. So, you can think of it as applying to a coordinate system where ω is along the z-axis and the observation point r is at angle θ = ψ to the z-axis. See equations 15 and 16 in the tutorial.
 
Last edited:
  • #12


I'll check my integration again to make sure I didn't make any arithmetic mistakes regarding that factor of 2. However going along with you, so we can now write [itex]\vec{{\gamma }} = (\bar{\gamma _{0x}}, \bar{\gamma _{0y}}, \bar{\gamma _{0z}}) = (0,-\frac{4M}{3R}\omega rsin\Psi , 0) = \frac{4M}{3R}(\vec{\omega} \times \vec{r})[/itex]? And that last expression is coordinate independent so do I have the freedom to then re - express it in coordinates it by shifting [itex]r[/itex] back to its original observation angle and shift the angular velocity back to the z - axis? Thanks!
 
  • #13


Yes. That's right.
 
  • #14


So I can rewrite it like this: [itex]\vec{r} = (x,y,z), \vec{\omega } = (0,0,\omega )[/itex] so that [itex]\vec{\gamma } = \frac{4M}{3R}(\vec{\omega} \times \vec{r}) = \frac{4M}{3R}(-\omega y,\omega x,0)[/itex] thus in particular [itex]\bar{\gamma} _{0x} = -\frac{4M}{3R}\omega y,\bar{\gamma} _{0y} = \frac{4M}{3R}\omega x ,\bar{\gamma} _{0z} = 0[/itex]. The vector potential for the gravito - magnetic field is defined as [itex]A_{i} = -\frac{1}{4}\bar{\gamma }_{0i}[/itex] therefore [itex]A_{x} = \frac{M}{3R}\omega y, A_{y} = -\frac{M}{3R}\omega x[/itex], [itex]A_{z} = 0[/itex]. The gravito - magnetic itself is [itex]B^{i} = \epsilon ^{ijk}\partial _{j}A_{k}[/itex] which gives [itex]B^{z} = -\frac{2M}{3R}\omega ,B^{x} = B^{y} = 0[/itex] so [itex]\vec{B} = -\frac{2M}{3R}\vec{\omega}[/itex]. I'm off by a sign from what the book says is the final answer but that isn't an issue I can just work back and see where I made the mistake. Can you just check if me re - expressing the angular velocity and field point position vector as I did above is a valid move? Thank you so much TSny!
 
  • #15


That looks correct to me.
 
  • #16


Ok one more thing then: I thought I understood it in the beginning but it seems I don't - why is it that aligning the arbitrary field point with the z - axis will cause the angular velocity, which was originally along the z - axis, to now be on only the x - z plane as opposed to some general direction?
 
  • #17


If I'm not mistaken, is it because, starting with the field point vector in an arbitrary direction and the angular velocity vector along the z - axis, we can first rotate about the z - axis until the field point vector is on the x - z plane, and this will still keep the angular velocity vector fixed along the z -axis because it is a rotation about the z - axis, and then rotate about the origin of the x - z plane until the field point vector is along the z - axis and now the angular velocity vector will lie at some angle on the x - z plane?
 
  • #18


That would be one way to get from the one coordinate system to the other. The main point is that we can choose the coordinate system anyway we wish. We just chose it so that ω would be in the x-z plane. We could just as well have chosen it such that ω was in the y-z plane. What made the problem easier was to choose the z-axis through the field point. That way, the angle between r and r' was just the spherical coordinate angle θ. Since the answer in this system could be expressed as a vector equation, the answer must be valid as a vector equation in any coordinate system. So, we write it out in any coordinate system we choose. We don't really need to worry about how to get from one system to the other.
 
  • #19


Yeah I was just about to ask why we couldn't have chosen to do it on the y - z plane and get the same answer. You answered that already however :). Well thank you so much, again, for all the help. Would you mind if I tried to summarize what we did just to make sure I got the procedure down?
 
  • #20


Not at all. Please do.
 
  • #21


Cool thanks! So we start off by re - orienting our axes, as described above, so that the field point vector is along the z - axis and the angular velocity is on the x - z plane (although we could have just as easily chosen the y - z plane). From here we expressed the components of the stress energy tensor in terms of this coordinate system, eliminated the components that would integrate to zero, and just kept the one surviving integral. We calculated the integral whichever way and found our result was [itex]\vec{\gamma } = -\frac{4M}{3R}\omega rsin\Psi \hat{y} [/itex] where [itex]\Psi [/itex] is the angle [itex]\vec{\omega }[/itex], which is on the x - z plane, makes with the z - axis.

Now this is the part I want to make sure I understand: we know from the way we oriented our coordinate system that the field point vector [itex]\vec{r}[/itex] is along the z - axis, and [itex]\vec{\omega }[/itex] is on the x - z plane, therefore [itex]-\omega rsin\Psi \hat{y}[/itex] is simply their cross product by the usual right hand rule argument (if [itex]\Psi [/itex] is negative then the right hand rule gives a positive y direction but [itex]sin\Psi [/itex] will be negative so that gives a negative y direction again and if [itex]\Psi [/itex] is positive then the right hand rule immediately gives a negative y direction correct?). Hence we can write [itex]-\omega rsin\Psi \hat{y} = \vec{\omega }\times \vec{{r}}[/itex] and [itex]\vec{\gamma } = \frac{4M}{3R}(\vec{\omega }\times \vec{{r}})[/itex]. This is of course coordinate indepedent so we have the freedom to re - express it any coordinate system we choose, in particular the original one we started with where the angular velocity vector was along the z - axis and the field point vector was off in some arbitrary direction.

Would that be accurate to say (especially that last paragraph)? Thank you TSny!
 

1. What is Poisson's equation for a slowly rotating spherical shell of mass?

Poisson's equation for a slowly rotating spherical shell of mass is a mathematical equation that describes the gravitational potential of a spherically symmetric mass distribution that is rotating at a slow rate. It takes into account both the mass and rotation of the shell to determine the gravitational potential at any point outside of the shell.

2. How is Poisson's equation derived for a slowly rotating spherical shell of mass?

Poisson's equation for a slowly rotating spherical shell of mass is derived from the general form of Poisson's equation, which relates the potential (Φ) and mass density (ρ) through the Laplacian operator (∇²). In the case of a slowly rotating spherical shell, the Laplacian operator is modified to take into account the effects of rotation, resulting in the specific form of the equation.

3. What are the assumptions made in solving Poisson's equation for a slowly rotating spherical shell of mass?

Some of the key assumptions made in solving Poisson's equation for a slowly rotating spherical shell of mass include: spherically symmetric mass distribution, slow rotation rate, and a vacuum outside of the shell. It is also assumed that the shell has a constant density and that the rotation is uniform throughout the shell.

4. What is the significance of solving Poisson's equation for a slowly rotating spherical shell of mass?

Solving Poisson's equation for a slowly rotating spherical shell of mass is important in understanding the gravitational potential and dynamics of rotating objects, such as planets or stars. It can also be used to model the potential of galaxies and other large-scale structures in the universe.

5. How is Poisson's equation for a slowly rotating spherical shell of mass used in practical applications?

Poisson's equation for a slowly rotating spherical shell of mass is used in many practical applications, such as in astrophysics, geophysics, and engineering. It is used to model the gravitational potential and dynamics of various systems, such as planets, stars, and satellites. It is also used in designing spacecraft trajectories and in studying the structure and evolution of galaxies.

Similar threads

  • Advanced Physics Homework Help
Replies
1
Views
757
  • Advanced Physics Homework Help
Replies
1
Views
1K
  • Advanced Physics Homework Help
Replies
1
Views
2K
  • Special and General Relativity
Replies
5
Views
223
  • Advanced Physics Homework Help
Replies
1
Views
716
  • Advanced Physics Homework Help
Replies
1
Views
249
  • Calculus and Beyond Homework Help
Replies
3
Views
474
  • Advanced Physics Homework Help
Replies
11
Views
1K
  • Advanced Physics Homework Help
Replies
6
Views
3K
  • Advanced Physics Homework Help
Replies
3
Views
1K
Back
Top