A filtering technique for system of reaction-diffusion equations

F. Dupros, M. Garbey, W. E. Fitzgibbon

Research output: Contribution to journalArticle

3 Scopus citations

Abstract

In this paper, we present a fast parallel solver designed for a system of reaction-convection-diffusion equations. A typical application is the large-scale computing of air quality models for which the main solver corresponds to reaction-diffusion-convection equations. Another potential application is the numerical simulation of population models where several colonies compete. Reaction-diffusion systems can be integrated in time by pointwise Newton iteration when all space-dependent terms are explicit in the time integration. A Newton-Jacobi iteration makes this scheme implicit. Such methods are easy to code and have scalable parallelism. However, they are numerically inefficient. An alternative method is to use operator splitting, decoupling the time integration of the reaction from the time integration of convection-diffusion. However, such methods may not be time accurate thanks to the stiffness of the reaction term and are complex to parallelize with good scalability. A second alternative is to use matrix-free Newton-Krylov methods. These techniques are particularly efficient provided that a good parallel preconditioner is customized to the application. The method is then not trivial to implement. We propose here a new family of fast, easy to code and numerically efficient reaction-diffusion solvers based on a filtering technique that stabilizes the explicit treatment of the diffusion terms. The scheme is completely explicit with respect to space, and the postprocessing to stabilize time stepping uses a simple FFT. We demonstrate the potential of this numerical scheme with two examples in air quality models that usually require the implicit treatment of diffusion terms and have compared our solution to classical schemes for two nonlinear reaction-diffusion problems. For general reaction-diffusion problems on tensorial product grids with regular space step, the filtering process can be applied as a black box postprocessing procedure. Further, we demonstrate on critical components of the algorithm the high potential of parallelism of our method on medium-scale parallel computers.

Original languageEnglish (US)
Pages (from-to)1-29
Number of pages29
JournalInternational Journal for Numerical Methods in Fluids
Volume52
Issue number1
DOIs
StatePublished - Sep 10 2006

Keywords

  • Air quality
  • Filtering
  • ODE
  • Reaction-diffusion

ASJC Scopus subject areas

  • Computational Mechanics
  • Mechanics of Materials
  • Mechanical Engineering
  • Computer Science Applications
  • Applied Mathematics

Fingerprint Dive into the research topics of 'A filtering technique for system of reaction-diffusion equations'. Together they form a unique fingerprint.

Cite this