Algebraic equation -- six degree polynomial

[gruba]

Homework Statement

Solve the equation for r,r>0,r<R.
\frac{-2\pi R^3}{3}-\frac{8\pi r^2\sqrt{R^2-r^2}}{3}+2r^2R+\frac{2\pi}{3}R^2\sqrt{R^2-r^2}=\frac{2\pi R^3}{3}

2. The attempt at a solution
After factoring,

-2R(2\pi R^2-3r^2)=2\pi\sqrt{R^2-r^2}(4r^2-R^2).

After squaring,

64\pi^2 r^6+(36R^4-96\pi^2 R^2)r^4+(32\pi^2 R^4-44\pi^2 R^4)r^2+12\pi^2 R^6=0.

r is unknown variable, so the last equation should be the equation of six degree polynomial.
How to solve it?

 

[Ray Vickson]

Homework Statement

Solve the equation for r,r>0,r<R.
\frac{-2\pi R^3}{3}-\frac{8\pi r^2\sqrt{R^2-r^2}}{3}+2r^2R+\frac{2\pi}{3}R^2\sqrt{R^2-r^2}=\frac{2\pi R^3}{3}

2. The attempt at a solution
After factoring,

-2R(2\pi R^2-3r^2)=2\pi\sqrt{R^2-r^2}(4r^2-R^2).

After squaring,

64\pi^2 r^6+(36R^4-96\pi^2 R^2)r^4+(32\pi^2 R^4-44\pi^2 R^4)r^2+12\pi^2 R^6=0.

r is unknown variable, so the last equation should be the equation of six degree polynomial.
How to solve it?

Try to solve for the new variable $S = \sqrt{R^2 - r^2}$, by converting your original equation into a cubic in $S$. This is solvable in principle using Cardan's formula for the solution of a cubic. Furthermore, $x = S/R$ obeys a cubic equation that does not contain the parameter $R$, so we can take $R = 1$ without loss of generality.

All solution values for $r$ are complex-valued, so your original equation does not have any real solutions in the region $-R \leq r \leq R$. Are you sure you have not made any errors in writing down the equation?

 

[haruspex]

The very first equation looks a bit unlikely to me. Are you sure it is correct? If not certain, please post the original problem as given to you.

 

[gruba]

The very first equation looks a bit unlikely to me. Are you sure it is correct? If not certain, please post the original problem as given to you.

Here is the original problem (the topic should be moved to Calculus Homework):

Circular tunnel of radius r is punctured through the center of a homogeneous sphere of radius R, (r<R).
1. What percentage of a sphere is lost?
2. What should be the value of r such that the sphere maintains exactly the half of it's volume?

This is a variation of the Napkin ring problem (en.wikipedia.org/wiki/Napkin_ring_problem).

1. Height of the spherical cap at the center is R-\sqrt{R^2-r^2}.
Volume of the one spherical cap that is removed is V_{cap}=\int_{R-\sqrt{R^2-r^2}}^R {(R^2-x^2)}\mathrm dx=\frac{\pi}{3}(r^2 - R^2)(\sqrt{R^2 - r^2} - 3R)
Remaining volume is V_{remaining}=\frac{4 \pi R^3}{3} - 2V_{cap} - 2 \pi r^2 \sqrt{R^2 - r^2}
Removed volume is V_{removed}=2(V_{cap}+\pi r^2\sqrt{R^2-r^2})=2\left(\frac{\pi}{3}(r^2 - R^2)(\sqrt{R^2 - r^2} - 3R)+\pi r^2\sqrt{R^2-r^2}\right)

Percentage of the sphere that is lost is x=\frac{V_{removed}\cdot 100}{V_{total}},V_{total}=\frac{4\pi R^3}{3}\Rightarrow x=\frac{150\left(\frac{\pi}{3}(r^2 - R^2)(\sqrt{R^2 - r^2} - 3R)+\pi r^2\sqrt{R^2-r^2}\right)}{\pi R^3}

2. To find the value of r such that the sphere maintains exactly the half of it's volume we solve the following equation for r:

