We give an introduction to a fuel cell model formulated as a nonsmooth differential-algebraic equation system with strong non-local coupling. A successful iterative method has been developed to solve this system. We talk about some improvements to make the method more robust and faster.