Le SPS (Spark plasma sintering) ou frittage flash est une technique innovante de compaction de poudre. Ce procédé fait intervenir le courant électrique pour chauffer l'échantillon en appliquant simultanément une pression. Grâce à la vitesse de chauffage, le procédé SPS apparaît comme étant une technologie prometteuse dans le secteur aéronautique servant à produire des matériaux denses à microstructure fine, composés par des intermétalliques difficiles à fondre, à former et à usiner avec les procédés conventionnels. Cependant, la fabrication de formes complexes est problématique à cause des hétérogénéités en densité qui peuvent apparaître lors de la compaction et qui proviennent de la distribution de la température et des contraintes dans la poudre compactée. La distribution du courant, de la température et des contraintes, ainsi que leurs différents effets, font l'objet d'une large étude, étant responsables de l'homogénéité de la microstructure. Une modélisation numérique 3D du procédé est réalisée, dans le cadre de la librairie CimLib. Elle englobe trois problèmes physiques fortement couplés : électrique, thermique et mécanique. Nous utilisons une approche monolithique qui consiste à résoudre une équation pour chaque problème sur un maillage unique représentant outils et poudre. Tout d'abord le couplage électrique-thermique est modélisé et les simulations numériques sont validées. Une loi de comportement viscoplastique compressible s'appuyant sur un modèle d'Abouaf est utilisée pour modéliser la densification de la poudre de TiAl. Ce modèle est validé par plusieurs cas tests de compaction de poudre dans un contexte lagrangien puis eulérien avant de passer à une simulation complète de couplage électrique-thermique-mécanique. Dans ce contexte monolithique, nous développons un modèle pour prendre en compte les effets du frottement entre la poudre et le moule. Enfin, la loi de comportement utilisée est identifiée pour la poudre intermétallique de TiAl. Le frittage par SPS d'échantillons de différentes tailles est simulé. Les résultats en termes de distribution de densité et déplacement sont validés grâce à une comparaison avec l'expérience. / Spark plasma sintering process (SPS) is a breakthrough technology for producing high quality sintered materials. An electric current is applied simultaneously with a vertical load to sinter the powder placed in a graphite mould. Joule effect leads to high heating rates which are favorable for enhancing the microstructure and physical properties. However, manufacturing complex shapes is problematic due to heterogeneities in density distribution that may appeari during compaction. For that reason, the development of a numerical model to predict sintering is necessary. The model should help controlling temperature and stress distributions, which are responsible for the microstructure homogeneity. A 3D numerical model is developed to ensure a predictive tool for SPS using CimLib, a code developed at CEMEF. The numerical model presents three physical problems strongly coupled: an electric problem, a thermal problem and a mechanical problem. A monolithic approach is used which consists in solving one equation for each problem using one unique mesh for tools and powder. First the electric thermal coupling is modeled and the numerical simulations are validated by comparison with commercial codes. A viscoplastic compressible law based on Abouaf model is implemented to model the densification of TiAl powder. This model is validated by comparing the numerical results of different compaction tests with analytic solutions using a Lagrangian and Eulerian framework. Then a fully coupled electric-thermal-mechanical simulation is carried out. In the monolithic framework, a model is developed to take into account friction effects between powder and mould. Finally, the parameters of the selected material law are identified for TiAl powder using our numerical model and SPS experiments. Sintering of different samples is then simulated. Results are compared with the experiments in terms of density distribution and displacement.

