Solution of overdetermined nonlinear pipe network - Сообщения
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 раз(а).
Fragments withal_nleqsolve() can be reduced to one line.
alnlegsolve(3).sm (635 КиБ) скачан 71 раз(а).
WroteHello,ndtma
Fragments withal_nleqsolve()can be reduced to one line.
Hello Ber7,
Thank you very much for your very valuable suggestion.
It works fine.
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.
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.
Best Regards
Carlos
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.
I always like your irrigation examples.
Hope crisis in your country resolves too.
Regards
Thank you for the sentiments and comments.
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.
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 раз(а).
WroteHello ndtmda
Here are the results of your example calculated
with a very old version of Pipe Flow Expert.
Thank you Carlos!
WroteHello Carlos
Can you please explain how the marked values of gauge pressures are calculated?
Thank you.
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
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.
WroteYou'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.
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
Best Regards
Carlos
-
Новые сообщения
-
Нет новых сообщений