A numerical comparison between degenerate parabolic and quasilinear hyperbolic models of cell movements under chemotaxis

We consider two models which were both designed to describe the movement of eukaryotic cells responding to chemical signals. Besides a common standard parabolic equation for the diffusion of a chemoattractant, like chemokines or growth factors, the two models differ for the equations describing the movement of cells. The first model is based on a quasilinear hyperbolic system with damping, the other one on a degenerate parabolic equation. The two models have the same stationary solutions, which may contain some regions with vacuum. We first explain in details how to discretize the quasilinear hyperbolic system through an upwinding technique, which uses an adapted reconstruction, which is able to deal with the transitions to vacuum. Then we concentrate on the analysis of asymptotic preserving properties of the scheme towards a discretization of the parabolic equation, obtained in the large time and large damping limit, in order to present a numerical comparison between the asymptotic behavior of these two models. Finally we perform an accurate numerical comparison of the two models in the time asymptotic regime, which shows that the respective solutions have a quite different behavior for large times.Comment: One sentence modified at the end of Section 4, p. 1

Stationary solutions with vacuum for a one-dimensional chemotaxis model with non-linear pressure

International audienceIn this article, we study a one-dimensional hyperbolic quasi-linear model of chemotaxis with a non-linear pressure and we consider its stationary solutions, in particular with vacuum regions. We study both cases of the system set on the whole line \Er and on a bounded interval with no-flux boundary conditions. In the case of the whole line \Er, we find only one stationary solution, up to a translation, formed by a positive density region (called bump) surrounded by two regions of vacuum. However, in the case of a bounded interval, an infinite of stationary solutions exists, where the number of bumps is limited by the length of the interval. We are able to compare the value of an energy of the system for these stationary solutions. Finally, we study the stability of these stationary solutions through numerical simulations

Density dependent diffusion models for the interaction of particle ensembles with boundaries

The transition from a microscopic model for the movement of many particles to a macroscopic continuum model for a density flow is studied. The microscopic model for the free flow is completely deterministic, described by an interaction potential that leads to a coherent motion where all particles move in the same direction with the same speed known as a flock. Interaction of the flock with boundaries, obstacles and other flocks leads to a temporary destruction of the coherent motion that macroscopically can be modeled through density dependent diffusion. The resulting macroscopic model is an advection-diffusion equation for the particle density whose diffusion coefficient is density dependent. Examples describing i) the interaction of material flow on a conveyor belt with an obstacle that redirects or restricts the material flow and ii) the interaction of flocks (of fish or birds) with boundaries and iii) the scattering of two flocks as they bounce off each other are discussed. In each case, the advection-diffusion equation is strictly hyperbolic before and after the interaction while the interaction phase is described by a parabolic equation. A numerical algorithm to solve the advection-diffusion equation through the transition is presented.Comment: 25 pages, 9 figure

Hyperbolic Techniques in Modelling, Analysis and Numerics

Several research areas are flourishing on the roots of the breakthroughs in conservation laws that took place in the last two decades. The meeting played a key role in providing contacts among the different branches that are currently developing. All the invitees shared the same common background that consists of the analytical and numerical techniques for nonlinear hyperbolic balance laws. However, their fields of applications and their levels of abstraction are very diverse. The workshop was the unique opportunity to share ideas about analytical issues like the fine-structure of singular solutions or the validity of entropy solution concepts. It turned out that generalized hyperbolic techniques are able to handle the challenges posed by new applications. The design of efficient structure preserving methods turned out to be the major line of development in numerical analysis

Differential Models, Numerical Simulations and Applications

This Special Issue includes 12 high-quality articles containing original research findings in the fields of differential and integro-differential models, numerical methods and efficient algorithms for parameter estimation in inverse problems, with applications to biology, biomedicine, land degradation, traffic flows problems, and manufacturing systems

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal