Dans ce manuscrit, nous présentons une méthode d’optimisation de forme qui se base sur des paramètres géométriques comme des longueurs, des angles, etc. Nous nous appuyons sur des techniques d’optimisation basées sur un gradient. La sensibilité de la fonction objectif par rapport à la position des noeuds du maillage nous est fournie par un solveur adjoint que l’on considère comme une boîte noire. Afin d’optimiser par rapport aux paramètres CAO, nous nous concentrons sur l’évaluation de la sensibilité de la position des noeuds par rapport à ces paramètres. Ainsi, nous proposons deux approches par différences finies. La première méthode s’appuie sur une projection harmonique afin de comparer dans un même espace le maillage initial et celui obtenu suite à la variation d’un paramètre CAO. Les développements présentés dans ce manuscrit permettent d’étendre l’application aux formes ayant plusieurs frontières comme les collecteurs d’échappement. Nous avons développé une méthode d’interpolation adaptée à cette comparaison. L’ensemble du processus a été automatisé et nous en montrons l’entière efficacité sur des applications industrielles en aérodynamique interne. La deuxième méthode se base directement sur les géométries CAO pour évaluer cette sensibilité. Nous utilisons la définition intrinsèque des patches dans l’espace paramétrique (u;v) pour effectuer cette comparaison. Grâce à l’utilisation des coordonnées exactes en tout point de la surface fournies par la CAO, nous évitons d’avoir recours à une interpolation afin d’avoir la meilleure précision de calcul possible. Cependant, contrairement à la première méthode, elle requiert d’identifier les correspondances entre les patches d’une forme à l’autre. Une application sur un cas académique a été faite en aérodynamique externe. La pertinence de la première méthode a été démontrée sur des cas représentatifs et multiobjectifs, ce qui permettrait de faciliter son déploiement et son utilisation dans un cadre industriel. Quant à la deuxième méthode, nous avons montré son fort potentiel. Cependant, des développements supplémentaires seraient nécessaires pour une application plus poussée. Du fait qu’elles sont indépendantes des solveurs mécaniques et du nombre de paramètres, ces méthodes réduisent considérablement les temps de développement des produits, notamment en permettant l’optimisation multiphysique en grande dimension. / In this manuscript, we present a shape optimization method based on CAD parameters such as lengths, angles, etc. We rely on gradient-based optimization techniques. The sensitivity of the objective function, with respect to the mesh nodes position, is provided by an adjoint solver considered here as a black box. To optimize with respect to CAD parameters, we focus on computing the sensitivity of the nodes positions with respect to these parameters. Thus, we propose two approaches based on finite differences. The first method uses a harmonic projection to compare in the same space the initial mesh and the one obtained after a change of the set of CAD parameters. The developments presented in this manuscript open up new doors like the application to shapes with multiple borders such as exhaust manifolds. We also developed an interpolation method suitable for this comparison. The entire process is automated, and we demonstrate the entire effectiveness on internal aerodynamics industrial applications. The second method is directly based on the CAD geometries to assess this sensitivity. To perform this comparison, we use the intrinsic definition of the patches in the parametric space (u;v). Through the use of the exact coordinates at any point on the surface provided by the CAD, we avoid using an interpolation to get the best calculation accuracy possible. However, unlike the first method, it requires to identify the correspondence between patches from one shape to another. An application on an external aerodynamics academic case was made. The relevance of the first method is demonstrated on a representative multi-objective case, which facilitate its deployment use in an industrial environment. Regarding the second method, we showed its great potential. However, further developments are needed to handle more advanced cases. Because they are independent of the mechanical solver and the number of parameters, these methods significantly reduce product development time, particularly by allowing large and multiphysics optimization.

Cette thèse s’inscrit dans le contexte de la simulation et du contrôle des vibrations de structures montées sur ressort qui peuvent apparaître sous l’effet de l’interaction avec l’écoulement de sillage instationnaire. Le contrôle de ce phénomène, appelé vibrations induites par vortex (VIV), est un enjeu critique dans l’optimisation de nombreux systèmes. Une méthode de frontières immergées (IBM) a été intégrée dans l’algorithme PISO du code OpenFOAM, dédié à la simulation d’écoulements fluides incompressibles. La méthode IBM permet une représentation précise de corps fixes ou en mouvement, tout en conservant des maillages structurés conduisant à des algorithmes plus précis et efficaces en termes de performances numériques. Pour calculer la divergence de l’équation de quantité de mouvement dans une boucle PISO et l’interpolation des flux, un calcul hybride orignal a été proposé avec une résolution analytique utilisant l’équation de la fonction noyau des quantités impliquant le terme force de l’IBM (quantités singulières). La méthode a été étendu au formalisme d’écoulements en régimes turbulents. Une loi de paroi a été intégrée permettant de modéliser la couche limite à grand nombre de Reynolds. Le travail de validation a été réalisé au regard des données expérimentales et numériques disponibles dans la littérature pour l’étude d’écoulements autour de cylindres et de sphères, sur une large gamme de nombres de Reynolds. Avec l’objectif de développer des lois de contrôle optimal pour le VIV, basées sur les mécanismes d’instabilité linéaire du système couplé dans le cadre de la théorie du contrôle, un solveur adjoint a été développé et validé. / This thesis is related to the simulation and the control of the vortex induced vibrations phenomenon (VIV), which can result from the fluid structure interactions between an unsteady wake and the body, when the shedding frequency in the wake is close to the natural frequency of the body. The control of VIV is a critical issue when optimizing many systems. An Immersed Boundaries Method (IBM) was implemented into the PISO algorithm as a new library of OpenFOAM, in order to perform reliable simulations of incompressible flows around bluff bodies.To compute the divergence of the momentum equation and the interpolation of the fluxes, an hybrid calculation with an analytical resolution of the quantities involving the force term (singular quantities) has been proposed. The mesh convergence of several errors was shown by means of a manufactured solution, allowing to analyze both the errors irelated to the discretization and to the IBM. The new algorithm was subsequently extended to the RANS and DDES formalism proposed in OpenFOAM for the simulation of turbulent flows. A wall law was integrated into theIBM method to model the boundary layers that develop around the bodies at large Reynolds numbers. Various 2D and 3D well-documented test cases of academic flows around fixed or moving solid bodies (cylinderand sphere) have been simulated and carefully validated against existing data from the literature in a large range of Reynolds numbers. With the objective of developing optimal control laws for VIV, based on the linear instability mechanisms of the coupled system within the framework of the control theory, a new adjoint solver was also developed and validated in OpenFOAM.