Dans le présent travail, nous étudions l'impact sur un obstacle rigide d'un écoulement transitoire à surface libre de fluide viscoplastique. Cette étude est conduite numériquement à l'aide de la méthode SPH (Smoothed Particle Hydrodynamics), en y intégrant le modèle rhéologique de Herschel-Bulkley. Le code employé est adapté à nos besoins et validé sur des cas test classiques. Les caractéristiques locales de l'écoulement à proximité de l'obstacle sont analysées et deux régimes d'impact sont mis en évidence en fonction de la pente d'écoulement. L'étude des pressions exercées sur l'obstacle, conduite spatialement et temporellement en fonction de ces régimes d'impact, nous permet de mettre en évidence les rôles respectifs des composantes gravitationnelle et cinétique de la pression. Nos résultats sont comparés systématiquement à des résultats expérimentaux issus de travaux précédents et sont cohérents avec ces derniers. Une étude comparative de nos écoulements de fluide viscoplastique avec des écoulements de matériau granulaires de propriétés similaires nous conduit à mettre en évidence des caractéristiques communes entre les deux matériaux. / In this work, we study the impact of a transient free-surface flow of viscoplastic fluid on a rigid obstacle. This study is conducted numerically using the SPH (Smoothed Particle Hydrodynamics) method, and the Herschel-Bulkley rheological model. The SPH code is adapted to our needs and validated on classic benchmarks. The local characteristics of the flow near the obstacle are analysed and two impact regimes are highlighted depending on the slope angle. By studying of the pressure exerted on the obstacle, both spatially and temporally, with regards to these impact regimes, we evidence the respective roles of the gravitational and kinetic components of the pressure. Our results are systematically compared with experimental data from a previous work and are shown to be consistent. A comparative study conducted on both our viscoplastic flows and flows of granular material of similar properties highlights common characteristics of the two materials.

Le travail présenté intègre l'utilisation des vibrations mécaniques dans les procédés de mise en forme à froid des matériaux de comportement viscoplastique. Les avantages de procédé de mise en forme assisté par vibrations se mesurent par une réduction de l'intensité de la force de mise en forme, une diminution de la contrainte d'écoulement, la réduction du frottement entre l'outil et lopin, une distribution uniforme de la déformation viscoplastique et une amélioration de l'état de surface de la pièce finale. A cet effet, un modèle analytique avec une seule source de vibration a été développé pour le procédé de forgeage. Les vibrations sinusoïdale et triangulaire ont été intégrées dans le modèle et il se trouve que les vibrations triangulaire est plus avantageux comparer à les vibrations sinusoïdal a cause du contrôle sur le rapport cyclique. On observe également que les paramètres principaux qui contrôlent la réduction de la force de forgeage sont : le rapport de vitesse définie par l'amplitude de la vitesse de vibration et la vitesse moyenne de l'outillage et la sensibilité à la vitesse de déformation. Le modèle analytique a été validé au moyen des simulations éléments finis réalisés avec le code de calcul FORGE2011 ® et les essais expérimentaux en dispositif expérimentale de pot vibrant a été conçu et les essais ont été réalisé sur les machines Lloyd LR30K ou ZWICK / Roell 1200. L'étude a ensuite été enrichi avec les plusieurs sources de vibrations. L'étude cinématique montre que les vibrations multiples générés une onde progressive. Cette nouvelle cinématique a été intégrée dans un modèle numérique sur le code de simulation par éléments finis FORGE2011 ® pour étudier le forgeage assisté par plusieurs sources des vibrations. L'analyse de résultats montre une réduction considérable de la force de forgeage et de la contrainte de l'écoulement ainsi qu'une modification de la loi de frottement générer par le changement de la direction du vecteur de vitesse.