Le forgeage à chaud est un procédé de formage des métaux utilisé pour former des matériaux qui sont difficiles à former à froid ainsi que pour réaliser des géométries complexes. La réduction de la limite d’élasticité à haute température et une augmentation subséquente de l’aptitude à la mise en forme constituent le principal mécanisme à l’origine du procédé. Les méthodes numériques constituent un moyen efficace de prédire les états de contrainte / déformation du produit à différentes étapes de la mise en forme. Bien que les méthodes classiques soient suffisamment précises pour fournir une représentation appropriée du procédé, elles ont tendance à être coûteuses en ressources informatiques. Cela limite leur utilisation dans des cas concret, en particulier pour des études d’optimisation du procédé. L’approche pseudo inverse (API), développée dans le contexte du forgeage à froid 2D axisymétrique, fournit une estimation rapide des champs de contrainte et de contrainte dans le produit final pour une forme initiale donnée. Dans ce travail, l’API est étendue pour inclure les effets thermiques et visco-plastiques dans le procédé de forgeage ainsi que dans le cas général 3D. Les résultats sont comparés aux codes commerciaux disponibles basés sur les approches classiques pour montrer l’efficacité et les limites de l’API. Les résultats obtenus indiquent que l’API est un outil assez efficace pouvant être utilisé à la fois pour des simulations 2D et 3D du forgeage à chaud. / Hot forging is a metal forming process used to form difficult-to-form materials as well as to achieve complex geometries. The reduction of yield stress at high temperatures and a subsequent increase in formability is the primary mechanism that drives the process. Numerical methods provide an efficient means to predict the material yield and the stress/strain states of the product at different stages of forming. Although classical methods are accurate enough to provide a suitable representation of the process, they tend to be computationally expensive. This limits its use in practical cases especially for process optimization. Pseudo Inverse Approach (PIA) developed in the context of 2D axisymmetric cold forming, provides a quick estimate of the stress and strain fields in the final product for a given initial shape. In this work, the PIA is extended to include the thermal and viscoplastic effects on the forging process as well as to the general 3D case. The results are compared with commercially available software, based on the classical approaches, to show the efficiency and the limitations of PIA. The results obtained indicate that PIA is a quite effective tool that can be used for both 2D and 3D simulations of hot forging.

L’acier 316L(N) est le matériau de référence pour les structures du circuit primaire des réacteurs nucléaires de quatrième génération, en raison de leur résistance mécanique à la température de fonctionnement, de l’ordre de 550°C. La thèse a permis de développer un modèle polycristallin, capable de prédire le comportement de ces aciers, basé sur la description du glissement viscoplastique des dislocation à haute température, de mise en œuvre simple et avec l’identification d’un nombre de paramètres matériau limité. La démarche de modélisation a été progressive. Lors de la première étape, nous avons proposé et validé une loi d'homogénéisation élasto-viscoplastique à champs moyens, grâce à de nombreux calculs par éléments finis, en considérant des durcissements plastique et des viscosités cristallines. Ensuite, un modèle de viscoplasticité cristalline, reposant sur les lois d’évolution des densités de différents types de dislocations, a été implémenté et les prédictions ont été validés en considérant un très grand nombre de résultats expérimentaux à faible. Le modèle a ensuite été enrichi afin de prendre en compte les mécanismes physiques supplémentaires observés à température élevée, comme la montée des dislocations, le vieillissement dynamique et l’apparition d’une structure de dislocation très hétérogène. Le modèle proposé nécessite uniquement l’ajustement de trois paramètres par identification inverse, utilisant seulement des essais de traction monotone avec saut de vitesse. Les prédictions du comportement mécanique en chargement uniaxial et cyclique sont également en bon accord avec les mesures expérimentales aux températures élevées. / The 316L(N) stainless steels is the reference material for the primary circuit structures of fourth-generation nuclear reactors. This alloy present high mechanical resistance at the operation temperature range of these reactors, of the order of 550 °C. This PhD allowed to develop a polycrystalline model based on the description of the viscoplastic dislocation slip at high temperatures, with straightforward implementation and with identification of a limited number of material parameters. The modeling process was progressive. In a first step, we proposed and validated a mean-field elastic-viscoplastic homogenization law, in comparison to numerous finite element calculations, considering crystalline plastic hardening and crystalline viscosity. Then, a model of crystalline viscoplasticity, based on the evolution laws of the different dislocations densities was implemented and the predictions were validated considering a very large number of experimental results at low temperature. The model was then enhanced to take into account the additional physical mechanisms observed at high temperature, such as dislocation climb, dynamic strain aging and the appearance of a very heterogeneous dislocation structure. The proposed model requires the adjustment of only three parameters by inverse identification, using only monotonic tensile tests at different strain rates. The mechanical behavior predictions in uniaxial and cyclic loading are also in good agreement with experimental measurements at high temperature.

