We propose a positivity preserving entropy decreasing finite volume scheme for nonlinear nonlocal equations with a gradient flow structure. These properties allow for accurate computations of stationary states and long-time asymptotics demonstrated by suitably chosen test cases in which these features of the scheme are essential. The proposed scheme is able to cope with non-smooth stationary states, different time scales including metastability, as well as concentrations and self-similar behavior induced by singular nonlocal kernels. We use the scheme to explore properties of these equations beyond their present theoretical knowledge.

Aquatic bacteria like Bacillus subtilis are heavier than water yet they are able to swim up an oxygen gradient and concentrate in a layer below the water surface, which will undergo Rayleigh-Taylor-type instabilities for sufficiently high concentrations. In the literature, a simplified chemotaxis-fluid system has been proposed as a model for bioconvection in modestly diluted cell suspensions. It couples a convective chemotaxis system for the oxygen-consuming and oxytactic bacteria with the incompressible Navier-Stokes equations subject to a gravitational force proportional to the relative surplus of the cell density compared to the water density. In this paper, we derive a high-resolution vorticity-based hybrid finite-volume finite-difference scheme, which allows us to investigate the nonlinear dynamics of a two-dimensional chemotaxis-fluid system with boundary conditions matching an experiment of Hillesdon et al. (Bull. Math. Biol., vol. 57, 1995, pp. 299-344). We present selected numerical examples, which illustrate (i) the formation of sinking plumes, (ii) the possible merging of neighbouring plumes and (iii) the convergence towards numerically stable stationary plumes. The examples with stable stationary plumes show how the surface-directed oxytaxis continuously feeds cells into a high-concentration layer near the surface, from where the fluid flow (recurring upwards in the space between the plumes) transports the cells into the plumes, where then gravity makes the cells sink and constitutes the driving force in maintaining the fluid convection and, thus, in shaping the plumes into (numerically) stable stationary states. Our numerical method is fully capable of solving the coupled chemotaxis-fluid system and enabling a full exploration of its dynamics, which cannot be done in a linearised framework.

No abstract

The paper is concerned with development of a new finite-volume method for a class of chemotaxis models and for a closely related haptotaxis model. In its simplest form, the chemotaxis model is described by a system of nonlinear PDEs: a convection-diffusion equation for the cell density coupled with a reaction-diffusion equation for the chemoattractant concentration. The first step in the derivation of the new method is made by adding an equation for the chemoattractant concentration gradient to the original system. We then show that the convective part of the resulting system is typically of a mixed hyperbolic-elliptic type and therefore straightforward numerical methods for the studied system may be unstable. The proposed method is based on the application of the second-order central-upwind scheme, originally developed for hyperbolic systems of conservation laws in Kurganov et al. (SIAM J Sci Comput 21:707-740, 2001), to the extended system of PDEs. We show that the proposed second-order scheme is positivity preserving, which is a very important stability property of the method. The scheme is applied to a number of two-dimensional problems including the most commonly used Keller-Segel chemotaxis model and its modern extensions as well as to a haptotaxis system modeling tumor invasion into surrounding healthy tissue. Our numerical results demonstrate high accuracy, stability, and robustness of the proposed scheme.

No abstract

No abstract

Shallow water models are widely used to describe and study free-surface water flow. Even though, in many practical applications the bottom friction does not have much influence on the solutions, the friction terms will play a significant role when the depth of the water is very small. In this paper, we study the Saint-Venant system of shallow water equations with friction terms and develop a well-balanced central-upwind scheme that is capable of exactly preserving its steady states. The scheme also preserves the positivity of the water depth. We test the designed scheme on a number of one-and two-dimensional examples that demonstrate robustness and high resolution of the proposed numerical approach. The data in the last numerical example correspond to the laboratory experiments reported in [L. Cea, M. Garrido, and J. Puertas, J. Hydrol, 382 (2010), pp. 88-102], designed to mimic the rain water drainage in urban areas containing houses. Since the rain water depth is typically several orders of magnitude smaller than the height of the houses, we develop a special technique, which helps to achieve a remarkable agreement between the numerical and experimental results.

We endow the nonlinear degenerate parabolic equation used to describe propagation of thermal waves in plasma or in a porous medium, with a mechanism for flux saturation intended to correct the nonphysical gradientflux relations at high gradients. We study both analytically and numerically the resulting equation: u t = [u n Q(g(u) x )] x , n > 0, where Q is a bounded increasing function. This model reveals that for n > 1 the motion of the front is controlled by the saturation mechanism and instead of the typical infinite gradients resulting from the linear flux-gradients relations, Q ∼ u x , we obtain a sharp, shock-like front, typically associated with nonlinear hyperbolic phenomena. We prove that if the initial support is compact, independently of the smoothness of the initial datum inside the support, a sharp front discontinuity forms in a finite time, and until then the front does not expand. Mathematics Subject Classification: 35K65, 35B65, 35B40