Le biogaz doit être purifié pour devenir un combustible renouvelable. De nombreux traitements actuels ne sont pas satisfaisants car, pour des raisons de coûts les procédés de séparation privilégiés aboutissent souvent au rejet direct ou indirect du sulfure d'hydrogène (H2S) à l’atmosphère, c’est le cas de la séparation à l’eau sous pression. Les objectifs de la thèse portent d’abord sur l’étude et la modélisation des méthodes connues de séparation de l'hydrogène sulfuré du méthane. Les concentrations typiques varient de 200 à 5000 ppm, et la séparation devra réduire la teneur résiduelle en H2S à moins de 1 ppm. Parallèlement seront étudiées les méthodes de traitement de H2S. Une fois la (ou les) méthode(s) de séparation sélectionnée(s), des essais de validation seront effectués sur un système traitant de l’ordre de 85 Nm3/h de méthane où seront injectées des quantités de H2S variant entre 1 et 100 ppm.Cette thèse requiert des modélisations réalistes sous Aspen Plus® ou sous un code équivalent pour établir a priori des efficacités de séparation selon différentes conditions opératoires et en prenant en compte le paramètre température. L’énergie dépensée pour la séparation effective sera un des critères fort de la comparaison, de même que l’économie de matière.Une approche système est indispensable pour étudier la rétroaction de la méthode de valorisation du H2S sur la ou les méthodes séparatives. A priori c’est aussi l’outil Aspen Plus® ou équivalent qui permettra cette approche système.L’étude du procédé sera menée selon la double approche modélisation et expérimentation, pour l’étude expérimentale des méthodes séparatives, l’échelle du banc sera semi-industrielle et le banc permettra d’étudier les méthodes de séparation jusqu’à -90°C. / Biogas must be purified for becoming a renewable fuel. At now, the most part of the purification techniques are not satisfactory because they imply hydrogen sulfides (H2S) rejection to the atmosphere. One example of these methods is the treatment with high pressure water. The first objective of the thesis is modeling the conventional methods for separating H2S from methane. Typical concentrations of H2S in methane vary from 200 to 5000 pm. Separation methods must decrease the concentration of H2S in methane to less than 1 ppm. At the same time, methods for H2S treatment will be studied.Once the most appropriated separation methods will be selected, some test will be carried out on a pilot plant capable of treating 85 Nm3/h of methane, where quantities of H2S ranging from 1 and 100 ppm will be injected. These tests will allow validating the modeling of the separation process. On the basis of the obtained results, a specific test bench will be conceived and constructed for validating the selected process.The thesis work requires simulating the separation process using the software Aspen Plus® or an equivalent one. The effectiveness of different operative conditions will be tested, varying also the parameter temperature. The energy necessary for the separation will be one of the most important criteria for the comparison, as well as the mass consumption of the different fluids involved in the process.A system approach is fundamental for evaluating the backward effect of the H2S valorization method on the separation techniques. The process simulator (Aspen Plus® or equivalent) will allow the system approach.The study will involve modeling and experimental parts. The experimental part will be carried out taking advantage of a semi-industrial size test bench, allowing studying the separation methods down to -90°C.

