Solution of overdetermined nonlinear pipe network

Solution of overdetermined nonlinear pipe network - Сообщения

#1 Опубликовано: 14.04.2022 07:52:18
NDTM Amarasekera

NDTM Amarasekera

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

Группа: User

At present, my country Sri Lanka is undergoing its worst crisis since its independence in 1948, with severe shortage of essential items like fuel, cooking gas, electricity, medicine, milk and food items. As a result, we are compelled to be at home, while thousands are protesting in the capital city of Colombo and other major cities.

So, during this period I spent some time on SMath, and attempted the attached example in the internet solved by Hardy Cross Method using EXCEL

Essentially it is an overdetermined nonlinear pipe network with 22 variables {pipe flows) and 23 equations (16 Nodes + 7 Loops). The solution using SMath is based on a snippet using al_nleqsolve() posted by Viacheslav N. Mezentsev (uni) in SMath Forum.

Overdetermined_Nonlinear_Pipe_Network.pdf (432 КиБ) скачан 101 раз(а).
Overdetermined_Nonlinear_Pipe_Network.sm (629 КиБ) скачан 104 раз(а).
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
4 пользователям понравился этот пост
sergio 14.04.2022 07:58:00, overlord 14.04.2022 09:39:00, Alvaro Diaz Falconi 14.04.2022 15:11:00, Davide Carpi 12.06.2022 09:03:00
#2 Опубликовано: 14.04.2022 09:40:58
Fridel Selitsky

Fridel Selitsky

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

Группа: User

Hello,ndtma
Fragments withal_nleqsolve() can be reduced to one line.
alnlegsolve(3).sm (635 КиБ) скачан 71 раз(а).
alglib3.png
1 пользователям понравился этот пост
NDTM Amarasekera 14.04.2022 10:06:00
#3 Опубликовано: 14.04.2022 10:03:58
NDTM Amarasekera

NDTM Amarasekera

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

Группа: User

Wrote

Hello,ndtma
Fragments withal_nleqsolve()can be reduced to one line.



Hello Ber7,

Thank you very much for your very valuable suggestion.
It works fine.
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#4 Опубликовано: 14.04.2022 20:33:11
Alvaro Diaz Falconi

Alvaro Diaz Falconi

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

Группа: User

Hi Nugegodage. I hope that the problems in your country can be solved, for the better, and that you and your family are well.

Thank you for this excellent example.

Your worksheet with a little different setup, using tables. No one new calc actually.

alnlegsolve v2.sm (904 КиБ) скачан 78 раз(а).

Best regards.
Alvaro.
1 пользователям понравился этот пост
NDTM Amarasekera 15.04.2022 06:34:00
#5 Опубликовано: 14.04.2022 23:09:43
CBG

CBG

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

Группа: User

Hello Nugegodage

I hope that you and your family are very well in these times of crisis in your country.

I have made your example in the Epanet program, the results are shown in the graph below.

For the network to be valid, Node A must have an input pressure equal to or greater than 129.4 mH2O,
and Node C equal to or greater than 150.28 mH2O.

Internet_Network.png

Best Regards

Carlos
1 пользователям понравился этот пост
NDTM Amarasekera 15.04.2022 02:01:00
#6 Опубликовано: 15.04.2022 00:33:25
NDTM Amarasekera

NDTM Amarasekera

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

Группа: User

Hi Jean, Alvaro, Carlos and others.
Thanks for expressing your sentiments and encouragement.
We are safe and hope the situation will settle as early as possible.

Many thanks to Alvaro and Carlos for your suggestions on the example.
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#7 Опубликовано: 17.04.2022 08:43:24
overlord

overlord

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

Группа: Moderator

Thanks for your sample.
I always like your irrigation examples.
Hope crisis in your country resolves too.

Regards
#8 Опубликовано: 17.04.2022 10:47:27
NDTM Amarasekera

NDTM Amarasekera

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

Группа: User

Hi overlord.
Thank you for the sentiments and comments.
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#9 Опубликовано: 26.05.2022 04:38:30
NDTM Amarasekera

NDTM Amarasekera

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

Группа: User

The subject revisited.
In my original post #1, the K values were simply computed using a given formula,
and head loss in each pipe computed accordingly.

In Method 1 of the attached example, friction factors (f) depend on the Reynold numbers (Re)
of each pipe and computed using a given formula involving Re.

In Method 2, NewtonRapson function of SMath is used for the numerical solution of the
implicit Colebrook-White equation to find f.

Many thanks to Ber7, Alvaro and Carlos for your earlier suggestions on the subject.

