I am trying to solve convection diffusion equation in a rectangular domain with some hole inside. The periphery of these hole are applied to no slip boundary condition. The flow has been initiated as a constant flow from a narrow inlet. Due to concentration gradient this could provide a flow pattern. But I am getting same contour plot for every time step. What is wrong about my code here? It would be very helpful for me if I can get any guidance from you.