We present a detailed study of the interactions among scalar and pseudoscalar mesons and glueballs within the framework of the generalized linear sigma model in the SU(3) flavor limit. The basis of our approach is to develop a global understanding of light scalar and pseudoscalar mesons (up to around 2 GeV) and thereby explore the underlying mixings among composite quark matter fields and glueballs. The chiral sector of the Lagrangian is formulated in terms of two chiral nonets representing quark-antiquarks and tetraquarks (in the leading order, this sector contains terms with eight or fewer number of quak or antiquark lines). The Lagrangian also contains a sector that represents scalar and pseudoscalar glueballs and their interactions with the matter chiral fields, in a manner that the axial and trace anomalies of QCD are exactly realized. In this construct, the model has two scalar octets and two pdeudoscalar octets (each a linear combination of two-and four-quark components), as well as three scalar SU(3) singlets and three pseudoscalar SU(3) singlets (each a linear combination of two-and four-quark components as well as a glueball component). With the inputs of the experimental masses of a0(980), a0(1450), and the masses of π(137), π(1300) and their decay constants, we perform an extensive numerical simulation to determine the boundaries of the parameter space of the model. We further incorporate experimental data on the mass spectrum of eta states as well as on several decay widths and decay ratios of f0 states to zoom in on the parameter space and make predictions for the substructure of pseudoscalar and scalar SU(3) octets and singlets as well as for the pseudoscalar and scalar glueball masses. We find the scalar and pseudoscalar glueball masses around 1.6 and 2.0 GeV, respectively.