Challenging differential equation for a physics problem( higher math needed)?
-
I thought out one problem in physics. But I have hard trouble setting up the equation. The problem is Lets say there is a elevator with a base area of A In a shaft at an height of H . The mass of the elevator is M. We assume that the shaft beneath the elevator is air tight. At beginning the air pressure in the shaft is po. As you may have thought the elevator starts a free fall. The gas in the shaft does not heat in the process. So po*Vo=p*V applies. What I want is a equation that shows how much has the elevator traveled in any given time. So what I thought was that. The acceleration at any given time is given by M*a=M*g-p*A where p is po*A*H=p*A*(H-x) p=po*H/(H-x) a=g - po*H*A/(M*(H-x)) So I also thought that a=x'' second derivative of x So I get x''=g - po*H*A/(M*(H-x)) Which must be solved to find a function x(t) 1)First question is if this logic is right 2)If not then how should this be done 3) I would also appreciate someone could solve this ( or if he or she gets another equation then solve that). 4) Finally is it possible to lose the assumption that the gas does not heat. I mean would the equation then be solvable. I hope there is someone out there who can solve this
-
Answer:
(I will be neglecting air resistance) There's one force (gravity) pulling the elevator downwards: Fg = M*g And the other one is the air pressure in the shaft: Fa = p*A Fres = M*a(t) = p(t)*A - M*g -> a(t) = p(t)*A/M - g p0*V0 = p(t)*V(t) -> p(t) = p0*V0/V(t) V0 = H*A ... V(t) = (H-b)*A b... what the elevator has travelled so far b = da(t)²/d²t a(t) = p0*V0/((H-b)*A) - g = p0*V0/((H - d²a(t)/dt²)*A) - g, which is about what you've come up with. Personally I think that this equation is really really hard to solve since it's a non-linear second-order differential equation where the differential is a reciprocal value. I've got no idea how to do this by hand...
ANZI at Yahoo! Answers Visit the source
Other answers
1/Your logic is right. Unfortunately, Eva messes it up completely by introducing a second derivative (!!??!) in place of the displacement itself. The second derivative is, as you correctly note,the acceleration which is the lhs of your final equation. 3/ I don't think it can be solved completely with elementary functions. However, you can get a first integral as follows: I rewrite the equation x"=g-k/(H-x) for brevity . Let's then multiply both sides by x' (the velocity) to get x'x" = gx'-k x'/(H-x) so that all three terms are time derivatives of easy to write functions. Integrating you get : 1/2x'^2 = gx +kLog(H-x) + C The integration constant C is easily found by observing that the initial velocity is 0, that is, when x=0, then x'=0. Hence: 0 = kLogH+C so that the equation for the velocity becomes: 1/2x'^2 = gx + kLog(1-x/H) Here you see that in the beginning (x close to 0) the velocity squared grows linearly with x, But later the log becomes strongly negative and the movement stops. The maximum displacement is given by x'=0, i.e. gx+kLog(1-x/H) = 0 which has, of course, the solution x=0 but also another one corresponding to x_max. This can only be solved numerically. You can find a graphical approximation by looking for the intersections of the two curves exp(-gx/k) and 1-x/H Unfortunately, I can't make a plot here. To go further, you could observe that the diffeq obtained is separable by taking the square root and writing it dx/sqrt(gx+kLog(1-x/H)) = dt and integrating and finally inverting. But that's very theoretical. I don't believe any simple primitive of the 1/sqrt(..) function exists and you would still have the problem to invert it. However, you might solve it numerically to get your results. 4/ Given that the whole story (elevator being set free, falling and being stopped by the air pressure) should be short, you might indeed think that the compression is adiabatic (no heat exchange during the process) rather than isothermal. In that case the air will indeed heat up. The pressure-volume relation becomes in this case pv^gamma = cte instead of pv=cte with gamma the ratio of constant pressure over constant volume heat capacities which is about 1.4 for air. Then you can follow the same course. Instead of having to integrate x'/(H-x) leading to a log, you'll have x'/(H-x)^gamma leading to a power law -(1-x/H)^(1-gamma) It might make the square root more manageable for analytical integration but I'm not sure and have no time to check.
jean-de-la-lune
Related Q & A:
- How to find the global equation of a dual curve?Best solution by Mathematics
- Is this problem a knapsack problem?Best solution by Programming Puzzles & Code Golf
- How to implement a physics engine?Best solution by Game Development
- How can I evaluate an equation within a String?Best solution by stackoverflow.com
- Do I have a learning disability in Math?Best solution by understood.org
Just Added Q & A:
- How many active mobile subscribers are there in China?Best solution by Quora
- How to find the right vacation?Best solution by bookit.com
- How To Make Your Own Primer?Best solution by thekrazycouponlady.com
- How do you get the domain & range?Best solution by ChaCha
- How do you open pop up blockers?Best solution by Yahoo! Answers
For every problem there is a solution! Proved by Solucija.
-
Got an issue and looking for advice?
-
Ask Solucija to search every corner of the Web for help.
-
Get workable solutions and helpful tips in a moment.
Just ask Solucija about an issue you face and immediately get a list of ready solutions, answers and tips from other Internet users. We always provide the most suitable and complete answer to your question at the top, along with a few good alternatives below.