V_{remaining}=\frac{1}{2}V_{total}
V_{remaining}=V_{total}-2\left(\frac{\pi}{3}(r^2-R^2)(\sqrt{R^2-r^2}-3R)-2\pi r^2\sqrt{R^2-r^2}\right),V_{total}=\frac{4\pi R^3}{3}
V_{remaining}=\frac{4\pi R^3}{3}-\frac{2r^2\sqrt{R^2-r^2}\pi}{3}+2R\pi r^2+\frac{2R^2\pi\sqrt{R^2-r^2}}{3}-2\pi R^3+4\pi r^2\sqrt{R^2-r^2}
\frac{4\pi R^3}{3}-\frac{2r^2\sqrt{R^2-r^2}\pi}{3}+2R\pi r^2+\frac{2R^2\pi\sqrt{R^2-r^2}}{3}-2\pi R^3+4\pi r^2\sqrt{R^2-r^2}=\frac{2\pi R^3}{3}\Rightarrow
-4\pi R^3+10r^2\sqrt{R^2-r^2}\pi+2R^2\sqrt{R^2-r^2}\pi+6R\pi r^2=0

Note: There are some mistakes in the equation of original post.

After dividing by two:

-2\pi R^3+5r^2\sqrt{R^2-r^2}\pi+R^2\sqrt{R^2-r^2}\pi+3R\pi r^2=0
5r^2\sqrt{R^2-r^2}\pi+R^2\sqrt{R^2-r^2}\pi=2\pi R^3-3R\pi r^2

After squaring:

25\pi^2 r^6+2\pi^2 r^4(8R^2-5)+\pi^2 R^2 r^2(10+11R^2)-3\pi^2 R^6=0

 

[haruspex]

The final equation looks more likely now in that the π all factor out, but every term should have order 6, e.g. r²R⁴.
Those 5 and 10 terms are not possible.
Why not make use of the napkin result? Think about the length of the cylindrical hole, measured along its side.

Edit: I suggest you check each equation by considering the extremes of r, namely 0 and R. E.g. take your expression for the volume of a cap and try r=R. It gives 0 instead of 2πR³/3.

 

[gruba]

I suggest you check each equation by considering the extremes of r, namely 0 and R. E.g. take your expression for the volume of a cap and try r=R. It gives 0 instead of 2πR³/3.

Do you think that the volume of a cap is not correct?

 

[gruba]

Volume of a cap is not correct. It should be
V_{cap}=\pi\cdot\frac{3R^3-3Rr^2-(R^2-r^2)^{3/2}}{3}

This gives that V_{removed}=2\left(\pi\left(\frac{3R^3-3Rr^2-(R^2-r^2)^{3/2}}{3}+r^2\sqrt{R^2-r^2}\right)\right)

1. Now x=\frac{150\left(\pi\left(\frac{3R^3-3Rr^2-(R^2-r^2)^{3/2}}{3}+r^2\sqrt{R^2-r^2}\right)\right)}{\pi R^3}

2. The equation for r is now
-3R^3+3Rr^2+(R^2-r^2)^{3/2}-3r^2\sqrt{R^2-r^2}=0

 

[haruspex]

Do you think that the volume of a cap is not correct?

Yes, I think it is wrong. Check the integration bounds.

 

[haruspex]

Your Vcap is still wrong. Try the test I suggested, what happens when r=R? What should it give there?

 

[gruba]

Your Vcap is still wrong. Try the test I suggested, what happens when r=R? What should it give there?

If r=R, then V_{cap}=\frac{2\pi R^3}{3}.

 

[gruba]

After some corrections, I am getting the equation -4r^6+12R^2r^4-12R^4r^2+3R^6=0.

 

[Ray Vickson]

After some corrections, I am getting the equation -4r^6+12R^2r^4-12R^4r^2+3R^6=0.

This is a cubic in $x = r^2$. Furthermore, if you set $r = Ru$ you get an equation for $u$ that does not involve $R$. In other words, you could have just set $R = 1$ right at the start, since both the questions you want to answer involve percentages of the sphere's volume, rather than actual numerical values.

The whole analysis would be much simpler if, instead of $r$ (the cylinder's radius) you use as your variable $h$ = height of the spherical cap. Of course, these are related by $r^2 = R^2 - (R-h)^2$. In terms of $h$ the volume of the removed portion is
V_{\text{removed}}= 2 \pi (R-h)[R^2 - (R-h)^2]+ \frac{2 \pi}{3} h^2 (3R-h),
using the formula for the spherical cap's volume from http://mathworld.wolfram.com/SphericalCap.html .
Setting $h = Rx$ you get the simple cubic equation
\frac{4}{3} x^3 - 4x^2 + 4x = \frac{2}{3}
or
4 x^3 - 12 x^2 + 12 x = 2

 

[haruspex]

After some corrections, I am getting the equation -4r^6+12R^2r^4-12R^4r^2+3R^6=0.

That looks right. Can you see how to 'complete the cube', so as to get an equation in the form (function of r and R)³=(function of R only)?
I note that Ray's equation is slightly different.

Edit: maybe not really different. I forgot Ray is not using x to represent r/R.