Comparison of two aspects of a PDE model for biological network formation

by   Clarissa Astuto, et al.

We compare the solutions of two systems of partial differential equations (PDE), seen as two different interpretations of the same model that describes formation of complex biological networks. Both approaches take into account the time evolution of the medium flowing through the network, and we compute the solution of an elliptic-parabolic PDE system for the conductivity vector m, the conductivity tensor ℂ and the pressure p. We use finite differences schemes in a uniform Cartesian grid in the spatially two-dimensional setting to solve the two systems, where the parabolic equation is solved by a semi-implicit scheme in time. Since the conductivity vector and tensor appear also in the Poisson equation for the pressure p, the elliptic equation depends implicitly on time. For this reason we compute the solution of three linear systems in the case of the conductivity vector m∈ℝ^2, and four linear systems in the case of the symmetric conductivity tensor ℂ∈ℝ^2× 2, at each time step. To accelerate the simulations, we make use of the Alternating Direction Implicit (ADI) method. The role of the parameters is important for obtaining detailed solutions. We provide numerous tests with various values of the parameters involved, to see the differences in the solutions of the two systems.


page 12

page 13

page 14

page 17

page 20


Asymmetry and condition number of an elliptic-parabolic system for biological network formation

We present results of numerical simulations of the tensor-valued ellipti...

Finite element discretization of a biological network formation system: a preliminary study

A finite element discretization is developed for the Cai-Hu model, descr...

Linearized Implicit Methods Based on a Single-Layer Neural Network: Application to Keller-Segel Models

This paper is concerned with numerical approximation of some two-dimensi...

A Multivariate Spline based Collocation Method for Numerical Solution of Partial Differential Equations

We propose a collocation method based on multivariate polynomial splines...

PARAOPT: A parareal algorithm for optimality systems

The time parallel solution of optimality systems arising in PDE constrai...

Spectrally accurate solutions to inhomogeneous elliptic PDE in smooth geometries using function intension

We present a spectrally accurate embedded boundary method for solving li...

Please sign up or login with your details

Forgot password? Click here to reset