Antonio Márquez

Antonio Márquez, Salim Meddahi, Francisco-Javier Sayas

The aim of this paper is to propose a systematic way to obtain convergent finite element schemes for the Darcy-Stokes flow problem by combining well-known mixed finite elements that are separately convergent for Darcy and Stokes problems. In the approach in which the Darcy problem is set in its natural $\mathbf{H}(\text{div})$ formulation and the Stokes problem is expressed in velocity-pressure form, the transmission condition ensuring global mass conservation becomes essential. As opposed to the strategy that handles weakly this transmission condition through a Lagrange multiplier, we impose here this restriction exactly in the space of global velocity field. Our analysis of the Galerkin discretization of the resulting problem reveals that, if the mixed finite element space used in the Darcy domain admits an $\mathbf{H}(\text{div})$-stable discrete lifting of the normal trace, then it can be combined with any stable Stokes mixed finite element of the same order to deliver a stable global method with quasi-optimal convergence rate. Finally, we present a series of numerical tests confirming our theoretical convergence estimates.

Antonio Márquez, Salim Meddahi, Francisco-Javier Sayas

We propose an efficient iterative method to solve the mixed Stokes-Dracy model for coupling fluid and porous media flow. The weak formulation of this problem leads to a coupled, indefinite, ill-conditioned and symmetric linear system of equations. We apply a decoupled preconditioning technique requiring only good solvers for the local mixed-Darcy and Stokes subproblems. We prove that the method is asymptotically optimal and confirm, with numerical experiments, that the performance of the preconditioners does not deteriorate on arbitrarily fine meshes.

Antonio Márquez, Salim Meddahi, Thanh Tran

The aim of this paper is to analyze a mixed discontinuous Galerkin discretization of the time-harmonic elasticity problem. The symmetry of the Cauchy stress tensor is imposed weakly, as in the traditional dual-mixed setting. We show that the discontinuous Galerkin scheme is well-posed and uniformly stable with respect to the mesh parameter $h$ and the Lamé coefficient $λ$. We also derive optimal a-priori error bounds in the energy norm. Several numerical tests are presented in order to illustrate the performance of the method and confirm the theoretical results.