A fast solver for systems of reaction-diffusion equations.
In this paper we present a fast algorithm for the numerical solution of systems of reaction-diffusion equations, {partial_derivative}{sub t} u + a {center_dot} {del}u = {Delta}u + f(x,t,u), and x element of {Omega} contained in R{sup 3}, t > 0. Here, u is a vector-valued function, u triple bond u(x,t) element of R{sup m} is large, and the corresponding system of ODEs, {partial_derivative}{sub t}u = F(x,t,u), is stiff. Typical examples arise in air pollution studies, where a is the given wind field and the nonlinear function F models the atmospheric chemistry. The time integration of Eq. (1) is best handled by the method of characteristics. The problem is thus reduced to designing for the reaction-diffusion part a fast solver that has good stability properties for the given time step and does not require the computation of the full Jacobi matrix. An operator-splitting technique, even a high-order one, combining a fast nonlinear ODE solver with an efficient solver for the diffusion operator is less effective when the reaction term is stiff. In fact, the classical Strang splitting method may underperform a first-order source splitting method. The algorithm we propose in this paper uses an a posteriori filtering technique to stabilize the computation of the diffusion term. The algorithm parallelizes well, because the solution of the large system of ODEs is done pointwise; however, the integration of the chemistry may lead to load-balancing problems. The Tchebycheff acceleration technique proposed in offers an alternative that complements the approach presented here. To facilitate the presentation, we limit the discussion to domains {Omega} that either admit a regular discretization grid or decompose into subdomains that admit regular discretization grids. We describe the algorithm for one-dimensional domains in Section 2 and for multidimensional domains in Section 3. Section 4 briefly outlines future work.
- Research Organization:
- Argonne National Lab., IL (US)
- Sponsoring Organization:
- US Department of Energy (US)
- DOE Contract Number:
- W-31-109-ENG-38
- OSTI ID:
- 786924
- Report Number(s):
- ANL/MCS/CP-104638; TRN: AH200133%%308
- Resource Relation:
- Conference: 13th International Conference on Domain Decomposition Methods, Lyon (FR), 10/09/2000--10/12/2000; Other Information: PBD: 20 Apr 2001
- Country of Publication:
- United States
- Language:
- English
Similar Records
Multigrid waveform relaxation on spatial finite element meshes
On an initial-boundary value problem for a class of nonlinear Schroedinger equations