La résolution analytique du problème de la stabilité d'une suspension solide (particules rigides de forme sphérique) dans le système de Taylor-Couette cylindrique a été menée. On s'est basé sur les travaux de Ali and Lueptow (2002) pour formuler les équations régissant la stabilité de l'écoulement dans le cadre d'une théorie linéaire. Ces équations sont valables dans l'approximation du faible espace annulaire et ont pour but la prévision de l'instabilité primaire. A cet effet, nous avons utilisé une méthode variationnelle telle que la méthode de Galerkin pour résoudre le problème aux valeurs propres conduisant à établir le diagramme de stabilité liée au nombre d'onde au voisinage de l'état critique du développement de la première instabilité. Tout d'abord, on a cherché à mettre au point les calculs dans le cas de l'écoulement de Taylor-Couette classique en se référant aux travaux de Ali and Lueptow (2002). Ensuite on a procédé à la résolution systématique des équations du mouvement et l'on cherche à déterminer le critère d'apparition des instabilités en présence de particules en suspension et l'on détermine simultanément les paramètres de couplage entre forces d'interaction liquide-solide. L'ensemble des travaux ainsi réalisés permettront de lever la contradiction fondamentale entre la théorie et l'expérience. L'étude expérimentale a permis d'analyser les effets de la concentration des particules en suspension (disques) et du rapport d'aspect radial ' sur l'apparition des instabilités dans le système de Taylor-Dean. Le dispositif expérimental utilisé est constitué d'un cylindre intérieur tournant et le cylindre extérieur est maintenu fixe. Le comportement rhéologique du fluide utilisé est viscoplastique obéissant au modèle de Herschel Bulkley. L'apparition des instabilités est examinée en utilisant une technique de visualisation. Pour une concentration donnée dans l'intervalle étudié, la nature des structures apparaissant dans le système d'écoulement dépend ', alors que pour une valeur donnée de ' dans l'intervalle étudié, la valeur du nombre de Taylor critique Tac dépend de la concentration des particules. Nous avons obtenu que le nombre de Taylor critique Tac correspondant au déclenchement de la première instabilité évolue non linéairement en fonction de '. De plus, nous avons examiné expérimentalement les effets de limitation axiale (effet de bords) sur le déclenchement des instabilités dans le système de Taylor-Dean. Les résultats obtenus montrent que les bords tournants n'affectent pas le type de structures qui apparaissent dans le système d'écoulement de Taylor-Dean. Cependant, ils influencent le seuil critique d'apparition des instabilités qui est marquée par des valeurs élevées du nombre de Taylor critique pour des bords tournants ce qui indique un effet stabilisant des bords mobiles.

Dans le présent travail, nous étudions l'impact sur un obstacle rigide d'un écoulement transitoire à surface libre de fluide viscoplastique. Cette étude est conduite numériquement à l'aide de la méthode SPH (Smoothed Particle Hydrodynamics), en y intégrant le modèle rhéologique de Herschel-Bulkley. Le code employé est adapté à nos besoins et validé sur des cas test classiques. Les caractéristiques locales de l'écoulement à proximité de l'obstacle sont analysées et deux régimes d'impact sont mis en évidence en fonction de la pente d'écoulement. L'étude des pressions exercées sur l'obstacle, conduite spatialement et temporellement en fonction de ces régimes d'impact, nous permet de mettre en évidence les rôles respectifs des composantes gravitationnelle et cinétique de la pression. Nos résultats sont comparés systématiquement à des résultats expérimentaux issus de travaux précédents et sont cohérents avec ces derniers. Une étude comparative de nos écoulements de fluide viscoplastique avec des écoulements de matériau granulaires de propriétés similaires nous conduit à mettre en évidence des caractéristiques communes entre les deux matériaux. / In this work, we study the impact of a transient free-surface flow of viscoplastic fluid on a rigid obstacle. This study is conducted numerically using the SPH (Smoothed Particle Hydrodynamics) method, and the Herschel-Bulkley rheological model. The SPH code is adapted to our needs and validated on classic benchmarks. The local characteristics of the flow near the obstacle are analysed and two impact regimes are highlighted depending on the slope angle. By studying of the pressure exerted on the obstacle, both spatially and temporally, with regards to these impact regimes, we evidence the respective roles of the gravitational and kinetic components of the pressure. Our results are systematically compared with experimental data from a previous work and are shown to be consistent. A comparative study conducted on both our viscoplastic flows and flows of granular material of similar properties highlights common characteristics of the two materials.

