Reaction-diffusion-chemotaxis systems have proven to be fairly
accurate mathematical models for many pattern formation problems in chemistry
and biology. These systems are important for computer simulations
of patterns, parameter estimations as well as analysis of the biological systems.
To solve reaction-diffusion-chemotaxis systems, efficient and reliable
numerical algorithms are essential for pattern generations. In this paper, a
general reaction-diffusion-chemotaxis system is considered for specific numerical
issues of pattern simulations. We propose a fully explicit discretization
combined with a variable optimal time step strategy for solving the reactiondiffusion-
chemotaxis system. Theorems about stability and convergence of the
algorithm are given to show that the algorithm is highly stable and efficient.
Numerical experiment results on a model problem are given for comparison
with other numerical methods. Simulations on two real biological experiments
will also be shown.
Citation: Chichia Chiu, Jui-Ling Yu. An optimal adaptive time-stepping scheme for solving reaction-diffusion-chemotaxis systems[J]. Mathematical Biosciences and Engineering, 2007, 4(2): 187-203. doi: 10.3934/mbe.2007.4.187
Related Papers:
[1] |
Fabrizio Clarelli, Roberto Natalini .
A pressure model of immune response to mycobacterium
tuberculosis infection in several space dimensions. Mathematical Biosciences and Engineering, 2010, 7(2): 277-300.
doi: 10.3934/mbe.2010.7.277
|
[2] |
Changwook Yoon, Sewoong Kim, Hyung Ju Hwang .
Global well-posedness and pattern formations of the immune system induced by chemotaxis. Mathematical Biosciences and Engineering, 2020, 17(4): 3426-3449.
doi: 10.3934/mbe.2020194
|
[3] |
Lin Zhang, Yongbin Ge, Xiaojia Yang .
High-accuracy positivity-preserving numerical method for Keller-Segel model. Mathematical Biosciences and Engineering, 2023, 20(5): 8601-8631.
doi: 10.3934/mbe.2023378
|
[4] |
Raimund Bürger, Gerardo Chowell, Elvis Gavilán, Pep Mulet, Luis M. Villada .
Numerical solution of a spatio-temporal predator-prey model with infected prey. Mathematical Biosciences and Engineering, 2019, 16(1): 438-473.
doi: 10.3934/mbe.2019021
|
[5] |
Fiona R. Macfarlane, Mark A. J. Chaplain, Tommaso Lorenzi .
A hybrid discrete-continuum approach to model Turing pattern formation. Mathematical Biosciences and Engineering, 2020, 17(6): 7442-7479.
doi: 10.3934/mbe.2020381
|
[6] |
Ali Moussaoui, Vitaly Volpert .
The impact of immune cell interactions on virus quasi-species formation. Mathematical Biosciences and Engineering, 2024, 21(11): 7530-7553.
doi: 10.3934/mbe.2024331
|
[7] |
Lin Zhang, Yongbin Ge, Zhi Wang .
Positivity-preserving high-order compact difference method for the Keller-Segel chemotaxis model. Mathematical Biosciences and Engineering, 2022, 19(7): 6764-6794.
doi: 10.3934/mbe.2022319
|
[8] |
Xiaomei Bao, Canrong Tian .
Turing patterns in a networked vegetation model. Mathematical Biosciences and Engineering, 2024, 21(11): 7601-7620.
doi: 10.3934/mbe.2024334
|
[9] |
M. B. A. Mansour .
Computation of traveling wave fronts for a nonlinear diffusion-advection model. Mathematical Biosciences and Engineering, 2009, 6(1): 83-91.
doi: 10.3934/mbe.2009.6.83
|
[10] |
Benedetto Bozzini, Deborah Lacitignola, Ivonne Sgura .
Morphological spatial patterns in a reaction diffusion model for metal growth. Mathematical Biosciences and Engineering, 2010, 7(2): 237-258.
doi: 10.3934/mbe.2010.7.237
|
Abstract
Reaction-diffusion-chemotaxis systems have proven to be fairly
accurate mathematical models for many pattern formation problems in chemistry
and biology. These systems are important for computer simulations
of patterns, parameter estimations as well as analysis of the biological systems.
To solve reaction-diffusion-chemotaxis systems, efficient and reliable
numerical algorithms are essential for pattern generations. In this paper, a
general reaction-diffusion-chemotaxis system is considered for specific numerical
issues of pattern simulations. We propose a fully explicit discretization
combined with a variable optimal time step strategy for solving the reactiondiffusion-
chemotaxis system. Theorems about stability and convergence of the
algorithm are given to show that the algorithm is highly stable and efficient.
Numerical experiment results on a model problem are given for comparison
with other numerical methods. Simulations on two real biological experiments
will also be shown.