A major instrumented geosynthetic reinforced approach embankment was constructed to 5.5 m elevation above ground, with prefabricated vertical drains, over a soft compressible clay deposit at Leneghan, Newcastle, Australia in May 1995. The field monitoring of settlements for over six years shows that the embankment manifests significant creep. The instrumentation, field performance and the finite element analyses for predicting the long-term performance of this embankment are described in this thesis. The maximum settlement of 1.1 m was observed one year after the completion of construction. However, the embankment continued to settle at a rate of 0.4 mm/day for the next 5 years. The horizontal displacements of 0.09-0.14 m at various locations and the maximum reinforcement strains of 0.67% were recorded. A numerical model was developed to perform a fully coupled large deformation elasto-viscoplastic finite element analysis for this performance prediction based on creep model proposed by Kutter and Sathialingam (1992). The foundation soil was modelled with creep material behaviour using six noded linear strain triangular elements. A well-documented case history ??? Sackville embankment, New Brunswick, Canada was analysed using this model as a benchmark problem and the model was found to predict all the behaviour characteristics reasonably well. The results obtained from finite element analysis using this model are shown to be in reasonable agreement with the observed performance of Leneghans embankment in terms of settlements, horizontal displacements, excess pore pressures and geosynthetic strains. But, the prediction of settlements was less than satisfactory beyond April 1999. Finite element analyses were performed to study the sensitivity of this embankment behaviour on the variation of hydraulic conductivity values and geosynthetic reinforcement properties. This sensitivity study indicated that the kv variation, the kh/kv ratio and the nominal values of geosynthetic properties adopted in the benchmark analysis are reasonable enough for the long-term behaviour prediction.

Αντικείμενο της παρούσας εργασίας είναι η πρόβλεψη τόσο της κίνησης, αλληλεπίδρασης και παραμόρφωσης δύο φυσαλίδων λόγω μεταβολής της πίεσης στο περιβάλλον ιξώδες υγρό, όσο και της ανοδικής κίνησης μιας φυσαλίδας λόγω άνωσης σε ένα Νευτωνικό ή ιξωδοπλαστικό ρευστό. Για τη μοντελοποίηση των αλληλεπιδρώντων φυσαλίδων, αναπτύχθηκε μιας νέα ελλειπτική μεθόδος κατασκευής του υπολογιστικού πλέγματος προκειμένου να αντιμετωπιστούν επιτυχώς τα ιδιάζοντα σημεία (πόλοι) των φυσαλίδων και οι μεγάλες παραμορφώσεις των διεπιφανειών τους. Με τη μέθοδο αυτή η πύκνωση του πλέγματος περιορίζεται μόνο στις περιοχές που είναι αναγκαίο, μειώνοντας έτσι το υπολογιστικό κόστος και αυξάνοντας την ακρίβεια των υπολογισμών. Για την επίλυση των παρακάτω προβλημάτων χρησιμοποιήθηκε η μέθοδος των μικτών πεπερασμένων στοιχείων κατά Galerkin. Στην περίπτωση των αλληλεπιδρώντων φυσαλίδων έχει εξετασθεί η επίδραση του σχετικού μεγέθους τους, της συχνότητας και του εύρους μεταβολής της επιβαλλόμενης πίεσης και πότε οδηγούν σε έλξη ή άπωση των φυσαλίδων. Στην περίπτωση ελκτικής δύναμης, ακολουθείται η κίνηση και η παραμόρφωσή τους μέχρι του σημείου που έρχονται σε επαφή, όπου αυτό είναι εφικτό. Για τη μελέτη του προβλήματος της φυσαλίδας που ανέρχεται λόγω άνωσης, υποθέτουμε αξονική συμμετρία και μόνιμη κατάσταση. Σύγκριση των προβλέψεών μας για το σχήμα των φυσαλίδων και το πεδίο ροής γύρω τους με προηγούμενα θεωρητικά και πειραματικά αποτελέσματα για Νευτωνικά ρευστά έδειξε άριστη συμφωνία. Στην περίπτωση του ιξωδοπλαστικού ρευστού εξετάστηκαν λεπτομερώς οι παραμορφώσεις των φυσαλίδων σαν συνάρτηση των αριθμών Bingham, Bond και Αρχιμήδη και υπολογίσθηκαν οι συνθήκες υπό τις οποίες είναι δυνατή η παγίδευση της φυσαλίδας μέσα σε αυτό. / The present study deals with the numerical simulation of the motion, interaction and deformation of two bubbles due to variation of the pressure of the ambient Newtonian fluid, and the buoyancy-driven rise of a bubble in a Newtonian or a viscoplastic fluid. A new elliptic mesh generation method is developed in order to deal with the singular points (poles) of the bubbles and the large deformations of their surface. This method permits us to increase the mesh resolution only in the regions that is necessary, decreasing thus the computational cost and increasing the precision of our calculations. The following problems are solved using the mixed finite element/Galerkin method. In the case of the interacting bubbles the effect of their relative size, the frequency and the width of the imposed pressure is examined as well as the conditions that lead in attraction or repulsion of the bubbles. In the case that attractive forces exist, the motion and the deformation of the bubbles followed up to the point that they come in contact, whenever this is possible. In order to study the problem of the bubble that rises due to buoyancy, axial symmetry and steady flow is assumed. Our results for the shape of the bubbles and the flow around them are in very good agreement with previous theoretical and experimental results for Newtonian fluids. The deformations of the bubbles rising in a viscoplastic material are also examined for various values of the Bingham, Bond and Archimedes numbers and the conditions under which entrapment of a bubble is possible are determined.