Le modèle de couche mince intégré sur l'épaisseur, Saint-Venant, utilisé classiquement pour simuler la propagation de laves torrentielles et coulées boueuses, repose sur plusieurs approximations concernant la forme des profils de vitesse en zones non-uniformes. Il est pourtant nécessaire d'utiliser ce type de modélisation, comme outil d'aide à la gestion des risques liés aux laves torrentielles. Nous proposons d'éprouver ses hypothèses, en observant une zone fortement non-uniforme, le front de coulées à surface libre et le champ de vitesse à l'intérieur de cette zone.En améliorant notre connaissance de l'évolution de la forme des profils de vitesse (de la dynamique interne) au front de coulées, nous cherchons à améliorer les modèles de couche mince. Cette thèse porte donc sur l'étude de la dynamique interne au front d'écoulements à surface libre de fluides newtoniens et viscoplastiques.Nous avons utilisé le dispositif du canal à fond mobile qui permet de générer des coulées stationnaires dans le référentiel de l'observateur au moyen d'un fond mobile remontant vers l'amont. Nous avons réalisé un travail technique sur ce canal et sur l'analyse des images pour pouvoir mesurer les champs de vitesse à haute résolution spatiale aux fronts de coulées à surface libre de fluides viscoplastiques. L'étude des fluides newtoniens a aussi été réalisée afin de valider les modèles et éprouver le dispositif expérimental.Nous avons comparé les résultats expérimentaux aux solutions théoriques de deux modèles de couche mince adaptés à la rhéologie de Herschel-Bulkley : le modèle classique de la lubrification, à la base du modèle de Saint-Venant et un modèle consistant à l'ordre 1 développé dans cette thèse. Le modèle consistant d'ordre 1 est la somme du modèle à l'ordre 0 (la lubrification) et de termes correctifs qui proviennent des contraintes normales et des termes d'inertie. Dans le cadre de notre configuration du fond mobile remontant vers l'amont, il est possible de déduire la forme du front en cherchant une solution de type onde progressive, sans passer par un modèle intégré dans l'épaisseur.Pour les fluides viscoplastiques, la structure classique du profil de vitesse, avec une zone cisaillée surmontée d'un plug non cisaillé est bien reconnaissable sur nos profils de vitesse en zone uniforme, et en zone faiblement variée. Mais à l'approche du front, cependant, la vitesse de surface augmente, les profils de vitesse expérimentaux deviennent cisaillés sur toute l'épaisseur, conduisant à la disparition du plug à proximité de la ligne de front.Le modèle de lubrification prédit l’existence d'un plug dans le front jusqu'à la ligne de contact, ce qui n'est pas observé expérimentalement. La vitesse de surface du modèle de lubrification augmente à l'approche du front, mais est largement sous-estimée par rapport à la vitesse de surface mesurée. Les vitesses de surface prédites par le modèle d'ordre 1 augmentent plus drastiquement au front, en meilleur accord avec les mesures que le modèle de lubrification. Pour certaines configurations expérimentales l'accord est même très bon. Remarquablement, le cisaillement des profils de vitesse à l'approche du front, observé expérimentalement, est aussi prédit par le modèle d'ordre 1.Les profils de vitesse présentent donc une évolution au front de coulées viscoplastiques en contradiction avec les hypothèses du modèle de Saint-Venant. Le modèle consistant d'ordre 1 permet d'améliorer les prédictions. Un modèle intégré dans l'épaisseur de type Saint-Venant basé sur les développements consistants d'ordre 1 est alors calculé, car il constitue l'étape nécessaire avant d'être intégré dans un outil de simulation opérationnel. / A depth-averaged model based on the thin-layer assumption, called Saint-Venant (Shallow-Water), is classically used to simulate the propagation and the spreading of debris and mud flows. It is based on several approximations concerning the shape of the velocity profile in non-uniform zones. We propose to test these hypotheses, examining a strongly non-uniform zone, the front of free-surface viscoplastic flows and the velocity field within this zone. By improving our knowledge about the internal dynamics in the front zone, we seek to improve the thin-layer models. This thesis therefore focuses on the study of the internal dynamics within the front of viscoplatic free-surface flows.We used the moving conveyor belt to generate stationary flows. We carried out a technical work on this set-up, and specific analysis of images obtained from the high-speed camera, in order to be able to measure velocity fields with a high resolution. The study of a Newtonian fluid was also carried out in order to validate the lubrication model and the experimental device.We compared experimental results to theoretical solutions of two thin-layer models taking into account the Herschel-Bulkley rheology: the classical model of lubrication, which is at the base of Saint-Venant model, and a consistent first-order model specifically developed in this thesis.The first-order model is equal to the zero-order model (lubrication), plus corrective terms derived from the normal stresses and inertia terms.In this study, for the purpose of comparison with our experimental results, we are interested in travelling-wave solutions. We are able to solve the shape of the front without using a depth-averaged model.Far from the front, experimental velocity profiles clearly display the characteristic 2-layer structure predicted by the lubrication solution, with constant values close to the free-surface (plug) and a sheared layer underneath. Closer to surge tip, the shape of experimental longitudinal velocity profilesthen begins to differ from the theoretical prediction. The 2-layer structure tends to disappear, and the profiles display shear across the whole depth ofthe flow. In this tip region, surface velocity also appears to increase faster than its theoretical counterpart. Surface velocity predicted by the first-order model increase more drastically in the tip region, in better agreement with the measurements than the lubrication model. The first-order model predicts a sheared velocity profile when approaching the front, as observed experimentally.The consistent first-order model then provides better predictions about internal dynamics than lubrication model. A depth-integrated model like Saint-Venant, based on consistent first-order developments is then calculated, as a first step before being integrated into an operational simulation tool.