La thermodiffusion, également appelé effet Soret, décrit le couplage entre les gradients de température et les flux massiques qui en résultent. Ce phénomène intervient dans de nombreux processus naturels et applications industrielles. En particulier, les réservoirs pétroliers sont sujets à ce phénomène impliquant des fluides multi constituants confinés dans une matrice poreuse et soumis à un gradient de température. Néanmoins, malgré beaucoup des progrès, il existe relativement peu de mesures fiables de ce phénomène et sa modélisation reste largement un problème ouvert. L’objectif principal de cette thèse s’inscrit dans ce cadre, à savoir développer une approche expérimental permettant de fournir des données de références sur la thermodiffusion notamment dans l’optique de quantifier l’effet de la pression sur cette dernière. Ainsi, durant cette thèse, nous avons développé une cellule de thermodiffusion en milieu libre qui permet d’étudier par shadowgraphie les fluctuations de non équilibre induites par effet Soret. L’appareil de mesure a ensuite été utilisé pour étudier deux mélanges binaires représentatifs de fluides pétroliers, à savoir le mélange équimassique tétraline/dodécane (en phase liquide) et le mélange dioxyde de carbone/méthane (en phases gaz et supercritique). En complément, des simulations de dynamique moléculaire ont été réalisées sur le mélange dioxyde de carbone/méthane. Par analyse dynamique des images de shadowgraphie nous avons pu déterminer les coefficients de diffusion et Soret en fonction de la pression pour le mélange tétraline/dodécane. Aux incertitudes près, nous observons une décroissance linéaire avec la pression pour ces coefficients. De plus nous avons observé l’effet du confinement de la cellule sur les fluctuations en très bon accord avec la théorie et les simulations. Pour le mélange dioxyde de carbone/méthane l’analyse dynamique a montré une cinétique difficilement accessible de par les limites physiques et informatiques du dispositif expérimental utilisé. L’analyse statique montre, quant à lui, une croissance rapide de l’amplitude des fluctuations avec la pression jusqu’à un seuil au-delà duquel elle décroît. Sur ce mélange les simulations de dynamique moléculaire ont montré un bon accord avec les prédictions théoriques. / Thermodiffusion, also called the Soret effect, describes the coupling between temperature gradient and resulting fluxes. This phenomenon is involved in a number of natural and industrial processes. In particular, multi components fluids in petroleum reservoirs are subjected to this phenomenon because of the geo-thermal gradient. Nevertheless, in spite of a lot of advances, there are few available data of this phenomenon and the establishment of a theoretical model, able to give a quantitative estimation of these transport coefficients whatever molecules in presence, is still an open question. The principal aim of this thesis is to develop an experimental approach allowing providing reference data on thermodiffusion as a function of the pressure. During this thesis, we developed a high pressure thermodiffusion cell in free medium, enabling us to study concentration non-equilibrium fluctuations induced by the Soret effect by means of shadowgraph optical technique. With this setup we investigated two binary mixtures representatives of petroleum fluids; namely the equimassic tetralin/dodecane mixture in liquid phase and the carbon dioxide/methane mixture in gaseous and super critical state. Furthermore, molecular dynamic simulations on the second mixture were performed. Using a dynamic image analysis, we have measured molecular diffusion and Soret coefficient for the tetralin/dodecane mixture. Within experimental uncertainties, we observed a linear decrease of these coefficients with the pressure. Furthermore, we were able to observe the effect of confinement (finite size effect induced by cell vertical boundary conditions) on fluctuation dynamics, in good agreement with calculations and simulations based on hydrodynamic fluctuation theory on similar solutal Rayleigh number. Concerning the carbon dioxide/methane mixture, the dynamic analysis revealed a kinetic too fast for our experimental apparatus. Conversely, static analysis revealed a rapid increase of the non-equilibrium fluctuation magnitude as a function of the pressure up to a threshold beyond which it decreases. On this mixture, performed molecular dynamic simulations provided results in good agreement with expected theoretical behaviour.

