This work presents a model where the Navier-Stokes equation is coupled to the advection-diffusion equation. The decomposition method is applied to the two-dimensional non-stationary Navier-Stokes equation and Duhamelβs Principle is used to obtain a series approximation of the solution. The concentration is represented by the Gaussian dispersion model, assuming constant diffusion coefficients