  Two phenomenological models describing an SU(N) quark-gluon plasma are
presented. The first is obtained from high temperature expansions of the free
energy of a massive gluon, while the second is derived by demanding color
neutrality over a certain length scale. Each model has a single free parameter,
exhibits behavior similar to lattice simulations over the range T_d - 5T_d, and
has the correct blackbody behavior for large temperatures. The N = 2
deconfinement transition is second order in both models, while N = 3,4, and 5
are first order. Both models appear to have a smooth large-N limit. For N >= 4,
it is shown that the trace of the Polyakov loop is insufficient to characterize
the phase structure; the free energy is best described using the eigenvalues of
the Polyakov loop. In both models, the confined phase is characterized by a
mutual repulsion of Polyakov loop eigenvalues that makes the Polyakov loop
expectation value zero. In the deconfined phase, the rotation of the
eigenvalues in the complex plane towards 1 is responsible for the approach to
the blackbody limit over the range T_d - 5T_d. The addition of massless quarks
in SU(3) breaks Z(3) symmetry weakly and eliminates the deconfining phase
transition. In contrast, a first-order phase transition persists with
sufficiently heavy quarks.
