In this work, we study a finite element method in order to simulate the motion of immersed rigid bodies. This method is of the fictitious domain type. The idea is to look for a smooth extension in the whole domain of the exact solution and to recover the optimal order obtain with a conformal mesh. This smooth extension is sought by minimizing a functional whose gradient is the solution of another fluid problem with a single layer distribution as a right hand side. We make the numerical analysis, in the scalar case, of the approximation of this distribution by a sum of Dirac masses. One of the advantage of this method is to be able to use fast solvers on cartesian mesh while recovering the optimal order of the error. Another advantage of this method is that the operators are not modified at all. Only the right hand side depends on the geometry of the original problem. We write a parallel C++ code in two and three dimensions that simulate fluid/rigid bodies flows with this method. We present the core blocks of this code to show how it works.
