An efficient computational method to approximate the solution of a general class of nonlinear reaction-diffusion systems in Cartesian grids is presented. The proposed scheme uses the Local Discontinuous Galerkin (LDG) method as spatial discretization and the symmetric Strang operator splitting as time marching scheme. As a result, not only diffusion and reaction are decoupled; but the reaction sub-step reduces to a set of completely independent small local nonlinear systems of ordinary differential equations. To reduce their computational cost a Lagrange interpolatory technique is used. Diffusion sub-steps are approximated by a directional splitting method derived from an algebraic decomposition of the discrete Laplacian operator, which amounts to solve sets of completely independent small block tridiagonal linear systems; one set per direction. Stability of the diffusive step is analyzed. Finally the efficiency and accuracy of the proposed scheme is numerically investigated on well known Turing pattern formation models in 2D/3D.

Castillo, P., Gomez, S. (2021). An interpolatory directional splitting-local discontinuous Galerkin method with application to pattern formation in 2D/3D. APPLIED MATHEMATICS AND COMPUTATION, 397(15 May 2021) [10.1016/j.amc.2021.125984].

An interpolatory directional splitting-local discontinuous Galerkin method with application to pattern formation in 2D/3D

Sergio Gomez
2021

Abstract

An efficient computational method to approximate the solution of a general class of nonlinear reaction-diffusion systems in Cartesian grids is presented. The proposed scheme uses the Local Discontinuous Galerkin (LDG) method as spatial discretization and the symmetric Strang operator splitting as time marching scheme. As a result, not only diffusion and reaction are decoupled; but the reaction sub-step reduces to a set of completely independent small local nonlinear systems of ordinary differential equations. To reduce their computational cost a Lagrange interpolatory technique is used. Diffusion sub-steps are approximated by a directional splitting method derived from an algebraic decomposition of the discrete Laplacian operator, which amounts to solve sets of completely independent small block tridiagonal linear systems; one set per direction. Stability of the diffusive step is analyzed. Finally the efficiency and accuracy of the proposed scheme is numerically investigated on well known Turing pattern formation models in 2D/3D.
Articolo in rivista - Articolo scientifico
Discontinuous Galerkin; Nonlinear reaction-diffusion systems; Splitting methods; Turing patterns;
English
2021
397
15 May 2021
125984
none
Castillo, P., Gomez, S. (2021). An interpolatory directional splitting-local discontinuous Galerkin method with application to pattern formation in 2D/3D. APPLIED MATHEMATICS AND COMPUTATION, 397(15 May 2021) [10.1016/j.amc.2021.125984].
File in questo prodotto:
Non ci sono file associati a questo prodotto.

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/10281/442744
Citazioni
  • Scopus 3
  • ???jsp.display-item.citation.isi??? 3
Social impact