L'Interaction Fluide-Structure (IFS) décrit une classe très générale de problème physique, ce qui explique la nécessité de développer une méthode numérique capable de simuler le problème FSI. Pour cette raison, un solveur IFS est développé qui peut traiter un écoulement de fluide incompressible en interaction avec des structures différente: élastique ou rigide. Dans cet aspect, le solveur peut couvrir une large gamme d'applications.La méthode proposée est développée dans le cadre d'une formulation monolithique dans un contexte Eulérien. Cette méthode consiste à considérer un seul maillage et résoudre un seul système d'équations avec des propriétés matérielles différentes. La fonction distance permet de définir la position et l'interface de tous les objets à l'intérieur du domaine et de fournir les propriétés physiques pour chaque sous-domaine. L'adaptation de maillage anisotrope basé sur la variation de la fonction distance est ensuite appliquée pour assurer une capture précise des discontinuités à l'interface fluide-solide.La formulation monolithique est assurée par l'ajout d'un tenseur supplémentaire dans les équations de Navier-Stokes. Ce tenseur provient de la présence de la structure dans le fluide. Le système est résolu en utilisant une méthode élément fini et stabilisé suivant la formulation variationnelle multiéchelle. Cette formulation consiste à décomposer les champs de vitesse et pression en grande et petite échelles. La particularité de l'approche proposée réside dans l'enrichissement du tenseur de l'extra contraint.La première application est la simulation IFS avec un corps rigide. Le corps rigide est décrit en imposant une valeur nul du tenseur des déformations, et le mouvement est obtenu par la résolution du mouvement de corps rigide. Nous évaluons le comportement et la précision de la formulation proposée dans la simulation des exemples 2D et 3D. Les résultats sont comparés avec la littérature et montrent que la méthode développée est stable et précise.La seconde application est la simulation IFS avec un corps élastique. Dans ce cas, une équation supplémentaire est ajoutée au système précédent qui permet de résoudre le champ de déplacement. Et la contrainte de rigidité est remplacée par la loi de comportement du corps élastique. La déformation et le mouvement du corps élastique sont réalisés en résolvant l'équation de convection de la Level-Set. Nous illustrons la flexibilité de la formulation proposée par des exemples 2D. / Numerical simulations of fluid-structure interaction (FSI) are of first interest in numerous industrial problems: aeronautics, heat treatments, aerodynamic, bioengineering... Because of the high complexity of such problems, analytical study is in general not sufficient to understand and solve them. FSI simulations are then nowadays the focus of numerous investigations, and various approaches are proposed to treat them. We propose in this thesis a novel monolithic approach to deal with the interaction between an incompressible fluid flow and rigid/ elastic material. This method consists in considering a single grid and solving one set of equations with different material properties. A distance function enables to define the position and the interface of any objects with complex shapes inside the volume and to provide heterogeneous physical properties for each subdomain. Different anisotropic mesh adaptation algorithms based on the variations of the distance function or on using error estimators are used to ensure an accurate capture of the discontinuities at the fluid-solid interface. The monolithic formulation is insured by adding an extra-stress tensor in the Navier-Stokes equations coming from the presence of the structure in the fluid. The system is then solved using a finite element Variational MultiScale (VMS) method, which consists of decomposition, for both the velocity and the pressure fields, into coarse/resolved scales and fine/unresolved scales. The distinctive feature of the proposed approach resides in the efficient enrichment of the extra constraint. In the first part of the thesis, we use the proposed approach to assess its accuracy and ability to deal with fluid-rigid interaction. The rigid body is prescribed under the constraint of imposing the nullity of the strain tensor, and its movement is achieved by solving the rigid body motion. Several test case, in 2D and 3D with simple and complex geometries are presented. Results are compared with existing ones in the literature showing good stability and accuracy on unstructured and adapted meshes. In the second, we present different routes and an extension of the approach to deal with elastic body. In this case, an additional equation is added to the previous system to solve the displacement field. And the rigidity constraint is replaced with a corresponding behaviour law of the material. The elastic deformation and motion are captured using a convected level-set method. We present several 2D numerical tests, which is considered as classical benchmarks in the literature, and discuss their results.