Στόχος της παρούσας διδακτορικής διατριβής είναι η μελέτη χρονικά μεταβαλλόμενων και ασταθών υπό συνθήκες ροών ιξωδοπλαστικών και ιξωδο-ελαστικών ρευστών λόγω της σημασίας που παρουσιάζουν τα υλικά αυτά σε διάφορες βιομηχανικές διεργασίες. Στα πλαίσια της εργασίας αυτής μελετήθηκε αρχικά το πρόβλημα της συμπίεσης ιξωδοπλαστικού ρευστού μεταξύ δύο δίσκων το οποίο είναι ένα βασικό πείραμα ρεολογικού χαρακτηρισμού μη Νευτωνικών υλικών. Η χρονικά μεταβαλλόμενη προσομοίωση επέτρεψε τον προσδιορισμό των σημαντικών διαφορών μεταξύ των δύο εκδοχών του συγκεκριμένου ρεολογικού πειράματος (δίσκοι κινούμενοι υπό σταθερή ταχύτητα ή υπό σταθερή δύναμη) που ήταν αδύνατος με τα προηγούμενα μοντέλα, ενώ παρουσιάζεται και πλήρης παραμετρική μελέτη της διεργασίας αυτής. Το υπόλοιπο μέρος της παρούσας εργασίας αφιερώθηκε στη μελέτη της διεργασίας εκβολής ιξωδοελαστικών υλικών η οποία συναντάται ευρέως στη βιομηχανία μορφοποίησης πολυμερών. Παρουσιάζεται μια πλήρης παραμετρική ανάλυση για την εκβολή ενός ιξωδοελαστικού ρευστού Phan-Tien & Tanner (PTT) από ένα δακτυλιοειδή αγωγό προκειμένου να εξετάσουμε την επίδραση της γεωμετρίας του αγωγού, και των ρεολογικών ιδιοτήτων του ρευστού. Επιπλέον ένα σημαντικό πρόβλημα που εμφανίζεται κατά την εκβολή από ένα κυλινδρικό, επίπεδο ή δακτυλιοειδή αγωγό είναι η εμφάνιση ασταθειών οι οποίες επηρεάζουν σημαντικά την ποιότητα του τελικού προϊόντος. Προκειμένου να διερευνηθεί ο μηχανισμός που προκαλεί την εμφάνιση αυτών των ασταθειών πραγματοποιήθηκε γραμμική ανάλυση ευστάθειας για τη αξονοσυμμετρική και καρτεσιανή δισδιάστατη ροή επικόλλησης-ολίσθησης (stick-slip flow) ενός ρευστού PTT. Γύρω από τη λύση μόνιμης κατάστασης πραγματοποιείται γραμμική ανάλυση ευστάθειας, και υπολογίζονται οι ιδιοτιμές του γενικευμένου προβλήματος με τη βοήθεια της μεθόδου Arnoldi. Έτσι προσδιορίζεται η επίδραση των ιδιοτήτων του υλικού στην ευστάθεια ή μη της ροής. / The purpose of this dissertation is to perform a study on time-dependent and conditionally unstable flows of viscoplastic and viscoelastic fluids because of their importance in various industrial processes. A typical experiment which is used widely for the rheological characterization of non-Newtonian fluids is squeeze flow which is the subject of the first chapter of this dissertation. The time-dependent simulation permitted the determination of the distinct differences between the two versions of this rheological experiment (disks moving with constant velocity or under constant force) which was impossible with the models used up to now, while a complete parametric analysis of this process is presented. The rest of this work focuses in the important problem of the extrusion process of a viscoelastic fluid, which is frequently encountered in the polymer industry. A complete parametric analysis is presented for the extrusion of a Phan-Thien Tanner (PTT) fluid from an annular die in order to examine the effect of the die geometry and the rheological properties of the fluid on this process. Moreover, it is widely known that during the extrusion from a cylindrical, planar or an annular die various flow instabilities may arise affecting significantly the quality of the final product. In order to investigate the mechanism, which causes the appearance of these instabilities, a linear stability analysis is performed for the cylindrical or planar stick-slip flow of a PTT fluid. The stability analysis is performed around the steady state solution and the eigenvalues of the generalized problem are calculated using the Arnoldi algorithm. With this method the effect of the various rheological properties of the fluid on the stability of the flow is determined.