Cette thèse se place dans le contexte de la fusion thermonucléaire pour la production d’énergie, dans le cadre de l’allumage rapide par faisceaux d’électrons chauds. Le travail présenté a pour but de caractériser la source de faisceaux d’électrons rapides, accélérés par lasers intenses (1019≠1020 W/cm2), et leur propagation dans des plasmas denses aussi bien à l’état solide quecomprimé.La première étude présentée avait pour but d’étudier la propagation d’électrons rapides, caractérisés par une densité de courant > 1011 A/cm2, dans des cibles d’aluminium chauffées à la température de Fermi par un choc plan contra-propagatif, qui les comprimait à deux fois la densité du solide. La géométrie de compression plane nous a permis de dissocier les pertesd’énergie dues aux effets résistifs et collisionnels, en comparant des cibles solides et comprimées de masses surfaciques identiques. Nous avons observé pour la première fois une augmentation des pertes d’énergie d’origine résistive dans les échantillons chauffés. La confrontation des données expérimentales avec les simulations, incluant une caractérisation complète de la source électronique, de l’état de compression des cibles et du transport d’électrons, a permis d’étudier l’évolution temporelle de la résistivité du matériau. Elle a notamment permis d’estimer que le pouvoir d’arrêt résistif dans les cibles tièdes et denses est d’amplitude comparable au pouvoir d’arrêt collisionnel.Dans la deuxième étude, nous avons analysé l’accélération et le transport d’électrons rapides produits lors de l’interaction d’un laser à haut contraste avec un cône de cuivre, enchâssé dans un bloc de carbone, et comprimé par un choc plan contra-propagatif. Un système d’imagerie X a permis de visualiser le couplage entre le faisceau laser intense et le cône à différents instants de la compression. Ce diagnostic, donnant accès à la distribution spatiale du faisceau d’électrons chauds, a montré une génération d’électrons dans tout le volume du cône pour des temps supérieurs au temps de débouché de choc au niveau de la pointe. Pour des temps antérieurs, l’interaction se produit à haut contraste, la source est restreinte au niveau de la pointe du cône, et la propagation collimatée des électrons vers l’intérieur de la cible est assurée par les champs magnétiques auto-générés. Ces conclusions ont été obtenues en confrontant les données expérimentales aux simulations.Une caractérisation hydrodynamique de la compression par choc de la cible a été effectuée à l’aide d’une technique de radiographie X, permettant de visualiser la propagation du front de choc dans la cible, sa collision avec la pointe du cône et son glissement le long des parois. Les mesures sont en accord avec des simulations hydrodynamiques. / The framework of this PhD thesis is the inertial confinement fusion for energy production, in the context of the electron fast ignition scheme. The work consists in a characterization of the transport mechanisms of fast electrons, driven by intense laser pulses (1019 ≠ 1020 W/cm2) inboth cold-solid and warm-dense matter.The first goal was to study the propagation of a fast electron beam, characterized by a current density > 1011 A/cm2, in aluminum targets initially heated close to the Fermi temperature by a counter-propagative planar shock. The planar compression geometry allowed us to discriminate the energy losses due to the resistive mechanisms from collisional ones by comparing solid and compressed targets of the same initial areal densities. We observed for the first time a significant increase of resistive energy losses in heated aluminum samples. The confrontation of the experimental data with the simulations, including a complete characterization of the electron source, of the target compression and of the fast electron transport, allowed us to study the time-evolution of the material resistivity. The estimated resistive electron stopping power in a warm-compressed target is of the same order as the collisional one.We studied the transport of the fast electrons generated in the interaction of a high-contrast laser pulse with a hollow copper cone, buried into a carbon layer, compressed by a counterpropagative planar shock. A X-ray imaging system allowed us to visualize the coupling of thelaser pulse with the cone at different moments of the compression. This diagnostic, giving access to the fast electron spatial distribution, showed a fast electron generation in the entire volume of the cone for late times of compression, after shock breakout from the inner cone tip. For earlier times, the interaction at a high-contrast ensured that the source was contained within the cone tip, and the fast electron beam was collimated into the target depth by self-generated magnetic fields. These conclusions were obtained by a confrontation of experimental data to simulation results.The hydrodynamic characterization of the shock-induced target compression was performed using a X-ray point projection radiography technique, allowing to visualize a propagation of the shock front into the target, its collision with the cone tip and its subsequent sliding along the cone walls. The measurements are in agreement with hydrodynamic simulations.

