Error by ODEs solution

#1 Опубликовано: 14.12.2023 01:35:18

55 сообщений из 634 понравились пользователям.

Группа: User

Here OK Rocket-1.sm (16 КиБ) скачан 135 раз(а).
Here not OK Rocket-2.sm (688 КиБ) скачан 118 раз(а).
About the problem Study-24-How-calculate-Rocket.doc (1 МиБ) скачан 138 раз(а).
from the book
https://www.routledge.com/25-Problems-for-STEM-Education/Ochkov/p/book/9780367345259

#2 Опубликовано: 14.12.2023 09:58:26

Alvaro Diaz Falconi

1 005 сообщений из 1 677 понравились пользователям.

Группа: User

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

Rocket-2.sm (688 КиБ) скачан 127 раз(а).

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 пользователям понравился этот пост

Valery Ochkov 14.12.2023 14:13:00

#3 Опубликовано: 14.12.2023 14:30:08

Valery Ochkov

55 сообщений из 634 понравились пользователям.

Группа: 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 КиБ) скачан 119 раз(а).

#4 Опубликовано: 14.12.2023 18:32:10

Alvaro Diaz Falconi

1 005 сообщений из 1 677 понравились пользователям.

Группа: 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 КиБ) скачан 133 раз(а).

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 Опубликовано: 14.12.2023 18:35:31

Valery Ochkov

55 сообщений из 634 понравились пользователям.

Группа: User

New question

Rocket-1.sm (734 КиБ) скачан 116 раз(а).

#6 Опубликовано: 15.12.2023 01:16:32

Alvaro Diaz Falconi

1 005 сообщений из 1 677 понравились пользователям.

Группа: 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 Опубликовано: 15.12.2023 02:39:15

Valery Ochkov

55 сообщений из 634 понравились пользователям.

Группа: 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 Опубликовано: 15.12.2023 03:25:33

Alvaro Diaz Falconi

1 005 сообщений из 1 677 понравились пользователям.

Группа: 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 КиБ) скачан 131 раз(а).

Best regards.
Alvaro.

#9 Опубликовано: 15.12.2023 04:04:15

Valery Ochkov

55 сообщений из 634 понравились пользователям.

Группа: User

Thanks, Alvaro.
And what about the 2-stage rocket?
Our first try
Rocket-2.sm (690 КиБ) скачан 126 раз(а).

#10 Опубликовано: 15.12.2023 04:30:05

Alvaro Diaz Falconi

1 005 сообщений из 1 677 понравились пользователям.

Группа: 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 Опубликовано: 15.12.2023 05:00:55

Valery Ochkov

55 сообщений из 634 понравились пользователям.

Группа: User

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

#12 Опубликовано: 15.12.2023 06:58:19

Martin Kraska

1 230 сообщений из 2 156 понравились пользователям.

Группа: 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 КиБ) скачан 137 раз(а).

Rakete2.sm (21 КиБ) скачан 143 раз(а).

Martin Kraska Pre-configured portable distribution of SMath Studio: https://en.smath.info/wiki/SMath%20with%20Plugins.ashx

#13 Опубликовано: 16.12.2023 15:47:47

Valery Ochkov

55 сообщений из 634 понравились пользователям.

Группа: 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 КиБ) скачан 118 раз(а).
Rakete-1.sm (16 КиБ) скачан 132 раз(а).
Rakete-2.sm (17 КиБ) скачан 129 раз(а).

1 пользователям понравился этот пост

Alvaro Diaz Falconi 17.12.2023 02:04:00

#14 Опубликовано: 17.12.2023 02:10:50

Alvaro Diaz Falconi

1 005 сообщений из 1 677 понравились пользователям.

Группа: 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 КиБ) скачан 118 раз(а).
Rakete-1.sm (16 КиБ) скачан 132 раз(а).
Rakete-2.sm (17 КиБ) скачан 129 раз(а).

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

Best regards.
Alvaro.

#15 Опубликовано: 17.12.2023 05:05:30

Jean Giraud

983 сообщений из 6 866 понравились пользователям.

Группа: 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 Опубликовано: 17.12.2023 05:49:57

Martin Kraska

1 230 сообщений из 2 156 понравились пользователям.

Группа: 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 КиБ) скачан 120 раз(а).
Rakete-1.sm (16 КиБ) скачан 132 раз(а).
Rakete-2.sm (17 КиБ) скачан 129 раз(а).

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

Martin Kraska Pre-configured portable distribution of SMath Studio: https://en.smath.info/wiki/SMath%20with%20Plugins.ashx

#17 Опубликовано: 17.12.2023 15:12:04

Valery Ochkov

55 сообщений из 634 понравились пользователям.

Группа: 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 Опубликовано: 17.12.2023 21:42:27

Jean Giraud

983 сообщений из 6 866 понравились пользователям.

Группа: 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 Опубликовано: 17.12.2023 22:33:09

Valery Ochkov

55 сообщений из 634 понравились пользователям.

Группа: User