A major instrumented geosynthetic reinforced approach embankment was constructed to 5.5 m elevation above ground, with prefabricated vertical drains, over a soft compressible clay deposit at Leneghan, Newcastle, Australia in May 1995. The field monitoring of settlements for over six years shows that the embankment manifests significant creep. The instrumentation, field performance and the finite element analyses for predicting the long-term performance of this embankment are described in this thesis. The maximum settlement of 1.1 m was observed one year after the completion of construction. However, the embankment continued to settle at a rate of 0.4 mm/day for the next 5 years. The horizontal displacements of 0.09-0.14 m at various locations and the maximum reinforcement strains of 0.67% were recorded. A numerical model was developed to perform a fully coupled large deformation elasto-viscoplastic finite element analysis for this performance prediction based on creep model proposed by Kutter and Sathialingam (1992). The foundation soil was modelled with creep material behaviour using six noded linear strain triangular elements. A well-documented case history ??? Sackville embankment, New Brunswick, Canada was analysed using this model as a benchmark problem and the model was found to predict all the behaviour characteristics reasonably well. The results obtained from finite element analysis using this model are shown to be in reasonable agreement with the observed performance of Leneghans embankment in terms of settlements, horizontal displacements, excess pore pressures and geosynthetic strains. But, the prediction of settlements was less than satisfactory beyond April 1999. Finite element analyses were performed to study the sensitivity of this embankment behaviour on the variation of hydraulic conductivity values and geosynthetic reinforcement properties. This sensitivity study indicated that the kv variation, the kh/kv ratio and the nominal values of geosynthetic properties adopted in the benchmark analysis are reasonable enough for the long-term behaviour prediction.

[pt] A superposição de um escoamento circular de Couette e um fluxo com gradiente de pressão axial, através de um espaço anular ocorre em muitas aplicações práticas, tais como: reatores químicos catalíticos, filtros, extratores líquido- líquido, mancais e o fluxo de retorno de lamas de perfuração entre a coluna de perfuração rotatória e a formação rochosa na perfuração de poços produtores de petróleo e gás. As linhas de corrente curvadas do fluxo circular de Couette podem causar uma instabilidade centrífuga que produz vórtices toroidais, conhecidos como vórtices de Taylor. A presença destes vórtices muda as características hidrodinâmicas e a transferência de calor no processo. Em conseqüência, é muito importante ser capaz de prever o aparecimento da instabilidade. A maioria das análises numéricas e experimentais disponíveis na literatura são para fluidos Newtonianos e viscoelásticos (soluções polimericas). Neste trabalho, o efeito das propriedades viscoplásticas de suspensões de altas concentrações neste tipo de escoamento e nas condições críticas para o aparecimento de vórtices são determinadas teoricamente através da solução das equações de conservação. As equações diferenciais foram integradas pelo método de elementos finitos-Galerkin e o sistema de equações algébricas não lineares resultante foi resolvido pelo método de Newton. / [en] The superposition of a circular Couette flow and a pressure- driven axial flow in an annulus occurs in many practical applications, such as catalytic chemical reactors, filtration devices, liquid-liquid extractors, journal bearings, and the return flow of drilling mud between the rotating drill string and the stationary wall in oil and gas well drilling. The curved streamlines of the circular Couette flow can cause a centrifugal instability leading to toroidal vortices, well known as Taylor vortices. The presence of these vortices changes the hydrodynamic and heat transfer characteristics of the process. Therefore, it is very important to be able to predict the onset of instability. Most of the available theoretical and experimental analyses are for Newtonian and viscoelastic (polymeric solutions) liquids. In this work, the effect of the viscoplastic properties of high concentration suspensions on the onset of the Taylor vortices are determined theoretically by solving the conservation equations and searching the critical conditions. The differential equations were solved by the Galerkin / finite element method and the resulting set of non-linear algebraic equations, by Newtons method.

