Cette thèse s'intéresse à la question récurrente qu'est la résolution d'un problème pour un grand nombre de configurations différentes. Malgré l'augmentation constante de la puissance de calcul que l'on connait aujourd'hui, le traitement direct d'un tel problème reste souvent hors de portée. La technique qui est développée ici est basée sur l'utilisation de la Proper Generalized Decomposition (PGD) dans le cadre de la méthode LATIN. On étudie tout d'abord la capacité de cette technique de réduction de modèle à résoudre un problème paramétré pour un espace de conception donné. Lors du traitement d'un tel problème, on génère une base réduite que l'on peut réutiliser et éventuellement enrichir en traitant un par un les problèmes correspondants aux jeux de paramètres étudiés. Le but devient alors de développer une stratégie, inspirée par la méthode " Reduced Basis ", afin d'explorer de façon rationnelle l'espace des paramètres. L'objectif étant de construire, avec le minimum de résolutions, une base réduite " complète " qui permet de résoudre tous les autres problèmes de l'espace de conception sans enrichir cette base. On commence dès lors par montrer l'existence d'une telle base complète en extrayant les informations les plus pertinentes des solutions PGD d'un problème pour tous les jeux de paramètres de l'espace de conception. On propose ensuite une stratégie rationnelle pour construire cette base complète sans la nécessité préalable de la résolution du problème pour tous les jeux de paramètres. Enfin, les performances de la méthode proposée sont illustrées sur plusieurs exemples, montrant des gains conséquents lorsque des études récurrentes doivent être menées.

Afin de tenir compte de l'amorçage et de la propagation des fissures en fatigue, il est nécessaire de connaître l'historique des différentes quantités mécaniques dans la zone d'amorçage. Cela nécessite une connaissance des états mécaniques cycliques stabilisés. D'un point de vue numérique, les simulations numériques d'agrégats polycristallin avec les méthodes de résolution classiques ne sont réalisées que pour quelques cycles. Ce travail présente le développement des méthodes accélérées de calcul pour réduire le temps de calcul de la méthode des Éléments Finis dans le cas des simulations d'agrégats polycristallins soumis à un chargement cyclique. La première idée est de maintenir la matrice de rigidité constante afin d'effectuer une seule factorisation. Un premier algorithme a été écrit dans ce sens avec une résolution incrémentale et non incrémentale. La deuxième proposition est d'utiliser une méthode de réduction dimensionnelle en espace/temps couplé à la méthode des éléments finis. La PGD a été choisie. Cette méthode permet de découpler les variables spatiales et temporelles et les champs de déplacement ne sont calculés que pour un certain nombre de modes. Deux possibilités d'utilisation de la méthode PGD dans le cadre de la plasticité ont été proposées. La troisième proposition consiste à utiliser la stratégie de réduction adaptative APR comme solveur afin de résoudre un modèle d'ordre réduit en termes de nombre de ddl. Une stratégie incrémentale d'amélioration de la qualité de la base pour un certain intervalle de temps choisi a été mis en place dans ce cadre. Quatre possibilités d'utilisation de la méthode APR ont été proposées. L'analyse des performances des différentes méthodes est effectuée tout d'abord sur un problème élasto-plastique classique tridimensionnel présentant un défaut sphérique, ensuite à l'échelle de la microstructure avec un calcul visco-élasto-plastique d'agrégats polycristallins tridimensionnels. Les analyses ont montré que les réponses mécaniques macroscopique et mésoscopique des méthodes de réduction de modèle sont très proches de la méthode incrémentale. Le temps des simulations a été réduit par les différentes méthodes. Les gains sont plus importants quand on augmente la taille des maillages et la non linéarité du problème. Toutefois, l'idée de garder la matrice de rigidité constante avec une résolution incrémentale a permis un gain de temps encore plus conséquent à l'échelle de la microstructure. / In order to take account of fatigue cracks initiation and growth, it is necessary to know the history of the various mechanical quantities in fatigue initiation site. This requires knowledge of the stabilized cyclic mechanical states. From a numerical approach, numerical simulation of polycrystalline aggregates with conventional resolution methods are only carried out for a few cycles. This work presents the development of accelerated numerical methods to reduce the computation time of the Finite Element method in the case of numerical simulation of polycrystalline aggregates under cyclic loading. The first idea is to keep a constant stiffness matrix during overall simulation in order to get just one single factorization to carry out. An algorithm has been proposed in this sense with an incremental and non incremental resolution. The second proposal is based on the use of a model reduction method coupled with the finite element method to solve space/time problem. The PGD has been selected. This method allows to decouple spatial and time coordinates and the displacement fields are computed for a certain number of modes. Two possibilities have been proposed to implement the PGD method in the context of plasticity. The third proposal is to use an a priori adaptative approach based on the use of APR strategy to solve a reduced order model in terms of number of degrees of freedom. An incremental adaptive strategie has been proposed in order to master the quality of the reduced-basis for a certain time steps. Four possibilities of using the APR method have been proposed. The applicability and the performance of the different methods have been first analyzed on a conventional three-dimensional elastoplastic problem with a spherical defect, then on the scale of the microstructure with numerical simulation of polycrystalline aggregates under cyclic elasto-visco-plastic loading. The analyzes have shown that the macroscopic and mesoscopic mechanical responses of the model reduction methods are very close to the incremental method. The simulation time has been reduced by the different methods. The time gains are more significant by increasing the size of the meshes and the non-linearity of the problem. Nevertheless, the idea of keeping a constant stiffness matrix with the incremental method has enabled more CPU time saving at the microstructural scale.

