The spectrum of two-gluon glueballs below 3 GeV is investigated in a potential model with dynamical gluon mass using variational method. The short distance potential is approximated by one-gluon exchange, while the long distance part is taken as a breakable string. The mass and size of the radial as well as orbital excitations up to principle quantum number n = 3 are evaluated. The predicted mass ratios are compared with experimental and lattice results.PACS numbers: 12.39.MK Quantum Chromodynamics (QCD) is widely accepted as the theory of strong interactions. It is generally believed that the gluon self-coupling in QCD implies the existence of bound states of confined gluons known as glueballs. The experimental discovery of these glueballs would be very important and would give further support to the theory of QCD. However, numerous technical difficulties have so far hampered our unequivocal identification of glueballs by experiment, largely because glueball states can mix strongly with nearby qq resonances. Nevertheless, the estimation of mass and size of pure gluon glueball states should still be pursued. This could guide experimental searches, as well as provide calibration for models of glueballs.Over the past 20 years there has been an on-going effort to obtain a nonperturbative form for the gluon propagator. Perhaps one of the most interesting result is that the gluon may have a dynamically generated mass [1]. The existence of a mass scale, or the absence of a pole at k 2 = 0, is natural if one assumes that gluons do not propagate to infinity, i.e., these propagators describe confined gluons. The concept of massive gluon has been widely used in independent field theoretic studies, and examples about the consequences of massive gluons can be found in the literature [2][3][4][5][6][7].In this paper, we focus on the calculation of two-gluon glueball systems and extend our previous work [8] on the estimation of the mass and size of low lying glueball states, using the variational method in the potential model of Cornwall and Soni [9,10]. The main feature of the present work is the consideration of radial as well as orbital excitations, up to principle quantum number n = 3.To exhibit both asymptotic freedom and the nonAbelian nature of QCD, gluon dynamics can be described as massive spin-one fields interacting through one-gluon exchange and a breakable string. At short distance the effective coupling constant of the gluon-gluon interaction becomes small and the interaction can be treated perturbatively. The short distance potential is approximated by one-gluon exchange and can be extracted from the tree-level Feynman amplitude of Fig. 1,where q is the momentum transfer of the system. At long distance, the non-Abelian nature of QCD implies gluon confinement via nonperturbative effects. These nonperturbative effects are implemented by introducing a string potential V str which is assumed to be spin independent,where β is related to the adjoint string tension K A viaIn the potential V str , the color screening of...