Dans ce travail de thèse, une méthode de pénalisation en volume pour la simulation d'écoulements magnétohydrodynamiques (MHD) en domaines confinés est présentée. Les équations incompressibles de la MHD résistives sont résolues par le truchement d'un solveur pseudo-spectral parallèlisé. La pénalisation en volume est une méthode de frontières immergées, caractérisée par une grande flexibilité dans le choix de la géométrie de l'écoulement. Dans le cas présent, elle permet d'utiliser des conditions aux limites non-périodiques dans un schéma pseudo-spectral Fourier. La méthode numérique est validée et sa convergence est quantifiée pour des écoulements hydrodynamiques et MHD, en deux et trois dimensions, en comparant les résultats numériques à ceux de la littérature et à des solutions analytiques. Dans un second temps, la génération spontanée de moment cinétique et magnétique est étudiée pour des écoulements MHD confinés 2D et 3D. L'influence du nombre de Reynolds et du rapport des énergies cinétique/magnétique est explorée, ainsi que les différences induites par les conditions aux limites. Le fait que l'axisymétrie des frontières résulte en un terme de pression non-nul est primordial pour engendrer de grandes valeurs du moment cinétique. L'exclusivité de cette auto-organisation aux écoulements 2D est étudiée en considérant la MHD 3D en présence d'un fort champ magnétique axial. La suite est consacrée à la simulation d'un fluide conducteur dans un cylindre avec un forçage magnétique axial et poloidal. En faisant varier l'amplitude du forçage poloidal, différents états dynamiques sont atteints. Enfin, l'effet du nombre de Prandtl sur le seuil des instabilitées est étudié. / A volume penalization method for the simulation of magnetohydrodynamic (MHD) flows in confined domains is presented. Incompressible resistive MHD equations are solved in 3D by means of a parallelized pseudo-spectral solver. The volume penalization technique is an immersed boundary method, characterized by a high flexibility in the choice of the geometry of the considered flow. In the present case, it allows the use of conditions different from periodic boundaries in a Fourier pseudo-spectral scheme. The numerical method is validated and its convergence is assessed for two- and three-dimensional hydrodynamical and MHD flows by comparing the numerical results with those of the literature or analytical solutions. Then, the spontaneous generation of kinetic and magnetic angular momentum is studied for confined 2D and 3D MHD flows. The influence of the Reynolds number and of the ratio of kinetic/magnetic energies is explored, as well as the differences induced by the boundary conditions. The fact that axisymmetric borders introduce a non-zero pressure term in the evolution equation of the angular momentum is essential to generate large values of the angular momentum. It is investigated whether this self-organization is exclusively observed in 2D flows by considering 3D MHD in the presence of a strong axial magnetic field. The last part is devoted to the simulation of a conducting fluid in a periodic cylinder with imposed axial and poloidal magnetic forcing, implying a resulting magnetic field. By varying the amplitude of the poloidal forcing, different dynamical states can be achieved. The effect of the Prandtl number on the threshold of the instabilities is then studied.

La vaporisation rapide du sodium liquide surchauffé est supposée être à l’origine des arrêts automatiques pour réactivité négative du réacteur Phénix.Un dispositif expérimental a été mis en œuvre pour reproduire la détente d'un gaz pressurisé, repoussant un liquide dans un canal de section rectangulaire très allongée.L’interface qui sépare les deux fluides, initialement plate, ondule du fait d'instabilités de Rayleigh-Taylor dont le caractère 2D est garanti par le rapport d'aspect de la section du canal. L’aire interfaciale augmente d'un facteur 50.L’expansion du gaz peut être divisée en deux phases principales : une phase dite « de Rayleigh-Taylor » (linéaire et non-linéaire) et une phase dite « à multi-structures » (transitionnelle et chaotique). La première est caractérisée par la dynamique de l'interface et l’aire interfaciale qui en résulte est proportionnelle à l’amplitude des ondulations. La deuxième est influencée par le comportement des structures liquides, dispersées dans la matrice gazeuse et l’aire interfaciale est alors proportionnelle au nombre de structures.La distribution de fraction volumique suggère un modèle d’écoulement composé de trois régions : une région où la frontière des bulles est clairement définie et régulière, une région compartimentée par des membranes liquides issues des frontières des bulles, une région diphasique formée de la queue de ces structures. L’analyse de sensibilité à la tension superficielle confirme que plus la tension est faible, plus les interfaces sont instables. Les ondes sont plus prononcées et plus de structures sont produites, ce qui conduit à une majoration du taux de production de l’aire interfaciale. / The sharp vaporization of superheated liquid sodium is investigated. It is suspected to be at the origin of the automatic shutdown for negative reactivity, occurred in the Phénix reactor at the end of the eighties.An experimental apparatus has been designed and operated to reproduce the expansion of overpressurized air, superposed to water in a narrow vertical rectangular section channel.When expansion begins, the initial flat interface separating the two fluids becomes corrugated under the development of two-dimensional Rayleigh-Taylor instabilities. The interface area increases significantly and becomes even 50 times larger than the initial value. Since the channel is very narrow, instabilities along the channel depth do not develop.The gas expansion in a narrow channel can be divided into two main phases: Rayleigh-Taylor (linear and non-linear) and multi-structures (transition and chaotic) phases. The former is characterized by the dynamic of corrugated profile and the interface area results proportional to the amplitude of corrugation The latter is influenced by the behavior of the liquid structures dispersed in gas matrix and the interface area is mainly proportional to the number of liquid structures.The distribution of volume fraction suggests a model of channel flow consisting of three regions: the regular profile of peaks, the spike region and the structures tails. The analysis of sensibility to surface tension confirms that, with a lower surface tension, the fluids configuration is more unstable. The interface corrugations are more pronounced and more structures are produced, leading to a higher increment of the interface area.

