Vivette Girault

Christine Bernardi, Martin Costabel, Monique Dauge et al.

The inf-sup constant for the divergence, or LBB constant, is explicitly known for only few domains. For other domains, upper and lower estimates are known. If more precise values are required, one can try to compute a numerical approximation. This involves, in general, approximation of the domain and then the computation of a discrete LBB constant that can be obtained from the numerical solution of an eigenvalue problem for the Stokes system. This eigenvalue problem does not fall into a class for which standard results about numerical approximations can be applied. Indeed, many reasonable finite element methods do not yield a convergent approximation. In this article, we show that under fairly weak conditions on the approximation of the domain, the LBB constant is an upper semi-continuous shape functional, and we give more restrictive sufficient conditions for its continuity with respect to the domain. For numerical approximations based on variational formulations of the Stokes eigenvalue problem, we also show upper semi-continuity under weak approximation properties, and we give stronger conditions that are sufficient for convergence of the discrete LBB constant towards the continuous LBB constant. Numerical examples show that our conditions are, while not quite optimal, not very far from necessary.

Andrea Bonito, Vivette Girault, Diane Guignard

We study a model describing the slow flow of a fluid through a deformable, porous, elastic solid undergoing small deformations. The stress-strain relationship of the solid incorporates nonlinear effects, formulated as a perturbation of the classical linear elasticity. To approximate the coupled system, we introduce a discrete scheme based on a first order semi-implicit time integration scheme combined with a standard finite element spatial discretization. We establish the existence and uniqueness of the discrete solution and derive a priori convergence estimates under the assumption that the nonlinear perturbations remain sufficiently small. Finally, we demonstrate the efficiency of the proposed scheme through numerical experiments that also highlight the nonlinear phenomena captured by the model.