Flow of Bingham plastics through straight, long tubes with non-circular cross-sections is studied by means of an analytical method that allows to model a wide spectrum of tube geometries. Shear stress normal to equal velocity lines, velocity field and plug zones are explored, in particular in a tube with equilateral triangular cross-section for small values of the Bingham number Bi, and they are compared with corresponding numerical solutions. We show that a circular plug is present at the center of the triangular tube cross-section, consistent with numerical simulations as well as with previous results in the literature, if calculations up to and including first order in the shape perturbation parameter ϵ are taken into account. However with the inclusion of the second order terms in the algorithm this structure is no longer present and no plug zone is predicted for the same pressure drop. We find that in that case normal shear stresses are always greater than the yield stress of the fluid. As a result, the central region becomes a pseudo-plug since it presents small but non-vanishing relative deformations and does not move as a rigid core. The energy dissipation function for the Bingham fluid flow is written in terms of natural coordinates. Its distribution depends only on the normal shear stress at any point with Bingham number as a parameter.