Michel Fournié

Sébastien Court, Michel Fournié, Alexei Lozinski

In the present work, we propose to extend to the Stokes problem a fictitious domain approach inspired by eXtended Finite Element Method and studied for Poisson problem in [Renard]. The method allows computations in domains whose boundaries do not match. A mixed finite element method is used for fluid flow. The interface between the fluid and the structure is localized by a level-set function. Dirichlet boundary conditions are taken into account using Lagrange multiplier. A stabilization term is introduced to improve the approximation of the normal trace of the Cauchy stress tensor at the interface and avoid the inf-sup condition between the spaces for velocity and the Lagrange multiplier. Convergence analysis is given and several numerical tests are performed to illustrate the capabilities of the method.

Michel Fournié, Alexei Lozinski

We study a fictitious domain approach with Lagrange multipliers to discretize Stokes equations on a mesh that does not fit the boundaries. A mixed finite element method is used for fluid flow. Several stabilization terms are added to improve the approximation of the normal trace of the stress tensor and to avoid the inf-sup conditions between the spaces of the velocity and the Lagrange multipliers. We generalize first an approach based on eXtended Finite Element Method due to Haslinger-Renard involving a Barbosa-Hughes stabilization and a robust reconstruction on the badly cut elements. Secondly, we adapt the approach due to Burman-Hansbo involving a stabilization only on the Lagrange multiplier. Multiple choices for the finite elements for velocity, pressure and multiplier are considered. Additional stabilization on pressure (Brezzi-Pitkäranta, Interior Penalty) is added, if needed. We prove the stability and the optimal convergence of several variants of these methods under appropriate assumptions. Finally, we perform numerical tests to illustrate the capabilities of the methods.

Sébastien Court, Michel Fournié

The paper extends a stabilized fictitious domain finite element method initially developed for the Stokes problem to the incompressible Navier-Stokes equations coupled with a moving solid. This method presents the advantage to predict an optimal approximation of the normal stress tensor at the interface. The dynamics of the solid is governed by the Newton's laws and the interface between the fluid and the structure is materialized by a level-set which cuts the elements of the mesh. An algorithm is proposed in order to treat the time evolution of the geometry and numerical results are presented on a classical benchmark of the motion of a disk falling in a channel.

Bertram Düring, Michel Fournié, Alain Rigal

We present new high-order Alternating Direction Implicit (ADI) schemes for the numerical solution of initial-boundary value problems for convection-diffusion equations with mixed derivative terms. Our approach is based on the unconditionally stable ADI scheme proposed by Hundsdorfer. Different numerical discretizations which lead to schemes which are fourth-order accurate in space and second-order accurate in time are discussed.