Year: 2008
Author: Hongyun Wang
International Journal of Numerical Analysis and Modeling, Vol. 5 (2008), Iss. 1 : pp. 1–23
Abstract
In mathematical studies of molecular motors, the stochastic motor motion is modeled using the Langevin equation. If we consider an ensemble of motors, the probability density is governed by the corresponding Fokker-Planck equation. Average quantities, such as, average velocity, effective diffusion and randomness parameter, can be calculated from the probability density. The WPE method was previously developed to solve Fokker-Planck equations (H. Wang, C. Peskin and T. Elston, J. Theo. Biol., Vol. 221, 491-511, 2003). The WPE method has the advantage of preserving detailed balance, which ensures that the numerical method still works even when the potential is discontinuous. Unfortunately, the accuracy of the WPE method drops to first order when the potential is discontinuous. Here we propose an improved version of the WPE method. The improved WPE method a) maintains the second order accuracy even when the potential is discontinuous, b) has got rid of a numerical singularity in the WPE method, and c) is as simple and easy to implement as the WPE method. Numerical examples are shown to demonstrate the robust performance of the improved WPE method.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/2008-IJNAM-794
International Journal of Numerical Analysis and Modeling, Vol. 5 (2008), Iss. 1 : pp. 1–23
Published online: 2008-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 23
Keywords: Fokker-Planck equation detailed balance numerical solutions.