Solution of overdetermined nonlinear pipe network

Solution of overdetermined nonlinear pipe network - Messages

#1 Posted: 4/14/2022 7:52:18 AM
NDTM Amarasekera

NDTM Amarasekera

130 likes in 352 posts.

Group: 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 KiB) downloaded 125 time(s).
Overdetermined_Nonlinear_Pipe_Network.sm (629 KiB) downloaded 134 time(s).
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
4 users liked this post
sergio 4/14/2022 7:58:00 AM, overlord 4/14/2022 9:39:00 AM, Alvaro Diaz Falconi 4/14/2022 3:11:00 PM, Davide Carpi 6/12/2022 9:03:00 AM
#2 Posted: 4/14/2022 9:40:58 AM
Fridel Selitsky

Fridel Selitsky

520 likes in 451 posts.

Group: User

Hello,ndtma
Fragments withal_nleqsolve() can be reduced to one line.
alnlegsolve(3).sm (635 KiB) downloaded 102 time(s).
alglib3.png
1 users liked this post
NDTM Amarasekera 4/14/2022 10:06:00 AM
#3 Posted: 4/14/2022 10:03:58 AM
NDTM Amarasekera

NDTM Amarasekera

130 likes in 352 posts.

Group: 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 Posted: 4/14/2022 8:33:11 PM
Alvaro Diaz Falconi

Alvaro Diaz Falconi

992 likes in 1675 posts.

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

Best regards.
Alvaro.
1 users liked this post
NDTM Amarasekera 4/15/2022 6:34:00 AM
#5 Posted: 4/14/2022 11:09:43 PM
CBG

CBG

77 likes in 312 posts.

Group: 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 users liked this post
NDTM Amarasekera 4/15/2022 2:01:00 AM
#6 Posted: 4/15/2022 12:33:25 AM
NDTM Amarasekera

NDTM Amarasekera

130 likes in 352 posts.

Group: 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 Posted: 4/17/2022 8:43:24 AM
overlord

overlord

554 likes in 1333 posts.

Group: Moderator

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

Regards
#8 Posted: 4/17/2022 10:47:27 AM
NDTM Amarasekera

NDTM Amarasekera

130 likes in 352 posts.

Group: 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 Posted: 5/26/2022 4:38:30 AM
NDTM Amarasekera

NDTM Amarasekera

130 likes in 352 posts.

Group: 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 KiB) downloaded 109 time(s).
Overdetermined_Pipe_Network_3 LOOPS-As per CHEGUIDE_al_nleqsolve.pdf (581 KiB) downloaded 78 time(s).
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#10 Posted: 5/29/2022 11:06:37 PM
CBG

CBG

77 likes in 312 posts.

Group: 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 Posted: 5/30/2022 5:19:26 AM
NDTM Amarasekera

NDTM Amarasekera

130 likes in 352 posts.

Group: 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 Posted: 6/6/2022 3:23:10 AM
NDTM Amarasekera

NDTM Amarasekera

130 likes in 352 posts.

Group: 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 Posted: 6/6/2022 11:31:51 PM
CBG

CBG

77 likes in 312 posts.

Group: 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 users liked this post
NDTM Amarasekera 6/7/2022 10:55:00 PM
#14 Posted: 6/7/2022 12:17:56 AM
NDTM Amarasekera

NDTM Amarasekera

130 likes in 352 posts.

Group: User

Thank you Carlos.
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
#15 Posted: 6/7/2022 5:20:35 PM
CBG

CBG

77 likes in 312 posts.

Group: 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 users liked this post
NDTM Amarasekera 6/7/2022 10:55:00 PM
#16 Posted: 6/7/2022 10:55:01 PM
NDTM Amarasekera

NDTM Amarasekera

130 likes in 352 posts.

Group: 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 Posted: 6/11/2022 12:26:18 AM
NDTM Amarasekera

NDTM Amarasekera

130 likes in 352 posts.

Group: 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 Posted: 6/11/2022 4:19:34 PM
CBG

CBG

77 likes in 312 posts.

Group: 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 users liked this post
Davide Carpi 6/12/2022 9:04:00 AM
#19 Posted: 6/11/2022 10:27:18 PM
NDTM Amarasekera

NDTM Amarasekera

130 likes in 352 posts.

Group: User

Thank you very much Carlos.
Look within!... The secret is inside you. Best Regards Eng. NDTM Amarasekera - Sri Lanka
  • New Posts New Posts
  • No New Posts No New Posts