We construct a hierarchy of loop equations for invariant circular ensembles.These are valid for general classes of potentials and for arbitrary inversetemperatures $ {\rm Re}\,\beta>0 $ and number of eigenvalues $ N $. Usingmatching arguments for the resolvent functions of linear statistics $f(\zeta)=(\zeta+z)/(\zeta-z) $ in a particular asymptotic regime, the globalregime, we systematically develop the corresponding large $ N $ expansion andapply this solution scheme to the Dyson circular ensemble. Currently we cancompute the second resolvent function to ten orders in this expansion and alsoits general Fourier coefficient or moment $ m_{k} $ to an equivalent length.The leading large $ N $, large $ k $, $ k/N $ fixed form of the moments can berelated to the small wave-number expansion of the structure function in thebulk, scaled Dyson circular ensemble, known from earlier work. From the momentexpansion we conjecture some exact partial fraction forms for the low $ k $moments. For all of the forgoing results we have made a comparison with theexactly soluble cases of $ \beta = 1,2,4 $, general $ N $ and even, positive $\beta $, $ N=2,3 $.