Error by ODEs solution

#1 Pubblicato: 14/12/2023 01:35:18

67 apprezzamenti in 713 post.

Gruppo: User

Here OK Rocket-1.sm (16,79 KiB) scaricato 1426 volta(e).
Here not OK Rocket-2.sm (688,11 KiB) scaricato 1457 volta(e).
About the problem Study-24-How-calculate-Rocket.doc (1,95 MiB) scaricato 845 volta(e).
from the book
https://www.routledge.com/25-Problems-for-STEM-Education/Ochkov/p/book/9780367345259

#2 Pubblicato: 14/12/2023 09:58:26

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Hi Valery. Just use cases() instead if() clauses.

Rocket-2.sm (688,31 KiB) scaricato 1490 volta(e).

But Newton's second law is not F=ma which is an expression valid only for systems where m is constant. The general expression given by Newton is F=d(mv)/dt or F=dp/dt with p=mv the linear momentum.

Best regards.
Alvaro.

1 utenti hanno apprezzato questo post

Valery Ochkov 14/12/2023 14:13:00

#3 Pubblicato: 14/12/2023 14:30:08

Valery Ochkov

67 apprezzamenti in 713 post.

Gruppo: User

Wrote
Hi Valery. Just use cases() instead if() clauses.

But Newton's second law is not F=ma which is an expression valid only for systems where m is constant. The general expression given by Newton is F=d(mv)/dt or F=dp/dt with p=mv the linear momentum.

Best regards.
Alvaro.

Thanks, Alvaro!
One questions.
How can we solve this problem with ODE as F=dp/dt, not F=m dv/dt?
Will we get others (numerical) answer?

Rocket-2-otherwise.sm (690,02 KiB) scaricato 1476 volta(e).

#4 Pubblicato: 14/12/2023 18:32:10

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

The speed of the rocket cannot be independent of the speed at which the propellant exits. Entering that information gives me this solution

Rocket-1.sm (20,07 KiB) scaricato 1432 volta(e).

Another way of looking at it, using the equation F = d(mv)/dt, we must take into account that in a rocket the reaction force is not applied to the same body as the action force, as the 3rd Law says, and therefore the temporal derivative of the product mv must be done by changing the reference frame, which is why it is F = m dv/dt + u dm/dt, where v is the speed of the rocket with respect to the inertial reference system attached to the earth, but u is the speed of the propellant with respect to the non-inertial reference system attached to the rocket

Best regards.
Alvaro.

#5 Pubblicato: 14/12/2023 18:35:31

Valery Ochkov

67 apprezzamenti in 713 post.

Gruppo: User

New question

Rocket-1.sm (734,95 KiB) scaricato 1477 volta(e).

#6 Pubblicato: 15/12/2023 01:16:32

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Wrote
New question...

Hi Valery. You can't take the temporal derivative directly. Please see my previous answer explaining why I think that.

Best regards.
Alvaro.

#7 Pubblicato: 15/12/2023 02:39:15

Valery Ochkov

67 apprezzamenti in 713 post.

Gruppo: User

Wrote

Wrote
New question...

Hi Valery. You can't take the temporal derivative directly. Please see my previous answer explaining why I think that.

Best regards.
Alvaro.

#8 Pubblicato: 15/12/2023 03:25:33

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Hi Valery. You're right encircling u. Its correct value is u = F/μ. Thus, the force over the rocket is F = - mg - u*dm/dt

Rocket-1.sm (23,73 KiB) scaricato 1469 volta(e).

Best regards.
Alvaro.

#9 Pubblicato: 15/12/2023 04:04:15

Valery Ochkov

67 apprezzamenti in 713 post.

Gruppo: User

Thanks, Alvaro.
And what about the 2-stage rocket?
Our first try
Rocket-2.sm (690,02 KiB) scaricato 1522 volta(e).

#10 Pubblicato: 15/12/2023 04:30:05

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Wrote
...
And what about the 2-stage rocket?
Our first try

Hi. Your original solution for rocket-1 is ok, and that implies that the movie's solution is also ok.

I was only pointing out the way to deduce it. Writing the relationship for several stages would imply asking for derivatives for functions with cases(), which cannot be done, because although functions can be integrated with cases() or using inside the Mathcad block, they definitely cannot be derived.

