  We review the pertinent features of the phase diagram of the zero-field
Blume-Capel model, focusing on the aspects of transition order, finite-size
scaling and universality. In particular, we employ a range of Monte Carlo
simulation methods to study the 2D spin-1 Blume-Capel model on the square
lattice to investigate the behavior in the vicinity of the first-order and
second-order regimes of the ferromagnet-paramagnet phase boundary,
respectively. To achieve high-precision results, we utilize a combination of
(i) a parallel version of the multicanonical algorithm and (ii) a hybrid
updating scheme combining Metropolis and generalized Wolff cluster moves. These
techniques are combined to study for the first time the correlation length of
the model, using its scaling in the regime of second-order transitions to
illustrate universality through the observed identity of the limiting value of
$\xi/L$ with the exactly known result for the Ising universality class.
