Algorithm for numerical solutions to the kinetic equation of a spatial population dynamics model with coalescence and repulsive jumps

DOI10.1007/S11075-020-00992-9MaRDI QIDQ2028046zbMATH OpenOpenAlexFDO

Authors Yuri Kozitsky, Krzysztof Pilorz, Igor P. Omelyan

Publication date 31 May 2021

Published in Numerical Algorithms (Search for Journal in Brave)

Copyright license Creative Commons Attribution 4.0 International

Full work available at URL https://arxiv.org/abs/2007.15969

zbMATH Keywords

coalescencehoppingintegro-differential equationsnumerical methodspopulation dynamicsspatial inhomogeneity

Mathematics Subject Classification ID

PDEs in connection with biology, chemistry and other natural sciences (35Q92) Integro-partial differential equations (35R09) Integro-partial differential equations (45K05) Jump processes on discrete state spaces (60J74) Jump processes on general state spaces (60J76) Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations (65L06) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Dynamic continuum models (systems of particles, etc.) in time-dependent statistical mechanics (82C21) Population dynamics (general) (92D25)

Abstract: An algorithm is proposed for finding numerical solutions of a kinetic equation that describes an infinite system of point articles placed in $mathbb{R}^d(dgeq1)$. The particles perform random jumps with pair wise repulsion, in the course of which they can also merge. The kinetic equation is an essentially nonlinear and nonlocal integro-differential equation, which can hardly be solved analytically. The derivation of the algorithm is based on the use of space-time discretization, boundary conditions, composite Simpson and trapezoidal rules, Runge-Kutta methods, adjustable system-size schemes, etc. The algorithm is then applied to the one-dimensional version of the equation with various initial conditions. It is shown that for special choices of the model parameters, the solutions may have unexpectable time behaviour. A numerical error analysis of the obtained results is also carried out.

Review text

The purpose of this paper is development of an efficient algorithm to solve the repulsion-jump coalescence kinetic equation in a high-dimensional space. Different techniques including time-space discretization (composite Simpson and trapezoidal rules, Runge-Kutta methods, and adjustable system-size schemes), periodic, Dirichlet, and asymptotic boundary conditions are combined for studying infinite systems on the basis of their finite samples. The algorithm is applied to population dynamics simulations of one-dimensional systems with various initial spatially inhomogeneous densities and forms of the jump, coalescence, and repulsion kernels. A comprehensive discussion of the obtained numerical results is provided for some specific choices of the model parameters and initial densities.

Recommendations

• An efficient stochastic algorithm for solving the Boltzmann equation
• Jumps and coalescence in the continuum: a numerical study
• Stochastic Lagrangian models and algorithms for spatially inhomogeneous Smoluchowski equation
• Approximation of coalescence integrals in population balance models with local mass conservation
• Solving the transport-coagulation problem in a two-dimensional spatial region

Cites work

• A kinetic equation for repulsive coalescing random jumps in continuum
• Coalescing Brownian flows: a new approach
• Flows, coalescence and noise.
• Fragmentation-Coagulation Models of Phytoplankton
• Jumps and coalescence in the continuum: a numerical study
• Kawasaki dynamics in continuum: micro- and mesoscopic descriptions
• Modified massive Arratia flow and Wasserstein diffusion
• Numerical recipes. The art of scientific computing.
• Phytoplankton Dynamics: from the Behavior of Cells to a Transport Equation
• Random jumps and coalescence in the continuum: evolution of states of an infinite particle system
• Spatially inhomogeneous population dynamics: beyond the mean field approximation
• The global evolution of states of a continuum Kawasaki model with repulsion
• The true self-repelling motion
• Vehicular traffic, crowds, and swarms: from kinetic theory and multiscale methods to applications and research perspectives

Cited in

(3)

• An interplay between attraction and repulsion in infinite populations
• Jumps and coalescence in the continuum: a numerical study
• A simple spatial model of population dynamics

This page was built for publication: Algorithm for numerical solutions to the kinetic equation of a spatial population dynamics model with coalescence and repulsive jumps

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2028046)