Dans le cadre de cette thèse, nous avons étudié l'interaction contenant-contenu entre un emballage à base de polyéthylène haute densité (PEHD) et l'acétate d'amyle considéré comme diffusant. Nous nous sommes plus précisément intéressés à l'impact du phénomène de sorption et de diffusion sur les propriétés mécaniques en sollicitations statique et de fluage de flacons en PEHD à différents stades du vieillissement physique. Une campagne expérimentale a été menée et une base de données issue d'essais menés sur des éprouvettes prélevées sur des flacons ainsi que d'essais de structure réalisés sur des flacons a été constituée. L'analyse des différents résultats expérimentaux a notamment mis en évidence la nature fickienne du phénomène de diffusion et la chute de la résistance à la compression verticale des flacons. Nous avons ensuite proposé la modélisation du comportement mécanique du PEHD avec prise en compte du transfert de masse induit par le phénomène de diffusion. Des modèles élastoviscoplastique et de fluage ont été proposés et implémentés dans le code de calcul en éléments finis ABAQUS. Les paramètres de ces modèles ont été identifiés à partir de sollicitations simples par la résolution d'un problème d'optimisation, et les modèles ainsi identifiés ont été validés par la simulation numérique d'essais réalisés sur des flacons. Les résultats numériques obtenus sont en très bonne adéquation avec les mesures expérimentales. / As part of this thesis, we studied the container-content interaction between a high-density polyethylene-based packaging (HDPE) and amyl acetate considered as diffusing agent. We were specifically interested in the impact of sorption and diffusion phenomenon on mechanical properties of HDPE bottles under static and creep conditions at different stages of physical aging. Several uniaxialtests have been performed on samples cut from bottles and structural tests have been performed on bottles. The analysis of such various experimental results has highlighted in particular the Fickean nature of the diffusion phenomenon and the fall of the vertical compression strength. We have then focused onthe modeling of the mechanical behavior of HDPE taking into account the mass transfer induced by the diffusion phenomenon. Elasto-viscoplastic and creep models have been proposed and implemented in the finite elements code ABAQUS. The parameters of these models have been identified from simple tests by solving an optimization problem, and thus identified models have been validated by numerical simulation of vertical compression of HDPE bottles. Numerical results are in very good agreement with experimental measurements.

