quasi 1-D flow solver