f_NR.jpg
Overdetermined_Pipe_Network_3 LOOPS-As per CHEGUIDE_al_nleqsolve.sm (286 КиБ) скачан 76 раз(а).
Overdetermined_Pipe_Network_3 LOOPS-As per CHEGUIDE_al_nleqsolve.pdf (581 КиБ) скачан 47 раз(а).
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#10 Опубликовано: 29.05.2022 23:06:37
CBG

CBG

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

Группа: User

Hello ndtmda

Here are the results of your example calculated
with a very old version of Pipe Flow Expert.

Overdetermined_pipe_network.png


Best Regards

Carlos
#11 Опубликовано: 30.05.2022 05:19:26
NDTM Amarasekera

NDTM Amarasekera

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

Группа: User

Wrote

Hello ndtmda

Here are the results of your example calculated
with a very old version of Pipe Flow Expert.



Thank you Carlos!
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#12 Опубликовано: 06.06.2022 03:23:10
NDTM Amarasekera

NDTM Amarasekera

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

Группа: User

Hello Carlos

Can you please explain how the marked values of gauge pressures are calculated?
Thank you.

bar.g.jpg
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#13 Опубликовано: 06.06.2022 23:31:51
CBG

CBG

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

Группа: User

Hello Ndtma


Wrote

Hello Carlos

Can you please explain how the marked values of gauge pressures are calculated?
Thank you.

bar.g.jpg




This pressure (0.9787 bar g) is caused by the water level in the tank, which is 10 meters.


The pressure of 0.6553 bar g is obtained by subtracting from the inlet pressure
of 0.9787 bar g. the losses in the pipe of 0.3234 bar g (3.304 m hd)

The system is not correctly resolved for the pressures, since nodes 3 and 8 have
pressures below atmospheric.

Best Regards

Carlos
1 пользователям понравился этот пост
NDTM Amarasekera 07.06.2022 22:55:00
#14 Опубликовано: 07.06.2022 00:17:56
NDTM Amarasekera

NDTM Amarasekera

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

Группа: User

Thank you Carlos.
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#15 Опубликовано: 07.06.2022 17:20:35
CBG

CBG

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

Группа: User

Hello Ndtma

Wrote

Thank you Carlos.



You're welcome.



Changing the diameter of tubes 9 and 10 from 0.25 m diameter to 0.5 m, the problem of negative pressures is solved.

Another way to solve it is by increasing the level below the tank from 0m to 20m. and with the liquid level at 10 m

Overdetermined_pipe_network_Solution.png


Best Regards

Carlos
1 пользователям понравился этот пост
NDTM Amarasekera 07.06.2022 22:55:00
#16 Опубликовано: 07.06.2022 22:55:01
NDTM Amarasekera

NDTM Amarasekera

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

Группа: User

Wrote


You're welcome.



Changing the diameter of tubes 9 and 10 from 0.25 m diameter to 0.5 m, the problem of negative pressures is solved.

Another way to solve it is by increasing the level below the tank from 0m to 20m. and with the liquid level at 10 m



Thank you Carlos. Your contributions are very much appreciated.
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#17 Опубликовано: 11.06.2022 00:26:18
NDTM Amarasekera

NDTM Amarasekera

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

Группа: User

Wrote

You're welcome.

Changing the diameter of tubes 9 and 10 from 0.25 m diameter to 0.5 m, the problem of negative pressures is solved.

Another way to solve it is by increasing the level below the tank from 0m to 20m. and with the liquid level at 10 m



Hello Carlos,

Can you please post the Pipe Data Table (with Pipe Number, Inner diameter, ..... Re, Friction Head loss)
for the case Changing the diameter of tubes 9 and 10 from 0.25 m diameter to 0.5 m..
as I get different results for this case.

Hope I am not troubling you.
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#18 Опубликовано: 11.06.2022 16:19:34
CBG

CBG

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

Группа: User

Hello Ndtma

Wrote



Hello Carlos,

Can you please post the Pipe Data Table (with Pipe Number, Inner diameter, ..... Re, Friction Head loss)
for the case Changing the diameter of tubes 9 and 10 from 0.25 m diameter to 0.5 m..
as I get different results for this case.

Hope I am not troubling you.




There is no problem.

Here is the requested table

Overdetermined_pipe_network_Solution_Change_Diam.png


Best Regards

Carlos
1 пользователям понравился этот пост
Davide Carpi 12.06.2022 09:04:00
#19 Опубликовано: 11.06.2022 22:27:18
NDTM Amarasekera

NDTM Amarasekera

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

Группа: User

Thank you very much Carlos.
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
  • Новые сообщения Новые сообщения
  • Нет новых сообщений Нет новых сообщений