- Journal Home
- Volume 36 - 2024
- Volume 35 - 2024
- Volume 34 - 2023
- Volume 33 - 2023
- Volume 32 - 2022
- Volume 31 - 2022
- Volume 30 - 2021
- Volume 29 - 2021
- Volume 28 - 2020
- Volume 27 - 2020
- Volume 26 - 2019
- Volume 25 - 2019
- Volume 24 - 2018
- Volume 23 - 2018
- Volume 22 - 2017
- Volume 21 - 2017
- Volume 20 - 2016
- Volume 19 - 2016
- Volume 18 - 2015
- Volume 17 - 2015
- Volume 16 - 2014
- Volume 15 - 2014
- Volume 14 - 2013
- Volume 13 - 2013
- Volume 12 - 2012
- Volume 11 - 2012
- Volume 10 - 2011
- Volume 9 - 2011
- Volume 8 - 2010
- Volume 7 - 2010
- Volume 6 - 2009
- Volume 5 - 2009
- Volume 4 - 2008
- Volume 3 - 2008
- Volume 2 - 2007
- Volume 1 - 2006
Commun. Comput. Phys., 34 (2023), pp. 132-172.
Published online: 2023-08
Cited by
- BibTex
- RIS
- TXT
Irreversible drift-diffusion processes are very common in biochemical reactions. They have a non-equilibrium stationary state (invariant measure) which does not satisfy detailed balance. For the corresponding Fokker-Planck equation on a closed manifold, using Voronoi tessellation, we propose two upwind finite volume schemes with or without the information of the invariant measure. Both schemes possess stochastic $Q$-matrix structures and can be decomposed as a gradient flow part and a Hamiltonian flow part, enabling us to prove unconditional stability, ergodicity and error estimates. Based on the two upwind schemes, several numerical examples – including sampling accelerated by a mixture flow, image transformations and simulations for stochastic model of chaotic system – are conducted. These two structure-preserving schemes also give a natural random walk approximation for a generic irreversible drift-diffusion process on a manifold. This makes them suitable for adapting to manifold-related computations that arise from high-dimensional molecular dynamics simulations.
}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.OA-2023-0021}, url = {http://global-sci.org/intro/article_detail/cicp/21883.html} }Irreversible drift-diffusion processes are very common in biochemical reactions. They have a non-equilibrium stationary state (invariant measure) which does not satisfy detailed balance. For the corresponding Fokker-Planck equation on a closed manifold, using Voronoi tessellation, we propose two upwind finite volume schemes with or without the information of the invariant measure. Both schemes possess stochastic $Q$-matrix structures and can be decomposed as a gradient flow part and a Hamiltonian flow part, enabling us to prove unconditional stability, ergodicity and error estimates. Based on the two upwind schemes, several numerical examples – including sampling accelerated by a mixture flow, image transformations and simulations for stochastic model of chaotic system – are conducted. These two structure-preserving schemes also give a natural random walk approximation for a generic irreversible drift-diffusion process on a manifold. This makes them suitable for adapting to manifold-related computations that arise from high-dimensional molecular dynamics simulations.