Numerical methods developed through the theory of variational inequalities are used to determine the plug flow regions in two classes of pipe flow problems in Bingham fluids. The first set, including a test problem in a concentric annulus, deals with flows in eccentric annuli in general, and the second one examines the flow in an L-shaped region.