The finite element method (FEM) is a widely-used numerical tool in the fields of structural dynamics, acoustics and electromagnetics. In this work, our goal is to develop robust FEM strategies for solving problems in the areas of acoustics, structures and electromagnetics, and then extend these strategies to solve multi-physics problems such as magnetohydrodynamics and structural acoustics. We now briefly describe the finite element strategies developed in each of the above domains. In the structural domain, we show that the trapezoidal rule, which is a special case of the Newmark family of algorithms, conserves linear and angular momenta and energy in the case of undamped linear elastodynamics problems, and an ‘energy-like measure’ in the case of undamped acoustic problems. These conservation properties, thus, provide a rational basis for using this algorithm. In linear elastodynamics variants of the trapezoidal rule that incorporate ‘high-frequency’ dissipation are often used, since the higher frequencies, which are not approximated properly by the standard displacement-based approach, often result in unphysical behavior. Instead of modifying the trapezoidal algorithm, we propose using a hybrid FEM framework for constructing the stiffness matrix. Hybrid finite elements, which are based on a two-field variational formulation involving displacement and stresses, are known to approximate the eigenvalues much more accurately than the standard displacement-based approach, thereby either bypassing or reducing the need for high-frequency dissipation. We show this by means of several examples, where we compare the numerical solutions obtained using the displacementbased and hybrid approaches against analytical solutions. We also present a monolithic formulation for the solution of structural acoustic problems based on the hybrid finite element approach. In the area of electromagnetics, since our goal is to ultimately couple the electromagnetic analysis with structural or fluid variables in a ‘monolithic’ framework, we focus on developing nodal finite elements rather than using ‘edge elements’. It is well-known that conventional nodal finite elements can give rise to spurious solutions, and that they cannot capture singularities when the domains are nonconvex and have sharp corners. The commonly used remedies of either adding a penalty term or using a potential formulation are unable to address these problems satisfactorily. In order to overcome this problem, we first develop several mixed finite elements in two and three dimensions which predict the eigenfrequencies (including their multiplicities) accurately, even for non-convex domains. In this proposed formulation, no ad-hoc terms are added as in the penalty formulation, and the improvement is achieved purely by an appropriate choice of the finite element spaces for the different variables. For inhomogeneous domains, ‘double noding’ is used to enforce the appropriate continuity conditions at an interface. Although the developed mixed FEM works very accurately for all 2D geometries and regular Cartesian 3D geometries, it has so far not yielded success for curved 3D geometries. Therefore, for 3D harmonic and transient analysis problems, we propose and use a modified form of the potential formulation that overcomes the disadvantages of the standard potential method, especially on non-convex domains. Electromagnetic radiation and scattering in an exterior domain traditionally involved imposing a suitable absorbing boundary condition (ABC) on the truncation boundary of the numerical domain to inhibit reflection from it. In this work, based on the Wilcox asymptotic expansion of the electric far-field, we propose an amplitude formulation within the framework of the nodal FEM, whereby the highly oscillatory radial part of the field is separated out a-priori so that the standard Lagrange interpolation functions have to capture a relatively gently varying function. Since these elements can be used in the immediate vicinity of the radiator or scatterer (with few exceptions which we enumerate), it is more effective compared to methods of imposing ABCs, especially for high-frequency problems. We show the effectiveness of the proposed formulation on a wide variety of radiation and scattering problems involving both conducting and dielectric bodies, and involving both convex and non-convex domains with sharp corners. The Time Domain Finite Element Method (TDFEM) has been used extensively to solve transient electromagnetic radiation and scattering problems. Although conservation of energy in electromagnetics is well-known, we show in this work that there are additional quantities that are also conserved in the absence of loading. We then show that the developed time-stepping strategy (which is closely related to the trapezoidal rule) mimics these continuum conservation properties either exactly or to a very good approximation. Thus, the developed numerical strategy can be said to be ‘unconditionally stable’ (from an energy perspective) allowing the use of arbitrarily large time-steps. We demonstrate the high accuracy and robustness of the developed method for solving both interior and exterior domain radiation problems, and for finding the scattered field from conducting and dielectric bodies. In the field of magneto-hydrodynamics, we develop a monolithic strategy based on a continuous velocity-pressure formulation that is known to satisfy the Babuska-Brezzi (BB) conditions. The magnetic field is interpolated in the same way as the velocity field, and the entire formulation is within a nodal finite element framework. Both transient and steady-state formulations are developed for two- and three-dimensional geometries. An exact linearization of the monolithic strategy ensures that rapid (quadratic) convergence is achieved within each time (or load) step, while the stable nature of the interpolations used ensure that no instabilities arise in the solution. Good agreement with analytical solutions, even with the use of very coarse meshes, shows the efficacy of the developed formulation.