Gravity Gradient

In standard kinematics problems involving gravity, we assume a constant gravitational field strength. Suppose instead that we have a fictitious scenario in which the gravity field strength varies as follows: g = g 0 + k y , g = g_0 + ky, where y y (in meters) is the distance above the ground. Here, the gravity is weakest at ground level and becomes stronger with increasing altitude (whilst always pointing toward the ground).

Consider the time t f t_f (in seconds) that it would take for an object to fall from an initial resting height y 0 y_0 to the ground. To 2 decimal places, what is the limiting value of t f t_f as y 0 y_0 approaches infinity?

Note: The quantity g g has units of m/s 2 \text{m/s}^2 , g 0 = 10 m/s 2 g_0 = 10 \text{ m/s}^2 , and k = 1 / s 2 k= 1/\text{s}^{2} .


The answer is 1.57.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

4 solutions

Mark Hennings
Feb 12, 2017

We have the differential equation v d v d y = a = g 0 y v \frac{dv}{dy} \; = \; a \; = \; -g_0 - y so that v 2 + ( g 0 + y ) 2 = ( g 0 + y 0 ) 2 v^2 + (g_0+y)^2 \; = \; (g_0 + y_0)^2 and hence d y d t = ( g 0 + y 0 ) 2 ( g 0 + y ) 2 \frac{dy}{dt} \; = \; -\sqrt{(g_0 + y_0)^2 - (g_0 + y)^2} Thus the time to fall to ground from height y 0 y_0 is t f = 0 y 0 d y ( g 0 + y 0 ) 2 ( g 0 + y ) 2 = 1 2 π sin 1 ( g 0 g 0 + y 0 ) t_f \; = \; \int_0^{y_0} \frac{dy}{\sqrt{(g_0 + y_0)^2 - (g_0 + y)^2}} \; = \; \tfrac12\pi - \sin^{-1}\big(\tfrac{g_0}{g_0 + y_0}\big) which tends to 1 2 π \boxed{\tfrac12\pi} as y 0 y_0 \to \infty .

This is a very rigorous and accurate method of calculating the time of fall. However, the time of fall for a very large initial height can be calculated directly by considering it as a simple harmonic oscillation. This is shown in the Brain's and Chirag's solution below.

Rohit Gupta - 4 years, 3 months ago

Log in to reply

Fine, if you are only interested in answering the question and obtaining the limiting value. My solution gives so much more.

Mark Hennings - 4 years, 3 months ago

Log in to reply

Yes, I agree with you. Your solution gives an accurate value for any height. It is just that if you want to calculate only the limiting value then the alternative approach can be taken.

Rohit Gupta - 4 years, 3 months ago

If you recognize the given differential equation as SHM, you can write the equation for y y given any time and initial height. I agree that you have explicitly shown more information, but you can derive everything from SHM as well, so please don't try to downplay the method (you seem to command a fair bit of respect here)

Brian Moehring - 4 years, 3 months ago

Log in to reply

@Brian Moehring By all means use SHM to obtain the explicit formula instead of solving the differential equation if you prefer. My point was merely to question why we should be content only to solve the limiting case of the problem when it is not hard to get more.

Mark Hennings - 4 years, 3 months ago

Log in to reply

@Mark Hennings All right. Personally, I just see a tradeoff in utility. Your solution, while more comprehensive, must be read rather than heard. My solution has less detail, but could reasonably be explained to someone orally. I don't really see one method as better than the other, but rather complementary.

Brian Moehring - 4 years, 3 months ago

I did it the diff eq way as well. I only asked for the limiting value because I was amused by the value it took. Might be fun to come up with a harder form that forces a diff eq approach.

Steven Chase - 4 years, 3 months ago

Although I like the spring analogy too

Steven Chase - 4 years, 3 months ago

@Mark Hennings How did you get v 2 + ( g 0 + y ) 2 = ( g 0 + y 0 ) 2 v^2 + (g_0+y)^2 \; = \; (g_0 + y_0)^2 from v d v d y = a = g 0 y v \frac{dv}{dy} \; = \; a \; = \; -g_0 - y

akshat rai - 3 years ago

Log in to reply

Integrate the differential equation with respect to y y and apply the initial condition that v = 0 v = 0 when y = y 0 y = y_0 .

