The work is devoted to the asymptotic and numerical analysis of the wave function propagating in two-dimensional quantum waveguides with confining potentials supported on slowly varying tubes. The leading term of the asymptotics of the wave function is determined by an adiabatic approach and the WKB approximation. Unlike other similar studies, in the present work we consider arbitrary bounded potentials and obtain exact solutions for the thresholds, and for the transverse modes in the form of power series of the spectral parameter. Our approach leads to an effective numerical method for the analysis of such quantum waveguides and for the tunnel effect observed in sections of the waveguide that shrink or widen too much. Several examples of interest show the applicability of the method.