Reservoir BC is working.
Calculate method in the problem class now is able to coordinate the pipes and boundar conditions with one flow rate only.
Junction was added but still not tested.
Calculate method in the problem class now is able to coordinate the pipes and boundar conditions with n pipes connected to it
Junction and Reservoir boundary conditions were added and are working.
Pipe regular MOC (no cavitation model) working.
Junction working.
Reservoir working.
Steady state calcilations working. This is based on Wylie's book, which changes the accoustic wave speed in order to have just two nodes per pipe so that calculations get faster