Mark Hennings - 3 years ago

Log in to reply

Thanks a ton sir.

akshat rai - 3 years ago
Brian Moehring
Feb 12, 2017

This describes a simple mass-spring system that obeys Hooke's law with spring constant k = m k=m and equilibrium position y = g 0 y=-g_0 .

As we stretch the spring further, the velocity of the mass will approach infinity near the equilibrium position, so the time it takes to travel from y = 0 y=0 to y = g 0 y=-g_0 will approach zero, and therefore we may instead ask how long it takes for the mass-spring system to reach its equilibrium position.

Since the angular velocity is ω = k / m = 1 \omega = \sqrt{k/m} = 1 , the time it takes to reach the equilibrium position is π / 2 ω = π 2 \frac{\pi/2}{\omega} = \frac{\pi}{2} .

This is a unique approach. Comparing this case with the simple harmonic motion simplifies the calculations a lot. However, you may want to note here, that for small initial heights this method will not work and the answer will depend on the initial height of drop.

Rohit Gupta - 4 years, 3 months ago

Log in to reply

My particular answer will not apply, but I think I make it clear in the whole second paragraph that I am dealing with the limiting case.

On the other hand, the first paragraph applies for ALL initial displacements, even small ones. Therefore, if you really want to compute the time for small heights, note that the equation for the SHM in the problem can written as y + g 0 = ( y 0 + g 0 ) cos t , y+g_0 = (y_0+g_0)\cos t, which we can use to solve for t t when y = 0 y=0 . This will give the time for any y 0 y_0 .

Brian Moehring - 4 years, 3 months ago

Swweeeet!!!!

Arijit ghosh Dastidar - 4 years, 2 months ago
Rohit Gupta
Feb 14, 2017

If we assume the motion to continue even after reaching the surface of the Earth (following the relation of acceleration given in the question), then the particle will go in simple harmonic motion.

Acceleration of the particle can be written as a = g 0 + y -a = g_0 + y .

a = Y a = -Y where, Y = g 0 + y Y = g_0 +y .

Comparing this equation with equation of simple harmonic motion, a = ω 2 y a = - \omega^2 y .

Here, ω \omega is the angular frequency.

ω = 1 rad s 1 \omega = 1 \, \text{rad s}^{-1} .

In an SHM, time taken to move from one extreme to the mean position is T 4 \dfrac{T}{4} .

Here, T T is the time period equals to 2 π ω \dfrac{2 \pi}{\omega} .

Therefore, the time taken to move from the extreme position to the mean position is 2 π 4 \dfrac{2 \pi}{4} .

The mean position is where the acceleration is zero, that will be at y = g 0 y= - g_0 .

If the object is dropped from a great height, then its speed near the Earth surface will be very large and the time to move from y = 0 y=0 to y = g 0 = 10 m y=-g_0 = -10 m will be very small. Therefore, the time to move from the infinity to the surface will be 2 π 4 \dfrac{2 \pi}{4} .

One fourth of the time period it is with angular frequency being 1 rad/s,the surface being the mean position,we can say time taken is T/4 where T is the total period,hence 2π/4=1.57,it is

I really like your thought, but this gives an answer independent of the initial height. So, are you suggesting that the answer won't depend on at which height the object is dropped?

Rohit Gupta - 4 years, 3 months ago

Log in to reply

If you look at the Mark's solution t = 1 2 π sin 1 ( g 0 g 0 + y 0 t= \frac12\pi - \sin^{-1}(\frac{g_0}{g_0 + y_0} This means that for small heights the answer depends on the initial height. But, If I consider your approach, it does not take care of the case of small heights. So, why your approach is not working for small heights?

Rohit Gupta - 4 years, 3 months ago

It certainly does,but since the height tends to infinity the restoring force dominates/overshoots the constant force/acceleration component which is 10 m/s^2,hence even if the constant component were 25 m/s^2 the answer would pretty much remain the same

Chirag Shyamsundar - 4 years, 3 months ago

As I said the restoring force is so dominant,actually the trajectory of the particle is y=cost+10 which is the superposition of a simple harmonic motion and a constant force motion

Chirag Shyamsundar - 4 years, 3 months ago

Log in to reply

Ok, I get you now.

Rohit Gupta - 4 years, 3 months ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...