Au cours des conflits récents, les véhicules ont été sujet à des attaques impliquant une grande quantité d'explosif enterré. En général, les planchers sont visés, et leurs déformations absorbent une partie de l'énergie. Le véhicule subit également une impulsion, transmise par la détonation, qui peut blesser gravement les occupants. L'intensité de l'impulsion dépend principalement de trois paramètres: le degré de saturation du sol, la nature du sol dans lequel l'explosif est enterré, et la profondeur d'enfouissement de la charge. Les simulations numériques doivent reproduire le processus intégral de l'explosion, incluant l'initiale interaction entre la charge et le sol, l'expansion des produits de détonation, la propagation de l'onde de choc jusqu'au plancher du véhicule, et la grande déformation des projections de sol. L'information mécanique transmise jusqu'au niveau du plancher doit être suffisamment précise car les effets sur le véhicule constituent l'objectif des simulations. Développer un modèle constitutif de sol est par conséquent délicat.%délicat. un défi. tâche difficile. Le modèle de sol est visco-plastique avec un cap écrouissable, une surface de rupture pour limiter les contraintes de cisaillement et un cut off contre les contraintes excessives de tension. La surface de charge est entièrement lisse afin que les incréments de déformation plastique soient continus. Un mécanisme de rigidification est présenté pour tenir compte de la compression de l'air et des très grandes contraintes dans l'environnement immédiat de la charge. Le taux d'humidité du sol est également pris en compte. Deux séries d'essais ont été réalisés à DGA TT pour servir de référence pour les simulations numériques. Une première campagne a mis en jeu de petites quantités d'explosif à plusieurs profondeurs d'enfouissement. Les pressions aériennes ont été mesurées à plusieurs hauteurs et comparées aux simulations. Un moyen d'essai a été utilisé lors de la seconde campagne d'essais pour mesurer à la fois la déformation d'une plaque représentant un plancher de véhicule, mais aussi l'impulsion transmise par l'explosion enfouie. Un certain nombre de taux d'humidité, profondeurs d'enfouissement, épaisseurs de plaques et gardes au sol a été étudié. Au cours des essais, le sol STANAG, composé de sable et de graviers, défini dans l'AEP 55 relatif au STANAG 4569, a été utilisé. Des simulations Eulériennes ont été réalisées et sont en accord avec les résultats d'essais. Le rôle des trois parties de la surface de charge, des mécanismes de rigidification élastique et d'écrouissage plastique a été évalué. La masse volumique initiale du sol et la position initiale du cap sont les paramètres ayant la plus grande influence sur les pressions aériennes. Quant à l'impulsion, elle est principalement contrôlée par le degré de saturation du sol. / In recent conflicts, vehicles have been facing underbelly attacks involving a large quantity of buried explosive. A part of the energy is absorbed by the deformation of the belly. Still the vehicle is subjected to the impulse transmitted by the detonation which may severely injure occupants. The intensity of the impulse is highly dependent on three main parameters which are the degree of saturation of the soil, the nature of the soil in which the explosive is buried and the depth of burial of the charge. Computer simulations should follow the complete process of the explosion, including the early interaction of the charge with the soil, the expansion of the detonation products, the propagation of the shock wave up to the vehicle floor, and the large deformation of the soil projections. The mechanical information transmitted up to the floor level should be sufficiently accurate because the impact on the vehicle is the key target of the simulations. Developing a constitutive soil model is therefore a challenging task. The soil model is visco-plastic with a hardening cap surface, a failure surface to limit shear stresses and a cut off surface against excessive tension. The entire yield surface is smooth so that the plastic strain increments are continuous. A stiffening mechanism is introduced to account for air compression and accommodate the huge stresses in the immediate neighborhood of the charge. The water content of the soil is also taken into account. Two sets of experiments were carried out at DGA TT to serve as a reference for computer simulations. A first campaign involved small quantity of explosive at several depths of burial. The aerial pressures were measured at various heights and compared to simulations. The second set of experiments used a test rig to measure both the deformation of a floor simulating plate and the impulse transmitted by the buried explosion. A range of soil water contents, depths of burial of the explosive, plate thicknesses and stand off distances has been explored. For all the experiments, the sandy gravel STANAG soil defined in AEP 55 of STANAG 4569 has been used. Eulerian simulations were run and fit with these experiments. The roles of the three parts of the yield surface, of the elastic stiffening and plastic hardening mechanisms could be assessed. The initial soil density and the initial cap position were found to have the largest influence on aerial pressures. As for the impulse transmitted to the metallic plate, it was observed to be mostly controlled by the degree of saturation of the soil.

