A numerical treatment of two-dimensional flow in a branching channel