Les techniques sur la récupération d'énergie qui visent à permettre aux réseaux de capteurs sans fil (Wireless Sensor Network, WSN) de devenir autonomes, sont reconnues comme des élément cruciaux pour répondre aux futurs besoins des objets connectés portés par l'internet des objets (Internet of Things, IoT). C’est dans ce contexte que les matériaux fonctionnels piézoélectriques et magnétostrictifs, qui peuvent être utilisés dans une large gamme de systèmes de récupération d'énergie, ont un regain d’intérêt au cours de ces dernières années. Cette thèse porte sur la modélisation multiphysique de ces deux matériaux fonctionnels avec la méthode éléments finis et par la réduction de modèle pour les systèmes qui en résultent, sur la base de la décomposition propre généralisée (Proper Generalized Decomposition, PGD). La modélisation de ces matériaux fonctionnels reste difficile bien que la recherche dans ce domaine a été l'objet de plusieurs études depuis des décennies. Une multitude de difficultés existent, parmi lesquelles les trois suivantes qui sont largement reconnues. La première difficulté résulte de la description mathématique des propriétés de ces matériaux qui est compliquée ; ce qui est particulièrement vrai pour les matériaux magnétostrictifs pour lesquels leurs propriétés dépendent de facteurs environnementaux externes tels que la température, la contrainte et le champ magnétique d’excitation. La deuxième difficulté résulte des effets de couplage entre les champs électromagnétiques, élastiques et thermiques qui doivent être considérés mutuellement, ce qui est au-delà de la capacité de la plupart des outils de simulation existants. La troisième difficulté vient du fait que les systèmes deviennent de plus en plus compacts pour être intégrés et/ou embarqués. Dans ce cas la modélisation multi-échelle est nécessaire, ce qui signifie que des modèles numériques tridimensionnels (3D) doivent être employés. Le travail présenté ici fournit des solutions pour répondre aux difficultés mentionnées. Une modélisation multiphysique sur la base des formes différentielles est d'abord établie. Dans cette modélisation, les quantités sont discrétisés en utilisant les éléments de Whitney appropriés. Après la discrétisation, le système est résolu en un bloc unique, ce qui évite les itérations entre les solutions physiques différentes tout en conduisant à des convergences rapides. La formulation prend en compte, la loi de comportement linéaire des matériaux piézoélectriques, et une loi de comportement non linéaire pour les matériaux magnétostrictifs basée sur le principe de l’énergie libre exprimé par le modèle (Discrete Energy-Averaged Model, DEAM). La mise en œuvre de notre formulation permet de décrire les comportements des matériaux fonctionnels piézoélectriques et magnétostrictifs à des coûts numériques raisonnables. Suite à cela, deux algorithmes basés sur la PGD pour la réduction de modèle sont proposés. Ces deux algorithmes ont permis de réduire considérablement le problème dimensionnel des modèles multiphysiques tout en en conservant de très bonnes précisions. Les algorithmes proposés fournissent également des moyens pour gérer le couplage avec la non-linéarité d’une manière efficiente. L’ensemble de nos modèles sont vérifiés et validés par des exemples représentatifs. / The energy harvesting technology that aims to enable wireless sensor networks (WSN) to be maintenance-free, is recognized as a crucial part for the next generation technology mega- trend: the Internet of Things (IoT). Piezoelectric and magnetostrictive materials, which can be used in a wide range of energy harvesting systems, have attracted more and more interests during the past few years. This thesis focuses on multiphysics finite element (FE) modeling of these two materials and performing model reduction for resultant systems, based on the Prop- er Generalized Decomposition (PGD). Modeling these materials remains challenging although research in this area has been under- going over decades. A multitude of difficulties exist, among which the following three issues are largely recognized. First, mathematically describing properties of these materials is com- plicated, which is particularly true for magnetostrictive materials because their properties depend on factors including temperature, stress and magnetic field. Second, coupling effects between electromagnetic, elastic, and thermal fields need to be considered, which is beyond the capability of most existing simulation tools. Third, as systems becoming highly integrated whole-scale simulations become necessary, which means three dimensional (3D) numerical models should be employed. 3D models, on the other hand, quickly turns intractable if not properly built. The work presented here provides solutions in respond to the above challenges. A differential forms based multiphysics FE framework is first established. Within this frame- work quantities are discreted using appropriate Whitney elements. After discretization, the system is solved as a single block, thus avoiding iterations between different physics solutions and leading to rapid convergences. Next, the linear piezoelectric, and a free energy based nonlinear magnetostrictive constitutive model called Discrete Energy Averaged Model (DE- AM) are incorporated into the framework. Our implementation describes underlying material behaviors at reasonable numerical costs. Eventually, two novel PGD based algorithms for model reduction are proposed. With our algorithms, problem size of multiphysics models can be significantly reduced while final results of very good accuracy are obtained. Our algo- rithms also provide means to handle coupling and nonlinearity conveniently. All our methodologies are demonstrated and verified via representative examples.

