Newtonian Gravity- Body 'falling' into the sun integration problems

[Romeo]

The problem is this:

Given sun of mass M and a body of mass m (M>>m) a distance r from the sun, find the time for the body to 'fall' into the sun (initially ignoring the radius of the sun).


Our first equation is therefore $$\ddot{r} = \frac {GM}{r^2}$$.

I am able to integrate this, giving:
$$\dot{r} = - {\sqrt{2GM}}{\sqrt{1/r - 1/R}}$$,

where R is the inital distance of the body from the sun. However, I am unable to integrate this again. I have shoved it into wolfram's integrator for an indicator of what to aim for, but cannot come close.

Any thoughts would be greatly appreciated.

Regards

Romeo

 

[whozum]

With G and M constant:

$$\int {\frac{GM}{r^2}dr = GM\int{\frac{1}{r^2}}dr$$

$$\frac{1}{r^2} = r^{-2}$$

Use the power rule for integration to find this integral. There are no square roots involved.

 

[Romeo]

Appreciated Woozum, but you mis-read my post. Maybe i should clarify further:

$$\dot{r} = \frac {dr}{dt},$$

and that i solved the first integral of this second order differential equation, but cannot, in particular, solve this:

$$\int \frac {1}{\sqrt{1/r - 1/R}} dr = -\sqrt{2GM}\int dt,$$

reminding that 'R' is the initial distance of the body from the sun and that 'r' is variable distance that you integrate to.

So, any ideas anyone?

Regards

Romeo

 

[BobG]

The problem is this:

Given sun of mass M and a body of mass m (M>>m) a distance r from the sun, find the time for the body to 'fall' into the sun (initially ignoring the radius of the sun).


Our first equation is therefore $$\ddot{r} = \frac {GM}{r^2}$$.

I am able to integrate this, giving:
$$\dot{r} = - {\sqrt{2GM}}{\sqrt{1/r - 1/R}}$$,

where R is the inital distance of the body from the sun. However, I am unable to integrate this again. I have shoved it into wolfram's integrator for an indicator of what to aim for, but cannot come close.

Any thoughts would be greatly appreciated.

Regards

Romeo

whozum's right, but, more importantly, how is that going to get your time?

Your first equation is your acceleration. To get your velocity, your integration should be:

$$\dot{r}=\int{a}\, dt$$

rdot is your velocity, which you can integrate wrt time to get your position.

That get's you an initial value that you can solve for time (just in case you didn't know the equation, already).

 

[Romeo]

whozum's right, but, more importantly, how is that going to get your time?

Your first equation is your acceleration. To get your velocity, your integration should be:

$$\dot{r}=\int{a}\, dt$$

rdot is your velocity, which you can integrate wrt time to get your position.

That get's you an initial value that you can solve for time (just in case you didn't know the equation, already).

BobG, I understand what you're intending to say, but haven't quite got what I'm trying to say:

I have the second order differential equation $$\ddot{r} = -GM /r^2$$ to solve, aiming to find an expression for 'r', which will then simply be be used to find an expression for the time it takes for a body to 'fall' into the sun from a general distance.

I neglected to show the work involved to solve the latter, but I am pretty certain that the expression for $$\dot{r}$$ is correct.

With that in mind, any ideas folks?

Regards

Romeo

 

[whozum]

Honestly don't know what your trying to do, but

How does [itex]\dot{r} = - {\sqrt{2GM}}{\sqrt{1/r - 1/R}} [/tex] stem from [itex]\ddot{r} = \frac {GM}{r^2} [/tex]<br /> <br /> For constant acceleration, <br /> $$x(t) = x_0 + v_0t + \frac{gt^2}{2} \mbox{ where } g = \frac{GM}{r^2} \mbox{ so then }$$<br /> <br /> $$x(t) = \frac{GMt^2}{2r^2}$$ with $$x_0 = v_0 t = 0$$<br /> <br /> and $$t = \sqrt{\frac{2xr^2}{GM}}$$[/itex][/itex]

 

[dextercioby]

Dog Maddit ! Do not double post !You got the same problem solved by Arildno in the Diff.Eq.forum.

Daniel.

 

[Romeo]

Honestly don't know what your trying to do, but

How does [itex]\dot{r} = - {\sqrt{2GM}}{\sqrt{1/r - 1/R}} [/tex] stem from [itex]\ddot{r} = \frac {GM}{r^2} [/tex]<br /> <br /> For constant acceleration, <br /> $$x(t) = x_0 + v_0t + \frac{gt^2}{2} \mbox{ where } g = \frac{GM}{r^2} \mbox{ so then }$$<br /> <br /> $$x(t) = \frac{GMt^2}{2r^2}$$ with $$x_0 = v_0 t = 0$$<br /> <br /> and $$t = \sqrt{\frac{2xr^2}{GM}}$$[/itex][/itex]

[itex]$<br /> Thank you anyway Whozum. I think maybe I should have stated that there will of course not be constant acceleration, hence your approach is invalid. However, there is now a solution in Diff Equ/ns, so check out what I was trying to do! <br /> <br /> Regards<br /> <br /> Romeo$[/itex]

 

[arildno]

A word of thanks could work wonders, you know.
Not a wholesale devotional of course, but an appreciative nod might be in order..

 

[Romeo]

Thank you all, apologies for the double post, all is accounted for in the Diff Equ forum. And again, thank you Arildno for the solution.

Regards

Romeo

 

[thinkRodrigue]

The problem can also be solved without performing any integration by applying Kepler's laws of planetary motion.

From Kepler's first law the orbit of any mass around the sun is an ellipse with the sun at one of the two foci. Let's call the semimajor axis of this ellipse 'a' and the semiminor axis of the elipse 'b'.

An object with zero initial velocity falling into the sun is the limiting case of this orbit as the minor axis of the ellipse (b) goes to 0. The distance f from the centre of the ellipse to the foci is f = (a^2 - b^2) ^(1/2). Obviously as b -> 0, f ->a. So the initial distance R from the object to the sun is R = 2a in this limit.

Kepler's third law gives the orbital period T in terms of the semimajor axis length.

T^2 = a^3 (4*Pi^2)/MG

Bearing in mind that on object falling into the Sun only completes half of its 'orbit' and using the relationship between a and R above you can solve the problem very simply.