Modeling mass transfer in the liquid boundary layer accounting for concentration polarization in pervaporation (PV) is particularly challenging since there is no practical way of experimentally determining solute concentration at the membrane surface. We have developed a computational fluid dynamics (CFD) approach to describe not only velocity distribution but also concentration profile in the liquid boundary layer of a slit membrane channel. The satisfactoriness of the numerical methodology used in CFD for obtaining concentration profiles were verified using a classic diffusion problem with its known analytical solution. The overall mass transfer coefficients from the numerical study were also compared with those from the experiment. (C) 2004 Elsevier Ltd. All rights reserved.