We present results of numerical simulations of the tensor-valued elliptic-parabolic PDE model for biological network formation. The numerical method is based on a nonlinear finite difference scheme on a uniform Cartesian grid in a two-dimensional (2D) domain. The focus is on the impact of different discretization methods and choices of regularization parameters on the symmetry of the numerical solution. In particular, we show that using the symmetric alternating direction implicit (ADI) method for time discretization helps preserve the symmetry of the solution, compared to the (non-symmetric) ADI method. Moreover, we study the effect of the regularization by the isotropic background permeability r> 0 , showing that the increased condition number of the elliptic problem due to decreasing value of r leads to loss of symmetry. We show that in this case, neither the use of the symmetric ADI method preserves the symmetry of the solution. Finally, we perform the numerical error analysis of our method making use of the Wasserstein distance.

Asymmetry and Condition Number of an Elliptic-Parabolic System for Biological Network Formation

Astuto C.
Investigation
;
Russo G.
Supervision
2023-01-01

Abstract

We present results of numerical simulations of the tensor-valued elliptic-parabolic PDE model for biological network formation. The numerical method is based on a nonlinear finite difference scheme on a uniform Cartesian grid in a two-dimensional (2D) domain. The focus is on the impact of different discretization methods and choices of regularization parameters on the symmetry of the numerical solution. In particular, we show that using the symmetric alternating direction implicit (ADI) method for time discretization helps preserve the symmetry of the solution, compared to the (non-symmetric) ADI method. Moreover, we study the effect of the regularization by the isotropic background permeability r> 0 , showing that the increased condition number of the elliptic problem due to decreasing value of r leads to loss of symmetry. We show that in this case, neither the use of the symmetric ADI method preserves the symmetry of the solution. Finally, we perform the numerical error analysis of our method making use of the Wasserstein distance.
2023
Asymmetry
Bionetwork formation
Cai-Hu model
Conditioning number
Finite-difference scheme
Leaf venation
Semi-implicit
Symmetric alternating direction implicit (ADI)
Wasserstein distance
File in questo prodotto:
File Dimensione Formato  
asymmetry.pdf

solo gestori archivio

Tipologia: Versione Editoriale (PDF)
Licenza: NON PUBBLICO - Accesso privato/ristretto
Dimensione 1.94 MB
Formato Adobe PDF
1.94 MB Adobe PDF   Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.11769/608051
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 0
  • ???jsp.display-item.citation.isi??? 1
social impact