Cette thèse est une étude expérimentale de la formation et de la maturation d’un anneau tourbillonnaire. Elle porte sur les écoulements faiblement inertiels (Reynolds : 5 à 500) en fluide newtonien puis non newtonien. Les anneaux sont générés par un système cylindre-piston. Ils sont analysés par visualisation et par vélocimétrie par images de particules (PIV). La dynamique en fluide newtonien à faible nombre de Reynolds se révèle plus complexe que celle à grands Reynolds avec l’apparition d’un anneau secondaire contrarotatif. Les résultats obtenus en fluide rhéofluidifiant montrent l’influence de l’indice de comportement ainsi que les zones de comportement rhéofluidifiant et newtonien. Les explorations en fluides viscoplastique et viscoélastique montrent la formation d’un ou plusieurs anneaux secondaires contrarotatifs, qui diffèrent dans leur formation et leur dynamique des anneaux observés en fluide newtonien et que l’on associe aux propriétés physiques intrinsèques du fluide. / The present work is an experimental study of the generation and the maturation of vortex rings, in order to characterize their structure and their global dynamics for small to moderate Reynolds numbers (between 5 and 500) in Newtonian and non-Newtonian fluids. The experimental set-up consists of a vertical cylindrical piston-tube system with the lower part immersed in a filled tank. Measurement campaigns have been carried out using dye visualization and Particle Image Velocimetry (PIV). A first part of the work is focussed on Newtonian fluid and allows the dynamics at low Reynolds numbers to be investigated qualitatively and quantitatively. This dynamics turns out to be more complex than the one classically observed at high Reynolds numbers, and is characterized by the production of a counter-rotating secondary vortex ring. The results obtained for shear thinning fluids show the influence of the power-law index on the development and the propagation of the ring. The computation of the shear rate field allows the results to be analyzed in terms of shear thinning and Newtonian regions. Finally a preliminary investigation for viscoplastic and viscoelastic fluids has been performed. In both cases, it is shown that one (for viscoplastic fluids) or several (for viscoelastic fluids) counter-rotating secondary vortex rings are generated, a phenomenon that can be associated with the intrinsic physical properties of the fluid. All these results provide several perspectives of study in the field of vortex rings dynamics in the weakly inertial regimes.

Le sujet de cette thèse porte sur quelques résultats de convergence pour les inéquations variationnelles avec applications dans l'étude des problèmes aux limites décrivant le contact entre un corps déformable et une fondation. La thèse est composée de deux parties. Dans la première partie, nous nous intéressons à l'analyse des inéquations quasivariationnelles, avec ou sans opérateurs de mémoire, dans un espace de Hilbert. Nous prouvons plusieurs résultats de convergence liés à la perturbation de l'ensemble des contraintes ainsi qu'à une méthode de pénalisation. Aussi, pour une classe d'inéquations quasivariationnelles avec opérateurs de mémoire nous étudions une formulation duale pour laquelle nous présentons des résultats d'existence, d'unicité et d'équivalence. La deuxième partie est consacrée à l'application de ces résultats abstraits dans l'étude de six problèmes de contact pour des matériaux élastiques, viscoélastiques et viscoplastiques, dans le cas statique ou quasistatique. Les lois de contact considérées sont la loi de Signorini, la loi de contact avec compliance normale et contrainte unilatérale et la loi de contact avec contrainte unilatérale et seuil critique. Enfin, nous étudions un nombre de problèmes de contrôle optimal associés aux certains modèles de contact. Pour ces problèmes nous obtenons des résultats d'existence et de convergence. / The topic of this thesis concerns some convergence results for variational inequalities with applications in the study of boundary value problems which describe the contact between a deformable body and a foundation. The thesis is divided into two parts. In the first part, we are interested in the analysis of quasivariational inequalities, with or without history-dependent operators, in Hilbert spaces. We prove some convergence results related to a perturbation of the set of constraints and a penalty method, as well. Moreover, for a class of history-dependent quasivariational inequalities we study a dual formulation for which we present existence, uniqueness and equivalence results. The second part is devoted to applications of these abstract results in the study of six contact problems with elastic, viscoelastic and viscoplastic materials, both in the static or quasistatic case. The contact conditions we consider are the Signorini condition, the normal compliance condition with unilateral constraint, the unilateral constraint condition with yield limit. Finally, we study a number of optimal control problems associated to some contact models. For these problems we provide existence and convergence results.

