Use triple integral to find center of mass

Click For Summary
SUMMARY

The discussion focuses on calculating the center of mass (CoM) of a uniform hemispherical shell with inner radius a and outer radius b using triple integrals in spherical coordinates. The participants clarify the use of Cartesian coordinates in the context of spherical coordinates and emphasize the need to treat the CoM as a 3D vector. The final integrals for the CoM are expressed as rCoM, θCoM, and φCoM, with the integrals correctly set up to account for the density ρ and the volume element in spherical coordinates.

PREREQUISITES
  • Understanding of triple integrals in calculus
  • Familiarity with spherical coordinates and their conversion to Cartesian coordinates
  • Knowledge of vector calculus and the concept of center of mass
  • Basic principles of density and volume in physics
NEXT STEPS
  • Study the derivation of the center of mass using triple integrals in spherical coordinates
  • Learn about the conversion formulas between Cartesian and spherical coordinates
  • Explore examples of calculating the center of mass for different geometric shapes
  • Investigate the implications of symmetry in center of mass calculations
USEFUL FOR

Students in physics and engineering, particularly those studying mechanics and calculus, as well as educators looking for examples of applying integrals to real-world problems involving center of mass.

BearY
Messages
53
Reaction score
8

Homework Statement


Find the centre of mass of a uniform hemispherical shell of inner radius a and outer radius b.

Homework Equations


##r_{CoM} = \sum \frac{m\vec{r}}{m}##

The Attempt at a Solution


Using ##x(r,\theta,\phi)## for coordinates,
$$x_{CoM}=\frac{\int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} \vec{x}\rho r^2\sin{\theta}drd\theta d\phi}{\int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} \rho r^2\sin{\theta}drd\theta d\phi}$$
My vector calculus is rusty, how do I handle the ##\vec{x}## in this integral?
 
Physics news on Phys.org
hvIt may help to first note that the LHS of the equation should be ##\vec x_{CoM}##. That is, it is a 3D vector that gives the location of the CoM.
Given that, treat the RHS as three separate integrals, one for each of the cartesian coordinates of the CoM: ##x,y## and ##z##.

##x_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} x(r,\theta,\phi) \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
##y_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} y(r,\theta,\phi) \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
##z_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} z(r,\theta,\phi) \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
where ##D## is the denominator in the formula in the OP.
Then use the formulas for converting from Cartesian to spherical coords to replace the ##x(r,\theta,\phi)## in the integrand by a function of ##r,\theta,\phi##. THen do the same for the integrals for ##y## and ##z##.
 
  • Like
Likes   Reactions: BearY
andrewkirk said:
hvIt may help to first note that the LHS of the equation should be ##\vec x_{CoM}##. That is, it is a 3D vector that gives the location of the CoM.
Given that, treat the RHS as three separate integrals, one for each of the cartesian coordinates of the CoM: ##x,y## and ##z##.

##x_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} x(r,\theta,\phi) \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
##y_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} y(r,\theta,\phi) \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
##z_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} z(r,\theta,\phi) \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
where ##D## is the denominator in the formula in the OP.
Then use the formulas for converting from Cartesian to spherical coords to replace the ##x(r,\theta,\phi)## in the integrand by a function of ##r,\theta,\phi##. THen do the same for the integrals for ##y## and ##z##.
I am confused about what Cartesian coordinates do here. I thought I am just looking for a point in the spherical coordinate?
 
BearY said:
I am confused about what Cartesian coordinates do here. I thought I am just looking for a point in the spherical coordinate?
Quite right. Sorry, I was distracted by the use of the character 'x' in the vector ##\vec x##.
So we can do it more simply as:

##r_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} r \,\rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
##\theta_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} \theta\, \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
##\phi_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} \phi\, \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
 
  • Like
Likes   Reactions: BearY
andrewkirk said:
Quite right. Sorry, I was distracted by the use of the character 'x' in the vector ##\vec x##.
So we can do it more simply as:

##r_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} r \,\rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
##\theta_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} \theta\, \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
##\phi_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} \phi\, \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
When I was studying vector calculus we only did parameterizing dimension x, y, and z with t and then integrate w.r.t dt only.
So without parameterization, we still integrate each dimension separately? That makes sense. I suppose I will have to redo vector calculus to know for sure.
 
andrewkirk said:
Quite right. Sorry, I was distracted by the use of the character 'x' in the vector ##\vec x##.
So we can do it more simply as:

##r_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} r \,\rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
##\theta_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} \theta\, \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
##\phi_{CoM} = \int_{0}^{\frac{\pi}{2}}\int_{0}^{2\pi}\int_{a}^{b} \phi\, \rho r^2\sin{\theta}\,dr\,d\theta\, d\phi\ /\ D##
I don't think this works because generally ##\langle \cos\theta \rangle \ne \cos\langle \theta \rangle##. Also, for a hemisphere of radius ##R##, the center of mass is on the z-axis at ##z=\frac {3}{16}{R}##, but the integral above with ##a=0## and ##b=R## gives ##r_{CoM} = \frac 38 R##.
 
  • Like
Likes   Reactions: BearY
@BearY: If that was a typical homework problem I would expect the student to consider the upper hemisphere at the origin, notice immediately that ##\bar x = \bar y = 0## from symmetry, and calculate$$
\bar z = \frac {\iiint_V z~dV}{\iiint_V 1~dV} =
\frac {\int_0^{2\pi}\int_0^{\frac \pi 2} \int_a^b \rho\cos\phi \rho^2\sin\phi~d\rho d\phi d\theta}
{\int_0^{2\pi}\int_0^{\frac \pi 2} \int_a^b 1 \rho^2\sin\phi~d\rho d\phi d\theta}$$
Notice the integral is conveniently worked in spherical coordinates, but you are calculating the center of mass ##(\bar x, \bar y, \bar z)##. And it generally is not true, for example, that ##\bar z = \bar \rho \cos\bar\phi## just because ##z =\rho \cos\phi##. You must set up the rectangular integral then change coordinates.
 
  • Like
Likes   Reactions: BearY
BearY said:
I am confused about what Cartesian coordinates do here. I thought I am just looking for a point in the spherical coordinate?

Your spherical coordinates ARE relative to a particular Cartesian coordinate system! In your convention---basically, the "physics" rather than "mathematics" convention---##\theta## is the azimuthal angle (between the ##z##-axis and the point ##(x,y,z)##) while ##\phi## is the "longitude" ( the counter-clockwise angle between the ##+x## axis and ##(x,y)##, as viewed from ##z > 0## looking down).
 
  • Like
Likes   Reactions: BearY

Similar threads

Replies
3
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 9 ·
Replies
9
Views
1K
  • · Replies 19 ·
Replies
19
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 14 ·
Replies
14
Views
2K
  • · Replies 9 ·
Replies
9
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
Replies
14
Views
2K