L’accroissement des besoins en énergie électrique pour les systèmes embarqués et leur augmentation de puissance nécessitent de concevoir des systèmes d’électronique de puissance toujours plus performants. Une solution d’avenir concerne la mise en œuvre de composants à base de diamant qui permettent l’augmentation conséquente des tensions et courants mis en jeux, mais aussi de la température maximale de jonction admissible. Le cadre de ces travaux est celui du projet de recherche Diamonix 2, qui concerne l’étude et l’élaboration d’un composant diamant fonctionnant à haute température. L’objectif du travail doctoral présenté ici est l’étude du packaging haute température de ce type de composant diamant. Plusieurs choix de matériaux et de techniques aptes à l’élaboration d’un assemblage de puce diamant sur un substrat métallisé ont été effectués. La caractérisation microstructurale et mécanique de trois types de jonctions ont été réalisées (refusion d’un alliage AuGe, frittage de nano pâtes d’argent et diffusion en phase solide d’indium dans des couches d’argent). Des essais mécaniques de cisaillement de divers assemblages ont permis d’évaluer le comportement thermomécanique des jonctions et des interfaces. Les essais de cisaillement ont servi à l’identification inverse des paramètres interfaciaux d’un modèle de zones cohésives, pour différents types d’interfaces. Des modèles éléments finis d’assemblage, incluant le comportement viscoplastique des jonctions et des lois d’endommagent des interfaces, ont servi à simuler le comportement thermomécanique du packaging d’un composant diamant. / The increase of electric power demand for embedded systems requires more efficient power electronics modules. A solution to reach this goal relates to the use of diamond-based components that allow high voltage, current density and the maximum allowable junction temperature. The framework is the same as that of the Diamonix 2 research project, which involves the elaboration and the study of a diamond-based die dedicated to high temperature environment. The purpose of the present work is to optimize and simulate the thermomechanical behavior of high temperature diamond die packaging. To reach this goal, the choice of materials that allow high temperature assemblies of diamond die/ceramic substrate was done (AuGe solder alloy, sintering of nano-silver paste, transient liquid phase bounding of indium in silver layers). Microstructural and mechanical characterization of the attachment and the diamond die/junction was realized. Nanoindentation and shear tests are performed for the mechanical characterization. Shear tests results carried out on the two assemblies have been used to identify the interfacial parameters of the bilinear cohesive zone model (CZM) for the diamond die/junction and ceramic substrate/junction interfaces. Finite element modelling of the diamond component packaging including viscoplastic behavior of the junctions and damage law of the interfaces of assemblies were built.

Les écoulements rampants, recirculants et instationnaires d’un fluide viscoplastique autour d’un cylindre ont été étudiés.Numériquement, les morphologies des écoulements, la localisation des zones rigides, les champs de contraintes et pression autour du cylindre ainsi que le coefficient de traînée, ont été déterminés sur un large domaine des nombres de Reynolds et d’Oldroyd.Expérimentalement, les fluides étudiés sont des gels de polymère Carbopol®. Le comportement élastoviscoplastique de ces gels a été modélisé par une loi d’Herschel-Bulkley adaptée. Le montage expérimental conçu et réalisé a été validé par l'étude de l'écoulement d'un fluide newtonien autour d'un cylindre et la mise en place d’une procédure adaptée pour les fluides à seuil.On a pu constater l'influence des conditions d’interface avec l’apparition d’une morphologie de lâchers de tourbillons simultanés et symétriques. / Creeping, recirculating and unsteady flows around a cylinder for yield stress fluid were studied.Numerically, morphologies of the flows, the location of the unyielded zones, the pressure and stress fields around the cylinder and the drag coefficient were determined over a wide range of Reynolds and Oldroyd numbers.Experimentally, fluids studied are polymer gels Carbopol®. The elastoviscoplastic behavior of these gels was modeled by a Herschel-Bulkley adapted law. The experimental setup was designed and validated by studying the flow of a newtonian fluid around a cylinder. An appropriate procedure for the viscoplastic fluid was implemented.We observe the influence of interface conditions with the appearance of a morphology showing simultaneous and symmetrical vortex shedding.