Salim Meddahi

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.

Felipe Lepe, Salim Meddahi, David Mora et al.

We introduce a discontinuous Galerkin method for the mixed formulation of the elasticity eigenproblem with reduced symmetry. The analysis of the resulting discrete eigenproblem does not fit in the standard spectral approximation framework since the underlying source operator is not compact and the scheme is nonconforming. We show that the proposed scheme provides a correct approximation of the spectrum and prove asymptotic error estimates for the eigenvalues and the eigenfunctions. Finally, we provide several numerical tests to illustrate the performance of the method and confirm the theoretical results.

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.

Felipe Lepe, Salim Meddahi, David Mora et al.

In this paper we analyze a finite element method for solving a quadratic eigenvalue problem derived from the acoustic vibration problem for a heterogeneous dissipative fluid. The problem is shown to be equivalent to the spectral problem for a noncompact operator and athorough spectral characterization is given. The numerical discretization of the problem is based on Raviart-Thomas finite elements. The method is proved to be free of spurious modes and to converge with optimal order. Finally, we report numerical tests which allow us to assess the performance of the method.

Ana Alonso Rodríguez, Salim Meddahi, Alberto Valli

We introduce and analyze a discontinuous Galerkin FEM/BEM method for a time-harmonic eddy current problem written in terms of the magnetic field. We use nonconforming Nédélec finite elements on a partition of the interior domain coupled with continuous boundary elements on the transmission interface. We prove quasi-optimal error estimates in the energy norm.

Norbert Heuer, Salim Meddahi

We present and analyze a discontinuous variant of the hp-version of the boundary element Galerkin method with quasi-uniform meshes. The model problem is that of the hypersingular integral operator on an (open or closed) polyhedral surface. We prove a quasi-optimal error estimate and conclude convergence orders which are quasi-optimal for the h-version with arbitrary degree and almost quasi-optimal for the p-version. Numerical results underline the theory.

Ana Alonso Rodríguez, Salim Meddahi, Alberto Valli

We introduce and analyze a discontinuous Galerkin method for a time-harmonic eddy current problem formulated in terms of the magnetic field. The scheme is obtained by putting together a DG method for the approximation of the vector field variable representing the magnetic field in the conductor and a DG method for the Laplace equation whose solution is a scalar magnetic potential in the insulator. The transmission conditions linking the two problems are taken into account weakly in the global discontinuous Galerkin scheme. We prove that the numerical method is uniformly stable and obtain quasi-optimal error estimates in the DG-energy norm.

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.