Discharge probability in GEM-based gaseous detectors has been numerically estimated using an axisymmetric hydrodynamic model.
Initial primary charge configurations in the drift region, obtained using Heed and Geant4, are found to have significant effect on the subsequent evolution of detector response.
Simulation of energy resolution has been performed to establish the capability of the hydrodynamic model to capture statistical nature of the experimental situation.
Finally, single and triple GEM configurations exposed to alpha sources have been simulated to estimate discharge probability which have been compared with available experimental data.
Despite the simplifying and drastic assumptions in the numerical model, the comparisons are encouraging.