Solving Partial Differential Algebraic Equations and Reactive Transport Models
Autor: | Jocelyne Erhel, Souhila Sabit, Caroline de Dieuleveult |
---|---|
Přispěvatelé: | Simulations and Algorithms on Grids for Environment (SAGE), Inria Rennes – Bretagne Atlantique, Institut National de Recherche en Informatique et en Automatique (Inria)-Institut National de Recherche en Informatique et en Automatique (Inria)-SYSTÈMES LARGE ÉCHELLE (IRISA-D1), Institut de Recherche en Informatique et Systèmes Aléatoires (IRISA), CentraleSupélec-Télécom Bretagne-Université de Rennes 1 (UR1), Université de Rennes (UNIV-RENNES)-Université de Rennes (UNIV-RENNES)-Institut National de Recherche en Informatique et en Automatique (Inria)-École normale supérieure - Rennes (ENS Rennes)-Université de Bretagne Sud (UBS)-Centre National de la Recherche Scientifique (CNRS)-Institut National des Sciences Appliquées - Rennes (INSA Rennes), Institut National des Sciences Appliquées (INSA)-Université de Rennes (UNIV-RENNES)-Institut National des Sciences Appliquées (INSA)-CentraleSupélec-Télécom Bretagne-Université de Rennes 1 (UR1), Institut National des Sciences Appliquées (INSA)-Université de Rennes (UNIV-RENNES)-Institut National des Sciences Appliquées (INSA)-Institut de Recherche en Informatique et Systèmes Aléatoires (IRISA), Université de Rennes (UNIV-RENNES)-Université de Rennes (UNIV-RENNES)-École normale supérieure - Rennes (ENS Rennes)-Université de Bretagne Sud (UBS)-Centre National de la Recherche Scientifique (CNRS)-Institut National des Sciences Appliquées - Rennes (INSA Rennes), Institut National des Sciences Appliquées (INSA)-Université de Rennes (UNIV-RENNES)-Institut National des Sciences Appliquées (INSA), Centre de Géosciences (GEOSCIENCES), MINES ParisTech - École nationale supérieure des mines de Paris, Université Paris sciences et lettres (PSL)-Université Paris sciences et lettres (PSL), ANDRA, ANR-12-MONU-0012,H2MNO4,Hydrogéologie Hétérogène avec un Modèle Numérique Original, Optimisé et Orienté Objets(2012), Université de Rennes (UR)-Institut National des Sciences Appliquées - Rennes (INSA Rennes), Institut National des Sciences Appliquées (INSA)-Institut National des Sciences Appliquées (INSA)-Université de Bretagne Sud (UBS)-École normale supérieure - Rennes (ENS Rennes)-Institut National de Recherche en Informatique et en Automatique (Inria)-Télécom Bretagne-CentraleSupélec-Centre National de la Recherche Scientifique (CNRS)-Université de Rennes (UR)-Institut National des Sciences Appliquées - Rennes (INSA Rennes), Institut National des Sciences Appliquées (INSA)-Institut National des Sciences Appliquées (INSA)-Université de Bretagne Sud (UBS)-École normale supérieure - Rennes (ENS Rennes)-Institut National de Recherche en Informatique et en Automatique (Inria)-Télécom Bretagne-CentraleSupélec-Centre National de la Recherche Scientifique (CNRS)-Institut de Recherche en Informatique et Systèmes Aléatoires (IRISA), Institut National des Sciences Appliquées (INSA)-Institut National des Sciences Appliquées (INSA)-Université de Bretagne Sud (UBS)-École normale supérieure - Rennes (ENS Rennes)-Télécom Bretagne-CentraleSupélec-Centre National de la Recherche Scientifique (CNRS), Mines Paris - PSL (École nationale supérieure des mines de Paris), contrat avec l'ANDRA, Topping, B.H.V. and Iványi, P. |
Jazyk: | angličtina |
Rok vydání: | 2013 |
Předmět: |
Mathematical optimization
Discretization Computer science [SDE.MCG]Environmental Sciences/Global Changes 0207 environmental engineering 010103 numerical & computational mathematics 02 engineering and technology 01 natural sciences L-stability symbols.namesake Matrix (mathematics) Multigrid method Simultaneous equations Distributed parameter system [SDU.STU.GC]Sciences of the Universe [physics]/Earth Sciences/Geochemistry Applied mathematics 0101 mathematics Algebraic number 020701 environmental engineering Newton's method ComputingMilieux_MISCELLANEOUS Partial differential equation Flux-corrected transport Independent equation Method of lines Delay differential equation Nonlinear systems of equations Nonlinear system Algebraic equation symbols Partial derivative [INFO.INFO-DC]Computer Science [cs]/Distributed Parallel and Cluster Computing [cs.DC] Differential algebraic geometry Differential algebraic equation [MATH.MATH-NA]Mathematics [math]/Numerical Analysis [math.NA] Numerical partial differential equations |
Zdroj: | The Third International Conference on Parallel, Distributed, Grid and Cloud Computing for Engineering The Third International Conference on Parallel, Distributed, Grid and Cloud Computing for Engineering, May 2013, PECS, Hungary Developments in Parallel, Distributed, Grid and Cloud Computing for Engineering Topping, B.H.V. and Iványi, P. Developments in Parallel, Distributed, Grid and Cloud Computing for Engineering, 31, Saxe-Coburg publications, pp.151-169, 2013, Computational Science, Engineering and Technology Series, ⟨10.4203/csets.31.7⟩ |
DOI: | 10.4203/csets.31.7⟩ |
Popis: | International audience; 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. |
Databáze: | OpenAIRE |
Externí odkaz: |