Solving Partial Differential Algebraic Equations and Reactive Transport Models
2013
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.
Keywords:
- Simultaneous equations
- Differential algebraic equation
- Mathematical optimization
- Flux-corrected transport
- Method of lines
- Computer science
- Multigrid method
- Nonlinear system
- Differential algebraic geometry
- Numerical partial differential equations
- Applied mathematics
- Delay differential equation
- Distributed computing
- Independent equation
- Correction
- Source
- Cite
- Save
- Machine Reading By IdeaReader
39
References
3
Citations
NaN
KQI