Les instabilités hydrodynamiques en cavité chauffée latéralement jouent un rôle important dans certains processus de fabrication de matériaux tels que le procédé de Bridgman horizontal. En effet, le fluide (métal liquide qui va se solidifier) est le siège d’une circulation thermoconvective due à l’existence d’un gradient de température horizontal qui est susceptible de devenir instationnaire via des instabilités oscillatoires. La connaissance et la maîtrise de ces instabilités sont donc primordiales afin de pouvoir améliorer la qualité des cristaux obtenus par cette technique. Dans cette thèse, nous nous sommes intéressés en premier aux instabilités affectant la circulation convective dans une cavité tridimensionnelle de dimensions 4×2×1. (longueur × largeur × hauteur). Grâce aux techniques numériques de continuation, nous avons pu obtenir les solutions stationnaires et oscillatoires, ainsi que leur stabilité, jusqu’à l’apparition de la quasi-périodicité en fonction du nombre de Grashof Gr et pour un nombre de Prandtl allant de 0 à 0,025. Ensuite, pour un éventuel contrôle des instabilités, nous nous sommes intéressés aux effets induits par la rotation de la cavité. Nous avons tout d’abord considéré un modèle monodimensionnel que nous avons développé durant cette thèse. Ce modèle analytique, bien que simplifié, est en très bon accord avec les observations en dynamique des écoulements atmosphériques (déviation des masses fluides vers la droite de la composante de vitesse dominante et vents thermiques). La stabilité linéaire de cet écoulement est ensuite effectuée en fonction du taux de rotation donné par le nombre de Taylor et du nombre de Grashof pour un nombre de Prandtl allant de 0 à 10. Nous avons pu montrer à travers ce modèle que la rotation possède un caractère stabilisant vis-à-vis de ce type d’écoulement. Enfin, nous nous sommes focalisés sur les effets de la rotation sur l’écoulement pleinement tridimensionnel dans la cavité de dimensions 4×2×1. Nous avons mis en évidence deux régimes d’écoulements : un régime dominé par la convection, où la circulation du fluide est déviée par la rotation dans la diagonale de la cavité, et un deuxième régime dominé par la rotation où la circulation du fluide est concentrée dans les couches limites dites d’Ekman et de Stewartson. Un très bon accord est observé entre le modèle analytique simplifié et la simulation numérique tridimensionnelle. / Hydrodynamic instabilities in laterally heated cavities play an important role in some material processing techniques such as the horizontal Bridgman process. Indeed, the fluid (liquid metal to be solidified) is the seat of a thermoconvective circulation due to the existence of a horizontal temperature gradient which is likely to become unsteady via oscillatory instabilities. The knowledge and the control of these instabilities are thus essential in order to be able to improve the quality of the crystals obtained by this technique. In this thesis, we are first interested in the instabilities of the convective circulation in a three-dimensional cavity of dimensions 4×2×1 (length × width × height). Thanks to the numerical continuation techniques, we were able to obtain the stationary and oscillatory solutions, as well as their stability, until the appearance of the quasi-periodicity according to the Grashof number Gr and for a Prandtl number Pr ranging from 0 to 0.025.Then, the effects induced by a rotation of the cavity around the vertical axis parallel to gravity (for a possible control of the instabilities) are studied and a one-dimensional model developed during this thesis was first considered. This analytical model, although simplified, is in very good agreement with the observations of the atmospheric flows (deviation of the fluid masses towards the right of the component of the dominant velocity and thermal winds). The linear stability of this flow as well as an energy analysis at the thresholds are then performed as a function of the rotation rate given by the Taylor number Ta and the Grashof number Gr for a Prandtl number Pr ranging from 0 to 10. Through this model, we have been able to show that the rotation has a stabilizing effect on this type of flow.We finally focused on the effects of this type of rotation on the steady fully threedimensional flow observed in the cavity 4×2×1 at low Grashof numbers.We have highlighted two flow regimes: a regime dominated by convection where the fluid circulation, deviated by the rotation, occurs in the diagonal of the cavity, and a second regime dominated by rotation where the fluid circulation is concentrated in the so-called Ekman and Stewartson boundary layers. A very good agreement is observed between the simplified analytical model and the three-dimensional numerical simulation.