Ce manuscrit présente trois expériences menées dans le cadre de la fusion nucléaire par confinement inertiel. La première expérience porte sur l'étude la propagation d'un faisceau d'électrons rapides dans un plasma pré-comprimé dans le cadre du schémas d'allumage rapide. Deux expériences sur la génération d'onde de choc dans des plasmas ont été menées dans le cadre du schéma d'allumage par choc. La première expérience a été consacre à l'étude le la propagation d'un faisceau d'électrons rapides dans une cible comprimée. L'implosion de la cible avec une géométrie cylindrique a été menée avec l'installation laser GEKKO XII (ILE Osaka, Japon). Le faisceau d'électrons rapides a été injecté en utilisant l'installation laser LFEX (_1019W/cm2). et sa propagation à travers le cylindre comprimé a été observée avec plusieurs diagnostics X. Cette expérience a démontré l'effet d'un guidage du faisceau d'électrons rapides résultant du champ magnétique auto-généré. Par ailleurs, les résultats de cette expérience ont été en bon accord avec deux des simulation numériques. Cette étude a fait l'objet de la publication Approach to the study of fast electron transport in cylindrically imploded targets, Laser and particle Beams, 33, 525-534, (2015). Deux autres expériences ont été réalisées pour l’étude de la propagation de chocs forts créés par laser dans un plasma. Celles-ci ont été réalisées avec différentes installations laser.Dans la première expérience avec le laser Gekko Xll, nous avons observé la création et la propagation de deux ondes de choc successives dans des plasmas d’ablation de CH et de Be. L’objectif de caractériser l’amplification d’un choc transmis par la collision des deux chocscontre propagatifs a été partiellement réalisé. La comparaison des résultats expérimentauxà des simulations hydrodynamiques a permis d’établir une amplification du choc d’unfacteur 2 en pression pour les conditions expérimentales réalisées. Les tirs sur une cible deBe a permis de développer et valider les outils d’exploitation de la propagation de deuxchocs par radiographie X. La deuxième expérience a été réalisée avec laser PHELIX de GSI(Darmstadt, Allemagne). Le but de cette expérience était d’étudier la génération de chocsforts qui ont été utilisés pour étudier l’équations d’état du carbone dans le domaine WDMpertinent pour le planétologie. Les conditions de pression et de densité pour le carboneont été obtenues en déduisant la pression et la vitesse de l’onde de choc des diagnosticschronométriques employés dans cette expérience. Des états du diamant en phase liquidemétallique pour des pressions de l’ordre de 7 Mbar et des températures de 15,000 degrésont été obtenus. / This manuscript presents three experiments conducted as part of a nuclear fusion byinertial confinement. The first experiment is the study of the fast electron beam propagationin a pre-compressed plasma in the fast ignition scheme. Two other experiments about theshock wave generation in plasmas were conducted in the ignition shock pattern.The first experiment was devoted to the study of the fast electron beam transport in a compressed target. The implosion of the target with a cylindrical geometry was carriedout with the GEKKO XII laser facility (ILE Osaka, Japan). The fast electron beam wasgenerated by the LFEX laser ( 1019W/cm2) and its propagation through the compressedcylinder was observed with several X-ray diagnostics. This experiment showed the guidingeffect of the electron beam resulting from self-generated magnetic fields. Furthermore, theresults of this experiment were in good agreement with numerical simulations. This studywas the subject of the publication Approach to the study of fast electron transport incylindrically imploded targets, Laser and Particle Beams, 33,525-534,(2015). Two other experiments were performed to study the propagation of strong shockscreated by lasers in a plasma. These were carried out with different laser systems. In the firstexperiment with the Gekko XII laser, we observed the creation and the propagation of twosuccessive shock waves in an ablation plasma in CH and Be. The objective to characterizethe amplification of a transmitted shock by the collision of two counter propagatingshocks has been partially realized. The comparison of the experimental results with thehydrodynamic simulations enabled us to confirm an amplification of the shock by a factor2 in pressure in the condition of this experiment. The shot with a Be target allowed todevelop and to validate the diagnostic method of X-ray radiography for shock propagation.The second experiment was performed with laser PHELIX GSI (Darmstadt, Germany).The purpose of this experiment was to study the generation of strong shocks. They wereapplied to study the equation of state of carbone in WDM state for the planetology. Thecondition of pressure and density for the carbon were obtained by deducting the pressureand the velocity of the shock wave chronometric diagnostics employed in this experiment.In this experiment, diamond was at the metallic liquid phase with the pressure of 7 Mbarand the temperature of 15,000 degrees

