Fully implicit local time-stepping methods for advection-diffusion problems in mixed formulations

DOI10.1016/J.CAMWA.2022.05.022MaRDI QIDQ6379234arXivFDO

Authors Thi-Thao-Phuong Hoang

Publication date 2 October 2021

Published in arXiv

Full work available at URL https://doi.org/10.1016/J.CAMWA.2022.05.022

Mathematics Subject Classification ID

PDEs in connection with fluid mechanics (35Q35) PDEs in connection with geophysics (35Q86) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Multigrid methods; domain decomposition for initial value and initial-boundary value problems involving PDEs (65M55) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Multigrid methods; domain decomposition for boundary value problems involving PDEs (65N55) Finite element methods applied to problems in fluid mechanics (76M10) Flows in porous media; filtration; seepage (76S05) Hydrology, hydrography, oceanography (86A05)

arXiv classification

cs.NAmath.NA

Abstract: This paper is concerned with numerical solution of transport problems in heterogeneous porous media. A semi-discrete continuous-in-time formulation of the linear advection-diffusion equation is obtained by using a mixed hybrid finite element method, in which the flux variable represents both the advective and diffusive flux, and the Lagrange multiplier arising from the hybridization is used for the discretization of the advective term. Based on global-in-time and nonoverlapping domain decomposition, we propose two implicit local time-stepping methods to solve the semi-discrete problem. The first method uses the time-dependent Steklov-Poincar'e type operator and the second uses the optimized Schwarz waveform relaxation (OSWR) with Robin transmission conditions. For each method, we formulate a space-time interface problem which is solved iteratively. Each iteration involves solving the subdomain problems independently and globally in time; thus, different time steps can be used in the subdomains. The convergence of the fully discrete OSWR algorithm with nonmatching time grids is proved. Numerical results for problems with various Pecl'et numbers and discontinuous coefficients, including a prototype for the simulation of the underground storage of nuclear waste, are presented to illustrate the performance of the proposed local time-stepping methods.

This page was built for publication: Fully implicit local time-stepping methods for advection-diffusion problems in mixed formulations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6379234)