One–dimensional blood flow models take the general form of nonlinear hyperbolic systems but differ in their formulation. One class of models considers the physically conserved quantities of mass and momentum, while another class describes mass and velocity. Further, the averaging process employed in the model derivation requires the specification of the axial velocity profile; this choice differentiates models within each class. Discrepancies among differing models have yet to be investigated. In this paper, we comment on some theoretical differences among models and systematically compare them for physiologically relevant vessel parameters, network topology, and boundary data. In particular, the effect of the velocity profile is investigated in the cases of both smooth and discontinuous solutions, and a recommendation for a physiological model is provided. The models are discretized by a class of Runge–Kutta discontinuous Galerkin methods.