Slot coating is a common method in the manufacturing of a wide variety of products. It belongs to a class of coating method known as premetered coating: in a steady state operation, the thickness of the coated liquid layer is set by the flow rate fed to the die and the speed of the substrate moving past, and is independent of other variables of the process. For many applications, the coating liquid is a particle suspension. The commom simplified approach is to study the flow as Newtonian and evaluate its viscosity based on the average particle concentration. However, experimental data shows that particle distribution is in fact non-uniform in shear flows. Moreover, particle distribution along the film thickness during the coating process will affect the final fim structure and consequently product performance. Hence, it is important to understand the fundamentals of coating process with particle suspension to better describe and predict the behavior of the flow and the particle distribution in the coated layer. The goal of the present work is to study different particle migration mechanisms in a suspension that may affect particle distribution in the coated film; such as diffusive mechanisms and sedimentation. The model presented takes into account the viscosity dependency and local particle concentration and surface-active particles, i.e. Marangoni effect. Two different approaches were used to study Marangoni flow. First, the bulk-interface transport is considered much faster than bulk diffusion such that the net flux is zero and surface concentration is equal bulk concentration. The second, more complete model, considers adsorption and desorption fluxes from the interface. The non-linear system of equations, with appropriate boundary conditions, is solved by Galerkin Finite Element Method and Newton s Method. The flow is two-dimensional and therefore two-dimensional elements are used to describe flow domain. Nonetheless, in order to evaluate particle concentration along the interface a one-dimensional element was created. The coupling between those two classes of elements of discussed. The results show that the particle distribution on the coated layer is a strong function of the film thickness and suspension properties, such as particle density and diffusion coefficients. The Marangoni flow associated with surface tension gradient due to particle concentration variation along the interface can change the recirculation pattern and particle concentration distribution. Furthermore, the results show that particle adsorption and desorption from interface have a strong effect on the particle distribution.