Cette thèse s’inscrit dans le cadre de l’étude du post-flambement local des grandes struc- tures composites raidies. La simulation du post- flambement par la méthode des éléments- finis est aujourd’hui limitée par le coût du calcul en particulier pour les grandes structures. Seules des zones restreintes peuvent être étudiées, en négligeant les interactions global/local. L’objectif de cette thèse est de proposer une stratégie de calcul performante pour la simula- tion du post-flambement local des grandes structures raidies à partir des connaissances sur le comportement mécanique des structures en post-flambement et d’un découpage naturel le long des raidisseurs favorable au calcul parallèle.Dans la littérature, les méthodes de réduction de modèle adaptative ont démontré leur ca- pacité à réduire le nombre d’inconnues tout en maîtrisant l’erreur d’approximation de la solution des problèmes non-linéaires. Par ailleurs, les méthodes de décomposition de do- maine avec localisation non-linéaire sont particulièrement adaptées au calcul parallèle en mécanique des structures en présence de non-linéarités locales.Les travaux de thèse portent dans un premier temps sur une stratégie de réduction de modèle adaptative spécifique au cas du post-flambement. Dans le cas d’un flambement local d’une grande structure raidie une combinaison avec une méthode de décomposition de domaine primale est ensuite proposée. Toutes ces stratégies sont implémentées dans un code de re- cherche programmé pendant la thèse. / This thesis is part of the study of local post-buckling of large stiffened composite struc- tures. The finite element simulation of structures subjected to post-bucking still faces com- putational limits, especially for large structures. Only restricted area may be studied for now, neglecting global/local interactions.The aim of the thesis is to propose an efficient computational strategy for local post-buckling analysis of large stiffened structures from knowledge on mechanical behavior of post-buckling structures and a natural partitionning along stiffeners conducive to parallel computation. In litterature, the adaptive model reduction solving techniques have demonstrated their abi- lity to drastically reduce the number of unknowns as well as to control the approximation error of solving non-linear problems. Furthermore, domain decomposition methods with a non-linear local step are suited to parallel computation in structural mechanics in the pre- sence of local non-linearities.Our work deals first with an adaptive model reduction strategy dedicated to post-buckling problems. In order to adress larger stiffened structures subjected to local post-buckling, like an aircraft fuselage, partitioning is then performed. The model reduction, as well as the adap- tive procedure are written in the framework of the primal domain decomposition method with a non-linear local step. These strategies are implemented in a research code developed for the purpose of the thesis.

