Numerical solution for two-dimensional flow in a branching channel using boundary-fitted coordinates