Abstract:In this paper, a decoupled characteristic Galerkin finite element procedure is provided for simulating growth-mediated autochemotactic pattern formation in self-propelling bacteria. In this procedure, a modified characteristic Galerkin method is established to solve the bacterial density equation, while the classical finite element procedure is considered for the self-secreted chemical density and polarization dynamics equations system. The convergence of this proposed method is considered under some regularity assumptions and the corresponding error estimate is derived. Numerical experiments are carried out to support the theoretical analysis. Furthermore, several new wave type pattern formations are found.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to simulate growth - mediated autochemotactic pattern formation in self - propelling bacteria. Specifically, the authors focus on the dynamic process of group behavior and pattern formation triggered by the chemicals secreted by the bacteria themselves during the bacterial reproduction process. This problem involves the interaction of multiple variables, including bacterial density \(\rho\), polarization \(p\) and autocrine chemical concentration \(c\). In order to describe the relationships between these variables more accurately, the authors introduced the classical logistic growth term to simulate the reproduction and death of bacteria, and used a series of partial differential equations (PDEs) to describe the changes of these variables over time and space:
\[
\begin{cases}
\frac{\partial \rho}{\partial t}=-\nabla\cdot(\rho v_0 p)+D_\rho\nabla^2\rho+\alpha\rho(1 - \frac{\rho}{K}),\\
\frac{\partial c}{\partial t}=D_c\nabla^2 c + k_0\rho - k_d c + k_a\nabla\cdot(\rho p),\\
\frac{\partial p}{\partial t}=-\gamma p+D_p\nabla^2 p+\beta\nabla c-\gamma_2|p|^2 p.
\end{cases}
\]
where:
- \(\rho\) is the bacterial density,
- \(c\) is the concentration of autocrine chemicals,
- \(p\) is the polarization vector,
- \(v_0\) is the self - propulsion speed of bacteria,
- \(D_\rho\), \(D_c\) and \(D_p\) are the diffusion coefficients of bacterial density, chemical concentration and polarization respectively,
- \(\alpha\) is the growth rate, \(K\) is the environmental carrying capacity,
- \(\beta\) is the chemotactic intensity, positive \(\beta\) indicates chemotactic attraction, and negative \(\beta\) indicates chemotactic repulsion,
- \(\gamma\) and \(\gamma_2\) are the polarization decay rate and saturation parameter respectively.
In order to better approximate the solutions of the above equations, the authors proposed a new modified characteristic Galerkin finite element method, which combines the advantages of the classical finite element method, can maintain mass conservation, and theoretically proved its convergence and error estimation. In addition, numerical experiments show that this method can effectively simulate the growth - mediated autochemotactic pattern formation process and discover several new waveform patterns.
In summary, this paper aims to study the growth - mediated autochemotactic pattern formation problem in self - propelling bacteria through numerical analysis methods, proposes a new numerical method and verifies its effectiveness and accuracy.