Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical approximation of a coagulation-Fragmentation Model for Animal Group Size Statistics
55
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We study numerically a coagulation-fragmentation model derived by Niwa and further elaborated by Degond et al., where a unique equilibrium distribution of group sizes is shown to exist in both cases of continuous and discrete group size distributions. We provide a numerical investigation of these equilibria using three different methods to approximate the equilibrium: a recursive algorithm based on the work of Ma et. al., a Newton method and the resolution of the time-dependent problem. All three schemes are validated by showing that they approximate the predicted small and large size asymptotic behaviour of the equilibrium accurately. The recursive algorithm is used to investigate the transition from discrete to continuous size distributions and the time evolution scheme is exploited to show uniform convergence to equilibrium in time and to determine convergence rates.
rate research
Read More
We translate a coagulation-framentation model, describing the dynamics of animal group size distributions, into a model for the population distribution and associate the blue{nonlinear} evolution equation with a Markov jump process of a type introduced in classic work of H.~McKean. In particular this formalizes a model suggested by H.-S. Niwa [J.~Theo.~Biol.~224 (2003)] with simple coagulation and fragmentation rates. Based on the jump process, we develop a numerical scheme that allows us to approximate the equilibrium for the Niwa model, validated by comparison to analytical results by Degond et al. [J.~Nonlinear Sci.~27 (2017)], and study the population and size distributions for more complicated rates. Furthermore, the simulations are used to describe statistical properties of the underlying jump process. We additionally discuss the relation of the jump process to models expressed in stochastic differential equations and demonstrate that such a connection is justified in the case of nearest-neighbour interactions, as opposed to global interactions as in the Niwa model.
Context. Grains in circumstellar disks are believed to grow by mutual collisions and subsequent sticking due to surface forces. Results of many fields of research involving circumstellar disks, such as radiative transfer calculations, disk chemistry, magneto-hydrodynamic simulations largely depend on the unknown grain size distribution. Aims. As detailed calculations of grain growth and fragmentation are both numerically challenging and computationally expensive, we aim to find simple recipes and analytical solutions for the grain size distribution in circumstellar disks for a scenario in which grain growth is limited by fragmentation and radial drift can be neglected. Methods. We generalize previous analytical work on self-similar steady-state grain distributions. Numerical simulations are carried out to identify under which conditions the grain size distributions can be understood in terms of a combination of power-law distributions. A physically motivated fitting formula for grain size distributions is derived using our analytical predictions and numerical simulations. Results. We find good agreement between analytical results and numerical solutions of the Smoluchowski equation for simple shapes of the kernel function. The results for more complicated and realistic cases can be fitted with a physically motivated black box recipe presented in this paper. Our results show that the shape of the dust distribution is mostly dominated by the gas surface density (not the dust-to-gas ratio), the turbulence strength and the temperature and does not obey an MRN type distribution.
We consider two simple models for the formation of polymers where at the initial time, each monomer has a certain number of potential links (called arms in the text) that are consumed when aggregations occur. Loosely speaking, this imposes restrictions on the number of aggregations. The dynamics of concentrations are governed by modifications of Smoluchowskis coagulation equations. Applying classical techniques based on generating functions, resolution of quasi-linear PDEs, and Lagrange inversion formula, we obtain explicit solutions to these non-linear systems of ODEs. We also discuss the asymptotic behavior of the solutions and point at some connexions with certain known solutions to Smoluchowskis coagulation equations with additive or multiplicative kernels.
Structured population models are a class of general evolution equations which are widely used in the study of biological systems. Many theoretical methods are available for establishing existence and stability of steady states of general evolution equations. However, except for very special cases, finding an analytical form of stationary solutions for evolution equations is a challenging task. In the present paper, we develop a numerical framework for computing approximations to stationary solutions of general evolution equations, which can also be used to produce existence and stability regions for steady states. In particular, we use the Trotter-Kato Theorem to approximate the infinitesimal generator of an evolution equation on a finite dimensional space, which in turn reduces the evolution equation into a system of ordinary differential equations. Consequently, we approximate and study the asymptotic behavior of stationary solutions. We illustrate the convergence of our numerical framework by applying it to a linear Sinko-Streifer structured population model for which the exact form of the steady state is known. To further illustrate the utility of our approach, we apply our framework to nonlinear population balance equation, which is an extension of well-known Smoluchowksi coagulation-fragmentation model to biological populations. We also demonstrate that our numerical framework can be used to gain insight about the theoretical stability of the stationary solutions of the evolution equations. Furthermore, the open source Python program that we have developed for our numerical simulations is freely available from our Github repository (github.com/MathBioCU).
We study the spectral statistics of the Dirac operator on a rose-shaped graph---a graph with a single vertex and all bonds connected at both ends to the vertex. We formulate a secular equation that generically determines the eigenvalues of the Dirac rose graph, which is seen to generalise the secular equation for a star graph with Neumann boundary conditions. We derive approximations to the spectral pair correlation function at large and small values of spectral spacings, in the limit as the number of bonds approaches infinity, and compare these predictions with results of numerical calculations. Our results represent the first example of intermediate statistics from the symplectic symmetry class.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
