We study two important properties of 2+1D QCD, namely confinement and Pseudoscalar glueball spectrum, using holographic approach. The confined state of the bounded quark-antiquark pair occurs in the self-coupling dominated nonperturbative regime, where the free gluons form the bound states, known as glueballs. The gauge theory corresponding to low energy decoupled geometry of isotropic non-supersymmetric D2 brane, which is again similar to the 2+1D YM theory, has been taken into account but in this case the coupling constant is found to vary with the energy scale. At BPS limit, this theory reduces to supersymmetric YM theory. We have considered NG action of a test string and calculate the potential of such confined state located on the boundary. The QCD flux tube tension for large quark-antiquark separation is observed to be a monotonically increasing function of running coupling. The mass spectrum of Pseudoscalar glueball is evaluated numerically from the fluctuations of the axion in the gravity theory using WKB approximation. This produces the mass to be related to the string tension and the levels of the first three energy states. The various results that we obtained quite match with those previously studied through the lattice approach.