Les suspensions colloïdales magnétiques sont utilisées dans de nombreuses applications pratiques. Toutefois, leurs propriétés structurales et dynamiques, en particulier en présence d'un champ appliqué, restent mal connues. Dans ce but, nous avons adopté une approche qui couple une modélisation par simulation numérique de dynamique brownienne à plusieurs techniques expérimentales sur des dispersions aqueuses préparées «sur mesure». L'emploi de la diffusion de rayonnement, en particulier de neutrons, a permis d'une part, l'élucidation de la structure locale et d'autre part, l'examen de la dynamique de translation à plusieurs échelles. L'effet d'un champ extérieur a été mesuré et évalué parallèlement par la simulation numérique. De plus, la dynamique de rotation a été sondée par la relaxation de biréfringence magnéto-optique. La confrontation des résultats expérimentaux et de simulation numérique complète la compréhension du système, soulignant le rôle des interactions dipolaires.

La thèse présente la modélisation, l'étude expérimentale et le contrôle actif des instabilités aéroacoustiques rencontrées en cavité sous écoulement turbulent à faible nombre de Mach. On propose une formulation problème de stabilité de l'interface fluide séparant deux écoulements uniformes de vitesse différente en integrant les effets acoustiques. Les modes d'instabilité de l'interface sont alors étudiés en fonction du nombre de Mach et de la configuration géométrique. Une maquette comportant une cavité de dimensions réglables est ensuite étudiée en soufflerie à l'aide de mesures de pression. Ces données valident en partie l'approche analytique adoptée. On conçoit alors un dispositif de contrôle des modes d'instabilité, appliqué en particulier dans le cas de leur couplage avec l'acoustique de la veine d'essais. Enfin, un système de vélocimétrie par images de particules synchronisé sur les modes d'oscillation permet de valider l'étude théorique et la stratégie de contrôle.

La plaine alluviale du Haut Sébaou, dépression synclinale à coeur miocène, d'orientation est*ouest est bordée au nord par les chaînons côtiers d'Afir et au sud par la chaîne calcaire du Djurdjura. Elle fait partie intégrante du Tell septentrional. Les faciès géologiques portent l'empreinte de l'orogenèse alpine. Les reconstitutions paléogéographiques révèlent l'antécédence au Pliocène du réseau hydrographique. Le remplissage alluvial est formé de six terrasses emboitées et étagées d'extension et d'épaisseur réduite réparties symétriquement sur les deux rives de l'oued. L'étude géophysique, complétée par les sondages mécaniques, définit un aquifère homogène de résistivités moyennes. L'étude hydroclimatologique confirme le type de climat méditerranéen du bassin versant. L'évaluation des paramètres: Pluviométrie/Evapotranspiration/ Ruissellement, a permis de proposer un essai de bilan hydrologique. L'étude piézométrique met en évidence la superposition des écoulements souterrains et superficiels: l'Oued Sebaou draine les eaux de la nappe alluviale. L'estimation des paramètres hydrodynamiques (T,S,K) au moyen d'essais de pompage a permis de sélectionner les zones propices à l'exploitation. L'utilisation conjointe des variations naturelles du niveau de la nappe, des paramètres hydrodynamiques et des débits de sortie au seuil du Belloua, nous a permis d'estimer les réserves régulatrices et de fixer un ordre de grandeur du débit d'exploitation rationnelle de l'aquifère. L'étude physico-chimique a mis en évidence deux types de faciès hydro-chimiques: bicarbonaté calcique-sulfaté et magnésien et bicarbonaté calcique-chloruré et magnésien. Les eaux ont une minéralisation moyenne homogène quoiqu'elles soient le siège d'échange de bases sous l'influence d'apports des versants. Elles présentent une bonne potabilité chimique, malgré leurs caractères d' eaux dures et incrustantes

