Adaptive space and time numerical simulation of reaction-diffusion models for intracellular calcium dynamics

详细信息	查看全文 | 推荐本文 |

• 作者：Ch. Nagaiah^a ; ^{nagaiah.chamakuri@uni-graz.at} ; S. Rü ; diger^b ; ^{sten.ruediger@gmail.com} ; G. Warnecke^c ; ^{warnecke@ovgu.de} ; M. Falcke^d ; ^{falcke@hmi.de}
• 关键词：Calcium dynamics ; Reaction&ndash ; diffusion system ; Hybrid stochastic and deterministic coupling ; Adaptive FEM ; Rosenbrock type methods
• 刊名：Applied Mathematics and Computation
• 出版年：2012
• 期刊代码：7_00963003
• 类别：cp
• 出版时间：15 June, 2012
• 卷：218
• 期：20
• 页码：10194-10210
• 文件大小：2787 K

摘要

Adaptivity in space and time for the numerical simulation of stochastic and deterministic equations for intracellular calcium dynamics is presented. The modeling of diffusion, reaction and membrane transport of calcium ions in cells leads to a system of reaction-diffusion equations. We describe the modulation of cytosolic and ER calcium concentrations close to the membrane of the cell organelle.

A conforming piecewise linear finite element method is used for the spatial discretization. Linearly implicit methods of Rosenbrock type are used for the time integration. We adopt a hybrid algorithm to solve the stochastic part. The space grid is adjusted to the strong localization of the calcium release following stochastic channel transitions. By automatically adapting the spatial meshes and time steps to the proper scales during the transition of channel states, the method accurately resolves the evolution of intracellular calcium concentrations as well as buffer concentrations. This article emphasizes adaptive and efficient hybrid numerical simulations in two space dimensions. The presented work establishes the basis for future simulations in a realistic 3D geometry.