This paper proposes a pseudo-potential describing slit flow in the formalis
m of lattice Monte Carlo simulation with the bond fluctuation algorithm as
the unique basic micro-relaxation mode. The main characteristics of slit fl
ow, such as the parabolic velocity profile and the pressure-flux relationsh
ip, are successfully reproduced in a three-dimensional self-avoiding multic
hain system. Both Newtonian and non-Newtonian regimes are revealed. The cha
in conformation and nonlinear rheological behavior are investigated. The si
mulation results agree with experimental measurements. This method call be
used to investigate the viscoelastic properties at different layers as well
as the global properties. Some peculiar phenomena in inhomogeneous flow ar
e found to be consistent with previous theoretical predictions by others.