Gardner Carl L, Jerome Joseph W, Eisenberg Robert S
Department of Mathematics, Arizona State University, Tempe, AZ 85287-1804, USA.
J Theor Biol. 2002 Dec 7;219(3):291-9. doi: 10.1006/jtbi.2002.3123.
Numerical methods are presented for simulating stochastic-in-time current pulses for an electrodiffusion model of the biological channel, with a fixed applied voltage across the channel. The electrodiffusion model consists of the parabolic advection-diffusion equation coupled either to Gauss' law or Poisson's equation, depending on the choice of boundary conditions. The TRBDF2 method is employed for the advection-diffusion equation. The rectangular wave shape of previously simulated traveling wave current pulses is preserved by the full set of partial differential equations for electrodiffusion.
