A three-dimensional and unsteady proton exchange membrane fuel cell (PEMFC) model with serpentine-parallel channels has been incorporated to simulate not only the fluid flow, heat transfer. species transport, electrochemical reaction, and current density distribution but also the behaviors of liquid water in the gas-liquid flow of the channels and porous media. Using this general model, the behaviors of liquid water were investigated by performing the motion, deformation, coalescence and detachment of water droplets inside the channels and the penetration of liquid through the porous media at different time instants. The results showed that: tracking the interface of liquid water in a reacting gas-liquid flow in PEMFC can be fulfilled by using volume-of-fluid (VOF) algorithm combined with solving the conservation equations of continuity, momentum, energy, species transport and electrochemistry: the presence of liquid water in the channels has a significant impact on the flow fields, e.g., the gas flow became unevenly distributed due to the blockage of liquid water where the high pressure would be suddenly built up and the reactant gas transport in the channels and porous media would be hindered by liquid water occupation. (C) 2008 Elsevier Ltd. All rights reserved.