Solving Partial Differential Algebraic Equations and Reactive Transport Models

Abstract : In some scientific applications, such as groundwater studies, several processes are represented by coupled models. For example, a density-driven flow model couples the flow equations with the transport of salt. A reactive transport model couples transport equations of pollutants with chemical equations. The coupled model can combine partial differential equations with algebraic equations, in a so-called PDAE system, which is in general nonlinear. A classical approach is to follow a method of lines, where space is first discretized, leading to a semi-discrete differential algebraic system (DAE). Then time is discretized by a scheme tuned for DAE, such that each time step requires solving a nonlinear system of equations. In some decoupled approaches, a fixed-point technique is used. However, a Newton method converges faster in general and is more efficient, even though each iteration is more CPU-intensive. In this talk, we deal with reactive transport models and show how a Newton method can be used efficiently. Numerical experiments illustrate the efficiency of a substitution technique. Moreover, it appears that using logarithms in the chemistry equations lead to ill conditioned matrices and increase the computational cost.
Type de document :
Communication dans un congrès
The Third International Conference on Parallel, Distributed, Grid and Cloud Computing for Engineering, May 2013, PECS, Hungary. 2013
Liste complète des métadonnées

https://hal.inria.fr/hal-00826660
Contributeur : Jocelyne Erhel <>
Soumis le : mardi 28 mai 2013 - 10:15:04
Dernière modification le : mardi 16 janvier 2018 - 15:54:11

Identifiants

  • HAL Id : hal-00826660, version 1

Citation

Jocelyne Erhel, Souhila Sabit, Caroline De Dieuleveult. Solving Partial Differential Algebraic Equations and Reactive Transport Models. The Third International Conference on Parallel, Distributed, Grid and Cloud Computing for Engineering, May 2013, PECS, Hungary. 2013. 〈hal-00826660〉

Partager

Métriques

Consultations de la notice

447