Solution of overdetermined nonlinear pipe network

Solution of overdetermined nonlinear pipe network - Messaggi

#1 Pubblicato: 14/04/2022 07:52:18
NDTM Amarasekera

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: 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,95 KiB) scaricato 993 volta(e).
Overdetermined_Nonlinear_Pipe_Network.sm (629,73 KiB) scaricato 1655 volta(e).
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
4 utenti hanno apprezzato questo post
sergio 14/04/2022 07:58:00, overlord 14/04/2022 09:39:00, Alvaro 14/04/2022 15:11:00, Davide Carpi 12/06/2022 09:03:00
#2 Pubblicato: 14/04/2022 09:40:58
Fridel Selitsky

Fridel Selitsky

520 apprezzamenti in 451 post.

Gruppo: User

Hello,ndtma
Fragments withal_nleqsolve() can be reduced to one line.
alnlegsolve(3).sm (635,89 KiB) scaricato 1599 volta(e).
alglib3.png
1 utenti hanno apprezzato questo post
NDTM Amarasekera 14/04/2022 10:06:00
#3 Pubblicato: 14/04/2022 10:03:58
NDTM Amarasekera

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: 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 Pubblicato: 14/04/2022 20:33:11
Alvaro

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: 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,97 KiB) scaricato 1611 volta(e).

Best regards.
Alvaro.
1 utenti hanno apprezzato questo post
NDTM Amarasekera 15/04/2022 06:34:00
#5 Pubblicato: 14/04/2022 23:09:43
CBG

CBG

78 apprezzamenti in 312 post.

Gruppo: 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 utenti hanno apprezzato questo post
NDTM Amarasekera 15/04/2022 02:01:00
#6 Pubblicato: 15/04/2022 00:33:25
NDTM Amarasekera

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: 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 Pubblicato: 17/04/2022 08:43:24
overlord

overlord

600 apprezzamenti in 1.381 post.

Gruppo: Moderator

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

Regards
#8 Pubblicato: 17/04/2022 10:47:27
NDTM Amarasekera

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: 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 Pubblicato: 26/05/2022 04:38:30
NDTM Amarasekera

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: 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,41 KiB) scaricato 1670 volta(e).
Overdetermined_Pipe_Network_3 LOOPS-As per CHEGUIDE_al_nleqsolve.pdf (581,61 KiB) scaricato 991 volta(e).
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#10 Pubblicato: 29/05/2022 23:06:37
CBG

CBG

78 apprezzamenti in 312 post.

Gruppo: 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 Pubblicato: 30/05/2022 05:19:26
NDTM Amarasekera

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: 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 Pubblicato: 06/06/2022 03:23:10
NDTM Amarasekera

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: 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 Pubblicato: 06/06/2022 23:31:51
CBG

CBG

78 apprezzamenti in 312 post.

Gruppo: 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 utenti hanno apprezzato questo post
NDTM Amarasekera 07/06/2022 22:55:00
#14 Pubblicato: 07/06/2022 00:17:56
NDTM Amarasekera

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: User

Thank you Carlos.
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#15 Pubblicato: 07/06/2022 17:20:35
CBG

CBG

78 apprezzamenti in 312 post.

Gruppo: 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 utenti hanno apprezzato questo post
NDTM Amarasekera 07/06/2022 22:55:00
#16 Pubblicato: 07/06/2022 22:55:01
NDTM Amarasekera

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: 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 Pubblicato: 11/06/2022 00:26:18
NDTM Amarasekera

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: 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 Pubblicato: 11/06/2022 16:19:34
CBG

CBG

78 apprezzamenti in 312 post.

Gruppo: 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 utenti hanno apprezzato questo post
Davide Carpi 12/06/2022 09:04:00
#19 Pubblicato: 11/06/2022 22:27:18
NDTM Amarasekera

NDTM Amarasekera

130 apprezzamenti in 352 post.

Gruppo: User

Thank you very much Carlos.
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
  • Nuovi post
  • Nessun nuovo post