A general 2D mathematical model was developed to simulate the purification of water from volatile organic compounds (VOCs) via vacuum membrane distillation (VMD) process in hollow fiber membrane contactors. The model was developed for hydrophobic membrane material conditions, taking into consideration axial and radial diffusion in the tube, membrane and compartments of the contactor and was simplified to the two-dimensional structure with a single porous membrane wall. The simulation has studied the mass and heat transfer of VMD system in the porous media, in which aqueous volatile organic solution was considered as an incompressible and steady fluid. Effect of the downstream pressure on the removal of 1, 1, 1-trichloroethane (TCA) was studied to validation of simulation results with experimental data that it was obtained from literature. The temperature, Reynolds number, and total mass flux (convective and diffusive) distribution of TCA are determined in the membrane module. (C) 2013 Society of Plastics Engineers