  We simulate the electrical response of multiple disjoint biological 3D cells
in the electropermeabilization process. Instead of solving the boundary value
problem in the volume, we reduce it to a system of boundary integrals equations
with nonlinear dynamics on the cell membranes via a coupling the local Multiple
Traces Formulation with a time semi-implicit scheme. Spatially, boundary
unknowns are approximated by spherical harmonics, thereby allowing for spectral
convergence rates for suitable time steps. Numerical results are provided to
validate our claims.