Cette thèse s'intéresse à la question récurrente qu'est la résolution d'un problème pour un grand nombre de configurations différentes. Malgré l'augmentation constante de la puissance de calcul que l'on connait aujourd'hui, le traitement direct d'un tel problème reste souvent hors de portée. La technique qui est développée ici est basée sur l'utilisation de la Proper Generalized Decomposition (PGD) dans le cadre de la méthode LATIN. On étudie tout d’abord la capacité de cette technique de réduction de modèle à résoudre un problème paramétré pour un espace de conception donné. Lors du traitement d’un tel problème, on génère une base réduite que l’on peut réutiliser et éventuellement enrichir en traitant un par un les problèmes correspondants aux jeux de paramètres étudiés. Le but devient alors de développer une stratégie, inspirée par la méthode « Reduced Basis », afin d’explorer de façon rationnelle l’espace des paramètres. L’objectif étant de construire, avec le minimum de résolutions, une base réduite « complète » qui permet de résoudre tous les autres problèmes de l’espace de conception sans enrichir cette base. On commence dès lors par montrer l’existence d’une telle base complète en extrayant les informations les plus pertinentes des solutions PGD d’un problème pour tous les jeux de paramètres de l’espace de conception. On propose ensuite une stratégie rationnelle pour construire cette base complète sans la nécessité préalable de la résolution du problème pour tous les jeux de paramètres. Enfin, les performances de la méthode proposée sont illustrées sur plusieurs exemples, montrant des gains conséquents lorsque des études récurrentes doivent être menées. / This thesis deals with the recurring question of the resolution of a problem for many different configu- rations, which can lead to highly expensive computations when using a direct treatment. The technique which is presented here is based on the use of Proper Generalized Decomposition (PGD) in the framework of the LATIN method. The feasibility of this model reduction technique approach is studied to compute the solution of a parametrized problem for a given space of parameters. For that purpose, a Reduced-Order Basis is generated, reused and eventually enriched, by treating, one-by-one, all the various parameter sets. The aim is to develop a strategy, inspired by the Reduced Basis method, to explore rationally the space of parameters. Then, the objective is to build, with the minimum of resolutions, a ‘‘complete’’ basis that enables to solve all the other problems without enriching the basis. We first exemplify the existence of a such complete basis by extracting the most relevant information from the PGD solutions of the problem for all the sets in the space of parameters. Secondly, we propose a rational strategy to build this complete basis without preliminary solving the problem for all the sets of parameters. Finally, the capabilities of the proposed method are illustrated through a variety of examples, showing substantial gains when recurrent studies need to be carried out.

Les liaisons boulonnées ou rivetées sont, en apparence, l’un des systèmes mécaniques parmi les plus simples qui soient. La fonction première de ces éléments est d’assurer un encastrement rigide entre les composants qu’elles assemblent. Le dimensionnement de ces pièces aux efforts statiques nominaux est un problème assez largement résolu. Néanmoins, ces composants suscitent l’intérêt des ingénieurs et chercheurs depuis plus de 50 ans. Ce paradoxe s’explique en partie par l’amortissement important que ces liaisons produisent au sein des structures. Les efforts subis par les liaisons d’assemblages provoquent le glissement partiel des interfaces de contact. Il résulte de ce glissement, une dissipation par frottement qui est une source majeure d’amortissement des structures aéronautiques. Le niveau vibratoire d’un système est directement lié à son amortissement, c’est à dire à sa capacité à dissiper ou à accumuler de l’énergie. Dans le domaine du transport, entre autres, les vibrations sont néfastes car elles nuisent au confort des usagers ou à l’intégrité des structures. Le bon dimensionnement des assemblages est donc de nature à améliorer le comportement vibratoire des systèmes mécaniques. Actuellement, les outils de calcul numériques permettent d’estimer assez précisément les modes et fréquences propres d’une structure a priori mais l’amortissement, c'est-à-dire le niveau vibratoire, reste une donnée mesurée, a posteriori au travers d’essais couteux. Ceci s’explique par le caractère multi-échelle et la complexité des problèmes de contact. Ainsi, le comportement dynamique des assemblages reste un sujet d’étude très privilégié. Les travaux de cette thèse cherchent à répondre de manière pratique au besoin de produire des modèles réduits de structures assemblées pour la dynamique. / At first sight, riveted and bolted connections seem to be one of the simplest mechanical systems possible. The primary function of these elements is to provide a rigid clamping between the components they assemble. The designing of these parts to nominal static stresses is already quite mastered. However, engineers are still bothered when it comes to modeling these components which have been studied for more than 50 years. This paradox can be explained by the large damping ratios joints generate in structures. The vibration level of a system is directly related to its damping ratio, i.e. its ability to dissipate or store energy. In the field of transport, among others, vibrations are unwanted because they affect the comfort of the users or the integrity of structures. A good designing of jointed structures is likely to improve the vibration behavior of mechanical systems. Nowadays, numerical calculations allow for the computing of the modes but the damping is still measured, a posteriori, through expensive tests. This is explained by the multi-scale nature and the complexity of contact physics. That is why the dynamic behavior of assembled structures remains an area of study for researchers. The work of this thesis aims to provide practical solutions for the identification and the designing of reduced order models for the dynamics of assembled structures.

