In this paper we study bulk viscosity in a thermal QCD model with large number of colors at two extreme limits: the very weak and the very strong 't Hooft couplings. The weak coupling scenario is based on kinetic theory, and one may go to the very strong coupling dynamics via an intermediate coupling regime. Although the former has a clear description in terms of kinetic theory, the intermediate coupling regime, which uses lattice results, suffers from usual technical challenges that render an explicit determination of bulk viscosity somewhat difficult. On the other hand, the very strong 't Hooft coupling dynamics may be studied using string theories at both weak and strong string couplings using gravity duals in type IIB as well as M-theory respectively. In type IIB we provide the precise fluctuation modes of the metric in the gravity dual responsible for bulk viscosity, compute the speed of sound in the medium and analyze the ratio of the bulk to shear viscosities. In M-theory, where we uplift the type IIA mirror dual of the UV complete type IIB model, we study and compare both the bulk viscosity and the sound speed by analyzing the quasi-normal modes in the system at strong IIA string coupling. By deriving the spectral function, we show the consistency of our results both for the actual values of the parameters involved as well for the bound on the ratio of bulk to shear viscosities. 65 5.1 The mirror type IIA model and its M-theory uplift 66 5.2 Quasi-normal modes, attenuation constant and the sound speed 70 5.3 The case with a vanishing bare resolution parameter 76 5.4 Shear viscosity, entropy and the bulk viscosity bound 80 6. Type IIA spectral function and the viscosity bound at strong coupling with non-zero flavors 85 6.1 Background gauge fluxes and perturbations on the flavor branes 86 6.2 Equation of motion for gauge field fluctuations 88 6.3 On-shell action and the strong coupling spectral function 95 6.4 The strong string coupling limit and pure classical supergravity 103 -1 -7. Conclusions and discussions 107 A. A Gauge-Invariant Combination of Scalar Modes of Metric Perturbations 112 A.1 The equation of motion for the fluctuation mode H tt 113 A.2 The equation of motion for the combined mode H s 115 A.3 The equations of motion for the remaining fluctuation modes 116 B. A derivation of the on-shell action and the Green's function 118C. Effective number of three-brane charges with background threeforms and the horizon radius 1201 Violation of this bound is seen in the presence of higher derivative terms, discussed first in [23]. In the absence of these terms, the KSS [22] bound continues to hold at strong 't Hooft coupling.2 The non-conformal string theory studied in [30] is different from what we consider here. In [30] it's the N = 2 * supersymmetric gauge theory obtained by a mass deformation of N = 4 Super Yang Mills theory. See also [31] and [32] for an even earlier study on bulk viscosity from first principles.