In this article, a probabilistic approach is applied to evaluate the impact of the GDL porosity uncertainty on the electrical performances of a PEMFC. The study is based on the use of a dynamical, symbolic, and acausal knowledge model. Some statistical distributions are introduced on the input model parameter (porosity) and the statistical distributions induced on the output parameters (cell voltage, resistance) are analyzed. The difference observed between the shapes of the input and output distributions (respectively some Gaussian and inverse Gamma distributions with a threshold phenomenon) is the result of strong nonlinearities linked with the integration of multiphase flow phenomena in the modelling (e.g. diffusion limit of the humidified air in the GDL). The study is also conducted for different conditions of temperature and pressure through a design of numerical experiments. One of the results obtained is that the variation coefficient related to the GDL porosity has, compared to the other parameters with their intervals of variation considered, little effect on the average output distributions. However, the dispersion introduced on the porosity impacts their shapes (e.g. significant effect on the standard deviation). (C) 2016 Hydrogen Energy Publications LLC. Published by Elsevier Ltd. All rights reserved.