Cette thèse est dédiée au développement de méthodes numériques pour la poromécanique, et à leur application dans un contexte de modélisation cardiaque.Elle est motivée par la prise en compte, dans les modèles de coeur humain, du réseau coronarien qui perfuse le myocarde, afin de mieux décrire les maladies vasculaires coronariennes.Nous appuyant sur des travaux existants, nous proposons un modèle de coeur perfusé, ainsi qu'une réduction 0D permettant de reproduire, à moindre coût de calcul, un cycle cardiaque réaliste avec masse et pression de perfusion. Le modèle mis au point nous permet de reproduire des phénomènes physiologiques auparavant inaccessibles dans les modèles, et d'une grande importance pour des applications cliniques, tels que la vasodilatation et les pathologies coronariennes.L'intégration d'un compartiment poreux pour représenter le myocarde perfusé dans les modèles 3D représente un défi technique d'un autre ordre. Nous inspirant des schémas en temps de type splitting établis en interaction fluide-structure pour modéliser les vaisseaux sanguins, nous proposons une discrétisation semi-implicite d'une formulation générale de poromécanique, satisfaisant un bilan d'énergie au niveau discret. Afin d'illustrer et valider notre démarche, l'environnement de calcul élément finis FreeFem++ nous permet de reproduire des cas tests classiques de gonflement et de drainage de milieux poreux en 2D, puis de vérifier le bilan énergétique discret.Enfin, motivés par l'étude de la discrétisation spatiale de notre problème, nous établissons dans un cadre linéaire un résultat de convergence totale du schéma sous conditions. Cela nous permet de proposer une méthode d'implémentation facile à mettre en oeuvre et présentant de bons résultats de stabilité. FreeFem++ nous permet à nouveau de valider nos résultats en illustrant les pathologies numériques associées à l'incompressibilité, et leur traitement efficace par les stratégies proposées, dans le cadre linéaire puis dans une situation plus générale. / This thesis is dedicated to the development of numerical methods for poromechanics, and to their application in a cardiac modeling context. It is motivated by the introduction into existing cardiac models of the coronary network that perfuses the myocardium, to better describe coronary vascular diseases.Drawing our inspiration from existing works, we propose a perfused heart model, and a 0D reduction allowing the cost-effective reproduction of a realistic cardiac cycle with perfusion mass and pressure. The model derived illustrates physiological phenomena inaccessible in former models, and with great clinical application potential, such as vasodilatation and coronary diseases.The integration of a porous compartment to represent the perfused myocardium within 3D models is more challenging. Relying on splitting time schemes established for fluid-structure interaction to model blood vessels, we propose a semi-implicit discretization of a general poromechanics formulation, satisfying a discrete energy balance. In order to illustrate and validate our approach, we reproduce in the finite element software FreeFem++ classical swelling and drainage 2D test cases, and we monitor the discrete energy balance.Finally, motivated by the study of spatial discretization aspects of our problem, we establish in a linear framework a conditional total convergence result. This enables us to propose a computational method easy to implement and presenting good stability results. FreeFem++ enables us again to validate our results, illustrating numerical pathologies associated with incompressibility, and their efficient treatment with the proposed strategies, first in a linear framework and then in a more general situation.

