Electrodiffusion Model Simulation of Ionic Channels: 1D Simulations

The drift-diffusion (Poisson-Nernst-Planck) model is applied to ionic channels in biological membranes plus surrounding solution baths. Simulations of the K channel in KCl solutions using the TRBDF2 method are presented which show significant boundary layers at the ends of the channel. The computed current-voltage curve for the K channel shows excellent agreement with experimental measurements.