Dans cette thèse, nous présentons une nouvelle méthode smoothed particle hydrodynamics (SPH) pour la résolution des équations de Navier-Stokes incompressibles, même en présence des forces singulières. Les termes de sources singulières sont traités d'une manière similaire à celle que l'on retrouve dans la méthode Immersed Boundary (IB) de Peskin (2002) ou de la méthode régularisée de Stokeslets (Cortez, 2001). Dans notre schéma numérique, nous mettons en oeuvre une méthode de projection sans pression de second ordre inspirée de Kim et Moin (1985). Ce schéma évite complètement les difficultés qui peuvent être rencontrées avec la prescription des conditions aux frontières de Neumann sur la pression. Nous présentons deux variantes de cette approche: l'une, Lagrangienne, qui est communément utilisée et l'autre, Eulerienne, car nous considérons simplement que les particules SPH sont des points de quadrature où les propriétés du fluide sont calculées, donc, ces points peuvent être laissés fixes dans le temps. Notre méthode SPH est d'abord testée à la résolution du problème de Poiseuille bidimensionnel entre deux plaques infinies et nous effectuons une analyse détaillée de l'erreur des calculs. Pour ce problème, les résultats sont similaires autant lorsque les particules SPH sont libres de se déplacer que lorsqu'elles sont fixes. Nous traitons, par ailleurs, du problème de la dynamique d'une membrane immergée dans un fluide visqueux et incompressible avec notre méthode SPH. La membrane est représentée par une spline cubique le long de laquelle la tension présente dans la membrane est calculée et transmise au fluide environnant. Les équations de Navier-Stokes, avec une force singulière issue de la membrane sont ensuite résolues pour déterminer la vitesse du fluide dans lequel est immergée la membrane. La vitesse du fluide, ainsi obtenue, est interpolée sur l'interface, afin de déterminer son déplacement. Nous discutons des avantages à maintenir les particules SPH fixes au lieu de les laisser libres de se déplacer. Nous appliquons ensuite notre méthode SPH à la simulation des écoulements confinés des solutions de polymères non dilués avec une interaction hydrodynamique et des forces d'exclusion de volume. Le point de départ de l'algorithme est le système couplé des équations de Langevin pour les polymères et le solvant (CLEPS) (voir par exemple Oono et Freed (1981) et Öttinger et Rabin (1989)) décrivant, dans le cas présent, les dynamiques microscopiques d'une solution de polymère en écoulement avec une représentation bille-ressort des macromolécules. Des tests numériques de certains écoulements dans des canaux bidimensionnels révèlent que l'utilisation de la méthode de projection d'ordre deux couplée à des points de quadrature SPH fixes conduit à un ordre de convergence de la vitesse qui est de deux et à une convergence d'ordre sensiblement égale à deux pour la pression, pourvu que la solution soit suffisamment lisse. Dans le cas des calculs à grandes échelles pour les altères et pour les chaînes de bille-ressort, un choix approprié du nombre de particules SPH en fonction du nombre des billes N permet, en l'absence des forces d'exclusion de volume, de montrer que le coût de notre algorithme est d'ordre O(N). Enfin, nous amorçons des calculs tridimensionnels avec notre modèle SPH. Dans cette optique, nous résolvons le problème de l'écoulement de Poiseuille tridimensionnel entre deux plaques parallèles infinies et le problème de l'écoulement de Poiseuille dans une conduite rectangulaire infiniment longue. De plus, nous simulons en dimension trois des écoulements confinés entre deux plaques infinies des solutions de polymères non diluées avec une interaction hydrodynamique et des forces d'exclusion de volume. / In this thesis we develop a new smoothed particle hydrodynamics (SPH) method suitable for solving the incompressible Navier-Stokes equations, even with singular forces. Singular source terms are handled in a manner similar to that in the immersed boundary (IB) method of Peskin (2002) or in the method of regularized Stokeslets (Cortez, 2001). The numerical scheme implements a second-order pressure-free projection method due to Kim and Moin (1985) and completely obviates the difficulties that may be faced in prescribing Neumann pressure boundary conditions. We present two variants of this approach, one Langrangian which is commonly used and one Eulerian, simply because we consider that the SPH particles are quadrature points on which the fluid properties are calculated, therefore, these points can be kept fixed in time. The proposed SPH method is first tested on the planar start-up Poiseuille problem and a detailed error analysis is performed. For this problem, the results are similar whether the SPH particles are free to move or fixed on a regular grid. Our hybrid SPH-IB method is then used to calculate the dynamics of a stretched immersed elastic membrane. The membrane is represented by a cubic spline along which the tension in the membrane is computed and transmitted to the surrounding fluid. The Navier-Stokes equations with singular force due to the membrane are then solved to determine the velocity of the fluid in which the membrane is immersed. The fluid velocity thus obtained is interpolated on the interface, to determine its displacement. We discuss the advantages, in this problem, of fixing the SPH particles, rather than allowing them to move with the fluid. A new coupled Brownian dynamics-SPH method for the computation of confined flows of non-dilute polymer solutions with full hydrodynamic interaction and excluded volume forces is next presented. The starting point for the algorithm is the system of coupled Langevin equations for polymer and solvent (CLEPS) (see Oono and Freed (1981) and Öttinger and Rabin (1989), for example) describing, in the present case, the microscopic dynamics of a flowing polymer solution with a bead-spring representation of the macromolecules. Numerical tests of some two-dimensional channel flows reveal that use of a second-order projection scheme coupled with fixed SPH quadrature points leads to second-order velocity convergence and almost second-order pressure convergence, provided that the solution is sufficiently smooth. In the case of large-scale dumbbell and bead-spring chain calculations, an appropriate scaling of the number of grid points as a function of the number of beads N ensures, in the absence of excluded volume forces, that the cost of our algorithm is O(N) flops. Finally, we begin calculations in three dimensions with our SPH model. To this end, we solve in three dimensions the problem of Poiseuille flow between two infinite and parallel plates and the problem of Poiseuille flow in a rectangular infinitely long duct. In addition, we carry out three dimensional computations of confined flows of non-dilute polymer solutions with full hydrodynamic interaction and excluded volume forces.

