Title :
Hybrid Finite Element Method for Describing the Electrical Response of Biological Cells to Applied Fields
Author :
Ying, Wenjun ; Henriquez, Craig S.
Author_Institution :
Dept. of Biomed. Eng., Duke Univ., Durham, NC
fDate :
4/1/2007 12:00:00 AM
Abstract :
A novel hybrid finite element method (FEM) for modeling the response of passive and active biological membranes to external stimuli is presented. The method is based on the differential equations that describe the conservation of electric flux and membrane currents. By introducing the electric flux through the cell membrane as an additional variable, the algorithm decouples the linear partial differential equation part from the nonlinear ordinary differential equation part that defines the membrane dynamics of interest. This conveniently results in two subproblems: a linear interface problem and a nonlinear initial value problem. The linear interface problem is solved with a hybrid FEM. The initial value problem is integrated by a standard ordinary differential equation solver such as the Euler and Runge-Kutta methods. During time integration, these two subproblems are solved alternatively. The algorithm can be used to model the interaction of stimuli with multiple cells of almost arbitrary geometries and complex ion-channel gating at the plasma membrane. Numerical experiments are presented demonstrating the uses of the method for modeling field stimulation and action potential propagation
Keywords :
Runge-Kutta methods; bioelectric potentials; biological effects of fields; biomembrane transport; cellular effects of radiation; finite element analysis; initial value problems; integration; linear differential equations; nonlinear differential equations; physiological models; Euler method; Runge-Kutta method; action potential propagation; biological cells; biological membranes; complex ion-channel gating; differential equations; electric flux; electrical response; hybrid finite element method; initial value problem; linear interface problem; linear partial differential equation; membrane currents; membrane dynamics; nonlinear initial value problem; nonlinear ordinary differential equation; plasma membrane; time integration; Biological cells; Biological system modeling; Biomembranes; Cells (biology); Differential equations; Finite element methods; Geometry; Partial differential equations; Plasmas; Solid modeling; Field stimulation; Laplace equation; hybrid finite element method; interface problem; transmembrane potential; Animals; Cell Membrane; Cell Physiology; Computer Simulation; Dose-Response Relationship, Radiation; Electric Stimulation; Electromagnetic Fields; Finite Element Analysis; Humans; Membrane Potentials; Models, Biological; Radiation Dosage;
Journal_Title :
Biomedical Engineering, IEEE Transactions on
DOI :
10.1109/TBME.2006.889172