Sunday, October 15, 2017

Calculating gravity when taking into account the change of gravitational force


This is a problem that has bothered me for a couple of weeks now, and I can't seem to wrap my head around it and understand it.


Let's say we have a planet with a mass of m. We also have an object of relatively small mass (so small that its gravitational field would not affect the planet), and we know that at time 0 it is at position h.


If we know the acceleration due to gravity (g), we can calculate the position y at time t:


$y = h - { g \; t^2 \over 2 } $


And we know that the acceleration due to gravity could be calculated:


$g = G { m \over y^2 }$


But since the acceleration depends on the position of the object relative to the planet, and the position of the object is dependent on the acceleration, this would obviously mean that the first formula won't work when we are working with displacement so large that it would change the acceleration significantly.


So, in essence, my question is, what formula can we use to calculate the position of an object, taking into consideration the change in acceleration due to gravity?




Answer



You are correct that the "normal" formula $y = h - gt^2/2$ doesn't work when the gravitational acceleration changes, so you need a different formula. The mathematical expressions are a little ugly, though. Steven laid the groundwork for this, but I'm going to point you to an earlier answer of mine where I did the calculation. The result comes out to be


$$t_f - t_i = \frac{1}{\sqrt{2G(m_1 + m_2)}}\biggl(\sqrt{r_i r_f(r_i - r_f)} + r_i^{3/2}\cos^{-1}\sqrt{\frac{r_f}{r_i}}\biggr)$$


$r_i$ and $t_i$ are the initial position (height) and time, respectively, and $r_f$ and $t_f$ are the corresponding final values. This equation is a little "backwards" in the sense that instead of expressing position as a function of time, it expresses time as a function of position. You can invert it to express position as a function of time, but you won't find a single nice function for it. You would have to do the inversion numerically, by plugging it into a computer, or by computing a power series or something like that.


No comments:

Post a Comment

classical mechanics - Moment of a force about a given axis (Torque) - Scalar or vectorial?

I am studying Statics and saw that: The moment of a force about a given axis (or Torque) is defined by the equation: $M_X = (\vec r \times \...