How Can I Solve the Rocket Motion Equation Numerically Using Maple?

AI Thread Summary
The discussion focuses on solving the rocket motion equation numerically using Maple, particularly addressing the equation dvr(t)/dt = - (dm/dt) / mr(t) * ve(t). The user is struggling with the formulation of ordinary differential equations (ODEs) and seeks assistance in deriving expressions for acceleration, velocity, and altitude over time. Concerns are raised about the units on both sides of the equation, indicating a potential error. Suggestions include analyzing the rocket's motion from its own reference frame to clarify the relationship between mass ejection and velocity change. The goal is to achieve a correct and solvable equation for numerical analysis.
Dulle
Messages
2
Reaction score
0
Hi there
I'm in desperate need for your help!

I need to solve the rocket motion equation numerically in Maple or some other program. I'm having trouble writing the ODs
in the program.
I have derived the rocket motion equation in a single dimension:

<br /> <br /> \frac {d\vec v_r(t)}{dt} = - \, \frac {\dot m_e(t)} {m_r(t)} \, \vec v_e(t)<br /> <br />

I assume that the change of mass is constant and the exhaust velocity is constant. Let's say that the rocket has an initial mass of 1 kg, an empty mass of 0.1 kg and the exhaust velocity 5 m/s.

So, i want to solve the equation numerically in Maple or some other program (Maple is preferred) and find the expression for acceleration as a function of time, velocity as a function of time and the altitude as a function of time.
 
Physics news on Phys.org
That equation doesn't look right to me. On the left you have dv/dt, which has units of acceleration, but the right side has units of velocity.

I would approach the problem by considering the rocket in its own reference frame. It ejects a mass dm, with velocity Ve. You can find how fast the rocket moves backwards because of this dm. This speed is of course the "dv", differential change in velocity, in any reference frame.

Once you get an equation, it should be easy to solve analytically.
 
I multiplied the values first without the error limit. Got 19.38. rounded it off to 2 significant figures since the given data has 2 significant figures. So = 19. For error I used the above formula. It comes out about 1.48. Now my question is. Should I write the answer as 19±1.5 (rounding 1.48 to 2 significant figures) OR should I write it as 19±1. So in short, should the error have same number of significant figures as the mean value or should it have the same number of decimal places as...
Thread 'Collision of a bullet on a rod-string system: query'
In this question, I have a question. I am NOT trying to solve it, but it is just a conceptual question. Consider the point on the rod, which connects the string and the rod. My question: just before and after the collision, is ANGULAR momentum CONSERVED about this point? Lets call the point which connects the string and rod as P. Why am I asking this? : it is clear from the scenario that the point of concern, which connects the string and the rod, moves in a circular path due to the string...
Thread 'A cylinder connected to a hanging mass'
Let's declare that for the cylinder, mass = M = 10 kg Radius = R = 4 m For the wall and the floor, Friction coeff = ##\mu## = 0.5 For the hanging mass, mass = m = 11 kg First, we divide the force according to their respective plane (x and y thing, correct me if I'm wrong) and according to which, cylinder or the hanging mass, they're working on. Force on the hanging mass $$mg - T = ma$$ Force(Cylinder) on y $$N_f + f_w - Mg = 0$$ Force(Cylinder) on x $$T + f_f - N_w = Ma$$ There's also...
Back
Top