Les pentes ont toujours été un sujet de discussion pour les ingénieurs en géotechnique. Leurs stabilités étant d’une importance majeure, puisque leur effondrement peut causer la perte d’inestimables et précieuses vies humaines. C’est l’une des raisons qui nous a poussé à étudier la sécurité des pentes. Dans ce projet, on a essayé de se concentrer sur les sols fins et plus précisément les argiles, pour voir l’effet de la variation de l’indice de plasticité sur la stabilité dynamique des pentes et des remblais construits sur des dépôt d’argile. On a aussi essayé d’étudier l’effet de la variation des zones sismiques au Québec (zone 3, 4 et 5) sur la stabilité des pentes face à des séismes de différentes caractéristiques. Alors, pour ce faire, plus que 2000 simulations numériques ont été réalisées sur un logiciel de calcul dynamique, FLAC, afin d’étudier le comportement de différentes pentes de 10 et 20 m d’épaisseur de dépôt (avec une hauteur de pente constante de 10 m) qui présentent une augmentation de cohésion de couche en sous-couche (delta) de 1,5, 2 et 2,5. La méthode pseudo-spectrale (récemment développée à l’université de Sherbrooke), qui se présente comme remplaçante à la méthode pseudostatique, a été utilisée pour comparer ses résultats avec ceux de la littérature et ensuite pouvoir la valider. L’effet de l’ajout d’une berme sur la stabilité dynamique a été également étudié.

On procède à la simulation numérique directe de l'interaction turbulence-sans-cisaillement/surface-plane, par une méthode pseudo-spectrale. L'originalité de la configuration tient au fait que la turbulence est entretenue à distance par un forçage aléatoire localisé au voisinage d'un plan parallèle à la surface. L'écoulement est donc statistiquement stationnaire et la couche de surface continûment alimentée par diffusion turbulente. Les évolutions des quantités statistiques et des bilans des tensions de Reynolds sont données et permettent de préciser la structure de la couche d'interaction. Une méthode originale d'analyse et de quantification des structures élémentaires participant au transfert d'énergie intercomposantes est proposée et testée. Cette grille d'analyse permet une meilleure compréhension du contenu de la corrélation pression-déformation en présence d'une surface imperméable.

La description magnétohydrodynamique est utilisée pour étudier les plasmas de fusion par confinement magnétique dans deux configurations: tokamak et reversed field pinch. Une méthode de Fourier pseudo-spectrale et une technique de pénalisation en volume sont employées pour résoudre les équations. La méthode de pénalisation permet d'introduire des conditions aux limites de Dirichlet et donc de faire varier facilement la géométrie considérée. Les simulations dans des géométries toroïdales de type tokamak montrent l'apparition spontanée de vitesses. Une importante composante toroïdale se développe si le système est peu dissipatif. Il est aussi montré que la brisure de symétrie dans la forme de la section du tore fait apparaitre un moment angulaire toroïdal. Pour le Reversed Field Pinch on montre l'émergence de structures hélicoïdales. La forme de ces structures varie en fonction des coefficients de transport ainsi que du paramètre de pincement du champ magnétique imposé. Pour compléter l'étude on compare les résultats du tore aux calculs dans un cylindre périodique. Les différences dans la dynamique des deux cas sont mises en avant. Finalement les simulations sont confrontées à des expériences et un meilleur accord est observé entre simulation et expérience pour la géométrie toroïdale que pour la géométrie cylindrique.

La description magnétohydrodynamique est utilisée pour étudier les plasmas de fusion par confinement magnétique dans deux configurations: tokamak et reversed field pinch. Une méthode de Fourier pseudo-spectrale et une technique de pénalisation en volume sont employées pour résoudre les équations. La méthode de pénalisation permet d’introduire des conditions aux limites de Dirichlet et donc de faire varier facilement la géométrie considérée. Les simulations dans des géométries toroïdales de type tokamak montrent l’apparition spontanée de vitesses. Une importante composante toroïdale se développe si le système est peu dissipatif. Il est aussi montré que la brisure de symétrie dans la forme de la section du tore fait apparaitre un moment angulaire toroïdal. Pour le Reversed Field Pinch on montre l’émergence de structures hélicoïdales. La forme de ces structures varie en fonction des coefficients de transport ainsi que du paramètre de pincement du champ magnétique imposé. Pour compléter l’étude on compare les résultats du tore aux calculs dans un cylindre périodique. Les différences dans la dynamique des deux cas sont mises en avant. Finalement les simulations sont confrontées à des expériences et un meilleur accord est observé entre simulation et expérience pour la géométrie toroïdale que pour la géométrie cylindrique. / A magnetohydrodynamic description is used to study magnetic fusion plasmas in two different configurations: tokamak and reversed field pinch. A Fourier pseudo-spectral method with a volume penalization technique are used to solve the system of equations. The penalization method is used to introduce Dirichlet boundary conditions and it al- lows to easily modify the consider geometry. The simulations of a tokamak configuration in a toroidal geometry show the spontaneous appearance of velocities. These velocities are dominated by their toroidal component if the system is little dissipative. It is also shown that the symmetry breaking of the cross section of the torus causes a toroidal angular momentum to develop. For the Reversed Field Pinch configuration we show the appearance of helical structures. The shape of these structures varies with the value of the transport coefficients and with the pinch ratio parameter of the imposed magnetic field. To complete the study, we compare the results of simulations obtained in toroidal and in periodic cylindrical geometries. The differences in the dynamics of these two cases are highlighted. Finally, simulations are compared to experimental data and a significant better agreement is observed between the simulation and the experiment for the toroidal geometry than for the cylindrical case.

