Keller-Segel chemotaxis model is described by a system of nonlinear partial differential equations:
a convection diffusion equation for the cell density coupled with a reaction-diffusion equation for
chemoattractant concentration. In this work, we study the phenomenon of Keller-Segel model coupled with
Boussinesq equations. The main objective of this work is to study the global existence and uniqueness and
boundedness of the weak solution for the problem, which is carried out by the Galerkin method