The purpose of this paper is twofold. First, we extend the applicability of Cattaneo's relaxation approach, one of the currently known relaxation approaches, to reformulate timedependent advection-diffusion-reaction equations, which may include stiff reactive terms, as hyperbolic balance laws with stiff source terms. We show that the approach works well for second-order partial differential equations. However, it does not for the higher-order case, as proved here for third-order partial differential equations, such as equations of the Korteweg-deVries-Burgers type. In this paper we restrict ourselves to selected second-order model equations. Second, we extend the applicability of existent high-order schemes to approximating this class of partial differential equations numerically. We start from linear problems, proceeding to nonlinear problems, and ending up with an application to an atherosclerosis model consisting of a system of time-dependent diffusionreaction equations. Several scalar examples are considered, and exact solutions are given for some of these. The proposed hyperbolization procedures turn out to give a generous stability range for the choice of the time step. This results in more efficient schemes than some methods with the parabolic restriction reported in the current literature. Implementations of our numerical schemes and convergence-rate assessments are carried out for methods of up to seventh order of accuracy in both space and time. © 2014 Society for Industrial and Applied Mathematics.

Advection-diffusion-reaction equations: Hyperbolization and high-order ader discretizations

Toro, Eleuterio Francisco;Montecinos Guzman, Gino Ignacio
2014

Abstract

The purpose of this paper is twofold. First, we extend the applicability of Cattaneo's relaxation approach, one of the currently known relaxation approaches, to reformulate timedependent advection-diffusion-reaction equations, which may include stiff reactive terms, as hyperbolic balance laws with stiff source terms. We show that the approach works well for second-order partial differential equations. However, it does not for the higher-order case, as proved here for third-order partial differential equations, such as equations of the Korteweg-deVries-Burgers type. In this paper we restrict ourselves to selected second-order model equations. Second, we extend the applicability of existent high-order schemes to approximating this class of partial differential equations numerically. We start from linear problems, proceeding to nonlinear problems, and ending up with an application to an atherosclerosis model consisting of a system of time-dependent diffusionreaction equations. Several scalar examples are considered, and exact solutions are given for some of these. The proposed hyperbolization procedures turn out to give a generous stability range for the choice of the time step. This results in more efficient schemes than some methods with the parabolic restriction reported in the current literature. Implementations of our numerical schemes and convergence-rate assessments are carried out for methods of up to seventh order of accuracy in both space and time. © 2014 Society for Industrial and Applied Mathematics.
5
Toro, Eleuterio Francisco; Montecinos Guzman, Gino Ignacio
File in questo prodotto:
File Dimensione Formato  
Toro VQR 3.pdf

Solo gestori archivio

Tipologia: Versione editoriale (Publisher’s layout)
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 568.82 kB
Formato Adobe PDF
568.82 kB 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: http://hdl.handle.net/11572/125876
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 42
  • ???jsp.display-item.citation.isi??? 33
social impact