L'approximation dite " anélastique " permet de filtrer les ondes acoustiques grâce à un développement asymptotique deséquations de Navier-Stokes, réduisant ainsi le pas en temps moyen, lors de la simulation numérique du développement d'instabilités hydrodynamiques. Ainsi, les équations anélastiques sont établies pour un mélange de deux fluides pour l'instabilité de Rayleigh-Taylor. La stabilité linéaire de l'écoulement est étudiée pour la première fois pour des fluides parfaits, par la méthode des modes normaux, dans le cadre de l'approximation anélastique. Le problème de Stokes issu des équations de Navier-Stokes sans les termes non linéaires (une partie de la poussée d'Archiméde est prise en compte) est défini ; l'éllipticité est démontrée, l'étude des modes propres et l'invariance liée à la pression sont détaillés. La méthode d'Uzawa est étendue à l'anélastique en mettant en évidence le découplage des vitesses en 3D, le cas particulier k = 0 et les modes parasites de pression. Le passage au multidomaine a permis d'établir les conditions de raccord (raccord Co de la pression sans condition aux limites physiques). Les algorithmes et l'implantation dans le code AMENOPHIS sont validés par les comparaisons de l'opérateur d'Uzawa développé en Fortran et à l'aide de Mathematica. De plus des résultats numériques ont été comparés à une expérience avec des fluides incompressibles. Finalement, une étude des solutions numériques obtenues avec les options anélastique et compressible a été menée. L'étude de l'influence de la stratification initiale des deux fluides sur le développement de l'instabilité de Rayleigh-Taylor est amorcée.

L’approximation dite « anélastique » permet de filtrer les ondes acoustiques grâce à un développement asymptotique deséquations de Navier-Stokes, réduisant ainsi le pas en temps moyen, lors de la simulation numérique du développement d’instabilités hydrodynamiques. Ainsi, les équations anélastiques sont établies pour un mélange de deux fluides pour l’instabilité de Rayleigh-Taylor. La stabilité linéaire de l’écoulement est étudiée pour la première fois pour des fluides parfaits, par la méthode des modes normaux, dans le cadre de l’approximation anélastique. Le problème de Stokes issu des équations de Navier-Stokes sans les termes non linéaires (une partie de la poussée d’Archiméde est prise en compte) est défini ; l’éllipticité est démontrée, l’étude des modes propres et l’invariance liée à la pression sont détaillés. La méthode d’Uzawa est étendue à l’anélastique en mettant en évidence le découplage des vitesses en 3D, le cas particulier k = 0 et les modes parasites de pression. Le passage au multidomaine a permis d’établir les conditions de raccord (raccord Co de la pression sans condition aux limites physiques). Les algorithmes et l’implantation dans le code AMENOPHIS sont validés par les comparaisons de l’opérateur d’Uzawa développé en Fortran et à l’aide de Mathematica. De plus des résultats numériques ont été comparés à une expérience avec des fluides incompressibles. Finalement, une étude des solutions numériques obtenues avec les options anélastique et compressible a été menée. L’étude de l’influence de la stratification initiale des deux fluides sur le développement de l’instabilité de Rayleigh-Taylor est amorcée. / The « anelastic » approximation allows us to filter the acoustic waves thanks to an asymptotic development of the Navier-Stokes equations, so increasing the averaged time step, during the numerical simulation of hydrodynamic instabilitiesdevelopment. So, the anelastic equations for a two fluid mixture in case of Rayleigh-Taylor instability are established.The linear stability of Rayleigh-Taylor flow is studied, for the first time, for perfect fluids in the anelastic approximation.We define the Stokes problem resulting from Navier-Stokes equations without the non linear terms (a part of the buoyancyis considered) ; the ellipticity is demonstrated, the eigenmodes and the invariance related to the pressure are detailed.The Uzawa’s method is extended to the anelastic approximation and shows the decoupling speeds in 3D, the particular casek = 0 and the spurius modes of pressure. Passing to multidomain allowed to establish the transmission conditions.The algorithms and the implementation in the existing program are validated by comparing the Uzawa’s operator inFortran and Mathematica langages, to an experiment with incompressible fluids and results from anelastic and compressiblenumerical simulations. The study of the influence of the initial stratification of both fluids on the development of the Rayleigh-Taylor instability is initiated.