Error by ODEs solution

#1 Posted: 12/14/2023 1:35:18 AM

55 likes in 633 posts.

Group: User

Here OK Rocket-1.sm (16 KiB) downloaded 26 time(s).
Here not OK Rocket-2.sm (688 KiB) downloaded 27 time(s).
About the problem Study-24-How-calculate-Rocket.doc (1 MiB) downloaded 42 time(s).
from the book
https://www.routledge.com/25-Problems-for-STEM-Education/Ochkov/p/book/9780367345259

#2 Posted: 12/14/2023 9:58:26 AM

Alvaro Diaz Falconi

992 likes in 1674 posts.

Group: User

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

Rocket-2.sm (688 KiB) downloaded 29 time(s).

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 users liked this post

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

#3 Posted: 12/14/2023 2:30:08 PM

Valery Ochkov

55 likes in 633 posts.

Group: 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 KiB) downloaded 27 time(s).

#4 Posted: 12/14/2023 6:32:10 PM

Alvaro Diaz Falconi

992 likes in 1674 posts.

Group: 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 KiB) downloaded 24 time(s).

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 Posted: 12/14/2023 6:35:31 PM

Valery Ochkov

55 likes in 633 posts.

Group: User

New question

Rocket-1.sm (734 KiB) downloaded 18 time(s).

#6 Posted: 12/15/2023 1:16:32 AM

Alvaro Diaz Falconi

992 likes in 1674 posts.

Group: 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 Posted: 12/15/2023 2:39:15 AM

Valery Ochkov

55 likes in 633 posts.

Group: 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 Posted: 12/15/2023 3:25:33 AM

Alvaro Diaz Falconi

992 likes in 1674 posts.

Group: 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 KiB) downloaded 24 time(s).

Best regards.
Alvaro.

#9 Posted: 12/15/2023 4:04:15 AM

Valery Ochkov

55 likes in 633 posts.

Group: User

Thanks, Alvaro.
And what about the 2-stage rocket?
Our first try
Rocket-2.sm (690 KiB) downloaded 27 time(s).

#10 Posted: 12/15/2023 4:30:05 AM

Alvaro Diaz Falconi

992 likes in 1674 posts.

Group: 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 Posted: 12/15/2023 5:00:55 AM

Valery Ochkov

55 likes in 633 posts.

Group: User

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

#12 Posted: 12/15/2023 6:58:19 AM

Martin Kraska

1222 likes in 2150 posts.

Group: 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 KiB) downloaded 30 time(s).

Rakete2.sm (21 KiB) downloaded 34 time(s).

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

#13 Posted: 12/16/2023 3:47:47 PM

Valery Ochkov

55 likes in 633 posts.

Group: 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 KiB) downloaded 26 time(s).
Rakete-1.sm (16 KiB) downloaded 27 time(s).
Rakete-2.sm (17 KiB) downloaded 29 time(s).

1 users liked this post

Alvaro Diaz Falconi 12/17/2023 2:04:00 AM

#14 Posted: 12/17/2023 2:10:50 AM

Alvaro Diaz Falconi

992 likes in 1674 posts.

Group: 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 KiB) downloaded 26 time(s).
Rakete-1.sm (16 KiB) downloaded 27 time(s).
Rakete-2.sm (17 KiB) downloaded 29 time(s).

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

Best regards.
Alvaro.

#15 Posted: 12/17/2023 5:05:30 AM

Jean Giraud

983 likes in 6866 posts.

Group: 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 Posted: 12/17/2023 5:49:57 AM

Martin Kraska

1222 likes in 2150 posts.

Group: 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 KiB) downloaded 20 time(s).
Rakete-1.sm (16 KiB) downloaded 27 time(s).
Rakete-2.sm (17 KiB) downloaded 29 time(s).

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 Posted: 12/17/2023 3:12:04 PM

Valery Ochkov

55 likes in 633 posts.

Group: 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 Posted: 12/17/2023 9:42:27 PM

Jean Giraud

983 likes in 6866 posts.

Group: 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 Posted: 12/17/2023 10:33:09 PM

Valery Ochkov

55 likes in 633 posts.

Group: User