Modeling for cardiac excitation propagation based on the Nernst-Planck equation and homogenization

The bidomain model is a commonly used mathematical model of the electrical properties of the cardiac muscle that takes into account the anisotropy of both the intracellular and extracellular spaces. However, the equations contain self-contradiction such that the update of ion concentrations does not consider intracellular or extracellular ion movements due to the gradient of electric potential and the membrane charge as capacitive currents in spite of the fact that those currents are taken into account in forming Kirchhoff's first law. To overcome this problem, we start with the Nernst-Planck equation, the ionic conservation law, and the electroneutrality condition at the cellular level, and by introducing a homogenization method and assuming uniformity of variables at the microscopic scale, we derive rational bidomain equations at the macroscopic level.