Une rivière alluviale s'écoule sur une épaisse couche de sédiments. Lorsqu'elle construit son lit, elle entraîne, transporte et dépose des sédiments, façonnant ainsi sa propre forme. Ainsi, le couplage entre l'écoulement et le transport sédimentaire régit la taille et la forme de la rivière. Dans cette thèse, nous étudions l'influence du transport sédimentaire sur la forme et la stabilité d'une rivière alluviale. Pour ce faire, nous reproduisons des rivières en laboratoire en laissant s'écouler un liquide visqueux sur un lit granulaire. L'aspect du chenal ainsi formé dépend des débits de liquide et de sédiment injectés en entrée. A l'aide de ces expériences, nous mettons en évidence les deux mécanismes qui contrôlent l'équilibre d'une rivière. D'abord, la gravité entraîne les grains vers le centre du chenal. Ce mécanisme érode continuellement les berges de la rivière, et tend donc à l'élargir. Cependant, les collisions d'un grain avec le lit dévient sa trajectoire dans la direction transverse à l'écoulement. Les grains se comportent ainsi comme des marcheurs aléatoires, qui, collectivement, diffusent vers les berges de la rivière. A l'équilibre, cette diffusion compense la gravité, et fixe ainsi la forme de la rivière. Lorsque la diffusion prend le dessus sur la gravité, elle peut induire une instabilité. En effet, si on perturbe un lit sédimentaire avec des stries longitudinales, le cisaillement fluide est plus faible là où l'écoulement est moins profond. Par conséquent, les grains diffusent depuis les creux de la perturbation vers ses crêtes. Cette rétroaction déstabilisante pourrait générer de nouveaux chenaux, et expliquer la formation des rivières en tresses. / An alluvial river builds its own bed with the sediment it transports. The channel bounds the flow, which in turns deforms the channel through erosion and deposition. This coupling between flow and sediment transport selects the shape and the size of the river. In this manuscript, we investigate it using laboratory experiments. The first ingredient of this coupling is gravity, which pulls the moving grains towards the center of the channel, thus continually eroding the banks. However, due to the roughness of the bed, the trajectory of a moving grain fluctuates across the stream. The bedload layer is therefore a collection of random walkers which diffuse towards the less active areas of the bed. In a river at equilibrium, this diffusion counteracts gravity to maintain the banks. When gravity and diffusion are out of balance, their interaction causes an instability. Indeed, if an initially flat bed of sediment is perturbed with longitudinal streaks, the flow-induced shear stress is weaker where the flow is shallower. Therefore, bedload diffusion induces a sediment flux towards the crests of the perturbation. This positive feedback induces an instability which can generate new channels. We suggest that this mechanism could initiate the braiding of alluvial rivers.

