The single mode relaxation time approximation has been demonstrated to greatly underestimate the lattice thermal conductivity of two-dimensional materials due to the collective effect of phonon normal scattering. Callaway's dual relaxation model represents a good approximation to the otherwise ab initio solution of the phonon Boltzmann equation. In this work we develop a discrete-ordinate-method (DOM) scheme for the numerical solution of the phonon Boltzmann equation under Callaway's model. Heat transport in a graphene ribbon with different geometries is modeled by our scheme, which produces results quite consistent with the available molecular dynamics, Monte Carlo simulations, and experimental measurements. Callaway's lattice thermal conductivity model with empirical boundary scattering rates is examined and shown to overestimate or underestimate the direct DOM solution. The length convergence of the lattice thermal conductivity of a rectangular graphene ribbon is explored and found to depend appreciably on the ribbon width, with a semiquantitative correlation provided between the convergence length and the width. Finally, we predict the existence of a phonon Knudsen minimum in a graphene ribbon only at a low system temperature and isotope concentration so that the average normal scattering rate is two orders of magnitude stronger than the intrinsic resistive one. The present work will promote not only the methodology for the solution of the phonon Boltzmann equation but also the theoretical modeling and experimental detection of hydrodynamic phonon transport in two-dimensional materials.