On iterative subdomain methods for the Stokes--Darcy problem

Iterative subdomain methods for the Stokes--Darcy problem that use Robin boundary conditions on the interface are reviewed. Their common underlying structure and their main differences are identified. In particular, it is clarified that there are different updating strategies for the interface conditions. For small values of fluid viscosity and hydraulic permeability, which are relevant in applications from geosciences, it is shown in numerical studies that only one of these updating strategies leads to an efficient numerical method, if this strategy is used in combination with appropriate parameters in the Robin boundary conditions. In particular, it is observed that the values of appropriate parameters are larger than those proposed so far. Not only the size but also the ratio of appropriate Robin parameters depends on the coefficients of the problem

A priori error analysis of a fully-mixed finite element method for a two-dimensional fluid-solid interaction problem

We introduce and analyze a fully-mixed finite element method for a fluid-solid interaction problem in 2D. The model consists of an elastic body which is subject to a given incident wave that travels in the fluid surrounding it. Actually, the fluid is supposed to occupy an annular region, and hence a Robin boundary condition imitating the behavior of the scattered field at infinity is imposed on its exterior boundary, which is located far from the obstacle. The media are governed by the elastodynamic and acoustic equations in time-harmonic regime, respectively, and the transmission conditions are given by the equilibrium of forces and the equality of the corresponding normal displacements. We first apply dual-mixed approaches in both domains, and then employ the governing equations to eliminate the displacement u of the solid and the pressure p of the fluid. In addition, since both transmission conditions become essential, they are enforced weakly by means of two suitable Lagrange multipliers. As a consequence, the Cauchy stress tensor and the rotation of the solid, together with the gradient of p and the traces of u and p on the boundary of the fluid, constitute the unknowns of the coupled problem. Next, we show that suitable decompositions of the spaces to which the stress and the gradient of p belong, allow the application of the Babuška–Brezzi theory and the Fredholm alternative for analyzing the solvability of the resulting continuous formulation. The unknowns of the solid and the fluid are then approximated by a conforming Galerkin scheme defined in terms of PEERS elements in the solid, Raviart–Thomas of lowest order in the fluid, and continuous piecewise linear functions on the boundary. Then, the analysis of the discrete method relies on a stable decomposition of the corresponding finite element spaces and also on a classical result on projection methods for Fredholm operators of index zero. Finally, some numerical results illustrating the theory are presented