L'explosion en supernova gravitationnelle représente le stade ultime de l'évolution des étoiles massives.La contraction du cœur de fer peut être suivie d'une gigantesque explosion qui donne naissance à une étoile à neutrons.La dynamique multi-dimensionnelle de la région interne, pendant les premières centaines de millisecondes, joue un rôle crucial sur le succès de l'explosion car des instabilités hydrodynamiques sont capables de briser la symétrie sphérique de l'effondrement.Les mouvements transverses et à grande échelle générés par deux instabilités, la convection induite par les neutrinos et l'instabilité du choc d'accrétion stationnaire (SASI),augmentent l'efficacité du chauffage de la matière par les neutrinos au point de déclencher une explosion asymétrique et d'impacter les conditions de naissance de l'étoile à neutrons. Dans cette thèse, les instabilités sont étudiées au moyen de simulations numériques de modèles simplifiés.Ces modèles permettent une vaste exploration de l'espace des paramètres et une meilleure compréhension physique des instabilités, généralement inaccessibles aux modèles réalistes.L'analyse du régime non-linéaire de SASI établit les conditions de formation d'un mode spiral et évalue sa capacité à redistribuer radialement le moment cinétique.L'effet de la rotation sur la dynamique du choc d'accrétion est également pris en compte.Si la rotation est suffisamment rapide, une instabilité de corotation se superpose à SASI et impacte grandement la dynamique.Les simulations permettent de mieux contraindre l'importance des modes non-axisymétriques dans le bilan de moment cinétique de l'effondrement du cœur de fer en étoile à neutrons.SASI pourrait sous certaines conditions accélérer ou ralentir la rotation du pulsar formé dans l'explosion.Enfin, une étude d'un modèle idéalisé de la région de chauffage est menée pour caractériser le déclenchement non-linéaire de la convection par des perturbations telles que celles produites par SASI ou les inhomogénéités de combustion pré-effondrement.L'analyse de la dimensionnalité sur le développement de la convection permet de discuter l'interprétation des modèles globaux et met en évidence les effets bénéfiques de la dynamique tridimensionnelle sur le déclenchement de l'explosion. / A core-collapse supernova represents the ultimate stage of the evolution of massive stars.The iron core contraction may be followed by a gigantic explosion which gives birth to a neutron star.The multidimensional dynamics of the innermost region, during the first hundreds milliseconds, plays a decisive role on the explosion success because hydrodynamical instabilities are able to break the spherical symmetry of the collapse.Large scale transverse motions generated by two instabilities, the neutrino-driven convection and the Standing Accretion Shock Instability (SASI),increase the heating efficiency up to the point of launching an asymmetric explosion and influencing the birth properties of the neutron star.In this thesis, hydrodynamical instabilities are studied using numerical simulations of simplified models.These models enable a wide exploration of the parameter space and a better physical understanding of the instabilities, generally inaccessible to realistic models.The non-linear regime of SASI is analysed to characterize the conditions under which a spiral mode prevails and to assess its ability to redistribute angular momentum radially.The influence of rotation on the shock dynamics is also addressed.For fast enough rotation rates, a corotation instability overlaps with SASI and greatly impacts the dynamics.The simulations enable to better constrain the effect of non-axisymmetric modes on the angular momentum budget of the iron core collapsing into a neutron star.SASI may under specific conditions spin up or down the pulsar born during the explosion.Finally, an idealised model of the heating region is studied to characterize the non-linear onsetof convection by perturbations such as those produced by SASI or pre-collapse combustion inhomogeneities. The dimensionality issue is examined to stress the beneficial consequences of the three-dimensional dynamics on the onset of the explosion.

Dans la configuration de Rayleigh-Bénard, on considère une fine couche de fluide placée entre deux parois horizontales, chauffée par le bas et refroidie par le haut. Cette couche peut être le siège d'une instabilité si le gradient thermique est suffisamment important : on passe alors de l'état conductif à l'état convectif et on parle de bifurcation primaire pour qualifier cette première transition. Cette mise en mouvement du fluide se fait de manière ordonnée : on constate ainsi l'émergence de différents motifs de convection comme des rouleaux, des carrés ou encore des hexagones. Ces structures vont à leur tour subir des instabilités qualifiées de secondaires qui vont limiter la gamme de nombres d'onde stables. On étudie ici théoriquement ces instabilités d'une part à proximité du seuil de la convection grâce à une approche faiblement non linéaire, d'autre part loin des conditions critiques grâce à une approche fortement non linéaire. Le fluide est rhéofluidifiant, ce qui correspond au comportement rhéologique le plus fréquemment rencontré, et est décrit par le modèle de Carreau. À proximité du seuil, on considère deux situations : la première correspond au cas où les plaques ont une conductivité finie, la seconde à celui d'un fluide thermodépendant. Dans chaque cas, l'influence du caractère rhéofluidifiant sur la nature du motif émergeant à la bifurcation primaire et sur les instabilités secondaires est mise en évidence. Pour étudier les motifs de convection loin des conditions critiques, on a recours à une procédure de continuation permettant de déterminer de proche en proche les caractéristiques de l'écoulement comme les champs de vitesse ou de température ainsi que le nombre de Nusselt. / In the Rayleigh-Bénard configuration, we consider a thin layer of fluid confined between two horizontal slabs which is heated from below and cooled from above. This layer undergoes an instability if the thermal gradient is strong enough: a transition from the conductive state to the convective state and called _ primary bifurcation _occurs. Moreover, it happens in an ordered way: we can notice the emergence of various convection patterns such as rolls, squares or hexagons. In their turn, these patterns undergo _ secondary instabilities _ that limit the range of stable wavenumbers. These instabilities are studied theoretically _firstly near the threshold thanks to a weakly nonlinear approach, and secondly far from critical conditions thanks to a strongly nonlinear approach. We consider a shear thinning fluid, the most common rheological behavior, which is described by the Carreau model. Near the threshold, two situations are considered: the first corresponds to finite conductivity plates, the second corresponds to a thermodependent fluid. In each case, the influence of the shear thinning effect on the nature of the pattern emerging at the primary bifurcation and on secondary instabilities is highlighted. To study the convection patterns far from the critical conditions, a continuation procedure is used to determine, step-by-step, the characteristics of the flow, such as the velocity or temperature fields and the Nusselt number.