So I guess the best solution is the one you gave at the beginning, simply replacing the if() with cases(). The deduction based on physical principles has already been done for one stage. Now it remains to apply the simplified result for several, which is exactly what the rocket-2 case does.

Best regards.
Alvaro.

#11 Pubblicato: 15/12/2023 05:00:55

Valery Ochkov

67 apprezzamenti in 713 post.

Gruppo: User

I think it is easier to solve two equations - for two stages, and then for one (last) stage of the rocket.

#12 Pubblicato: 15/12/2023 06:58:19

Martin Kraska

1.361 apprezzamenti in 2.275 post.

Gruppo: Moderator

Analytical reference solution for h(t) (from a mechanics textbook) and with a single second order DE. in() would be even better without braces (operator form)

Rakete.sm (21,1 KiB) scaricato 1448 volta(e).

Rakete2.sm (21,45 KiB) scaricato 1488 volta(e).

Technische Mechanik mit SMath Studio: https://link.springer.com/book/10.1007/978-3-658-50592-9

#13 Pubblicato: 16/12/2023 15:47:47

Valery Ochkov

67 apprezzamenti in 713 post.

Gruppo: User

I would like to publish an article about Rockets in the Journal "Physic in school".
Martin, Rasonar!
Would you like to be coauthors?
See please attach
Polet na Lunu.docx (585,34 KiB) scaricato 787 volta(e).
Rakete-1.sm (16,31 KiB) scaricato 1447 volta(e).
Rakete-2.sm (17,18 KiB) scaricato 1452 volta(e).

1 utenti hanno apprezzato questo post

Alvaro 17/12/2023 02:04:00

#14 Pubblicato: 17/12/2023 02:10:50

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Wrote
I would like to publish an article about Rockets in the Journal "Physic in school".
Martin, Rasonar!
Would you like to be coauthors?
See please attach
Polet na Lunu.docx (585,34 KiB) scaricato 787 volta(e).
Rakete-1.sm (16,31 KiB) scaricato 1447 volta(e).
Rakete-2.sm (17,18 KiB) scaricato 1452 volta(e).

Hi Valery. Sure. It will be a great experience.

Best regards.
Alvaro.

#15 Pubblicato: 17/12/2023 05:05:30

Jean Giraud

983 apprezzamenti in 6.866 post.

Gruppo: User

Wrote
I would like to publish an article about Rockets in the Journal "Physic in school".

The book about Holistic is the rocket model ... works fine in Smath 30405.
Once data will be collected clean, it may be interesting to collect as a function.
A good candidate may be Thiele.
The rocket system does not follow a single DE law, Like ISA.
Thiele fits ISA form - 1km ... 100 km
Jean.

#16 Pubblicato: 17/12/2023 05:49:57

Martin Kraska

1.361 apprezzamenti in 2.275 post.

Gruppo: Moderator

Wrote
I would like to publish an article about Rockets in the Journal "Physic in school".
Martin, Rasonar!
Would you like to be coauthors?
See please attach
Polet na Lunu.docx (510,06 KiB) scaricato 831 volta(e).
Rakete-1.sm (16,31 KiB) scaricato 1447 volta(e).
Rakete-2.sm (17,18 KiB) scaricato 1452 volta(e).

Of course. If I remember correctly, some time there were examples for three-body problems earth-moon-spacecraft. Can't find it just now but that could make another co-author.

Edit: https://en.smath.com/forum/yaf_postsm76986_Mathcad-Toolbox.aspx#post76986

Technische Mechanik mit SMath Studio: https://link.springer.com/book/10.1007/978-3-658-50592-9

#17 Pubblicato: 17/12/2023 15:12:04

Valery Ochkov

67 apprezzamenti in 713 post.

Gruppo: User

> Of course. If I remember correctly, some time there were examples for three-body problems

See see please solutions from Slava

three bodies problem

#18 Pubblicato: 17/12/2023 21:42:27

Jean Giraud

983 apprezzamenti in 6.866 post.

Gruppo: User

Wrote
A good candidate may be Thiele.

I tried Cholesky ... Quite good.
For this 6 points Rocket from the book, Holistic works fine.
Unfortunately, we don't have an exportable formula.
Re-Conciliate the project ... Thiele
Jean

#19 Pubblicato: 17/12/2023 22:33:09

Valery Ochkov

67 apprezzamenti in 713 post.

Gruppo: User