A numerical method for the solution of mass and energy conservation equations in a network of control volumes and links