Les surfaces cristallines jouent un rôle essentiel dans grand nombre d'applications des nanotechnologies, une surface étant souvent un support nécessaire pour des nano-objets. La surface support peut être à son tour nano-structurée, et intervenir ainsi directement dans les propriétés du matériau réalisé. Ainsi, parmi les défis à relever pour maîtriser la morphologie de la surface et donc contrôler les nanostructures, il y a la compréhension des mécanismes moteurs à l'origine des instabilités de surface. Ce travail de thèse constitue une étude théorique accompagnée par des simulations Monte-Carlo de la dynamique des instabilités de croissance. C'est une contribution à la compréhension des mécanismes microscopiques gouvernant ces instabilités. Deux démarches ont été suivies dans cette étude : la première est une étude statistique qui consiste à analyser la morphologie de la surface instable, puis de proposer des modèles numériques (barrière Ehrlich-Schwoebel, impuretés) et des équations analytiques capables de décrire ces instabilités et remonter enfin aux mécanismes atomistiques sous-jacents. La deuxième, est une étude de la dynamique des marches qui a pour but de prévoir puis de déterminer les mécanismes atomistiques exacts gouvernants le mouvement de ces marches et les échanges de matière entre elles à l'aide de mesures de la distribution de largeur des terrasses (TWD). Les résultats de ce travail ont abouti à la mise en évidence de nouveaux mécanismes de nanostructuration spontanée [impuretés], ainsi que sur le développement et la validation numérique de nouveaux outils mathématiques [approche de Fokker-Planck] pour décrire les fluctuations des séparations entre les marches (TWDs) sur les surfaces loin de l'équilibre.

Les satellites artificiels baignent dans un environnement radiatif hostile qui conditionne en partie leur fiabilité et leur durée de vie en opération : les ceintures de Van Allen. Afin de les protéger, il est nécessaire de caractériser la dynamique des électrons énergétiques piégés dans les ceintures radiatives. Elle est déterminée essentiellement par les interactions entre les électrons énergétiques et les ondes électromagnétiques existantes. <br /><br />Le travail de cette thèse a consisté à concevoir un schéma numérique original pour la résolution du système d'équations modélisant ces interactions : les équations de Vlasov-Maxwell relativistes. Notre choix s'est orienté vers des méthodes d'intégration directe. Nous proposons trois nouvelles méthodes spectrales pour discrétiser en impulsion les équations : une méthode de Galerkin et deux méthodes de type collocation. Ces approches sont basées sur des fonctions de Hermite qui ont la particularité de dépendre d'un facteur d'échelle permettant d'obtenir une bonne résolution en vitesse. <br /><br />Nous présentons dans ce manuscript les calculs conduisant à la discrétisation et à la résolution du problème de Vlasov-Poisson monodimensionnel ainsi que les résultats numériques obtenus. Puis nous étudions les extensions possibles des méthodes au problème complet relativiste. Afin de réduire les temps de calcul, une parallélisation et une optimisation des algorithmes ont été mises en \oe uvre. Enfin, les calculs de validation du code 1Dx-3Dv, à partir d'instabilités de types Weibel et whistlers, à une ou deux espèces d'électrons, sont détaillés.

