Free fall trajectory in non constant gravitational field

Click For Summary
SUMMARY

The discussion focuses on deriving the free fall trajectory of an object in a non-constant gravitational field, specifically above Earth where gravitational acceleration (g) varies with distance. Participants emphasize the need to apply Newton's law of gravitation, represented as F = -GmM/r², instead of assuming g is constant. The integration process involves using the method of quadrature to relate velocity and distance, ultimately leading to an expression for time as a function of distance, which results in an elliptic integral that cannot be expressed in elementary functions.

PREREQUISITES
  • Understanding of Newton's law of gravitation
  • Familiarity with calculus, particularly integration techniques
  • Knowledge of differential equations and their applications
  • Basic concepts of elliptic integrals
NEXT STEPS
  • Explore the derivation of elliptic integrals and their applications in physics
  • Study advanced integration techniques, particularly quadrature methods
  • Learn about gravitational fields and their variations with altitude
  • Investigate numerical methods for solving non-elementary integrals
USEFUL FOR

Physicists, mathematicians, and engineers involved in gravitational studies, trajectory analysis, and advanced calculus applications.

Creator
Messages
566
Reaction score
6
Rusty on the math...
I am working on a problem and need the derivation of a free fall trajectory for an object at a distance above Earth where the change in acceleration is not negligible. How do I integrate the distance/ acceleration formula taking into account the change in g to get the total distance traveled.
The method and final formula...?
I'm looking for total distance after time t, trying to see how it scales as what power of t.
thanks.
 
Physics news on Phys.org
You need to know what the function of g with respect to position is, and then plug it into your integral instead of assuming it is a constant.
 
Let me get this right. The object is tracing a trajectory where 'g' changes according to the distance from the Earth?
 
siddharth23 said:
Let me get this right. The object is tracing a trajectory where 'g' changes according to the distance from the Earth?

Over long enough range, g is not constant as the Earth is not of uniform density. Although I suspect what the OP was interested in was the decreasing effect of gravity with altitude. In which case you would substitute Newton's law of gravitation in for the constant g in your integral g(x)=G M1 M2 / x^2.
 
You possibly can't determine the density at every point. Has to be g(x). Like QuatumPion said, put it in the equation in terms of x and integrate, In this case, it'll be in terms of 'y' as 'x' would be the horizontal distance traveled by the projectile.
 
siddharth23 said:
You possibly can't determine the density at every point. Has to be g(x). Like QuatumPion said, put it in the equation in terms of x and integrate, In this case, it'll be in terms of 'y' as 'x' would be the horizontal distance traveled by the projectile.

Probably, but even for the first case if you have a gravity map and know the trajectory/inclination than you could make a function of g with respect to position, although that would be fairly complicated (depends on orbital parameters, rotation of earth, etc).
 
Creator said:
Rusty on the math...
I am working on a problem and need the derivation of a free fall trajectory for an object at a distance above Earth where the change in acceleration is not negligible. How do I integrate the distance/ acceleration formula taking into account the change in g to get the total distance traveled.
The method and final formula...?
I'm looking for total distance after time t, trying to see how it scales as what power of t.
thanks.
Sidhartha23 and QuantumPion have argued back and forth but we haven't heard back from Creator so I am going to look at the specific case Creator asked about: An object falling from such a great height the we must use the general law for gravity: F= -GmM/r^2, rather than treating F as a constant, still satisfies "Force= mass times acceleration". Here "m" is the mass of the object, M is the mass of the Earth so we have ma= -GmM/r^2 or
a= \frac{dv}{dt}= -\frac{GM}{r^2}

We would like to integrate that but we need to integrate with respect to t and have "r" on the right. To fix that, we can use a method called "quadrature". By the chain rule, dv/dt= (dv/dr)(dr/dt)= v(dv/dr) becaues v, the speed, is dr/dt. Now we have
v\frac{dv}{dr}= -\frac{GM}{r^2}
which we can write, in "differential form" as
vdv= -\frac{GM}{r^2}dr= -GMr^{-2}dr

Integrating both sides,
\frac{1}{2}v^2= GMr^{-1}+ C
assuming that the object falls from a standstill at height R, v= 0, so 0= GMR^{-1}+ C, C= -GMR^{-1}.

From that,
v= \frac{dr}{dt}= \sqrt{2GMr^{-1}- 2GMR^{-1}}= \sqrt{2GM}\sqrt{r^{-1}- R^{-1}}

We can write that as
\frac{dr}{\sqrt{2GM}\sqrt{r^{-1}+ R^{-1}}}= dt

Integrating both sides of that will give t as a function of r which can then solved for r as a function of t. However, the left side cannot be integrated in terms of "elementary functions". It is an "elliptic integral" (the name, unsurprisingly, connected with the fact that planets move in ellipses around the sun). There used to be multi-volumes of tables of elliptic integrals.
 
Last edited by a moderator:
Thanks for all the quick responses...
You are correct Siddarth23 and Quantum Pion...about my question.
Creator
 
HallsofIvy said:
Sidhartha23 and QuantumPion have argued back and forth but we haven't heard back from Creator so I am going to look at the specific case Creator asked about: An object falling from such a great height the we must use the general law for gravity: F= -GmM/r^2, rather than treating F as a constant, still satisfies "Force= mass times acceleration". Here "m" is the mass of the object, M is the mass of the Earth so we have ma= -GmM/r^2 or
a= \frac{dv}{dt}= -\frac{GM}{r^2}

We would like to integrate that but we need to integrate with respect to t and have "r" on the right. To fix that, we can use a method called "quadrature". By the chain rule, dv/dt= (dv/dr)(dr/dt)= v(dv/dr) becaues v, the speed, is dr/dt. Now we have
v\frac{dv}{dr}= -\frac{GM}{r^2}
which we can write, in "differential form" as
vdv= -\frac{GM}{r^2}dr= -GMr^{-2}dr

Integrating both sides,
\frac{1}{2}v^2= GMr^{-1}+ C
assuming that the object falls from a standstill at height R, v= 0, so 0= GMR^{-1}+ C, C= -GMR^{-1}.

From that,
v= \frac{dr}{dt}= \sqrt{2GMr^{-1}- 2GMR^{-1}}= \sqrt{2GM}\sqrt{r^{-1}- R^{-1}

We can write that as
\frac{dr}{\sqrt}2GM}\sqrt{r^{-1}+ R^{-1}= dt

Integrating both sides of that will give t as a function of r which can then solved for r as a function of t. However, the left side cannot be integrated in terms of "elementary functions". It is an "elliptic integral" (the name, unsurprisingly, connected with the fact that planets move in ellipses around the sun). There used to be multi-volumes of tables of elliptic integrals.
Thank you Halls ...that is exactly what I was looking for...however, I only followed you up to this eqn. (4th one)...
\frac{1}{2}v^2= GMr^{-1}+ C

And your next step (5) I don't understand since I thought the v here was final velocity...

Your last two eqns. remained in ITEX form and did not print correctly...could you give the last two again. Maybe its my browser...mot too famailiar with mathjax.
Thanks;
Creator
 
Last edited:

Similar threads

  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 13 ·
Replies
13
Views
3K
  • · Replies 25 ·
Replies
25
Views
3K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 11 ·
Replies
11
Views
2K
  • · Replies 7 ·
Replies
7
Views
10K