“…Details on the computation of these integrals are discussed in detail for axisymmetric and 2D geometries. 13,15 Discretized forms of these three equations are solved simultaneously for r(s) on electrodes and c(c) on interfaces. Equations (7) to (9) invoke superposition of sources that include free charge, r, on electrodes (S 0 ); polarization charge, c, and trapped charge, k, on material and physical interfaces (C 0 ); and mobile and trapped space charge, q in the volume (V 0 ).…”