Nous proposons dans ce manuscrit une nouvelle méthode pour prendre en compte l’effet du gradient en Fretting-fatigue. Les champs mécaniques présents à proximité du front de contact sont décrits à travers des facteurs d’intensité non locaux. L’objectif est d’aboutir à une description du champ de vitesse sous la forme d’une somme de termes exprimés chacun comme le produit d’un facteur d’intensité (Is, Ia, Ic), qui dépend des chargements macroscopiques appliqués à l’ensemble et d’une fonction de forme (ds, da, dc), qui est liée à la géométrie locale du contact. Cette description est obtenue à travers un processus non intrusif de post-processing des résultats obtenus avec des calculs à éléments finis. De plus, elle a été pensée pour être implémentée dans un contexte industriel. En pratique, pour chaque chargement macroscopique et pour chaque géométrie, il est possible de calculer un ensemble de facteurs d’intensité non locaux qui permettent de décrire les champs mécaniques locaux près du front de contact. Cette description non locale a l’avantage d’être (i) indépendante de la géométrie du contact employé et (ii) utilisable dans des modèles à éléments finis utilisés dans l’industrie qui sont caractérisés par des maillages plus grossiers par rapport à ceux utilisés pour étudier le fretting-fatigue dans des milieux académiques. Une étude est menée pour vérifier que les facteurs d’intensité non locaux peuvent être utilisés pour transposer les résultats expérimentaux d’une géométrie à une autre. / In this manuscript a new method to describe the stress gradient effect in fretting-fatigue is proposed. It is based on the description of the mechanical fields arising close to the contact edges through nonlocal intensity factors. For this purpose, the kinetic field around the contact ends is partitioned into a summation of multiple terms, each one expressed as the product between intensity factors, Is, Ia, Ic, depending on the macroscopic loads applied to the mechanical assembly, and spatial reference fields, ds, da, dc, depending on the local geometry of the part. This description is obtained through nonintrusive post-processing of FE computation and is conceived in order to be easily implementable in the industrial context. As a matter of fact, for any given macroscopic load and geometry, a set of nonlocal intensity factors is computed that permits to characterize the mechanical fields close to the contact edges. Such nonlocal description has the advantage of being (i) geometry independent so that the nonlocal intensity factors can be used to compare laboratory test with real-scale industrial assembly, (ii) applicable to industrial FE models usually characterized by rougher meshes compared to the ones used to describe fretting-fatigue in the academic context. The procedure is applied to fretting-fatigue test data in order to verify whether the nonlocal intensity factors can be used to transpose experimental results to different contact geometries from the one in which they have been obtained.

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

La dynamique d'instabilité d'un écoulement laminaire décollé est étudiée expérimentalement et son contrôle par le biais de la simulation numérique. La configuration étudiée est une couche limite laminaire décollée au dessus d'une géométrie de type bosse.Pour une certaine gamme de paramètres, l'écoulement de recirculation en aval de la bosse est caractérisé par un battement basse fréquence. L'étude expérimentale de cette dynamique a permis de retrouver les différents régimes d'instabilité mis a jour par voie numérique. Ces résultats prouvent notamment que les instabilités basse fréquence, dont l'existence a été surtout mise en évidence dans des configurations d'écoulements compressibles, sont un phénomène générique pour des bulles de recirculations allongées. Le contrôle du battement basse fréquence est ensuite étudié par voie numérique suivant deux approches complémentaires. Un asservissement en boucle fermée de la dynamique de perturbation linéaire est tout d'abord proposé. Les modes d'instabilité linéaires sont utilisés afin de construire des modèles réduits de la dynamique de perturbation. Cette réduction de modèle donne lieu à des estimateurs de faible dimension capables d'estimer la dynamique et de la contrôler. Ainsi la dynamique d'instabilité linéaire peut être supprimée en couplant le système de Navier-Stokes linéarisé avec le contrôleur.Le contrôle de la dynamique non linéaire est ensuite étudié en utilisant une méthode d'optimisation Lagrangienne. Cette méthode permet de calculer les lois de contrôle à partir de la dynamique non linéaire des équations de Navier-Stokes. / The dynamics and control of a separated boundary-layer flow have been investigated. Separation is triggered by a bump mounted on a flat plate and the transition dynamics has been investigated experimentally. For a certain parameter range, the recirculation region is subject to self-sustained low-frequency oscillations, and results from the numerical simulation for the same geometry are recovered. These results show that low frequency oscillations, observed mainly in compressible flow regimes, are inherent to elongated recirculation bubbles.The control of this low-frequency instability has been investigated using modern control theory based on two complementary approaches. Feedback control of the linear perturbation dynamics is first considered. Global instability modes are used to build reduced-order estimators. This model reduction gives rise to low-dimensional compensators capable of controlling the unstable dynamics. Once coupled to the unstable linearised Navier-Stokes system, the compensator is seen to succesfully control the unstable dynamics. The control of the nonlinear dynamics is then investigated using adjoint-based optimisation procedures. This method is used to compute control laws based on a complete knowledge of the nonlinear dynamics. Although the low-frequency instability is clearly attenuated, it seems difficult to control the flow towards its steady state, using only a few blowing/suction actuators localized on the wall.