The boundary conditions for this simulation assumed that no current was flowing into or
out of the control volume. Therefore, insulating Neumann boundary conditions were assigned to
every face of the volume surrounding the heart. Initial conditions characterize an initial uniform
potential of distribution of 1 V throughout a section of the model for the activation variable,
while the adjacent sections remain at zero. For the inhibitor variable, the adjacent sections have a
value of 0.3 V.
Values for F, a, 0, y and 6 were obtained from literature from an FEM model done by
Filippi et al. 29 to be 0.01, 0.1, 0.5, 1 and 0 respectively, which are standard values used in simple
FitzHugh-Nagumo models. Preliminary data from these models is presented.