Cette thèse vise à modéliser les instabilités hydrodynamiques générant des détachements tourbillonnaires pariétaux (ou VSP) responsables des Oscillations De Pression dans les moteurs à propergol solide longs et segmentés par interaction avec l’acoustique du moteur. Ces instabilités sont modélisées en tant que modes de stabilité linéaire globaux de l’écoulement d’un conduit à parois débitantes. En supposant que les structures pariétales émergent d’une perturbation de l’écoulement de base, des modes discrets et indépendants du maillage utilisé sont calculés. Dans ce but, une discrétisation par collocation spectrale multi-domaine est implémentée dans un solveur parallèle afin de s’affranchir de la croissance polynomiale des fonctions propres et de la présence de couches limites. Les valeurs propres ainsi calculées dépendent explicitement des frontières du domaine, à savoir la position de la perturbation et celle de la sortie, et sont ensuite validées par simulation numérique directe. On montre alors qu’elles permettent bien de décrire la réponse à une perturbation initiale de l’écoulement modifié par une rupture de débit pariétale. Ensuite, la simulation d’une réponse forcée par l’acoustique se fait sous forme de structures tourbillonnaires dont les fréquences discrètes sont en accord avec celles des modes de stabilité. Ces structures sont réfléchies en ondes de pression de même fréquences remontant l’écoulement. Finalement, la simulation numérique et la théorie de la stabilité permettent de montrer que le VSP, dont la réponse est linéaire vis-à-vis d’un forçage compressible comme l’acoustique, est le phénomène moteur des Oscillations De Pression. / The current work deals with the modeling of the hydrodynamic instabilities that play a major role in the triggering of the Pressure Oscillations occurring in large segmented solid rocket motors. These instabilities are responsible for the emergence of Parietal Vortex Shedding (PVS) and they interact with the boosters acoustics. They are first modeled as eigenmodes of the internal steady flowfield of a cylindrical duct with sidewall injection within the global linear stability theory framework. Assuming that the related parietal structures emerge from a baseflow disturbance, discrete meshindependant eigenmodes are computed. In this purpose, a multi-domain spectral collocation technique is implemented in a parallel solver to tackle numerical issues such as the eigenfunctions polynomial axial amplification and the existence of boundary layers. The resulting eigenvalues explicitly depend on the location of the boundaries, namely those of the baseflow disturbance and the duct exit, and are then validated by performing Direct Numerical Simulations. First, they successfully describe flow response to an initial disturbance with sidewall velocity injection break. Then, the simulated forced response to acoustics consists in vortical structures wihich discrete frequencies that are in good agreement with those of the eigenmodes. These structures are reflected into upstream pressure waves with identical frequencies. Finally, the PVS, which response to a compressible forcing such as the acoustic one is linear, is understood as the driving phenomenon of the Pressure Oscillations thanks to both numerical simulation and stability theory.

