A combined mixed finite element ADI scheme for solving Richards' equation with mixed derivatives on irregular grids

Authors:
A. Arrarás;L. Portero;J. C. Jorge
Affiliations:
Dpto. de Ingeniería Matemática e Informática, Universidad Pública de Navarra, Campus de Arrosadía s/n, 31006, Pamplona (Navarra), Spain;Dpto. de Ingeniería Matemática e Informática, Universidad Pública de Navarra, Campus de Arrosadía s/n, 31006, Pamplona (Navarra), Spain;Dpto. de Ingeniería Matemática e Informática, Universidad Pública de Navarra, Campus de Arrosadía s/n, 31006, Pamplona (Navarra), Spain
Venue:
Applied Numerical Mathematics
Year:
2009

Citing 5
Cited 2

An alternating direction implicit scheme for parabolic equations with mixed derivative and convective terms

Journal of Computational Physics
Mixed Finite Elements for Elliptic Problems with Tensor Coefficients as Cell-Centered Finite Differences

SIAM Journal on Numerical Analysis
Enhanced Cell-Centered Finite Differences for Elliptic Equations on General Geometry

SIAM Journal on Scientific Computing
Stability results for linearly implicit fractional step discretizations of non-linear time dependent parabolic problems

Applied Numerical Mathematics
Stability of ADI schemes applied to convection-diffusion equations with mixed derivative terms

Applied Numerical Mathematics

Contractivity of domain decomposition splitting methods for nonlinear parabolic problems

Journal of Computational and Applied Mathematics
Variable step-size fractional step Runge-Kutta methods for time-dependent partial differential equations

Applied Numerical Mathematics

Quantified Score

Hi-index	0.00

Visualization

Abstract

We propose and analyze an efficient numerical method for solving semilinear parabolic problems with mixed derivative terms on non-rectangular domains. The spatial semidiscretization process is based on an expanded mixed finite element scheme which, combined with suitable quadrature rules, is converted into a cell-centered finite difference scheme. This choice preserves the asymptotic accuracy and local conservation of mass of the method, while substantially reducing the computational cost of the totally discrete scheme. To obtain it, an alternating direction implicit scheme is used for the integration in time. The resulting numerical algorithm involves sets of uncoupled tridiagonal systems which can be solved in parallel. We set out some theoretical results of unconditional convergence (of second order in space and first order in time) for our method. Finally, a numerical experiment is shown in order to illustrate the theoretical results.