Les gaz ionisés contenant des particules solides, appelés plasmas poussiéreux, sont principalement utilisés pour le dépôt de couches minces et la synthèse de nanoparticules aux propriétés maîtrisées.<br />L'objectif de cette thèse était de synthétiser des poussières de taille et/ou composition chimique connues en vue d'applications en microélectronique, photovoltaïque (nanocristaux de silicium) et astrophysique (tholins, ISD grains).<br />La localisation et l'étude du démarrage de la phase d'agrégation des nanocristaux, par l'étude paramétrique d'une instabilité, ont été réalisées dans des plasmas d'Argon/Silane. Puis, une étude de la forme et du comportement du nuage de poussières a été menée dans deux décharges différentes. Elles ont mis en évidence la forte influence de la croissance des poussières sur le plasma et sur le nuage, et d'optimiser nos paramètres de dépôt pour les nanocristaux de silicium non-agglomérés. Une étude paramétrique des caractéristiques électriques du plasma (courant, tension et densités électronique et ionique) en Azote/Méthane a permis de cibler les meilleurs paramètres pour synthétiser des tholins (aérosols analogues de Titan). Par ailleurs, l'ellipsométrie in-situ de Mie-Rayleigh en chimie Azote/Acétylène nous a donné certaines informations sur les analogues de poussières interstellaires synthétisés (structure, indice optique). Finalement, l'effet des basses températures de gaz a été exploré afin d'augmenter la taille des nanocristaux de silicium. Différentes hypothèses (chimiques et thermodynamiques) sont discutées dans le but d'expliquer les effets observés : accélération de la cinétique de croissance et augmentation de la taille des nanocristaux.

Les moteurs à propergol solide sont parfois le siège d'instabilités aéroacoustiques résultant d'un couplage entre l'hydrodynamique des gaz brûlés et les modes acoustiques de la chambre de combustion. Ces instabilités se traduisent par de fortes Oscillations de Pression (ODP) dans la chambre de combustion du moteur. Ces ODP entrainent des vibrations de la structure, qui si elles venaient à dépasser certains niveaux pourraient nuire à la charge utile. Au vu du coût d'un essai, il est important de disposer d'outils permettant de prédire l'apparition de ces instabilités au moment de la conception. L'objectif de cette thèse est en premier lieu la mise au point d'une chaîne de couplage permettant d'évaluer l'impact des interactions fluide-structure sur l'amplitude des oscillations aéroacoustiques présentes au sein du propulseur. Une attention particulière est portée à l'algorithme de couplage entre les solveurs fluide et solide afin d'assurer une bonne conservation de l'énergie à l'interface fluide-structure, point clé dans l'étude d'instabilités. La chaîne numérique ainsi conçue est appliquée à une configuration réduite du moteur à propergol solide d'Ariane 5 dans le cadre de deux études. La première porte sur l'impact des vibrations de la structure sur les d'instabilités aéroacoustiques. L'effet d'un croisement de fréquences des modes propres longitudinaux de la structure et un des modes acoustiques de la chambre de combustion est traité. La seconde étude s'intéresse à l'effet des battements des protections thermiques du propulseur dans l'écoulement. Une structuration de l'écoulement et un net renforcement des ODP sont mis en évidence. / Large solid propellant rocket motors may be subjected to aero-acoustic instabilities arising from a coupling between the burnt gas flow and the acoustic eigenmodes of the combustion chamber. These instabilities lead to large pressure oscillations in the combustion chamber. These pressure oscillations cause vibrations which might jeopardize the payload if they happen to be larger than a certain threshold. Given the size and cost of any single firing test or launch, it is of first importance to rely on numerical tools able to predict these instabilities so that they can be avoided at the design level. The first purpose of this thesis is to build a numerical tool in order to evaluate how the coupling of the fluid flow and the whole structure of the motor influences the amplitude of the aeroacoustic oscillations living inside the rocket. A particular attention was paid to the coupling algorithm between the fluid and the solid solvers in order to ensure the best energy conservation through the interface.The numerical chain is applied to a sub-scaled configuration of Ariane 5 solid rocket motor in two studies. The first relates to the impact of vibration of the structure on aeroacoustic instabilities. The effect of a crossover frequency between the longitudinal modes of the structure and the acoustic modes of the combustion chamber is assessed. The second study examines the effect of thermal protection oscillations in the flow. An increased of the flow organisation and a significant strengthening of pressure oscillations are highlighted.