Systems of reaction–diffusion partial differential equations (RD-PDEs) are widely applied for modeling life science and physico-chemical phenomena. In particular, the coupling between diffusion and nonlinear kinetics can lead to the so-called Turing instability, giving rise to a variety of spatial patterns (like labyrinths, spots, stripes, etc.) attained as steady state solutions for large time intervals. To capture the morphological peculiarities of the pattern itself, a very fine space discretization may be required, limiting the use of standard (vector-based) ODE solvers in time because of excessive computational costs. By exploiting the structure of the diffusion matrix, we show that matrix-based versions of time integrators, such as Implicit–Explicit (IMEX) and exponential schemes, allow for much finer problem discretizations. We illustrate our findings by numerically solving the Schnakenberg model, prototype of RD-PDE systems with Turing pattern solutions, and the DIB-morphochemical model describing metal growth during battery charging processes.

D'Autilia M.C., Sgura I., Simoncini V. (2020). Matrix-oriented discretization methods for reaction–diffusion PDEs: Comparisons and applications. COMPUTERS & MATHEMATICS WITH APPLICATIONS, 79(7), 2067-2085 [10.1016/j.camwa.2019.10.020].

Matrix-oriented discretization methods for reaction–diffusion PDEs: Comparisons and applications

Simoncini V.
2020

Abstract

Systems of reaction–diffusion partial differential equations (RD-PDEs) are widely applied for modeling life science and physico-chemical phenomena. In particular, the coupling between diffusion and nonlinear kinetics can lead to the so-called Turing instability, giving rise to a variety of spatial patterns (like labyrinths, spots, stripes, etc.) attained as steady state solutions for large time intervals. To capture the morphological peculiarities of the pattern itself, a very fine space discretization may be required, limiting the use of standard (vector-based) ODE solvers in time because of excessive computational costs. By exploiting the structure of the diffusion matrix, we show that matrix-based versions of time integrators, such as Implicit–Explicit (IMEX) and exponential schemes, allow for much finer problem discretizations. We illustrate our findings by numerically solving the Schnakenberg model, prototype of RD-PDE systems with Turing pattern solutions, and the DIB-morphochemical model describing metal growth during battery charging processes.
2020
D'Autilia M.C., Sgura I., Simoncini V. (2020). Matrix-oriented discretization methods for reaction–diffusion PDEs: Comparisons and applications. COMPUTERS & MATHEMATICS WITH APPLICATIONS, 79(7), 2067-2085 [10.1016/j.camwa.2019.10.020].
D'Autilia M.C.; Sgura I.; Simoncini V.
File in questo prodotto:
File Dimensione Formato  
DAutilia_Sgura_Simoncini_CAMWA_R2.pdf

accesso aperto

Tipo: Postprint
Licenza: Licenza per Accesso Aperto. Creative Commons Attribuzione - Non commerciale - Non opere derivate (CCBYNCND)
Dimensione 48.12 MB
Formato Adobe PDF
48.12 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/11585/784042
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 18
  • ???jsp.display-item.citation.isi??? 16
social impact