A hybrid numerical method to compute erythrocyte TMP in low-frequency electric fields

This paper presents a coupling method of the finite element method and the boundary element method to compute the transmembrane potential (TMP) of an erythrocyte in a low-frequency electric field. We compute an in vitro erythrocyte's TMP induced by external electric fields by this hybrid method. It takes advantage of the homogeneous characteristics from both intracellular region and extracellular region. Moreover, we may use a fine three-dimensional (3-D) mesh around the thin membrane and avoid 3-D meshes in other regions. Numerical results of a spherical cell show that the hybrid method is accurate. The computed threshold of the applied electric field for membrane electric breakdown agrees well with those experimental results. Numerical results can also guide us to locate the maximum induced TMP on the erythrocyte membrane in various electric fields. Some further applications of the hybrid method are also discussed.