La résolution analytique du problème de la stabilité d’une suspension solide (particules rigides de forme sphérique) dans le système de Taylor-Couette cylindrique a été menée. On s’est basé sur les travaux de Ali and Lueptow (2002) pour formuler les équations régissant la stabilité de l’écoulement dans le cadre d’une théorie linéaire. Ces équations sont valables dans l’approximation du faible espace annulaire et ont pour but la prévision de l’instabilité primaire. A cet effet, nous avons utilisé une méthode variationnelle telle que la méthode de Galerkin pour résoudre le problème aux valeurs propres conduisant à établir le diagramme de stabilité liée au nombre d’onde au voisinage de l’état critique du développement de la première instabilité. Tout d’abord, on a cherché à mettre au point les calculs dans le cas de l’écoulement de Taylor-Couette classique en se référant aux travaux de Ali and Lueptow (2002). Ensuite on a procédé à la résolution systématique des équations du mouvement et l’on cherche à déterminer le critère d’apparition des instabilités en présence de particules en suspension et l’on détermine simultanément les paramètres de couplage entre forces d’interaction liquide-solide. L’ensemble des travaux ainsi réalisés permettront de lever la contradiction fondamentale entre la théorie et l’expérience. L’étude expérimentale a permis d’analyser les effets de la concentration des particules en suspension (disques) et du rapport d’aspect radial ’ sur l’apparition des instabilités dans le système de Taylor-Dean. Le dispositif expérimental utilisé est constitué d’un cylindre intérieur tournant et le cylindre extérieur est maintenu fixe. Le comportement rhéologique du fluide utilisé est viscoplastique obéissant au modèle de Herschel Bulkley. L’apparition des instabilités est examinée en utilisant une technique de visualisation. Pour une concentration donnée dans l’intervalle étudié, la nature des structures apparaissant dans le système d’écoulement dépend ’, alors que pour une valeur donnée de ’ dans l’intervalle étudié, la valeur du nombre de Taylor critique Tac dépend de la concentration des particules. Nous avons obtenu que le nombre de Taylor critique Tac correspondant au déclenchement de la première instabilité évolue non linéairement en fonction de ’. De plus, nous avons examiné expérimentalement les effets de limitation axiale (effet de bords) sur le déclenchement des instabilités dans le système de Taylor-Dean. Les résultats obtenus montrent que les bords tournants n’affectent pas le type de structures qui apparaissent dans le système d’écoulement de Taylor-Dean. Cependant, ils influencent le seuil critique d’apparition des instabilités qui est marquée par des valeurs élevées du nombre de Taylor critique pour des bords tournants ce qui indique un effet stabilisant des bords mobiles. / The analytical solution of the stability problem of a solid suspension (rigid spherical particles) in the system of cylindrical Taylor-Couette was conducted. We are based on the work of Ali and Lueptow (2002) to formulate the equations governing the stability of the flow in a linear theory. These equations are valid in the approximation of small gap configuration and aim to predict the primary instability. For this purpose, we used a variational method such as the Galerkin method to solve the eigenvalue problem leading to establish the stability diagram related to the wave number in the vicinity of the critical state of development of the first instability. First, we develop the calculations in the case of Taylor-Couette flow with reference to classic work of Ali and Lueptow (2002). Then, we carried out a systematic solution of the equations of motion and we search to determine the criterion of onset of instabilities in the presence of suspended particles and coupling parameters are simultaneously determined from liquid-solid interaction force. All work carried out and will remove the fundamental contradiction between theory and experiment. The experimental study has analyzed the effect of the concentration of suspended particles (disks) and radial aspect ratio ' on the occurrence of instabilities in the Taylor-Dean flow system. The experimental device used consists of a rotating inner cylinder and the outer cylinder is stationary. The rheological behavior of the fluid is viscoplastic obeying to Herschel Bulkley model. The onset of instability is examined using a visualization technique. For a given concentration in the range studied, the nature of the structure appearing in the flow system depends on ', while for a given value of ' in the range studied, the value of the critical Taylor number Tac depends on the particle concentration. We obtain that the critical Taylor number Tac corresponding to the onset of the first instability evolves nonlinearly versus '. In addition, we examined experimentally the effect of axial limitation (endwall effects) on the onset of instabilities in the Taylor-Dean flow system. The results show that the rotating endwalls do not affect the type of structures that appear in the Taylor-Dean flow system. However, they influence the threshold of appearance of instabilities which is characterized by high values of the critical Taylor number for rotating endwalls indicating a stabilizing effect of the rotating endwalls.