Dans le but de réduire le bruit de jet, source principale de nuisance sonore au décollage d'un avion, une compréhension fine des mécanismes aéroacoustiques mis en jeu est nécessaire. Les structures cohérentes de grande échelle se développant dans la couche de mélange d'un jet semblent responsables d'une part importante du bruit observé en champ lointain, surtout dans les basses fréquences. Une approche permettant d'étudier ces structures turbulentes est fournie par la théorie de stabilité, notamment au moyen des équations de stabilité parabolisées (PSE). L'étude de ces ondes d'instabilité est alors complémentaire d'autres approches (LES ou expériences), puisqu'elle permet de mettre en évidence la nature et la dynamique de ces structures, également présentes dans les résultats de simulations ou de mesures.Au cours de ces travaux de thèse, nous nous sommes intéressés aux structures cohérentes se développant dans des jets à double flux étudiés au cours du projet européen CoJeN (Coaxial Jet Noise). En particulier, nous avons exploité une base de données issues de mesures de fluctuations de pression réalisées en champ proche et en champ lointain de ces jets. Nous avons alors pu comparer les résultats de notre modélisation PSE à ces mesures en périphérie immédiate du jet, confirmant ainsi la pertinence d’un tel modèle, même dans des configurations aussi complexes. De plus, le calcul du rayonnement acoustique en champ lointain engendré par les fluctuations de pression modélisées nous a permis de faire des comparaisons directes avec les niveaux et les directivités mesurés. Nous avons ainsi pu mettre en évidence quantitativement la contribution de ces structures turbulentes de grande échelle au bruit total rayonné par le jet. / Increasingly stringent aircraft noise regulations require the development of innovative noise reduction strategies. Jet noise is a dominant acoustic component during take-off and a fine understanding of the underlying aeroacoustics mechanisms is then necessary. Large-scale coherent structures that develop in the mixing layer of jets appear to be the dominant acoustic source responsible for the lowfrequency far-field noise observed at low emission angles. A stability analysis based on the parabolized stability equations (PSE) is a suitable tool for studying these coherent structures, revealing the nature and the dynamics of the fluctuations obtained by simulations or experiments. The present work is focused on coherent structures developing in the two mixing layers of dualstream jets studied in the course of the European project CoJeN (Coaxial Jet Noise). In particular, pressure fluctuations measurements acquired in the near and far fields of two coaxial jets have been thoroughly analyzed. A direct comparison of these experimental results with linear PSE calculations has been performed in the vicinity of the jets, referred to as the linear-hydrodynamic region, confirming the relevance of the approach even in such complex industrial configurations. Furthermore, the acoustic projection to the far-field of the wavepackets issued by this model and calibrated in the near-field allows a direct comparison of the acoustic levels and directivity with far field sound measurements. A quantitative assessment of the contribution of the instability waves to the total jet noise measured has therefore been obtained.