Numerical simulations based on Monte Carlo dynamics are used to investigate the resistivity behavior of granular superconductors containing a random distribution of pi junctions, as in superconducting materials with d-wave symmetry. The presence of pi junctions leads to quenched in circulating currents (chiralities) and to chiral glass behavior at low temperatures, even without an external magnetic field. An XY spin glass model in the phase representation is used to determine the current-voltage characteristics and critical exponents of the resistivity transition. In two dimensions, the linear resistivity is nonzero at finite temperatures and the dynamic scaling analysis of the nonlinear resistivity is consistent with a phase transition at zero temperature. In three dimensions, we find a transition at finite temperatures below which the linear resistivity vanishes and the corresponding critical exponents are determined from the scaling analysis. The results are in good agreement with Langevin simulations in the phase representation. The dynamic exponent z is significantly different from previous results obtained in the vortex representation.
